Abstract:
A system and method for obtaining the contact thermal resistance for a sample and pan without a priori knowledge of the properties of either sample, by measuring the reversing heat capacity of a sample and its pan (or an empty pan on the reference side of the DSC) at a long period and a short period during a quasi-isothermal MDSC experiment and then finding the value of contact thermal resistance that makes the short and long period heat capacities match. Several different methods may be used to find the contact thermal resistance using quasi-isothermal MDSC with a long and a short period. Two methods are direct calculation methods that use the results from an MDSC experiment used with model equations to calculate the contact thermal resistance. A third method is another direct calculation method, based upon the phase angle between the heat flow and temperature signals. A fourth and fifth method use curve fitting of the apparent heat capacity for multiple values of pan contact thermal resistance.

Description:
[0001]    The present application incorporates by reference the following patent applications: U.S. patent application Ser. No. 09/533,949 (the “&#39;949 application”), ENTITLED “Heal Flux Differential Scanning Calorimeter,” filed Mar. 23, 2000; U.S. patent application Ser. No. 09/643,870 (the “&#39;870 application”), entitled “Heat Flux Differential Scanning Calorimeter,” filed Aug. 23, 2000; U.S. patent application Ser. No. 09/643,869 (the “&#39;869 application”), entitled “Power Compensation Differential Scanning Calorimeter,” filed Aug. 23, 2000; U.S. patent application Ser. No. 09/767,303 (the “&#39;903 application”) entitled “Differential Scanning Calorimeter,” filed Jan. 26, 2001; and U.S. patent application Ser. No. 09/769,313 (the “&#39;313” application), filed Jan. 26, 2001, entitled “Modulated Differential Scanning Calorimeter.” 
     
    
     
       BACKGROUND  
         [0002]    1. Field of the Invention  
           [0003]    The present invention relates generally to differential scanning calorimeters (DSCs), and more specifically to systems and methods for calibrating the contact thermal resistance between pans and sensors in DSCs.  
           [0004]    2. Background of the Invention  
           [0005]    Differential Scanning Calorimeters measure the heat flow to a sample as the sample temperature is varied in a controlled manner. There are two basic types of DSC, heat flux and power compensation. Brief descriptions of the two types of DSC are included below. A detailed description of the construction and theory of DSCs is disclosed in “Differential Scanning Calorimetry an Introduction for Practitioners”, G. Höhne, W. Hemminger and H. -J. Flammersheim (Springer-Verlag, 1996).  
           [0006]    Heat flux DSCs include a sensor to measure heat flow to a sample to be analyzed. The sensor has a sample position and a reference position. The sensor is installed in an oven whose temperature is varied dynamically according to a desired temperature program. As the oven is heated or cooled, the temperature difference between the sample and reference positions of the sensor is measured. This temperature difference is assumed to be proportional to the heat flow to the sample.  
           [0007]    Power compensation DSCs include a sample and a reference holder installed in a constant temperature enclosure. Each of the holders has a heater and a temperature sensor. The average of the sample and reference holder temperatures is used to control temperature, which follows the desired temperature program. In addition, differential power proportional to the temperature difference between the holders is added to the average power to the sample holder and subtracted from the average power to the reference holder in an effort to reduce the temperature difference between sample and reference holders to zero. The differential power is assumed to be proportional to the sample heat flow and is obtained by measuring the temperature difference between the sample and reference holders. In commercial power compensation DSCs, the difference between sample and reference temperature is generally not zero because a proportional controller is used to control the differential power.  
           [0008]    In both heat flux and power compensation DSCs, a sample to be analyzed is loaded into a pan and placed on the sample position of the DSC. An inert reference material may be loaded into a pan and placed on the reference position of the DSC, although usually the reference pan is empty. The temperature program for conventional DSCs typically includes combinations of constant heating (or cooling) rate and constant temperature segments. Modulated DSC (“MDSC”), which is described in U.S. Pat. No. 5,224,775 (the “&#39;775 patent”), and which is incorporated by reference herein, uses a temperature program in which periodic temperature oscillations are superposed on the constant heating rate and constant temperature segments. The experimental result is the sample heat flow versus temperature or time. The heat flow signal is the result of heat flow to or from the sample due to its specific heat and as a result of transitions occurring in the sample.  
           [0009]    For both heat flux and power compensation DSCs, a temperature difference is created between the sample and reference positions of the DSC during the dynamic portion of the DSC experiment. In heat flux DSCs, the temperature difference results principally from the combination of three differential heat flows: the difference between the sample and reference heat flow, the difference between sample and reference sensor heat flow and the difference between sample and reference pan heat flow. In power compensation DSCs, the temperature difference results principally from the combination of three differential heat flows plus the differential power supplied to the sample holders: the difference between the sample and reference heat flow, the difference between sample and reference holder heat flow and the difference between sample and reference, pan heat flow. The heat flow difference between the sample and reference consists of heat flow due to the heat capacity difference between the sample and reference, the heat flow of a transition, or the difference in heating rate that occurs during am MDSC experiment. The heat flow difference between the sample and reference sections of the DSC is the result of thermal resistance and capacitance imbalances in the sensor or between the holders and the difference in heating rate that occurs between the sample and reference sections of the DSC during a transition or during an MDSC experiment. Similarly, the heat flow difference between the sample and reference pans is the result of mass differences between the pans and the difference in heating rate that occurs during a sample transition or during a MDSC experiment.  
           [0010]    In conventional heat flux DSCs the sensor imbalance and pan imbalance are assumed to be insignificant and the differences in heating rates are ignored. In conventional power compensation DSCs the holder imbalance and pan imbalance are assumed to be insignificant and the differences in heating rates are ignored. When the balance assumptions are satisfied and the sample heating rate is the same as the reference heating rate, the temperature difference is proportional to the sample heat flow and the differential temperature gives an accurate measure of the sample heat flow. The sample heat flow is only proportional to the measured temperature difference between sample and reference when the heating rate of the sample and reference are identical, the sensor is perfectly symmetrical, and the pan masses are identical. Proportionality of sample heat flow to temperature difference for a balanced sensor and pans occurs only during portions of the experiment when the instrument is operating at a constant heating rate, the sample is changing temperature at the same rate as the instrument and there are no transitions occurring in the sample. During Modulated DSC experiments, the heating rates of the sample and reference are generally not the same and the difference between measured sample and reference temperatures is not proportional to the sample heat flow.  
           [0011]    Thus, the sample heat flow from a conventional DSC is not the actual sample heat flow, but includes the effects of imbalances and differences in heating rates; in other words the DSC sample heat flow measurement is smeared. For many DSC experiments, the smeared sample heat flow is a sufficiently accurate result. For example, when the desired experimental result is the total energy of the transition, such as the heat of fusion of a melt, the total peak area is integrated over a suitable baseline and the result from a conventional DSC is sufficiently accurate. If however, partial integration of the peak area is required (for example, in the study of reaction kinetics), the smeared sample heat flow of conventional DSC should not be used. Another example of when the conventional DSC result is inadequate is when two or more transitions in a sample occur within a small temperature interval. In that case, the transitions may be poorly separated in prior art DSCs because of the smearing effects.  
           [0012]    During a transition, the heat flow to the sample increases or decreases from the pre-transition value depending upon whether the transition is exothermic or endothermic and whether the DSC is being heated or cooled. The change in sample heat flow causes the heating rate of the sample to be different from that of the DSC and as a consequence, the sample pan and sensor heating rates become different from the programmed heating rate.  
           [0013]    U.S. patent applications Ser. Nos. 09/533,949 and 09/643,870, incorporated by reference above, disclose a heat flux DSC that uses a four term heat flow equation to account for sensor imbalances and differences in heating rate between the sample and reference sections of the sensor. The four term DSC heat flow equation derived in the &#39;949 application is:  
       q   =       Δ                     T   0     ·     (         R   r     -     R   s           R   r     ·     R   s         )         -       Δ                 T       R   r       +       (       C   r     -     C   s       )     ·            T   s            τ         -       C   r     ·            Δ                   T          τ                                 
 
           [0014]    The first term accounts for the effect of the difference between the sample sensor thermal resistance and the reference sensor thermal resistance. The second term is the conventional DSC heat flow. The third term accounts for the effect of the difference between the sensor sample thermal capacitance and the sensor reference thermal capacitance. The fourth term accounts for the effect of the difference between the heating rates of the sample and reference sides of the DSC.  
           [0015]    U.S. patent application Ser. No. 09/643,869, incorporated by reference above, discloses a power compensation DSC that uses a five term heat flow equation to account for sample and reference holder imbalances and differences in heating rate between the sample and reference holders. The five term power compensation DSC heat flow equation derived in the &#39;869 application is:  
       q   =       Δ                 p     +     Δ                     T   0     ·     (         R   r     -     R   s           R   r     ·     R   s         )         -       Δ                 T       R   r       +       (       C   r     -     C   s       )     ·            T   s            τ         -       C   r     ·            Δ                   T          τ                                 
 
           [0016]    The first term is the difference in power supplied to the sample position versus the power supplied to the reference position. The second term accounts for differences between the thermal resistances of the sample and reference holders. The third term accounts for the heat flow that results from the difference in temperature between the sample and reference holders and is the conventional power compensation DSC heat flow. The fourth term is the heat flow resulting from imbalances in thermal capacitance between the sample and reference holders. The fifth term reflects heat flow resulting from differences in heating rate between the sample and reference holders.  
           [0017]    The contact thermal resistances between the sample and reference pans and the DSC sensor are needed to practice the improved methods of DSC and MDSC disclosed in U.S. patent applications Ser. Nos. 09/767,903 and 09/769,313, incorporated by reference above. Two methods for determining the contact thermal resistance are described in these applications. The first method uses the slope of the onset of a first order transition (typically a melt) to find the contact thermal resistance. The second method applies quasi-isothermal MDSC to a sample of known heat capacity. These two methods can be used to obtain typical values of the contact thermal resistances for use in a contact thermal resistance model equation that uses the thermal conductivity of the pans, sensor, DSC purge gas and three geometric parameters to calculate the contact thermal resistance. Using the typical values found by the above two methods, the geometric parameters may be determined and the model equation may be used in heat flow measurements.  
           [0018]    However, there is an uncertainty in the heat flow measurement resulting from the use of typical values of contact thermal resistance, because the actual contact thermal resistance during an experiment almost never matches the typical value exactly. Ideally the contact thermal resistance should be determined for each pan and used in the DSC experiment that used that pan. The two methods disclosed in the &#39;903 and 313 applications may not be used to determine the contact thermal resistance for any given experiment because they require a priori knowledge of sample properties. Either the properties of a transition in the sample in the case of the onset slope method or the sample heat capacity in the MDSC method must be known. Because these properties are the properties that one wishes to determine from DSC experiments, they are not available to determine contact thermal resistance.  
           [0019]    References  
           [0020]    The following references may be consulted for additional background information:  
           [0021]    A. Boller, Y. Jin and B. Wunderlich, “Heat Capacity Measurement by Modulated DSC at Constant Temperature”, Journal of Thermal Analysis V42 (1994) 307-330.  
           [0022]    B. Wunderlich, “Temperature-Modulated Calorimetry of Polymers with Single and Multiple Frequencies to Determine Heat Capacities as well as Reversible and Irreversible Transition Parameters”,  Material Characterization by Dynamic and Modulated Thermal Analysis Techniques,  ASTM STP 1402, A. Riga and L. Judovits, Eds., ASTM, Conshohocken, Pa., 2001.  
           [0023]    B. Wunderlich, R. Androsch, M. Pyda and Y. K. Kwon, “Heat Capacity by Multi-Frequencies Saw-Tooth Modulation”, Submitted to Thermochimica Acta, September 1999.  
           [0024]    R. L. Danley and P. A. Caulfield, “Baseline Improvements Obtained by a New Heat Flow Measurement Technique”, Proceedings 29 th  N. American Thermal Analysis Society, 2001.  
           [0025]    R. L. Danley and P. A. Caulfield, “DSC Resolution and Dynamic Response Improvements Obtained by a New Heat Flow Measurement Technique”, Proceedings 29 th  N. American Thermal Analysis Society, 2001.  
           [0026]    Definitions  
           [0027]    C p  is the heat capacity;  
           [0028]    T p  is the sample or reference pan temperature;  
           [0029]    Ts is the temperature of the sample sensor;  
           [0030]    Tr is the temperature of the reference sensor;  
           [0031]    q is the sample or reference heat flow;  
           [0032]    q1 is the amplitude of the first heat flow signal;  
           [0033]    q2 is the amplitude of the second heat flow signal;  
           [0034]    q s  and q c  are the magnitudes of the sine and cosine components  
           [0035]    qs1 is the magnitude of the sine component of the first heat flow signal;  
           [0036]    qs2 is the magnitude of the sine component of the second heat flow signal;  
           [0037]    qc1 is the magnitude of the cosine component of the first heat flow signal;  
           [0038]    q c2  is the magnitude of the cosine component of the second heat flow signal;  
           [0039]    T s  and T c  are the magnitudes of the sine and cosine components of the temperature signal, respectively;  
           [0040]    Ts1 is the magnitude of the sine component of the first temperature signal;  
           [0041]    Tc1 is the magnitude of the cosine component of the first temperature signal;  
           [0042]    Ts2 is the magnitude of the sine component of the second temperature signal;  
           [0043]    T c2  is the magnitude of the cosine component of the second temperature signal;  
           [0044]    ω 1  is the angular frequency corresponding to the first period;  
           [0045]    ω 2  is the angular frequency corresponding to the second period;  
           [0046]    R p  is the contact thermal resistance;  
           [0047]    R n  is the nominal value of the contact thermal resistance;  
           [0048]    θ is the product of the angular frequency ω and the elapsed time t of the experiment; and  
           [0049]    Q is the magnitude of the heat flow.  
         SUMMARY OF THE INVENTION  
         [0050]    The present invention is a system and method for obtaining the contact thermal resistance for a sample and pan without a priori knowledge of the properties of either the sample or the pan. The method of the present invention involves measuring the reversing heat capacity of a sample and its pan (or an empty pan on the reference side of the DSC) at a long period and a short period during a quasi-isothermal MDSC experiment and finding the value of contact thermal resistance that makes the short and long period heat capacities match. Several different methods may be used to find the contact thermal resistance using quasi-isothermal MDSC with a long and a short period.  
           [0051]    The first two methods are direct calculation methods that use the results from an MDSC experiment used with model equations to calculate the contact thermal resistance. Both methods start with the measured, or apparent, heat capacity for a sample and its pan, which is given by:  
         C   p     =       q   _       ω                     T   _     p                               
 
           [0052]    Where q is the sample or reference heat flow and T p  is the sample or reference pan temperature and the overbar indicates the amplitude of the signal. The pan temperature is found from:  
             T   p   =T−qR   p ,  
           [0053]    where T is the temperature of the sample or reference sensor and R p  is the sample or reference contact thermal resistance. Both the measured temperatures and the measured heat flows are sinusoidal signals with both sine and cosine components:  
             q=q   s  sin θ+q c  cos θ 
             T=T   s  sin θ+T c  cos θ 
           [0054]    Coefficients q s  and q c  are the magnitudes of the sine and cosine components of the heat flow signal. Coefficients T s  and T c  are the magnitudes of the sine and cosine components of the temperature signal, and θ is the product of the angular frequency ω and the elapsed time t of the experiment. Substituting into the above equation for the pan temperature and collecting sine and cosine components yields:  
             T   p =( T   s   −q   s   R   p )sin θ+( T   c   −q   c   R   p )cos θ 
           [0055]    The magnitude of the pan temperature is the square root of the sum of the squares of the sine and cosine components, which is substituted into the heat capacity equation giving:  
         C   p     =       q   _       ω                         (       T   s     -       q   s          R   p         )     2     +       (       T   c     -       q   c          R   p         )     2                                   
 
           [0056]    The magnitudes of the sine and cosine components of the sensor temperature and the heat flow may be found using the deconvolution technique currently employed in MDSC as disclosed in the &#39;775 patent or any other technique for separating sine and cosine components of a signal. The first two direct calculation methods diverge at this point. The first method can be used to complete the calculation when the longer of the two periods is so long that the dependence on period essentially vanishes (i.e., vanishes for all practical purposes). In that case, the heat capacity may be found accurately regardless of the value of contact thermal resistance because the sensor and pan temperatures become essentially identical. The value of the heat capacity measured during the long period is substituted for C p , and the sine and cosine magnitudes of the short period data are used on the right hand side of the equation. The result is a quadratic equation, which may be solved to find R p .  
         R   p     =         -   b     ±         b   2     -     4      a                 c             2      a                             
 
           [0057]    Where:  
           
         a=q 
         s 
         2 
         +q 
         c 
         2  
       
             b=− 2( T   s   q   s   +T   c   q   c )        c   =       T   s   2     +     T   c   2     -     (         q   s   2     +     q   c   2         ω                   C   p         )                               
           [0058]    A decision must be made as to which of the two roots is correct. The contact thermal resistance must be real and positive. If both roots are positive, the smaller of two positive roots is the correct root. If one root is negative, then the positive root is the correct root.  
           [0059]    The second direct calculation method (which is the preferred method; it is essentially a more general version of the first method) can be used to complete the calculation by using two equations for C p  as above, one each with the sine and cosine components of the temperature and heat flow signals from the long and short period measurements. The heat capacities are eliminated with the result being a quadratic equation in the contact thermal resistance R p . In that case, the coefficients of the quadratic equation are:  
       a   =       q   s2   2     +     q   c2   2     -         (         ω   1          q   2           ω   2          q   1         )     2          (       q   s1   2     +     q   c1   2       )                 b   =     -     2   [         q   s2          T   s2       +       q   c2          T   c2       -         (         ω   1          q   2           ω   2          q   1         )     2          (         q   s1          T   s1       +       q   c1          T   c1         )         ]                 c   =       T   s2   2     +     T   c2   2     -         (         ω   1          q   2           ω   2          q   1         )     2          (       T   s1   2     +     T   c1   2       )           ,                         
 
           [0060]    where the subscripts 1 and 2 indicate the results from the long and short period modulations, respectively, and q is the amplitude of the modulated heat flow (not the sine and cosine components).  
           [0061]    When the long period part of the calibration experiment is very long, both the first and second methods give the same result. In principle, the second method allows any two periods to be used and does not require one of the periods to be so long that the period dependence vanishes. However, the two periods should be chosen far enough apart so that there will be substantial differences in the magnitudes of the signals. The contact thermal resistances can be tested by running a quasi-isothermal heat capacity experiment at multiple frequencies. If the contact thermal resistance is correct, the sample heat capacity and the sample and reference apparent heat capacities will be essentially independent of modulation period.  
           [0062]    A third method is also a direct calculation method. This third method uses the phase angle between the heat flow signal and sensor temperatures, and begins by assuming that the measured heat flow is a cosine function. This assumption is not necessary, but simplifies the calculations. For quasi-isothermal MDSC experiments, when the sample is not reacting or undergoing a transition, the heat flow will be sinusoidal for a sinusoidal temperature modulation:  
           q=Q cos ωt,  
           [0063]    where Q is the magnitude of the heat flow. The heat flowing to the sample and its pan is just the sensible heat associated with the specific heat capacity of the sample and pan:  
       q   =       C   p                 T   p            t                               
 
           [0064]    Setting the equations equal gives the sample pan temperature:  
         T   p     =       Q       C   p        ω          sin                 ω                 t                           
 
           [0065]    As mentioned above, the equation for pan temperature in terms of sensor temperature and the measured heat flow is:  
           
         T 
         p 
         =T−qR 
         p  
       
           [0066]    Substituting the measured heat flow and the pan temperature in terms of the magnitude of the modulated heat flow into the equation for pan temperature and solving for measured temperature yields:  
       T   =         QR   p        cos                 ω                 t     +       Q     ω                   C   p            sin                 ω                 t                             
 
           [0067]    The tangent of the phase angle between the measured heat flow and measured temperature is the ratio of the magnitude of the sine and cosine components of the measured temperature:  
         tan                 φ     =     1     ω                   C   p          R   p                               
 
           [0068]    Thus, the pan contact thermal resistance is:  
         R   p     =     1     ω                   C   p        tan                 φ                             
 
           [0069]    The phase angle between the sensor temperature and the heat flow can be found using the deconvolution method of the original &#39;775 patent, or any other convenient technique that gives the phase angle. The phase angle is measured during the short period part of the calibration experiment and the heat capacity is found from the long period part of the calibration experiment. This calculation always yields a single positive valued contact thermal resistance.  
           [0070]    A fourth and a fifth method for implementing the present invention use curve fitting of the apparent heat capacity for multiple values of pan contact thermal resistance at a short period to determine the contact thermal resistance. The calibration experiments for both of these methods use a long period and a short period quasi-isothermal modulation. In the fourth method, the long period modulation gives values of the sample and reference heat capacities that are assumed to be independent of the pan contact thermal resistance and, hence, independent of period. These values of heat capacity may be considered to be the target values. During the short period part of the calibration experiment, the sensor temperature and the heat flow are separated into their sine and cosine components. Using those signal components, several values (at least four if the function to be fitted is a quadratic, because it has three variables; more than four would be preferable) of contact thermal resistance are used to calculate values of the apparent heat capacity C′ p , using the following equation:  
         C   p   ′     =       q   _       ω              (       T   s     -       q   s          R   p         )     2     +       (       T   c     -       q   c          R   p         )     2                                   
 
           [0071]    In practice, the number of points used depends upon the fitting method and the function to be fitted. Since, in this method calculates these points from one data set, there is no reason to limit the number of points, except in the sense that increasing the number of points would increase the computational time. In any event, if the number of variables in the function is n, then at least n+1 points are needed, with better results being obtained with a greater number of points. In a preferred embodiment, the resulting data is fitted to a quadratic equation:  
           
         C′ 
         p 
         =aR 
         p 
         2 
         +bR 
         p 
         +c  
       
           [0072]    Any suitable routine for curve fitting may be used. For example the “method of least squares” is a suitable curve fitting routine. The contact thermal resistance is calculated from the quadratic formula:  
         R   p     =         -   b     ±         b   2     -     4        a        (     c   -     C   p       )                 2      a                             
 
           [0073]    where C p  is the heat capacity measured at the long period modulation.  
           [0074]    This equation gives two positive roots. The smaller root is the correct value. Occasionally, a real solution cannot be found. When that situation occurs, the maximum value of the pan contact thermal resistance from the quadratic fit is used. The maximum value for the quadratic fit is:  
         R     p                 max       =         -   b       2      a       .                           
 
           [0075]    It has been found that the maximum value typically results in values for the pan contact thermal resistance that yield multiple frequency apparent heat capacities that vary only a few percent, which is a substantial improvement over the methods used in the &#39;903 or &#39;313 applications, which use typical values for the pan contact thermal resistance.  
           [0076]    The fifth method is similar to the fourth method, but uses two arbitrary (although substantially different) modulation periods. The apparent heat capacity versus contact thermal resistance is fitted for the results from both periods. The heat capacities are eliminated between the two equations, yielding a quadratic equation, which is solved to find the contact thermal resistance. The quadratic equation is:  
             R   p   2 ( a   1   −a   2 )+ R   p ( b   1   −b   2 )+ c   1   −c   2 =0  
           [0077]    Where, a, b and c are the coefficients of the fitted data and the subscripts 1 and 2 indicate the coefficient of each of the two different periods. The contact thermal resistance is given by:  
         R   p     =         -     (       b   1     -     b   2       )       ±           (       b   1     -     b   2       )     2     -     4        (       a   1     -     a   2       )          (       c   1     -     c   2       )               2        (       a   1     -     a   2       )                               
 
           [0078]    The correct root must be chosen. It is either the positive real root, or the smaller of two real positive roots. In the case where a complex conjugate pair of roots is found, the contact thermal resistance is taken to be the maximum value of the contact thermal resistance equation:  
         R   p     =       -     (       b   1     -     b   2       )         2        (       a   1     -     a   2       )                               
 
           [0079]    In practice, methods one through five may be applied in one of two ways. In the first way, one or more quasi-isothermal calibration experiments are done at one or more temperatures to determine the pan contact thermal resistance at each of the calibration temperatures. Using this data, the DSC or MDSC experiment is performed using these values of contact thermal resistance to measure the heat flow. When multiple temperatures are used, a curve may be fitted to the data, or interpolation may be used between points. This first way is limited to samples that either do not have transitions within the range of the contact thermal resistance calibration or that have fully reversible transitions.  
           [0080]    The second way for using the pan contact thermal resistance calibration is to use it to adjust the typical values of pan contact thermal resistance used for a given pan type. In this implementation, the contact thermal resistance calibration is performed at the beginning of a thermal program and the contact thermal resistance obtained is compared to the typical value for that pan type and purge gas at the calibration temperature. The ratio of the calibration value to the typical value is used to scale the typical value during the experiment that follows. Several additional calibration steps may be inserted at desired points within the thermal program.  
       
    
    
     BRIEF DESCRIPTION OF THE DRAWINGS  
       [0081]    [0081]FIG. 1 is a schematic drawing of an exemplary heat flux differential scanning calorimeter adapted to practice the present invention.  
         [0082]    [0082]FIG. 2 is a plot showing a comparison of the typical results obtained using the present invention to calibrate the pan contact thermal resistance to the results obtained without calibrating the pan contact thermal resistance. 
     
    
     DETAILED DESCRIPTION OF THE INVENTION  
       [0083]    Apparatus  
         [0084]    [0084]FIG. 1 is a schematic diagram of an exemplary heat flux differential scanning calorimeter adapted to practice the present invention. The differential scanning calorimeter comprises three main components, DSC cell  100 , DSC module  200  and computer  300 . DSC cell  100  comprises a sensor assembly  101  with sample position  102 S and a reference position  102 R. A sample within a sample pan  103 S and a reference within a reference pan  103 R are placed on the sample and reference positions. [Heat is exchanged between each of the pans and its sensor position by a sample thermal contact thermal resistance  104 S and a reference thermal contact thermal resistance  104 R.] Sensor assembly  101  is installed within an enclosure  105  heated by heating element  106 . The temperature of the DSC cell is controlled via power control  205 , according to instructions received from embedded microprocessor  230 . DSC module  200  includes T 0 , ΔT and ΔT 0  amplifiers  220   a,    220   b  and  220   c,  respectively, which receive input from thermocouples  210 ,  211  and  212  as shown in FIG. 1. The output signals of the T 0 , ΔT and ΔT 0  amplifiers are converted from analog to digital signals by A/D converters  221   a,    221   b  and  221   c.  The output of the A/D converters is provided to embedded microprocessor  230 . Embedded microprocessor  230  comprises thermocouple lookup application  231 , sensor coefficient application  232 , contact thermal resistance model  233 , heat flow calculation  234 , temperature control application  235  and temporary data storage  236 .  
         [0085]    Thermocouple Lookup  231  is a program resident in embedded microprocessor  230  that converts the digital signal representing the output signal of the T0 thermocouple to a temperature. The temperature of the ends of the T0 thermocouple wire is measured by a thermistor, and that temperature is converted to the equivalent voltage of a thermocouple at that temperature. The equivalent thermocouple voltage is summed with the output voltage of the T0 thermocouple. The resultant reference junction compensated voltage is converted to temperature by using a table of temperature versus voltage based on NIST monograph  175 .  
         [0086]    Sensor Coefficients  232  is a program resident in embedded microprocessor  230  that determines the sensor coefficients (R s , R r , C s , C r ) used in the heat flow calculation. The temperature of the DSC cell as indicated by the T0 thermocouple is used to determine the appropriate value for each of the coefficients. Sensor coefficients are generated using the calibration procedures disclosed in the &#39;949 and &#39;870 applications and saved in the module in tabular form. The program supplies the sensor coefficients to heat flow calculator  234 .  
         [0087]    Contact Resistance Model  233  is a program resident in the embedded microprocessor that calculates the pan contact thermal resistance using the thermal contact thermal resistance model equation described above.  
         [0088]    Heat Flow Calculator  234  is a program resident in embedded microprocessor  230  that calculates the heat flow using the methods described in the &#39;949 and &#39;870 applications. The ΔT, ΔT0 and T0 signals are used as inputs. Sensor coefficients required by the program are supplied by sensor coefficient program  232  and pan contact thermal resistances needed by the program are supplied by contact thermal resistance model program  233 .  
         [0089]    Temperature Control  235  is a program resident in embedded microprocessor  230  that determines the power to be supplied to DSC heater  106 . It implements a proportional plus integral plus derivative control algorithm.  
         [0090]    Temporary Data Storage  236  is non-volatile RAM within module  200  that stores the results of an experiment during the experiment.  
         [0091]    Embedded microprocessor  230  is in communication over, e.g., an ethernet network  301 , with computer  300 . Computer  300  comprises instrument control interface  302 , data analysis module  303 , data storage module  304 .  
         [0092]    Instrument Control Interface  302  is a program resident in computer  300  that provides the user interface to module  200 . It is used to program the thermal method for the experiment to select any options and to control the instrument, e.g., start and stop experiments, select purge gas flow rates, and select instrument mode: MDSC or standard DSC.  
         [0093]    Data Analysis  303  is a program resident in computer  300  that is used to display and process the results of the experiment. The user may select the signals to be displayed and display options such as axis scaling and selection of the abcissa. Analysis of the results may also be performed, such as integration of the area of a peak to determine the enthalpy of a transition.  
         [0094]    Data Storage  304  is a non-volatile storage of the experimental results, e.g., a hard disk drive.  
         [0095]    The following are step-by-step instructions for carrying out the first three methods described above in the summary of the invention. Although each of the methods is described herein as selecting parameters, such as a calibration temperature, a temperature amplitude, and a modulation period, the present invention is not limited to the manner by which these parameters are selected. For example, there are many different ways of defining a temperature program that is characterized by a calibration temperature, temperature amplitude and a modulation period (or frequency). Any way of defining a temperature program that is characterized by a calibration temperature, temperature amplitude and period (or frequency) will involve selecting a calibration temperature, temperature amplitude and period (or frequency).  
         [0096]    Procedure  
         [0097]    Method 1 (this Method can be used when the Dependence of the Measured Heat Capacity on Period Vanishes for the Longer Period)  
         [0098]    The first method described above can be carried out according to the following exemplary steps:  
         [0099]    1. Select the calibration temperature, the temperature amplitude (typically, 0.5° C.), a long modulation period (typically, 100 seconds), and a short modulation period (typically, 30 seconds).  
         [0100]    2. Allow the DSC cell to equilibrate at the selected calibration temperature.  
         [0101]    3. Hold the temperature at the selected calibration temperature for an isothermal hold period of 5 minutes, to allow the signals to stabilize.  
         [0102]    4. Modulate the temperature of the DSC cell at the long period until the temperature and heat flow amplitudes stabilize.  
         [0103]    5. Calculate the long period apparent heat capacity using the following equation:  
         C   p     =         q   _       ω                     T   _     p         .                           
 
         [0104]    6. Modulate the temperature of the DSC cell at the short period until the temperature and heat flow amplitudes stabilize.  
         [0105]    7. Separate the resulting temperature and heat flow amplitudes into their sine and cosine amplitudes using any convenient mathematical technique. (The preferred method uses the technique disclosed in the &#39;775 patent.)  
         [0106]    8. Use the apparent heat capacity from step 5 and the temperature and heat flow sine and cosine amplitudes from step 7 to calculate the pan contact thermal resistance using the quadratic formula:  
         R   p     =         -   b     ±         b   2     -     4      ac             2      a                             
 
         [0107]    where:  
         
       a=q 
       s 
       2 
       +q 
       c 
       2  
     
           b=− 2( T   s   q   s   +T   c   q   c )  
         [0108]    [0108]       c   =       T   s   2     +     T   c   2     -     (         q   s   2     +     q   c   2         ω                   C   p         )                               
         [0109]    9. Select the smallest positive root as the contact thermal resistance unless complex roots are found. Complex roots occur when:  
           b   2 −4 ac &lt;0.  
         [0110]    in which case, the contact thermal resistance is taken to be the real part:  
         R   p     =         -   b       2      a       .                           
 
         [0111]    10. Calculate the contact thermal resistance correction factor to be applied to the contact thermal resistance function:  
         PF   =     Rp   Rn       ,                         
 
         [0112]    where R p  is the contact thermal resistance determined from step 9 and R n  is the nominal value of the contact thermal resistance function at the calibration temperature as described in the &#39;903 and &#39;313 applications.  
         [0113]    11. Continue with the DSC experiment using PF to adjust the value of contact thermal resistance used in calculating the heat flow.  
         [0114]    Steps 1 through 11 may be repeated several times during a DSC or an MDSC experiment at temperatures selected by the user.  
         [0115]    Note: For DSC experiments, only the sample pan contact thermal resistance is calculated and it is used to calibrate both the sample and reference contact thermal resistances. The contact thermal resistance of the reference pan does not affect the heat flow measurement because the sample and reference calorimeters are independent. For MDSC experiments, both the sample and reference pan contact thermal resistances are calibrated. This requires that steps 5 through 11 are applied to both the sample and reference temperature and heat flow signals.  
         [0116]    Method 2 (Two Arbitrary Periods)  
         [0117]    The second method described above can be carried out according to the following exemplary steps:  
         [0118]    1. Select the calibration temperature, the temperature amplitude (typically, 0.5° C.), a long modulation period (typically, 100 seconds), and a short modulation period (typically, 30 seconds).  
         [0119]    2. Allow the DSC cell to equilibrate at the selected calibration temperature.  
         [0120]    3. Hold the DSC cell at the selected calibration temperature for an isothermal hold period of 5 minutes, to allow the signals to stabilize.  
         [0121]    4. Modulate the temperature of the DSC cell at the long period until the temperature and heat flow amplitudes stabilize.  
         [0122]    5. Separate the resulting temperature and heat flow amplitudes into their sine and cosine amplitudes using any convenient mathematical technique. Store the component values. (The preferred method uses the technique disclosed in the &#39;775 patent.)  
         [0123]    6. Modulate the temperature of the DSC cell at the short period until the temperature and heat flow amplitudes stabilize.  
         [0124]    7. Separate the resulting temperature and heat flow amplitudes into their sine and cosine amplitudes using any convenient mathematical technique. Store the component values. (The preferred method uses the technique disclosed in the &#39;775 patent.)  
         [0125]    8. Use the temperature and heat flow sine and cosine amplitudes from steps 5 and 7 to calculate the pan contact thermal resistance using the quadratic formula:  
           R   p     =         -   b     ±         b   2     -     4      ac             2      a         ,     
          where   :     
                a   =       q   s2   2     +     q   c2   2     -         (         ω   1          q   2           ω   2          q   1         )     2          (       q   s1   2     +     q   c1   2       )                     b   =     -     2        [         q   s2          T   s2       +       q   c2          T   c2       -         (         ω   1          q   2           ω   2          q   1         )     2          (         q   s1          T   s1       +       q   c1          T   c1         )         ]                     c   =       T   s2   2     +     T   c2   2     -         (         ω   1          q   2           ω   2          q   1         )     2            (       T   s1   2     +     T   c1   2       )     .                                         
 
         [0126]    9. Select the smallest positive root as the contact thermal resistance unless complex roots are found. Complex roots occur when:  
           b   2 −4 ac &lt;0  
         [0127]    in which case, the contact thermal resistance is taken to be the real part of the complex roots:  
         R   p     =         -   b       2      a       .                           
 
         [0128]    10. Calculate the contact thermal resistance correction factor to be applied to the contact thermal resistance function:  
         PF   =     Rp   Rn       ,                         
 
         [0129]    where R p  is the contact thermal resistance determined from step 9 and R n  is the nominal value of the contact thermal resistance function at the calibration temperature as described in the &#39;903 and &#39;31 3 applications.  
         [0130]    11. Continue with the DSC experiment using PF to adjust the value of contact thermal resistance used in calculating the heat flow.  
         [0131]    Steps 1 through 11 may be repeated several times during a DSC or an MDSC experiment at temperatures selected by the user.  
         [0132]    Note: For DSC experiments, only the sample pan contact thermal resistance is calculated and it is used to calibrate both the sample and reference contact thermal resistances. The contact thermal resistance of the reference pan does not affect the heat flow measurement because the sample and reference calorimeters are independent. For MDSC experiments, the sample and reference pan contact thermal resistances are calibrated. This requires that steps 5 through 11 are applied to both the sample and reference temperature and heat flow signals.  
         [0133]    Method 3 (Phase Angle Method)  
         [0134]    The third method described above can be carried out according to the following exemplary steps:  
         [0135]    1. Select the calibration temperature, the temperature amplitude (typically, 0.5° C.), a long modulation period (typically, 100 seconds), and a short modulation period (typically, 30 seconds).  
         [0136]    2. Allow the DSC cell to equilibrate at the selected calibration temperature.  
         [0137]    3. Hold the DSC cell at the selected calibration temperature for an isothermal hold period of 5 minutes, to allow the signals to stabilize.  
         [0138]    4. Modulate the temperature of the DSC cell at the long period until the temperature and heat flow amplitudes stabilize.  
         [0139]    5. Calculate the long period apparent heat capacity using the following equation:  
         C   p   ′     =       q   _       ω                   T   p                               
 
         [0140]    6. Modulate the temperature of the DSC cell at the short period until the temperature and heat flow amplitudes stabilize.  
         [0141]    7. Separate the resulting temperature and heat flow amplitudes into their sine and cosine amplitudes using any convenient mathematical technique. (The preferred method uses the technique disclosed in the &#39;775 patent.)  
         [0142]    8. Calculate the phase angle between the heat flow and the sensor temperature using the short period components of heat flow and sensor temperature:  
           tan                 φ     =         q   s       q   c       -       T   s       T   c           ,                         
 
         [0143]    where the subscripts s and c indicate the amplitudes of the sine and cosine components.  
         [0144]    9. Use the apparent heat capacity from step 5 and the phase angle from step 8 to calculate the pan contact thermal resistance according to:  
         R   p     =       1     ω                   C   p        tan                 φ       .                           
 
         [0145]    10. Calculate the contact thermal resistance correction factor to be applied to the contact thermal resistance function:  
         PF   =     Rp   Rn       ,                         
 
         [0146]    where R p  is the contact thermal resistance determined from step 9 and R n  is the nominal value of the contact thermal resistance function at the calibration temperature as described in the &#39;903 and &#39;313 applications.  
         [0147]    11. Continue with the DSC experiment using PF to adjust the value of contact thermal resistance used in calculating the heat flow.  
         [0148]    Steps 1 through 11 may be repeated several times during a DSC or an MDSC experiment at temperatures selected by the user.  
         [0149]    Note: For DSC experiments, only the sample pan contact thermal resistance is calculated and it is used to calibrate both the sample and reference contact thermal resistances. The contact thermal resistance of the reference pan does not affect the heat flow measurement because the sample and reference calorimeters are independent. For MDSC experiments, the sample and reference pan contact thermal resistances are calibrated. This requires that steps 5 through 11 are applied to both the sample and reference temperature and heat flow signals.  
         [0150]    Method 4 (This Method can be used when the Dependence of the Measured Heat Capacity on Period Vanishes for the Longer Period)  
         [0151]    The fourth method described above can be carried out according to the following exemplary steps:  
         [0152]    1. Select the calibration temperature, the temperature amplitude (typically, 0.5° C.), a long modulation period (typically, 100 seconds), and a short modulation period (typically, 30 seconds).  
         [0153]    2. Allow the DSC cell to equilibrate at the selected calibration temperature.  
         [0154]    3. Hold the temperature at the selected calibration temperature for an isothermal hold period of 5 minutes, to allow the signals to stabilize.  
         [0155]    4. Modulate the temperature of the DSC cell at the long period until the temperature and heat flow amplitudes stabilize.  
         [0156]    5. Calculate the long period apparent heat capacity using the following equation:  
         C   p     =       q   _       ω                   T   p                               
 
         [0157]    6. Modulate the temperature of the DSC cell at the short period until the temperature and heat flow amplitudes stabilize.  
         [0158]    7. Separate the resulting temperature and heat flow amplitudes into their sine and cosine amplitudes using any convenient mathematical technique. (The preferred method uses the technique disclosed in the &#39;775 patent.)  
         [0159]    8. Calculate the heat capacity over a range of contact thermal resistances using the temperature and heat flow component amplitudes from step 7 and the following equation:  
         C   p     =       q   _       ω              (       T   s     -       q   s          R   p         )     2     +       (       T   c     -       q   c          R   p         )     2                                   
 
         [0160]    9. Fit a quadratic polynomial to the data generated in step 8:  
         
       C 
       p 
       =aR 
       p 
       2 
       +bR 
       p 
       +c  
     
         [0161]    10. Solve the quadratic to find the value of R p  that corresponds to the heat capacity from step 5.  
         R   p     =         -   b     ±         b   2     -     4        a        (     c   -     C   p       )                 2      a                             
 
         [0162]    11. Select the smallest positive root as the contact thermal resistance unless complex roots are found. Complex roots occur when:  
           b   2 −4 a ( c−C   p )&lt;0  
         [0163]    in which case, the contact thermal resistance is taken to be the real part:  
         R   p     =       -   b       2      a                             
 
         [0164]    12. Calculate the contact thermal resistance correction factor to be applied to the contact thermal resistance function:  
         P                 F     =       R                 p       R                 n                             
 
         [0165]    where R p  is the contact thermal resistance determined from step 10 and R n  is the nominal value of the contact thermal resistance function at the calibration temperature as described in the &#39;903 and &#39;313 applications.  
         [0166]    13. Continue with the DSC experiment using PF to adjust the value of contact thermal resistance used in calculating the heat flow.  
         [0167]    Steps 1 through 12 may be repeated several times during a DSC or an MDSC experiment at temperatures selected by the user.  
         [0168]    Note: For DSC experiments, only the sample pan contact thermal resistance is calculated and it is used to calibrate both the sample and reference contact thermal resistances. The contact thermal resistance of the reference pan does not affect the heat flow measurement because the sample and reference calorimeters are independent. For MDSC experiments, both the sample and reference pan contact thermal resistances are calibrated. This requires that steps 5 through 12 be applied to both the sample and reference temperature and heat flow signals.  
         [0169]    Method 5 (Two Arbitrary Periods)  
         [0170]    The fifth method described above can be carried out according to the following exemplary steps:  
         [0171]    1. Select the calibration temperature, the temperature amplitude (typically, 0.5° C.), a long modulation period (typically, 100 seconds), and a short modulation period (typically, 30 seconds).  
         [0172]    2. Allow the I)SC cell to equilibrate at selected calibration temperature.  
         [0173]    3. Hold the DSC cell at the selected calibration temperature for an isothermal hold period of 5 minutes, to allow the signals to stabilize.  
         [0174]    4. Modulate the temperature of the DSC cell at the long period until the temperature and heat flow amplitudes stabilize.  
         [0175]    5. Separate the resulting temperature and heat flow amplitudes into their sine and cosine amplitudes using any convenient mathematical technique. Store the component values. (The preferred method uses the technique disclosed in the &#39;775 patent.)  
         [0176]    6. Calculate the heat capacity over a range of contact thermal resistances using the temperature and heat flow component amplitudes from step 5 and the following equation:  
         C   p     =       q   _       ω              (       T   s     -       q   s          R   p         )     2     +       (       T   c     -       q   c          R   p         )     2                                   
 
         [0177]    7. Fit a quadratic polynomial to the data generated in step 6.  
         
       C 
       p 
       =a 
       1 
       R 
       p 
       2 
       +b 
       1 
       R 
       p 
       +c 
       1  
     
         [0178]    8. Modulate the temperature of the DSC cell at the short period until the temperature and heat flow amplitudes stabilize.  
         [0179]    9. Separate the resulting temperature and heat flow amplitudes into their sine and cosine amplitudes using any convenient mathematical technique. Store the component values. (The preferred method uses the technique disclosed in the &#39;775 patent.)  
         [0180]    10. Calculate the heat capacity over a range of contact thermal resistances using the temperature and heat flow component amplitudes from step 9 and the following equation:  
         C   p     =       q   _       ω              (       T   s     -       q   s          R   p         )     2     +       (       T   c     -       q   c          R   p         )     2                                   
 
         [0181]    11. Fit a quadratic polynomial to the data generated in step 10.  
         
       C 
       p 
       =a 
       2 
       R 
       p 
       2 
       +b 
       2 
       R 
       p 
       +c 
       2  
     
         [0182]    12. Using the coefficients of the polynomials from steps 7 and 11, solve the quadratic equation below to find R p .  
         R   p     =         -     (       b   1     -     b   2       )       ±           (       b   1     -     b   2       )     2     -     4        (       a   1     -     a   2       )          (       c   1     -     c   2       )               2        (       a   1     -     a   2       )                               
 
         [0183]    13. Select the smallest positive root as the contact thermal resistance unless complex roots are found. Complex roots occur when:  
         ( b   1   −b   2 ) 2 −4( a   1   −a   2 )( c   1   −c   2 )&lt;0,  
         [0184]    in which case, the contact thermal resistance is taken to be the real part:  
         R   p     =       -     (       b   1     -     b   2       )         2        (       a   1     -     a   2       )                               
 
         [0185]    14. Calculate the contact thermal resistance correction factor to be applied to the contact thermal resistance function:  
         P                 F     =       R                 p       R                 n                             
 
         [0186]    where R p  is the contact thermal resistance determined from step 10 and R n  is the nominal value of the contact thermal resistance function at the calibration temperature as described in the &#39;903 and &#39;313 applications.  
         [0187]    15. Continue with the DSC experiment using PF to adjust the value of contact thermal resistance used in calculating the heat flow.  
         [0188]    Steps 1 through 14 may be repeated several times during a DSC or an MDSC experiment at temperatures selected by the user.  
         [0189]    Note: For DSC experiments, only the sample pan contact thermal resistance is calculated and it is used to calibrate both the sample and reference contact thermal resistances. The contact thermal resistance of the reference pan does not affect the heat flow measurement because the sample and reference calorimeters are independent. For MDSC experiments, both the sample and reference pan contact thermal resistances are calibrated. This requires that steps 5 through 12 are applied to both the sample and reference temperature and heat flow signals.  
         [0190]    Results  
         [0191]    [0191]FIG. 2 shows the result of calibration of pan contact thermal resistances using the first method described above. The sample used was 26.25 mg of alumina powder in a nominal 25 mg aluminum pan, with an empty reference pan. A quasi-isothermal MDSC thermal method was executed at 100° C. using the MDSC method described in the &#39;313 application, with periods of 20, 40, 60, 80 and 100 seconds. Two sets of results are shown: the sample heat capacity using the pan contact thermal resistance method of the present invention (C p  Res Cal) and the sample heat capacity without pan contact thermal resistance calibration (C p  Res Uncal). The reduction of variation of C p  with period using the present invention is clear. Thus, using the first method to calibrate the pan contact thermal resistance, the variation of heat capacity was reduced from 2.02 mJ/° C. to 0.28 mJ/° C. However, these results may be considered illustrative of the improvement expected of all of the disclosed methods.  
         [0192]    The foregoing disclosure of the preferred embodiments of the present invention has been presented for purposes of illustration and description. It is not intended to be exhaustive or to limit the invention to the precise forms disclosed. Many variations and modifications of the embodiments described herein will be obvious to one of ordinary skill in the art in light of the above disclosure. For example, selecting a modulation frequency inherently selects a modulation period, and there are many different ways by which a modulation amplitude could be selected. The present invention is not limited to the way in which a temperature amplitude, modulation period, modulation frequency, or calibration temperature or other parameter may be selected, or to the exact equations disclosed herein. The scope of the invention is to be defined only by the claims appended hereto, and by their equivalents.  
         [0193]    Further, in describing representative embodiments of the present invention, the specification may have presented the method and/or process of the present invention as a particular sequence of steps. However, to the extent that the method or process does not rely on the particular order of steps set forth herein, the method or process should not be limited to the particular sequence of steps described. As one of ordinary skill in the art would appreciate, other sequences of steps may be possible. Therefore, the particular order of the steps set forth in the specification should not be construed as limitations on the claims. In addition, the claims directed to the method and/or process of the present invention should not be limited to the performance of their steps in the order written, and one skilled in the art can readily appreciate that the sequences may be varied and still remain within the spirit and scope of the present invention.