Abstract:
An analytical model simulates the propagation of radiation through a coated continuous web where layer thickness and refractive index, as variables, determine the speed and direction of transmitted radiation. The model predicts characteristics of transmitted radiation based on characteristics of incident radiation and initially assigned values for layer thicknesses. Coating thickness(s) are ascertained in a process whereby incident radiation of known characteristics is directed onto a coated web and thereafter, actual measurements of transmitted radiation are compared to predicted characteristics. Using a fitting algorithm, the assigned thickness(es) of the layer(s) of the model are adjusted and the process repeated until the actual and predicted values are within desired limits at which time, the assigned thickness(es) represent the measured calipers. Radiation measurements are obtained using terahertz time-domain spectroscopy apparatus that creates a continuous reference whereby a sample pulses&#39; phase and amplitude can be tracked and corrected.

Description:
REFERENCE TO RELATED APPLICATIONS 
       [0001]    This application claims priority under 35 U.S.C. §119(e) to co-pending application No. 61/596,595 filed Feb. 8, 2012, which is incorporated herein by reference 
     
    
     FIELD OF THE INVENTION 
       [0002]    The present invention generally relates to techniques for determining the thickness, basis weight and related properties of coatings and film layers that employ terahertz spectroscopy measurements. 
       BACKGROUND OF THE INVENTION 
       [0003]    Industrial production of continuous web materials often requires measurement of caliper or basis weight (weight per area unit) for testing and control purposes. For example, in the case of manufacturing coatings, co-extruded plastics, and related multi-layer products, it is often necessary to measure the caliper or basis weight of each of the layers and/or coats. Standard means for determining the caliper of continuous web coatings include, for instance, differential weigh techniques that rely on nuclear or x-ray measurements. This non-contacting scheme employs dual scanners, each equipped with a sensor, that measure the web weight before and after the coating procedure, with the difference being the coating weight. An analogous dual scanner method is the contacting, differential caliper measurement where subtraction of the web caliper before and after the coating procedure yields the coating weight. Finally, infrared radiation detection uses a distinctive IR signature of the coating to detect and measure coating weight. 
         [0004]    The utility of these prior art caliper measurement techniques is restricted by a number of drawbacks. Two-scanner differential systems are prohibitively expensive and require excessive space for accommodating both scanners. Furthermore, the technique itself is not reliable for detecting relatively thin coatings that are applied onto a thick base material. Indeed, dual-scanner measurements often yield calculated coating thicknesses with significant errors. Variability of the base layer is another source of error in dual-scanner systems. To compensate for this variability, dual-scanner systems often include same-spot software that compares measurement results obtained by both scanners from the same spots on the moving web. Unfortunately, the scanners do not remain perfectly aligned and their inability to follow the same spot renders the measurement useless especially when the scanners are measuring a traveling web whose thickness is not uniform in the machine direction. It has been demonstrated that conventional dual scanner systems when measuring a thin polymer coating on a baseboard can exhibit an error that is larger than the measured thickness of the coating. Finally, IR sensing is sometimes impossible to implement due to the high IR attenuation in the measured material or lack of distinctive IR signatures that would allow for differentiation of the coating material. The art is in need of an accurate and repeatable technique for measuring the thickness and related properties of coatings and films that are formed on continuous, traveling non-uniform webs. 
       SUMMARY OF THE INVENTION 
       [0005]    The present invention is based in part on the development of an analytical model that describes the propagation of terahertz radiation pulse through a coated web and especially through the coating and base substrate layers of a continuous web sample material. As radiation travels through the layers, its direction(s) and speed(s) of travel are influenced by, among other things, the thicknesses of the layers and their refractive indices. Accordingly, the mathematical relationships used in the model incorporate layer thickness and layer refractive index as variables to predict how incident radiation behaves as it is transmitted through the coated web. Given that incident and transmitted radiation can be measured and the refractive indices can be initially ascertained or measured, the thickness(s) of the coating(s) in the sample are ascertained in an iterative process whereby incident electromagnetic radiation of known characteristics is directed to a coated web and thereafter, actual measurements are compared to predicted transmitted radiation characteristics. The assigned thickness(es) of the layer(s) of the model are adjusted and the process repeated until the actual measurement and predicted transmitted electromagnetic radiation value are within desired limits at which point, the assigned thickness(es) represent the measured caliper values. 
         [0006]    The caliper measurement is particularly accurate and precise when implemented with terahertz time-domain spectroscopy (THz-TDS) that creates a continuous set of reference pulses whereby a sample pulses&#39; phase and amplitude can be tracked and corrected. The invention can be implemented with an apparatus that preferably employs two sensor heads and is equipped with means for measuring the z-distance between the sensor heads. In operation, a sample being analyzed, such as a travelling web having non-uniform thickness, is positioned between the two sensor heads. By measuring the z-distance and THz signals in real time, the THz signals can be corrected to enhance the spectrometer measurement precision. The non-contact caliper coating measurements are practically insensitive to variations in the basis weight, thickness or other physical parameters of the base substrate such as paper on which the coating is applied. 
         [0007]    Accordingly, in one aspect, the invention is directed to a method of measuring the thickness of a coating on a coated continuous web that includes the steps of: 
         [0008]    (a) developing an analytical model of the coated continuous, web having one or more coating layers and a substrate layer wherein the model expresses the transmission of an electromagnetic field through the coated continuous web as a function of (i) the thickness of each coating layer, (ii) the thickness of the substrate layer, (iii) the refractive indices of each coating layer and of the substrate which values are initially assigned, and (iv) angle of incident of terahertz radiation onto the coated continuous web; 
         [0009]    (b) establishing a sample measurement for the coated continuous web sample by employing a time-domain terahertz spectroscopy apparatus, that includes an emitter that generates pulses of radiation and a detector that receives pulses of radiation, by a technique which includes the steps of: 
         [0010]    (i) positioning a beam splitter along an optical path between the emitter and the detector to form reference radiation pulses and sample radiation pulses that are directed to the sample and the reference radiation pulses do not interact with the sample; 
         [0011]    (ii) positioning a first mirror to reflect sample radiation pulses to the detector; and 
         [0012]    (iii) positioning a second mirror to reflect reference radiation pulses to the detector and using the reflected radiation pulses to monitor fluctuations in phase and amplitude of the reference radiation pulses; 
         [0013]    (c) establishing a reference measurement for the coated continuous web sample with the time-domain spectroscopy by executing steps (b)(i), (b)(ii) and (b)(iii) without placement of any sample for measurement; and 
         [0014]    (d) comparing the sample measurement to the reference measurement and iteratively assigning values for the thickness of the coating and the thickness of the substrate layers based on the model defined in step (a) so as to achieve a best-fit relationship between the sample measurement and the reference measurement to thereby determine the thickness of the coating. 
         [0015]    In another aspect, the invention is directed to an apparatus for analyzing a coated continuous web sample, having one or more coating layers and a substrate layer, to measure the thickness of a coating on the continuous web, that includes: 
         [0016]    means for establishing sample measurement data for the coated continuous web sample that employs a time-domain terahertz spectroscopy device; 
         [0017]    means for establishing reference measurement data for the coated continuous web sample that employs the time-domain terahertz spectroscopy device; and 
         [0018]    electronic data processing means programmed and operatively connected for processing the sample, measurement data and the reference measurement data wherein the data processing means, is programmed to compare the sample measurement data to the reference measurement data and to iteratively assign values for the thickness of the coating and the thickness of the substrate layer based on an analytical model of the coated continuous web so as to achieve a best fit relationship between the sample measurement data and the reference measurement data to thereby determine the thickness of the coating. 
     
    
     
       BRIEF DESCRIPTION OF THE DRAWINGS 
         [0019]      FIG. 1  is an analytical model for a coated board depicted as a four layer optical system; 
           [0020]      FIG. 2  is a representation of an exemplary THz or near THz-TDS transmission-based sheet measurement system; 
           [0021]      FIG. 3  depicts a device for generating continuous referencing (CR) in time-domain spectroscopy; 
           [0022]      FIG. 4  shows a representative continuous CR trace generated by a beam splitter; 
           [0023]      FIG. 5  compares the inter-pulse jitter between two sets of pulses; and 
           [0024]      FIG. 6  are graphs that show the delays between two sets of pulses and corrected delays thereto. 
       
    
    
     DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 
       [0025]      FIG. 1  depicts a model representing the coated board  4  as a multilayer optical system with four layers consisting of air  5 , a coating material  6 , a base material or board  7  and finally air  8 . The layers have respective refractive indices denoted n air , n c  and n b . As an electro-magnetic field E of a terahertz pulse propagates through the system, it encounters three different boundaries between the layers denoted as 1, 2, and 3, respectively, that modify the electric and magnetic fields. 
         [0026]    To describe the THz electro-magnetic field interaction in this system, the ABCD matrix analysis ( Optics , Eugene Hecht, 2d edition, Addison-Wesley Publishing Co., Inc., 1987, p. 373-378) was used and it was assumed that the THz beam can be viewed as a plane wave. Note that the formalism for Gaussian beam can be adapted. Thus, at a boundary, the relation between the electric (E) and magnetic (H) fields is given by the characteristic matrix M, defined by layer i and i+1: 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       [ 
                       
                         
                           
                             
                               E 
                               i 
                             
                           
                         
                         
                           
                             
                               H 
                               i 
                             
                           
                         
                       
                       ] 
                     
                     = 
                     
                       
                         M 
                         i 
                       
                        
                       
                         [ 
                         
                           
                             
                               
                                 E 
                                 
                                   i 
                                   + 
                                   1 
                                 
                               
                             
                           
                           
                             
                               
                                 H 
                                 
                                   i 
                                   + 
                                   1 
                                 
                               
                             
                           
                         
                         ] 
                       
                     
                   
                    
                   
                       
                   
                 
               
               
                 
                   Eq 
                   . 
                   
                       
                   
                    
                   1 
                 
               
             
           
         
       
     
         [0000]    M i  is a 2 by 2 matrix define by the wave number k, the thickness of the layer d, and the angle of the field in the layer as respect to the boundary (θ i+1 ): 
         [0000]    
       
         
           
             
               
                 
                   
                     M 
                     i 
                   
                   = 
                   
                     [ 
                     
                       
                         
                           
                             cos 
                              
                             
                               ( 
                               kd 
                               ) 
                             
                           
                         
                         
                           
                             
                               jsin 
                                
                               
                                 ( 
                                 kd 
                                 ) 
                               
                             
                             
                               
                                 Y 
                                 
                                   i 
                                   + 
                                   1 
                                 
                               
                                
                               
                                 ( 
                                 
                                   θ 
                                   
                                     i 
                                     + 
                                     1 
                                   
                                 
                                 ) 
                               
                             
                           
                         
                       
                       
                         
                           
                             j 
                              
                             
                                 
                             
                              
                             
                               sin 
                                
                               
                                 ( 
                                 kd 
                                 ) 
                               
                             
                              
                             
                               
                                 Y 
                                 1 
                               
                                
                               
                                 ( 
                                 
                                   θ 
                                   
                                     i 
                                     + 
                                     1 
                                   
                                 
                                 ) 
                               
                             
                           
                         
                         
                           
                             cos 
                              
                             
                               ( 
                               kd 
                               ) 
                             
                           
                         
                       
                     
                     ] 
                   
                 
               
               
                 
                   Eq 
                   . 
                   
                       
                   
                    
                   2 
                 
               
             
           
         
       
     
         [0000]    Where j is the imaginary number (to avoid confusion with the index i) and Y i+1 (θ i+1 ′, a function defined by 
         [0000]    
       
         
           
             
               
                 
                   
                     
                       Y 
                       
                         i 
                         + 
                         1 
                       
                     
                      
                     
                       ( 
                       
                         θ 
                         
                           i 
                           + 
                           1 
                         
                       
                       ) 
                     
                   
                   = 
                   
                     
                       
                         
                           ɛ 
                           o 
                         
                         
                           μ 
                           o 
                         
                       
                     
                      
                     
                       n 
                       
                         i 
                         + 
                         1 
                       
                     
                      
                     
                         
                     
                      
                     
                       cos 
                        
                       
                         ( 
                         
                           θ 
                           
                             i 
                             + 
                             1 
                           
                         
                         ) 
                       
                     
                   
                 
               
               
                 
                   Eq 
                   . 
                   
                       
                   
                    
                   3 
                 
               
             
           
         
       
     
         [0000]    When the electric field is not in the plane of incidence, n i +1 is the refractive index of the layer i+1. So, in the case of a multilayer system, the characteristic matrix M of the system is given by the product between all characteristic matrix at each boundaries M i : 
         [0000]    
       
         
           
             
               
                 
                   M 
                   = 
                   
                     
                       
                         M 
                         i 
                       
                        
                       
                         M 
                         
                           i 
                           + 
                           1 
                         
                       
                        
                       
                         M 
                         
                           i 
                           + 
                           2 
                         
                       
                        
                       
                           
                       
                        
                       … 
                     
                     = 
                     
                       [ 
                       
                         
                           
                             
                               m 
                               11 
                             
                           
                           
                             
                               m 
                               12 
                             
                           
                         
                         
                           
                             
                               m 
                               21 
                             
                           
                           
                             
                               m 
                               22 
                             
                           
                         
                       
                       ] 
                     
                   
                 
               
               
                 
                   Eq 
                   . 
                   
                       
                   
                    
                   4 
                 
               
             
           
         
       
     
         [0000]    where the m ij  are the matrix elements. In time-domain terahertz measurement, the THz electric field without any sample or coated board E o  is preferably measured first. Hereafter, E o  will be referred as the incident electric field E i  in the system. Because an application of the model is for representation of the transmitted electric field E t , the relation between the E i  and E t  is given by the transfer function I as follows: 
         [0000]    
       
         
           
             
               
                 
                   t 
                   = 
                   
                     
                       
                         E 
                         o 
                       
                       
                         E 
                         t 
                       
                     
                     = 
                     
                       
                         2 
                          
                         
                           
                             Y 
                             o 
                           
                            
                           
                             ( 
                             
                               θ 
                               o 
                             
                             ) 
                           
                         
                       
                       
                         
                           
                             
                               Y 
                               o 
                             
                              
                             
                               ( 
                               
                                 θ 
                                 o 
                               
                               ) 
                             
                           
                            
                           
                             m 
                             
                               1 
                               1 
                             
                           
                         
                         + 
                         
                           
                             
                               Y 
                               o 
                             
                              
                             
                               ( 
                               
                                 θ 
                                 o 
                               
                               ) 
                             
                           
                            
                           
                             
                               Y 
                               t 
                             
                              
                             
                               ( 
                               
                                 θ 
                                 t 
                               
                               ) 
                             
                           
                            
                           
                             m 
                             
                               1 
                               2 
                             
                           
                         
                         + 
                         
                           m 
                           
                             2 
                             1 
                           
                         
                         + 
                         
                           
                             
                               Y 
                               t 
                             
                              
                             
                               ( 
                               
                                 θ 
                                 t 
                               
                               ) 
                             
                           
                            
                           
                             
                               m 
                               2 
                             
                             2 
                           
                         
                       
                     
                   
                 
               
               
                 
                   Eq 
                   . 
                   
                       
                   
                    
                   5 
                 
               
             
           
         
       
     
         [0000]    The indices o and t for the Y function (Eq. 3) correspond respectively to the initial and final layers where the THz beam propagates. In this case, both layers are air. In a preferred specific application to optimize sensor sensitivity, the THz pulse goes through the coated board twice so that the final transfer function is equation 5 squared. If the board thickness is unknown, preferably a sufficient contrast between the coating and the base material is established in order for the technique to yield good results. The minimal contrast depends on the sensor sensitivity. In one implementation of the model and apparatus, the coating refractive index n c  is preferably kept constant with frequency while the board refractive index n b  is assumed to be like paper (see Applied Optics, vol. 48, no. 33, p. 6541-6546). These refractive indices are used in the equations 2 and 3. Because of the subtlety of the THz time-domain technique in detecting the path length difference to the reference pulse, the phase kd in equation 2 should be adjusted accordingly. Instead of using n c  or n b  alone to define the wave number, (n c −1) or (n b −1) is used which is precisely the path length difference to vacuum or air. With the present invention, by incorporating n c , n b , Eo and Et in a non-linear least squares fitting algorithm, the coating caliper can be extracted accurately. Note that although the incident radiation is shown in  FIG. 1  to be normal to the coating, it is understood that the angle can be an oblique angle and the analytical model adjusted accordingly. Indeed, the angle can be measured initially and the appropriate model that is already programmed into a processor is selected to calculate the coating thickness. Once the thickness is determined, the basis weight and related properties can also be calculated using standard techniques with density data. The caliper coating measurement technique is particularly suited for gauging continuous webs where each layer in the coating or substrate layer is fabricated from paper, board, plastic, polymer, natural rubber, metal, natural fiber and/or synthetic fiber. 
         [0027]      FIG. 2  illustrates a continuous real time referencing generating device  10  that is positioned between an emitter  12  and a detector  14  of a time-domain spectrometer. Device  10  is particularly suited for establishing sample and reference measurements for a coated sample. The device includes a beam splitter  16  which splits incidental light  18  from emitter  12  into reflected light  22  and transmitted light  20 . The transmitted light  20  is reflected back as light  21  by flat mirror  30  through beam splitter  16  before being received by detector  14 . Reflected light  22  is transmitted through a sample  40  and thereafter reflected as light  23  by flat mirror  32  back toward the beam splitter  16  where it is transmitted into detector  14 . In the case where this device  10  augments a terahertz spectrometer, the emitter and detector can be, for example, photoconductive antennae. For terahertz radiation, a preferred beam splitter comprises a thick piece of high resistivity (&gt;10,000 O-cm) silicon slab. The thickness of the beam splitter is selected to be large enough so that reflections from the back surface thereof lie outside the measurement window. If a thinner beam splitter is used, the pulse shape will look different since multiple reflections from the backside are included in the window. In principal either configuration can be employed with the present invention. A thicker beam splitter is preferred because it is less prone to vibrate. A high resistivity silicon slab is particularly suited for use as a beam splitter as it has a high refractive index with acceptable absorption in terahertz frequencies. Alternatively, sapphire or polyethylene could be used as the beam splitter although they have higher absorption than that of silicon. 
         [0028]    In a preferred embodiment, device  10  includes an upper head or module  51  that houses emitter  12 , detector  14  and mirror  30  and a lower head  52  that houses mirror  32 . The gap or displacement distance “z” between sensor heads  51 ,  52  through which the coated non-uniform continuous web sample  40  traverses can vary particularly when the dual sensor heads are in motion as part of a scanner. To account for this z “wander” and therefore the varying length of radiation traveled, the gap separation can be continuously measured. Dynamic measurements can be achieved with conventional devices, such as, for example, a displacement sensor  26 , which employs inductive or magnetic measuring devices. Distance signals from z measurements can be communicated to a processor that calculates the time delay based on the z measurements and generates time delay signals. While the terahertz time-domain spectroscopy-based device  10  is illustrated in the transmission mode, it is understood that the invention is also applicable in the reflection mode. 
         [0029]    In operation, laser pulses are directed to emitter  12  and detector  14 , respectively. (A continuous wave photoconductor can be used instead.) Initially, a terahertz radiation pulse is generated by the photoconductive antenna emitter and is incident on silicon beam splitter  16  at a 45 degrees angle resulting in two pulses traveling in perpendicular directions. One pulse is transmitted through sample  40  whereas the other pulse is used to track the fluctuations in time and amplitude. Since the two pulses are generated at the same time and position, their phases and amplitudes correlate very strongly and this correlation is used to correct measurement errors. This method is particularly suited for fiber optic based time domain terahertz spectrometers where the bending of fibers and temperature fluctuations could have strong adverse effects in precision in amplitude and phase measurements. Only one detector is required although multiple detectors can be employed. 
         [0030]    The continuous real time referencing generating device  10 , configured as depicted in  FIG. 2 , but without separate heads  51  and  52 , was incorporated into a THz or near THz-TDS transmission-based sheet measurement system  100  as shown in  FIG. 3 . It has been demonstrated that with this THz-TDS system  100 , jitter limits the measurement precision of the system. It is expected that implementation of the constant referencing technique will significantly enhance its precision. System  100  comprises a near THz or THz generator including at least one pulsed laser source (e.g. femtosecond (fs) laser)  105  for emitting a beam of optical pulses. A beam splitter  106  splits the optical beam into two beams, a reflected beam  102  and a transmitted beam  103 . The reflected beam  102  is directed to reflective delay comprising optics  108  including a “dummy” delay stage. The purpose of the dummy delay is to make both the source (transmitter) and receiver (detector) arms of the THz-TDS system  100  have nominally equal optical path lengths; this results in the source and receiver fs-pulses being derived from the same original fs-pulse. The intention of the dummy delay is to minimize noise. However this dummy delay is not generally required and the THz-TDS system  100  can be also generally operated without it in certain applications. 
         [0031]    The transmitted beam  103  is directed via mirrors  104  to delay comprising optics  109  shown as linear delay stage. The delay comprising optics  108  and  109  are configured to make the optical path length of the reflected beam  102  to the detector  110  be nominally equal to the optical path length of the transmitted beam  103  to the near THz or THz transmitter  111 . 
         [0032]    The near THz or THz transmitter  111  includes a transmit antenna operable for emitting THz or near THz radiation pulses having a frequency between 0.05 THz and 50 THz at a sample location on the sheet material. The THz transmission antenna will generally have a bias voltage applied to it (not shown), which can be modulated be used for the utilization of a lock-in detection scheme in system  100 . The detector  110  includes a receiver antenna operable to receive near THz or THz radiation transmitted by the location on the sample  130  that is irradiated by the incident radiation. Embodiments of the invention are not limited to using photoconductive antennae as described herein. For example, other methods such as optical rectification using crystals such as Zinc Telluride (ZnTe) may also be used. The THz detector  110  is coupled to the receive antenna and is also coupled to receive the delayed optical pulses from delay comprising optics  108  for synchronously detecting the THz or near THz radiation transmitted by the sheet material sample  130 . The detector  110  generates electrical detection signals. Although  FIG. 3  shows a THz focused beam interacting with the sample  130 , it is understood that other optical geometries such as a collimated geometry can also be used. 
         [0033]    A signal processing system  125  is coupled to the detector  110  to receive the electrical detection signals. The signal processing system  125  comprises memory  126  Signal processing system  125  also includes processing electronics  128  that generally includes a transimpedance (current to voltage)-amplifier, filter and analog to digital (A/D) converter. A processor  127  receives processed electrical signal (amplified, filtered and converted to a digital signal) from processing electronics  128 . The processor  127  can combine a signal associated with the transmitted pulse together with the calibration data and a reference signal pulse to determine at least one property of the sheet material sample selected from the moisture content, basis weight and caliper. 
         [0034]    With respect the reference pulses and to techniques for applying the correction, processor  127  in signal processing system  125  can compare the continuous set of reference pulses with sample pulses and execute the necessary corrections for time jitter and amplitude variations. In particular, process  127  calculates the sample and reference measurement data for the coated web sample and is programmed to compare the sample measurement data to the reference measurement data. By iteratively assigning values for the thickness of the coating and the thickness of the substrate layer based on the analytical model of the coated continuous web a best fit relationship between the sample measurement data and the reference measurement data is achieved to calculate the caliper of the coating. 
         [0035]    The detection electronics generally utilizes a lockin detection scheme (not shown) in which a modulated bias voltage is applied to the transmitter antenna. This modulation signal is also fed to lock-in detection electronics that increases the signal-to-noise ratio of the system as well as minimizes the effect of any background signal. A mechanical chopper can be used in the source femtosecond laser beam  105  to realize the lock-in detection—in the case a DC bias voltage is applied to the antenna. 
         [0036]    In implementing the inventive technique for measuring caliper coatings, the refractive indices of the measured materials (base and coating) are determined from published data or they can be measured with a conventional device. In addition, it may be useful to add another sensor to measure a specific parameter of the base layer (or any other layer). For instance, an infrared gauge can be use to measure the moisture of the base material to correct the refractive index. 
         [0037]    Thereafter, the reference measurement is taken wherein no sample in used in device  10  for generating continuous referencing shown in  FIG. 2 . The reference signal may be corrected using the Z-sensor signal. The coated, traveling continuous web is positioned between the gap of the device and sample measurements are taken. Measure signals from the detector are captured and saved for future analysis. The measure signal may be corrected using the Z-sensor signal. The fitting algorithm is applied in which the delay and amplitude of the measuring signal versus the reference signal and refractive indexes of materials are used in the algorithm to determine the thickness of the base and coating material. Thickness of the coating material is supplied as a result of the procedure. 
         [0038]    With THz-TDS transmission-based sheet measurement system  100  of  FIG. 3 , for each pulse that is generated by the emitter of the terahertz spectrometer, the beam splitter generates two pulses one that passes through the sample of interest and the other pulse traveling through air.  FIG. 4  depicts two representative pulses that are formed by the beam splitter. Each pulse has a characteristic profile and peak and, the left pulse is designated the “left” pulse and the right one the “right” pulse. An arbitrary point somewhere in the middle to separate the two can be selected. In the above analysis, the “right” pulse is used both as a reference (i.e., where the pulse travels through air) and as a sample pulse (i.e., where the pulse travels through a sample). These are two separate measurements conducted one after the other. 
         [0039]    One method of implementing the continuous referencing procedure is to initially collect a set of reference pulses. In this set, the two pulses (“left” and “right”) travel through air. Thereafter, a sample of interest is inserted in one of the arms of the spectrometer (i.e., the right arm) and a new data set is collected. All analysis can now be performed using the reference and sample “right” pulse measurements. The left pulse is used to trace the changes in time and amplitude of the right pulse. Since any changes in the “left” pulse will most likely manifest in the “right” pulse, the “left” pulse could be monitored and used to correct the “right” pulse and thereby reduce the effects of jitter. 
         [0040]      FIG. 5  shows the inter-pulse time jitter between the “left” and “right” pulses for 100 different measured pulses. The “index” designation on the x-axis refers to the experiment measurement number (1-100) so that each point on the graph corresponds to a set of left and right pulses. The plot shows the pair wise time delay between the left and right pulses. This correlation plot further establishes that the two pulses move synchronously which permits the correction of time jitter in one of the pulses. The comparison demonstrates that the standard deviation of jitter could be reduced by a factor of 2. 
         [0041]    For fiber optics based terahertz time-domain spectrometers where precise phase measurements are particularly critical, the continuous referencing procedure can compensate for the adverse effects of temperature and other factors that otherwise would interfere with phase measurements. In particular, for a fiber optics based spectrometer, a polarization and/or path length change occurs as a result of fibers being twisted or stretched that occurs with movement of the sensors as they scan the sample. The top plot of  FIG. 6  is the inter-pulse time jitter between the “left” and “right” pulses for 200 different measured pulses as the fiber connections of the terahertz time-domain spectrometer were manually twisted to simulate the effects of movement for a scanning type system typically used in a paper machine. The plot shows that, the polarization change affects the delay of the pulses. The lower plot of  FIG. 6  demonstrates that the inventive technique can be employed to track and correct these delays as the two pulses are affected similarly and one could be used to track the changes in the other. 
         [0042]    Once the continuous set of reference pulses is obtained, the sample pulses&#39; phase and amplitude can be tracked and corrected, as necessary. For example, time jitter can be corrected using the reference pulses by adding or subtracting a delay to the sample pulse. To correct for amplitude variations with the reference pulses, scaling of the sample amplitude by multiplicative factor can be used to increase or decrease its intensity. 
         [0043]    The foregoing has described the principles, preferred embodiments and modes of operation of the present invention. However, the invention should not be construed as being limited to the particular embodiments discussed. Thus, the above-described embodiments should be regarded as illustrative rather than restrictive, and it should be appreciated that variations may be made in those embodiments by workers skilled in the art without departing from the scope of the present invention as defined by the following claims.