Abstract:
The concentration of optically active substances is determined by measuring the polarization of light passed through the substance using a light source, a polarizer and a modulator operated at a predetermined modulation frequency f F , as well as a measuring cell, an analyzer and a detector. The output signal from the detector 16 is alternately integrated during each half period of the modulation frequency to form voltage values U I+  and U I- . These voltage values are stored and the quotient Q of these stored voltage values is formed. Subsequently, an adjustable constant voltage is subtracted from the quotient. The measuring arrangement can have a miniature design for implantation. No mechanically moving parts and only a single detector are used. The measurement is accomplished by a simple, miniature linear ray path having high resolution.

Description:
FIELD OF THE INVENTION 
     The invention relates to polarimeters for measuring the concentration of optically active substances and, more particularly, polarimeters which can be used to measure the glucose concentration in body liquids. The apparatus measures the polarization of polarized light fed to a detector through the substance to be analyzed. The polarized light passes through a modulator which is operated at a predetermined modulation frequency. From the modulator, the light passes through a measuring cell, which contains the substance, and an analyzer before reaching the detector. 
     BACKGROUND OF THE INVENTION 
     A prior known polarimeter operates in accordance with the principle of automatic optical null balance. The prior polarimeter contains two line radiators, a mercury and a sodium vapor lamp, as well as five optical filters arranged on a filter wheel. A tilting mirror is coupled to the filter wheel shaft. The mirror allows the radiation of the corresponding light source to be inserted into the ray path of the polarimeter at the same time the measuring wavelength is selected. Further, a quartz-iodine and a deuterium lamp are provided to furnish continuous radiation. A grid monochromator is provided to select the desired wavelength. 
     In this known design, the monochromatic light travels through the polarizer, the cell containing the sample, and the analyzer before arriving at a photo multiplier. The polarizer and analyzer are formed by rotatably arranged &#34;glan&#34; prisms made of calcite. The polarizer, including the respective plane of polarization of the transmitted light, vibrates about the optical longitudinal axis of the system at an excursion of about ±0.7 radians at 50 Hz. The 50 Hz signal is generated in the photo multiplier, in the unbalanced state of the system. The 50 Hz signal is amplified and fed, with the correct sign, to the power input of a servo motor. The servo motor is mechanically coupled to the analyzer for rotating the analyzer until the 50 Hz signal becomes zero. 
     An optically active sample inserted into the ray path causes rotation of the plane of polarization. The analyzer is turned by means of the servo system into the new balanced position. The angular difference between the new and the original balanced position corresponds to the optical capability of rotation of the sample. This type of apparatus, such as the Perkin Elmer model 241 MC, requires mechanical movement of the polarizer, a mechanically acting modulator and a motor rotated analyzer thus making the device unsuitable for miniaturization and/or implantation. 
     A further prior known method for determining the blood glucose using polarimetry operates by feeding polarized light via a Faraday modulator, operated at a predetermined modulation frequency, through the sample and an analyzer to the detector. In this known method, a further detector is provided because the ray is divided into both measuring and reference signals. Subsequently, the quotient is formed from the two output signals. However, methods using ray division, such as that shown in German Patent DE-OS No. 2944113, require suitable outlays for the optical system. Particularly, the two detectors must meet synchronization requirements. 
     There is therefore a need for a method and apparatus for measuring the concentration of optically active substances while allowing miniaturization and consequently, implantation into a living body. 
     SUMMARY OF THE INVENTION 
     The present invention solves this problem through the use of clocked, integrated signal processing in conjunction with a modulator and a single ray path to obtain a resolution of 5×(10 -4 ) for an angle of rotation α M  of the plane of polarization. The apparatus of the present invention further operates without the use of any moving parts. For example, a relative resolution of 4% is obtained for a cell length of 2.7 cm and a resolution for a glucose concentration C of 34 mg/l (referred to the physiological standard value C 0  =800 mg/l). 
    
    
     BRIEF DESCRIPTION OF THE DRAWINGS 
     FIG. 1 schematically illustrates the design and the circuit for carrying out the method of the present invention. 
     FIGS. 2-2H are diagrams serving to explain the operation of FIG. 1. 
     FIG. 3 shows a preferred embodiment of a polarimeter according to the invention. 
     FIG. 4 is a diagram showing the photo current waveform as a function of the measured angle and the modulation angle. 
     FIG. 5 schematically illustrates an embodiment of the two memory contents in the processor of the present invention. 
    
    
     DETAILED DESCRIPTION 
     Referring to FIG. 1, there is shown an apparatus for carrying out the method of determining the glucose concentration in body liquids, particularly in blood serum. A light source 2, preferably a light-emitting diode LED, emits light which passes through a polarizer 4, and a modulator 6 such as a Faraday modulator. This modulator 6 includes a coil 7 and a crystal 8. The apparatus further includes a lens 10 which may, for example, be made of glass or plastic, a measuring cell 12 for containing a liquid substance to be measured, an analyzer 14 and a detector 16. In the measuring cell 12, the plane of vibration of the light is rotated by an angle α M  because of the optically active substance (not shown). The light is at least approximately linearly polarized. The analyzer is adjusted to a permeability minimum with the modulator 6 disconnected and in absence of any optically active substance in the measuring cell 12. 
     An oscillator 20 controls the light source 2, a frequency divider 21 and a counting circuit 22 with its output signal U 20  having a frequency f 0 . The modulator 6 and the counting circuit 22 are controlled by the modulator frequency f F  of the output signal U 21  from the frequency divider 21. One output signal U 22  from the counting circuit 22 is sent as a logical control signal for an integrator 26. An output signal U 16  from the detector 16 is fed, preferably via an amplifier 25 as input to the integrator 26. Two further output signals U 23  and U 24  from the counting circuit 22 are fed into one input of memories 27 and 28, respectively. The memories 27, 28 also have a second input which receives the output signal U 26  from the integrator 26. The memories 27 and 28 are preferably of the sample-and-hold type. The output signals U 27  and U 28  from the memories 27 and 28 are fed into a divider 30. The divider 30 further has an additional summing input. This additional input is provided with a given constant voltage U K  from a constant-voltage source 29. 
     In operation of the embodiment of FIG. 1, the measuring cell 12 receives a liquid substance that is to be measured. It should be noted that the invention can also be used for measuring a solid substance, for example, an organic substance having asymmetric carbon bonds, or any general optically active crystal such as quartz, cinnabar, or one of the sulfuric acid double salts from the alkali metals. The solid substance would be placed in the measuring cell 12 for carrying out the measurement. 
     Referring to FIG. 2, there is shown a diagram plotting the output voltage U 20  from the oscillator 20 as a function of the time t. The oscillator 20 supplies a constant output frequency f O . 
     According to the diagram of FIG. 2A, there is shown the output signal from the frequency divider 21. The output signal U 21  has a modulation frequency f F  which is preferably at least one, and in particular, about two orders of magnitude smaller than the oscillator frequency f 0 . 
     The photo current U 16  from the detector 16 is shown in the diagram of FIG. 2B. The diagram indicates the noise of the photo current, or the voltage U 16  proportional to the noise of the photo current, as well as the different mean amplitudes of the photo currents during both phases of the Faraday current. Using a vanishing angle of rotation α M  =0 within the measuring cell 12 allows both amplitudes to have the same value. 
     After amplification of the photo current U 16  and its integration, an output signal from the integrator 26 is obtained as shown in the diagram of FIG. 2D. The output signal U 26  is reset, at times t 3  and t 7  by resetting pulses U 22  having a length τ R  from t 3  to t 4  and from t 7  to t 8  as shown in FIG. 2C. 
     According to the diagrams of FIGS. 2E and 2F, at times t 1  to t 2  and t 4  to t 6 , the control pulses U 23  and U 24  from the counting circuit 22 alternately switch the memories 27 and 28 into the sampling state. 
     According to the diagrams of FIGS. 2G and 2H, the output voltages U 27  and U 28  from the memories 27 and 28 change only during corresponding sample pulses τ S  + and τ S  -. The difference or the quotient Q from output voltages U 27  and U 28  furnish, after subtracting the constant volta U K , the measuring signal U A . The measuring angle U A  is proportional to the angle of rotation α M  of the plane of polarization and thereby, to the concentration of the optically active substance contained in the measuring cell 12. 
     This embodiment of the measuring apparatus produces a resolution for the angle of rotation of the plane of polarization given by the equation: 
     
         Δα.sub.M =5×(10.sup.-4) 
    
     Referring to FIG. 3, there is shown another embodiment for an implantable glucose measuring polarimeter. In this embodiment, the light source 2, polarizer 4 and Faraday crystal 8 are cemented to each other. Imaging of the radiating crystal surface from the light source 2 on the detector 16 is accomplished exclusively by the lens action obtained by the convex back surface 9 of the Faraday crystal 8. The transmission of the surface 9 is preferably enhanced by a λ/4 layer of silicon oxide SiO or a more highly oxidized mix of SiO x  with 1&lt;x&lt;2. The measuring cell 12, having a length L, for example, L=27 mm, is closed off on its input side with a window 32. The window 32 is preferably a neutral glass window. The output side of the measuring cell 12 is closed off by the analyzer 14. Two holes 34, 35 serve to allow the measuring liquid to flow through the cell 12. The liquid typically may be a filtered blood serum. Two pairs of electric terminals 36, 37 and 38, 39 are provided for coupling to the light source 2 and the detector 16, respectively. The terminals used for the Faraday coil 7 are not shown in the figure for simplification. The indicated ray bundle is limited by the maximum usable angle opening 2 ν max  for the radiation emitted by the light source 2. 
     In the embodiment of FIG. 3, in which the light source 2, the polarizer 4 and the Faraday crystal 8 are cemented together by immersion layers, practically no reflection losses are obtained. Further, a mechanically strong assembly results. 
     A relatively small, inexpensive and stable light source 2, preferably a light-emitting diode LED, can be provided which furnishes monochromatic radiation which can be modulated at high frequency. 
     Because the light source 2 can be imaged on the detector 16 through the use of lenses, the maximum amplitude of the photo current and therefore, a favorable signal-to-noise ratio, are obtained. The convex back surface 9 of the Faraday crystal 8 produces an immersion lens action with respect to the light source 2. This allows the utilization of the high index of refraction, for example, n=3.4 for gallium phosphide GaP when using a wavelength λ=575 mm (yellow). By this means, the cost of the lenses for imagining the light source 2 on the detector can be saved altogether or at least partially reduced. A semiconductor laser diode is also suitable for use as the light source because of its high radiation output and the maximum measuring sensitivity achievable thereby. 
     Preferably suited for the optical medium of the Faraday crystal 8 in the modulator 6 is an A III  B V  semiconductor crystal. These compounds have a large Verdet constant V, i.e., a strong Faraday effect is obtained with a relatively small magnetic field in the Faraday coil 7. An especially well suited compound is gallium phosphide having a Verdet constant V=4.5×10 3  min/(T×cm) at a wavelength λ=665 nm and a constant V=6.7×10 3  min/(T×cm) for a wavelength λ=560 nm. Preferably suited for the detector 16 is a silicon photo diode or a silicon photo cell. 
     It is preferable to use a garnet, which can be produced in large crystals according to the Czochralski method, IG for the optical medium of the Faraday crystal 8. The garnet should have great purity and a large Faraday constant. One such garnet has a gadolinium-gallium base. These garnets are transparent throughout the entire visible spectral range. Also highly suited are iron garnets of the type M 3  Fe 4  O 12 , where M is one of the rare earth metals, particularly ytterbium, thulium, yttrium, erbium or holmium. These rare earth metals have a very large Faraday constant below the magnetic saturation, for example, V=4×10 4  min/(T×cm) for the yttrium--iron garnet (Y 3  Fe 5  O 12 ) for a wavelength of λ=1 μm. Therefore, very short Faraday crystals 8 less than 1 mm in length can be used. The coil 7 surrounding the crystals 8 can thus be very small which again decreases the amount of electric energy required for generating the magnetic field. This is particularly important for the implantable glucose sensor, since only very-low-power supplies can be employed. 
     In conjunction with iron garnets, an infrared light-emitting diode of gallium-indium-arsenide-phosphide Ga x  In 1-x  As y  P 1-y  is preferably used as the light source. This LED has an emission wavelength in the range 0.9 μm≦λ≦1.4 μm and, more preferably, in the range of 1.0 μm≦λ≦1.11 μm. In this range, the radiation penetration depth in water is larger than 60 mm, so that cell lengths of this order of magnitude can be employed. Thereby, correspondingly low concentrations of the optically active substance, particularly glucose, can be measured. In conjunction with iron garnets, a germanium photo diode or a germanium photocell can be employed as the detector 16, especially for wavelengths in the range 1.0 μm≦λ≦1.7 μm. 
     The Faraday coil&#39;s 7 control voltage U 21 , having the modulation frequency f F , is a symmetrical squarewave voltage as seen from the diagram of FIG. 2A. With this control voltage U 21 , the largest possible detector signal U 16 , for a given amplitude of the coil current, is obtained. The oscillator frequency f 0  becomes approximately 100×f F  and the modulation frequency f F  is chosen ≲500 Hz. 
     The oscillator frequency f 0 , representing the clock frequency f L  of the light source 2, is preferably chosen to be substantially larger than the modulating frequency f F , as seen from the diagrams of FIGS. 2 and 2A. The modulation frequency f F  is preferably chosen to be in the range 50 Hz≦f F  ≦500 Hz. The clock frequency f L  provides an increased light yield in conjunction with an LED light source 2 by utilizing the nonlinear emission-operating current characteristic of the LED, especially in the visible range of the spectrum. Thereby, an improved signal-to-noise ratio is obtained and thus a higher measuring sensitivity. Using an oscillator frequency f 0  as the clock frequency f L  for the light source 2, the exact same number of radiation pulses from the light source 2 is obtained in each integration period without any additional circuitry cost. The integration period has a duration of 1/(2×f F ). This provides a measuring signal which is a highly stable measuring signal. 
     The starting times t 1  and t 5  for the sample pulse time intervals τ s+   and τ s-  of the sample-and-hold memories 27 and 28 as well as the starting times t 3  and t 7  of the resetting time intervals τ R  for the integrator 26 are determined by counting out f O  and f L  periods during each period f F . Thereby, very accurate and equal integration times for each of the integral values to be stored in the two memories 27 and 28 are obtained. 
     In the method of the present invention, an output signal is obtained having low noise and drift U A  =(U 27  /U 28 )×U N  -U K  =Q-U K . The output signal U A  is zero by a suitable choice of the constant voltage U K  at a zero or background angle of rotation of the polarization plane in the measuring cell 12. The output signal is thus proportional as a measure for the desired concentration of the optically active substance to the angle of rotation Δα M  caused by this substance. Therefore, U A  ˜Δα M . U N  is the constant internal normalizing voltage of the divider 30. 
     The primary measurement signal is the photo current i ph  generated in the detector 16. The principal passage directions of the polarizer 4 and the analyzer 14 are perpendicular to each other, i.e., they are adjusted for minimum transmission. The direction of vibration corresponding to the situation when the polarized radiation strikes the analyzer 14 is characterized by the angle of rotation α=0 for the plane of vibration. The plane of polarization of the radiation leaving the polarizer 4 is alternately rotated the maximum angles ±α F  of the same magnitude, by the a-c current of the frequency f F  flowing through the Faraday coil 7. In the cell 12, the optically active substance causes a further rotation by the angle α M  proportional to the concentration C of the substance. Thereby, the photo current i ph  oscillates with the frequency f F  between a minimum value depending on the degree of polarization and the alternately generated extreme values as shown by Eq. (1). 
     
         i.sub.ph ±=sin.sup.2 (±α.sub.F +α.sub.M)+i.sub.u ×cos.sup.2 (±α.sub.F ±α.sub.M)=i.sub.u +(i.sub.o -i.sub.u)×sin.sup.2 (±α.sub.F +.sub.M).    Eq. (1) 
    
     With an alternating squarewave current in the Faraday coil 7, only the values i ph+  and i ph-  are alternately generated. 
     An ideal case has completely linear polarization of the radiation reaching the analyzer 14. In this case, the &#34;depolarization component&#34; i u  =0. When the Faraday modulator 6 is turned off and the optically active substance in the measuring cell 12 is absent, i.e., with α F  =α M  =0, then the photo current i ph  disappears. 
     For the general case of incomplete polarization 0&lt;i u  &lt;i o , the formation of the total current U 16 , having different amplitudes i ph+  and i ph- , is shown in FIG. 4. The depolarization component i u  is shown exaggerated such that i u  =i o  /4. With commercially available polarizers, i u  can be less than i o  /100. If the ratio: 
     
         ε=iu/io=(1-p)                                      Eq. (2) 
    
     of Eq. (2) is introduced as a measure for the degree of polarization p of the radiation reaching the detector 16, then Eq. (1) for the two photo currents which are generated by the Faraday coil 7, assume the following form: 
     
         i.sub.ph ±=i.sub.o [ε+(1-ε)×sin.sup.2 (±α.sub.F +α.sub.M)].                      Eq. (3) 
    
     The photo current i ph  is amplified in the amplifier 25 to form the voltage U v  =v×i ph  and is subsequently integrated by integrator 26 to form the two voltages ##EQU1## where τ J  is the effective integration time. τ J  is fixed by the oscillator 20, the frequency divider 21 and the counting circuit 22 and is equal to the time difference t 5  -t 4  between the end of a resetting pulse U 22  at time t 4  and the next following storage pulse U 24  at time t 5 . The integrator 26 acts in a manner known per se as a mean value amplifier having adjustable gain. 
     The integral values U I  ± of alternately generated voltage U 26  are stored separately during the time intervals τ s  + and τ S  - defined by the voltage pulses U 23  and U 24 . The changing contents of the memory and the additive constant voltage U K  from the constant voltage source 29, allows the output voltage (U A ) to be finally generated in the divider 30 as follows: ##EQU2## The bars over i ph+  and i ph-  in Eq. (5) indicate that each of the two photo currents are averaged over the corresponding period τ J . U N , as mentioned above, is the internal normalizing voltage of the divider 30, for example, 10 V. If the constant voltage U K  is set equal to the internal normalizing voltage U N , then one obtains: ##EQU3## With the photo current calculated in accordance with Eq. (3), or the voltage U 16  proportional thereto, the output signal normalized to the voltage U N  is as follows: ##EQU4## Similar to the integration of Eq. (4), the variables α M , α F  and ε on the right side of Eq. (7) are interpreted as mean values over an integration period τ J . Thereby, the output voltage U A  and the normalized output signal y are also mean values having a corresponding reduced noise component. For a further reduction of the noise bandwidth, the divider 30 is preferably chosen to have a time constant τ Q  &gt;1/f F . 
     For the case of small rotations in the plane of polarization (this is particularly important for determining glucose concentration) and a small Faraday rotation sin α F  &lt;&lt;1,  for example, α F  =2°, sin α F  =0.035, then Eq. (7) simplifies to: ##EQU5## In addition, if the radiation is linearly polarized to a sufficient degree, such that: ##EQU6## then the fractional term in Eq. (8) can be neglected and the output signal becomes proportional to the angle of rotation α M  and therefore also to the concentration C of the optically active substance to be measured as shown by Eqs. (10) and (11), respectively. 
     
         y=4×α.sub.M /α.sub.F                     Eq. (10) 
    
     
         U.sub.A =4×(α.sub.M /α.sub.F)×U.sub.NEq. (11) 
    
     According to the different possible algebraic representations of Eq. (5), as for example, shown in Eq. (5A) ##EQU7## the output voltage U A  can also be determined from the two output signals U 27  and U 28  in a manner known per se by a subtraction and subsequent division. Accurate adjustment of the null point for compensating small deviations between the photo currents i ph+  and i ph-  or between the output signals U 27   and U 28  only requires the addition of a small offset voltage γU γU A  =U N  -U K  when adjusting the analyzer 14 for minimum transmission. Therefore, the output circuitry of the two sample and hold memories 27 and 28 in FIG. 1 can also be represented by an equivalent circuit diagram as shown in FIG. 5. 
     In the circuit of FIG. 5, the difference voltage U 31  =(U 27  -U 29 )/2 is formed in the subtractor 40. The subtractor 40 can include, for example, an inverter 31 and two resistors 42 and 43. Subsequently, the difference voltage U 31  is divided by the output signal U 28  in a divider 30, taking into consideration the factor 2. The constant-voltage source 44 only needs to supply the small offset voltage γU A , which for an ideal adjustment equals zero. 
     In the embodiment according to FIG. 1, an analog signal processor is provided. The functional elements, shown in the lower part of FIG. 1 however, can partially be replaced, for example, by a microprocessor. The microprocessor would function as the oscillator 20, the frequency divider 21, the counting circuit 22 as well as the integrator 26, the sample-and-hold memories 27 and 28 and the divider 30 having the constant-voltage source 29. 
     Further the embodiment according to FIG. 1 provides a Faraday modulator as the modulator 6. In some cases it may be advantageous to use a Kerr modulator in which an artificial birefringence is generated. Thereby, a rotation of the plane of polarization by the electric field strength also occurs.