|
1.IntroductionThere have been extensive efforts to develop optical techniques to measure blood glucose concentration for diabetes care.1, 2, 3 There is also the application of in situ real-time monitoring of glucose in cell culture, microassay arrays, and bioreactors. In both applications, accuracy, reliability, and speed are important. This work is about a laser-based photometric technique that can address these issues. This work aims primarily for the latter application, but the technique is potentially applicable for blood glucose. A broad review of minimally or noninvasive blood glucose monitoring has been given by Klonoff.4 More recently, Khalil5 and Cote6 specifically reviewed optical spectroscopy, including absorption, polarimetry, Raman, and fluorescence. Absorption spectroscopy in particular, for both the near-IR (NIR) and mid-IR (MIR) spectral region, has been explored extensively by many research groups and commercial concerns over many years. 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24 For the MIR “fingerprint” region, Heise 7 in early work7 analyzed blood glucose in vitro using Fourier transform infrared (FTIR) spectroscopy. Bhandare 9 also used FTIR to investigate the spectral interference of other blood constituents on the glucose signature. In more recent works, 10, 11, 12, 13 various groups focused on using only a few key MIR wavelengths rather than the entire spectrum to determine blood glucose. Specifically, Shen 11 found that spectroscopy at just two wavenumbers, 1082 and , was sufficient. Rhiel 14 employed FTIR to monitor not only glucose but also lactose and ammonia in cell culture in real time. Opposite to absorption, MIR radiation was also used to measure glucose based on Kirchoff’s law of thermal emission.15 Some other works are based on laser photoacoustic effects.25, 26 There are generally three aspects in these studies: the spectroscopy aspect is concerned with spectral features and processing algorithms; the photometry aspect is concerned with optical system engineering and measurements; and the deployment aspect is concerned with application. Photometry is often undervalued and underdeveloped compared to spectroscopy but essential for real-life applications. The MIR works, 7, 8, 9, 10, 11, 12, 13, 14, 15 for example, demonstrated the spectroscopic merit, but the photometry was simply lacking. The thermal-source FTIR instrument employed in these works is useful as a research instrument but not necessarily practical or cost effective. Thermal sources have low power spectral density, low radiance, high relative intensity noise (RIN), and large thermal noise. These impose limitations on dynamic range and speed for accuracy and real-time monitoring. In comparison, the NIR approach does not have any intrinsic spectroscopic advantage, but benefits greatly from more advanced photometric technologies. Recent progress can help exploit the MIR region. This work describes the use of advanced MIR semiconductor lasers for in vitro aqueous glucose measurements. Lasers have well-established photometric advantages including high spectral intensity, brightness, high modulation bandwidth, low RIN, and noise. They are also well suited for applications such as microfluidic or microarray bioassay, as the beam can be focused to a small volume. For practicality, they are robust, compact, and potentially low cost. Advanced MIR semiconductor lasers have been demonstrated for tissue imaging,27 and application for conjunctiva glucose measurement is being developed.24 A potential issue with MIR lasers is the lack of large spectral coverage, but some studies have shown that only a few wavelengths may be needed.11, 13, 24 Also, recent progress in tunable MIR lasers 27, 28, 29, 30, 31 promises broad wavelength coverage. Photometry is crucial to the issue of accuracy. There is often some ambiguity on what contributes to glucose measurement errors: how much it is due to spectroscopic problems (interference, distorted spectra, etc.), and how much is due to photometric errors. Glucose errors were usually determined by comparing with wet chemistry results, but not analytically linked to instrument performance. As a basic principle of system engineering, a system should fully account for all errors, and establish an intrinsic measurement confidence level based on the causal and deterministic relation with the instrument noises and errors. This is a focus of this work. The experiments were performed for glucose in water and (IL) solution. It is chemically much simpler than blood, since the objective is to establish the baseline photometric merit without the issue of chemical complexity. It is hoped that it may eventually motivate blood applications with many-wavelength systems to handle the chemical complexity. The rest of the work is organized as follows. Section 2 describes the experimental approach. Results and discussion are presented in Sec. 3; and Sec. 4 presents a summary and conclusion. 2.Experimental approach2.1.Overall Concept2.1.1.Regression ApproachGiven a glucose solution with concentration , which also contains an interferent substance such as IL of concentration , a system can determine these concentrations via the equation: where is the difference between the absorption coefficient of the solution and that of pure water , and coefficients are quantities to be obtained from calibration. For glucose and IL solution, it is expected that in the MIR owing to Rayleigh’s scattering law and the latter can be neglected. In Eq. 1, it is seen that higher-order terms are negligible; and it is convenient to use the matrix format:where the notation indicates that it is only an inferred estimate of the true value . In more complex solutions with substances, would represent values of concentration.The calibration for is obtained from measured for solutions of known and based on the relation: where superscript “cal” indicates that these are calibrated concentrations, and is the specific absorption coefficient. As it turns out, all higher-order nonlinear terms were found insignificant and the relation becomes:where:where is the number of calibrated samples. With linear regression, the estimate for is:The relation between and is simply: and Eq. 2 becomes:2.1.2.Uncertainty and confidence level analysisThe critical issue and an objective of this work is the confidence level of . Mathematically, it is a function of two independent random variables and . The error of is based on the photometric system: laser, detector, electronic noises, and other measurement apparatus errors. is a random variable not in the sense that it varies from one measurement to the next, but that it has an intrinsic uncertainty associated from imperfect calibration. Ideally, if is very accurately calibrated, the uncertainty of was caused only by photometric measurement error. In this case, error can be characterized by a normal distribution . Then, is also a normal distribution function with an estimated variance: Realistically however, there is no perfect calibration and this is a key issue for optical glucose technologies. Calibration may have to be done frequently, as it can be unique to an individual or apparatus and can change over time. The expected variance of from Eq. 7 is:which also depends on photometric error, though it is not necessarily the same as that in Eq. 9. Photometry thus has a two-fold impact on error, and the quality of calibration is critical.With uncertainty taken into account, the uncertainty of can then be complex, and variance can be a useless concept. For example, consider the simplest case when is a normal-distribution scalar with mean and standard of deviation ; the probability density function (PDF) of of Eq. 8 is then: which has an infinite variance owing to the well-known Cauchy tail at large . Instead of variance, a rigorous measure is the confidence interval, e.g., the interval given by equationwhich has a 99% probability of containing the true value of glucose concentration.If uncertainties are small: , a cut-off can be imposed to keep only the Gaussian component and an approximated variance can be used: where . The uncertainty of is degraded from Eq. 9 by an extra term , where is mean square. The simple and intuitively obvious implication of Eq. 12 is that for any system, there is no need to measure with much less error than that of and vice versa. More importantly, is essentially a spectrum, and arbitrary applications of techniques such as spectral smoothing and other digital filtering, while producing appealing-looking calibration spectra, violates its integrity and may induce spectral distortion that affects the confidence of . A merit and objective of laser-based photometry is to obtain and with high confidence that is causally and rigorously linked to the system performance.2.2.Experimental Setup2.2.1.Optical configuration and experimental conditionThe experiments were set up to measure sample absorption in the transmission mode, using a tunable MIR semiconductor lasers around and a reference laser at . A schematic diagram is shown in Fig. 1 . The sample is a thin liquid film held between two windows by surface tension. One window is fixed, and the other is movable with to a stepping-motor-driven micrometer. The sample thickness was varied 50 to . The laser beam was focused to a spot diameter. The sample volume was , which is in the range for microassay or microfluidic applications. The sample temperature was not measured, but the room temperature was controlled at . Owing to very low laser power (compare Sec. 2.2.3), the instantaneous temperature rise due to laser adiabatic heating is calculated to be , and the steady-state average temperature rise is (based on water-specific heat of and conductivity of ). Thus, the laser-induced temperature rise is insignificant compared to the room temperature fluctuation. There was no effort to study temperature-dependence absorption. 2.2.2.Sample preparationThe anhydrous glucose (Avocado Research Chemicals, Heysham, England) was used to prepare the glucose deionized water (DI) water solution samples. Glucose concentrations ranged from 0 to . The solution of 10% (Baxter, Deerfield, Illinois) was diluted by DI water and mixed with stock glucose solution. The uncertainty of glucose concentration was calculated to be based on the errors of glucose weight and water volume measurement. Intralipid concentration estimated uncertainty was , based on errors of measured liquid volume for mixing. 2.2.3.Laser and detectorThe MIR tunable laser was operated in single mode, covering from 9.5 to (1020 to ) as illustrated in Fig. 2 , with a superposition of glucose and water absorption spectrum. The other laser has a fixed wavelength at , and serves as a reference. The lasers were operated in pulse mode with repetition rate from 25 to ; the pulse width ranged from 200 to . The pulse peak power was from 5 to . The duty cycle was 0.5 to 1%, and the average power was only 25 to . A MCT photodiode (Kolmar Technologies) detector with a transimpedance amplifier (TIA) was used to measure the transmitted power. 2.2.4.Data acquisition, signal processing, and system noiseAn automated data acquisition system was used. The computer generated a sweep of sample thickness and the transmission was measured. The signal was digitized at sampling rate. The overall photometric system noise was determined from three sources: laser noise, detector noise, and electronic noise. The laser RIN power spectral density was measured to be , which appear to be white-noise like without component down to . The combined detector and electronic noise equivalent power spectral density was measured to be . 3.Experimental ResultsThe experimental results include: 1. demonstration of large dynamic range and sample thickness for accurate measurement of absorption coefficients , 2. calibration for specific absorption coefficient and its accuracy, 3. the effect of IL as an interferent substance, and 4. application to the case of fixed and unknown sample optical thickness. 3.1.System Dynamic Range and Measurements of Absorption CoefficientFigure 3 shows a typical transmittance that varies over 4 orders of magnitude as a function of sample thickness, which is essential for the accuracy and precision of obtained as the slope of the curve. The lower curve was for a glucose solution, and the upper curves were for pure water as a reference. The full signal range was actually over 6 orders of magnitude,32 but only results for thickness were shown. This signal dynamic range and sample thickness is much larger than those in FTIR-based works. Speed is also a figure of merit. The time for a typical scan was a few seconds, limited by the computer control and acquisition time. The actual laser on-time was only and its total energy was for the entire scan. It appears in Fig. 3 that can be obtained via Beer’s law: , where is the transmittance and the sample optical thickness. However, Beer’s law is not necessarily correct, especially for laser measurements. The reason is that a focused laser beam contains many transverse components, and the effective path length is longer than the actual. Proper analysis is needed and the exact theory is given in the Appendix in Sec. 5. Deviation from Beer’s law was indeed observed for a tightly focused spot and small sample thickness in this work. Calculation was performed to ensure that the results here were obtained for sufficiently large spot size and sample thickness to correctly infer the absorption coefficient. An example of measurement is the upper curve of Fig. 3, which actually consists of five single-scan results but are indistinguishable on this scale. The slopes yield ranging from 530.17 to , with regression standard errors (SE) varying from 0.12 to , which is to in relative terms, and the regression ranges from 0.999996 to 0.999999. More important is the residual results shown in Fig. 4 . Ideally, they should be consistent with photometric error as discussed in Eq. 9. Large nonrandom or systematic residuals might indicate the presence of uncontrolled variables. The residuals in Fig. 4 appear random with magnitude and behavior consistent with expected photometric system noise, which is plotted as the dashed curve in Fig. 4. It is based on the transmittance signal expressed as: where and are the detector noise and laser RIN, respectively. At small sample thickness, the transmitted signal is strong and the noise is dominated by laser RIN, since and the noise iswhere is the laser RIN power spectral density mentioned in Sec. 2.2.4, and is the bandwidth. At large thickness, the reverse is true: the signal is weak and dominated by detector noise. The dashed curve and data are in excellent agreement, verifying the causal relation between measurement accuracy and the basic photometry performance.However, more extensive measurements with repeated scanning indicated another source of error unrelated to photometry. It was the mechanical variation of the sample thickness control assembly. This mechanical error was intermittent but sufficiently frequent to degrade precision to to . It could be traced to the scan-to-scan variation of the micrometer, independently observed via interferometry. This variation was actually within the manufacturer specification of repeatability. Another potential source of errors is the thermal expansion of the assembly, which was built on a -long steel frame with a thermal expansion of . If repeated scanning heated the frame, or room temperature varied by , this could also contribute to errors, although evidence suggested it was less likely than the former. Observation of this mechanical error actually underscores the importance of establishing a baseline photometric accuracy. This type of micrometer is usually adequate for many applications, but the reason for its error to be detected was because of the high sensitivity and accuracy of photometric measurements. An implication is that the full advantage of laser accuracy can be exploited for cases without moving parts or with high mechanical precision, such as microfiuidic channels or optical waveguide probes. 3.2.Specific Absorption Coefficient of Glucose (without )Calibration for specific glucose absorption coefficient was obtained for solutions without IL to establish the baseline glucose accuracy without interferent. It was derived from the differential absorption as described by Eqs. 3, 4, 5, 6, 7. A challenge is that is a just small difference between two large absorption slopes, as seen in Fig. 3. Figure 5a shows as a function of from 0 to at , 9.63, and . The resulting regression statistics based on Eq. 3 are listed in Table 1 . There are two main results: the validity of the linear model, and the accuracy of calibration. Table 1Regression results of differential absorption coefficients of glucose in water.
The linear model discussed in Sec. 2.1 is clearly valid, as part A of Table 1 shows that the nonlinear terms ( -dependency) for both the 9.63- and results are insignificant, as indicated in column “ ” (which indicates the probability for the null hypothesis that the nonlinear coefficient= ). The linear regression results are shown in part B of Table 1, with and for 9.63 and , respectively. Their uncertainties are given in the SE column, which are and , and in terms of relative unit, 0.76 and 0.62%. Ideally, SE should be determined only by photometric error. Unfortunately, the mechanical error dominated these results. As evidenced in Fig. 5b, the root-mean-squared (RMS) residual error of is for , which is consistent with mechanical errors for ten scan averages. The photometric error contribution should be .Without the mechanical error, the expected relative SE of should have been and . If this calibration is used to infer glucose concentration in unknown samples based on the simplified Eq. 8: , the confidence level can be rigorously established based on Eq. 11, and a calculation result is shown in Fig. 6 . The upper solid curves represent the confidence level for an actual worst-case mechanical error, and the lower dashed curves are for the hypothetical photometric-only error. For both, the calculation was for a single-scan error of . The results show that given the mechanical error, glucose could be determined within with 90% confidence and with 99% confidence. For hypothetical photometry-only errors, the results would be and for 90 and 99% confidence, respectively. It is also clear in part B in Table 1 that the linear coefficient for is within the SE, suggesting that glucose has no significant absorption at this wavelength. Thus, it can be used as a reference for samples of fixed but unknown thickness. This is further discussed in Sec. 3.4. Also evident in Fig. 5b is that the RMS error for is smaller than that of 9.5 to by a factor proportional to at these wavelengths, which are at and at . This behavior indicates that error is not a constant-amplitude additive term, but is consistent with the mechanical error, which is proportional to . 3.3.Effect of InterferenceA matrix of samples with glucose concentration ranging from 0 to and IL concentrations from 0 to 8% were measured at and . The differential absorption coefficients are listed in Table 2 and plotted in Fig. 7 . Each measurement resulted from a ten-scan average. The symbols are experimental data, and the planes are the linear regression fits. It is clear that two substances have nearly orthogonal absorption behaviors at these wavelengths. Glucose shows strong absorption at , while IL has little absorption; the reverse is true at . Interestingly, IL has a negative absorption coefficient at both wavelengths. This can be interpreted that the IL volume that displaces water has a lower absorption coefficient than water. Similar negative effects were reported for blood at 4.6 and in a previous work.32 Table 2Water-substracted absorption coefficients (cm−1) for various glucose and Intralipid® concentrations (absorption of pure water in parentheses).
The analysis according to Eq. 7 was applied to obtain matrix of Eq. 6. The main results include: 1. the validity of the linear model in Eq. 4; 2. calibration matrix ; and 3. the effect of IL on the precision and accuracy of glucose determination. 3.3.1.Probability of nonlinear termsThe regression results for nonlinear terms are shown in Table 3 , part A. For , a -test clearly shows that various nonlinear terms can be neglected. However, the results for indicates a probable quadratic term for with . This quadratic term persisted even when glucose is eliminated from the regression model, as shown in the last two rows of Table 3, part A. Further analysis, however, suggests that this nonlinear term may not be as probable as suggested by the statistics. The regression residuals for data of different glucose concentrations are not random, but systematically form clusters with respect to each value of IL concentration. A possibility for this peculiar result is that each set of samples with the same IL but different glucose concentrations were measured under similar room temperature and mechanical conditions. But sets with different IL concentrations were measured at different times and conditions, and thus had different errors. This systematic error would generate false quadratic behavior only for IL but none for glucose. Also, for both wavelengths, the cross-correlation terms are insignificant, indicating that their absorptions are simply additive to each other. Table 3Regression results for glucose and Intralipid® absorptions.
3.3.2.Linear coefficientsWith the nonlinear terms excluded, the linear regression results are shown in Table 3, part B. For both, the -ratio statistics indicated that the linear model alone was quite sufficient to account for the data. The coefficients and along with the contour of 95% joint-confidence ellipse (based on -ratio distribution) and individual 84% confidential intervals (based on -distribution) are plotted in Fig. 8 . For , the value of is comparable to the value 0.05697 for no-IL obtained in Sec. 3.2. For , is clearly not significant. 3.3.3.Inferred glucose and intralipid concentrationWith calibrated coefficients , the reverse inference of glucose and IL concentration can be obtained via Eq. 8. There was no need to prepare unknown samples given the available extensive number of samples. The results are shown in Fig. 9 . The filled circles represent known values, the open squares with center dots represent inferred values. The residuals of all 24 concentrations are shown in the scatter plot of Fig. 10 . It appears from Figs. 9 and 10 that the residuals of IL concentration do not seem to be random around the expected values. All inferred IL concentrations in Fig. 9 are either higher or lower than expected values, forming two distinctive clusters in Fig. 10. This reflects the systematic error discussed before for the IL data. More important are the glucose results. The RMS of all residual glucose errors shown in Fig. 10 is , indicated by the two thick dashed lines. Figure 11 plots the glucose RMS error for each IL. Here, the result in Sec. 3.2 is just a subset with . These RMS errors are consistent with the mechanical errors, and there is no obvious correlation between glucose RMS and , indicating that IL has no effect on glucose error. The result is hardly surprising, considering that the two substances have orthogonal absorption behaviors at the two wavelengths. Ellipsoid confidence contours can be calculated for inferred glucose and IL concentrations; however, they are not relevant with regard to IL, owing to the IL systematic errors, and IL is also not of interest. Therefore, the bivariate PDF can be integrated over to yield PDF for only, and numerical approximation can be used to obtain the confidence interval. The result is shown in Fig. 12 . It is remarkably similar to the result without IL in Fig. 6, confirming that IL does not interfere with glucose accuracy. Comparison of Figs. 12 and 10 shows that all 24 glucose data points fall within the calculated 99% confidence level, and 23 out of 24 (96%) fall within the calculated 90% confidence level. 3.4.Approach for Fixed and Unknown ThicknessIn applications such as monitoring a microfluidic channel, microassay array, or diffuse scattering from tissues, the sample effective optical thickness is usually fixed or unknown. The system must infer the concentrations without varying the thickness. For these cases, the algorithm is slightly modified from that in Sec. 2.1 by treating as an additional unknown variable. The details are described in Appendix B in Sec. 6. The essential result is similar to Eq. 8: where is modified from by including , is modified by including water absorption, and a key difference is to use instead of , where is the optical signal that can be transmission or backscattering.It is sufficient to consider an example to illustrate the principle, and this is summarized in Table 4 . The first two rows show the numerical values of , and , which are obtained from the transmittance data of an aqueous glucose sample set at a fixed assembly thickness. The last two rows show inferred value and , to be compared with known values and , respectively. There is good agreement for , but is a bit off, which is likely due to some unaccounted optical loss (Appendix B in Sec. 6.). A crucial aspect of this result is the choice of wavelengths. A requirement of Eq. 14 is that is not badly conditioned, which in principle is no different from Eq. 8. In practice, however, it is a more critical requirement, since contains , which provides spectral contrast to infer . Since is large, the ’s must be chosen such that has variation much larger than their uncertainty. In the second row of Table 4, the small-font numbers for indicate its uncertainty. Their small uncertainty values ensure that is meaningful, owing to the large water absorption difference between 5.4 and . In comparison, any wavelength set in the laser tuning range 9.5 to was not as useful, since water absorption varies too slowly over this range. Both this and the IL results suggest that large wavelength diversity can be advantageous for accuracy and confidence, as discussed in Appendix B Sec. 6. 4.Summary and ConclusionThis work describes the use of mid-IR semiconductor lasers to measure aqueous glucose concentration and examines the issue of accuracy. Compared with many previous works 7, 8, 9, 10, 11, 12, 13, 14 that focused mainly on spectroscopy, this work aims at photometric engineering to demonstrate the merit of the laser-based technique as opposed to others such as FTIR. Specifically, the work focuses on the laser properties of high signal-to-noise ratio, large dynamic range, and speed. The laser capability to focus onto a small sample volume is also suitable for microfluidic and microarray bioassay applications. The approach measures the absorption coefficients of aqueous glucose and solutions of various concentrations for different wavelengths (5.4, 9.54, and ). Linear regression is applied to obtain the specific absorption coefficients for each solute, which are then applied to infer their concentrations and compare with known values. The important issue is accuracy, and the errors are shown to be causally and deterministically related to the basic instrument noises and system errors. There are no unaccounted sources of error. This causal relation allows inferring measurement confidence level based on rigorous system understanding instead of phenomenological assessment. The results show that glucose concentration could be determined and with 90 and 99% confidence, respectively, in a single measurement. These values improve to and with ten measurements. These uncertainties were not limited by photometry, but by a larger mechanical error of the sample thickness control assembly. This mechanical error is not fundamental, and in applications without this problem, the expected uncertainty due to photometry-only error would be and for 90 and 99% confidence level, respectively, for a single measurement. up to 8% concentration had no effect on the glucose measurement. The choice of wavelengths is also evidently important. The large difference between 5.4 and offers large spectral contrast between glucose, IL, and water that contributes to the measurement confidence. In general, for the laser method, fewer ’s entails a simpler system, whereas FTIR offers the entire spectrum. However, an insight from these results is that with regard to accuracy, measurements with many wavelengths but noisy is not necessarily better than those with fewer wavelengths but low noise. Other recent MIR works 10, 11, 12, 13, 24 even with FTIR also suggest that a few wavelengths might be sufficient. The laser approach can thus be applicable and practical. 5.Appendix A: Deviation from Beer’s Law for Focused Laser BeamIn transmission or diffuse backscattering, the laser light coherence property may require the full optical solution for accurate interpretation of the result. For an absorptive thin film with parallel interfaces, i.e., a lossy etalon, the transmittance is: where TE and TM denote the transverse electric and magnetic polarization, respectively, and is the sample thickness. For all other quantities, the subscripts 1 and 2 correspond to the window and sample, respectively. is the complex dielectric constant; with , and is the transverse wave number of the beam parallel to the interface.For small , the paraxial approximation can be used: where is the beam propagation angle in the medium. Substituting the imaginary component in of Eq. 15 results in the absorption term:which indicates that the beam experiences a larger loss factor associated with a longer path length. A focused laser beam is a superposition of many plane waves with nonvanishing transverse , therefore it experiences a higher loss than a normal-incident plane wave.In this work, the index of refraction of the window is very close to that of water, and there is little difference between TE and TM transmissions. Thus the net beam transmittance is: where , are the power spectral density functions of transverse spatial frequencies , . For an - Gaussian beam:With paraxial approximation, applying Eq. 19 with Eq. 20 yields:From Eq. 21, Beer’s law would be obtained for large beam waists: , . But for small beam waists, Eq. 22 shows a deviation by the term , which would lead to substantial error of if not included.In this work, the semiconductor laser beam is not Gaussian and has - asymmetry. A numerical calculation for a modeled laser beam produces results shown in Figs. 13, 14, 15 . Figure 13 shows that the transmittance is superficially similar to Beer’s law on a large scale. But as evidenced in Fig. 14, the residuals from the exponential fit clearly deviate from Beer’s law, even if the interference fringes are ignored. The deviation is worse for smaller spot size as expected. The “effective” absorption coefficient deviates from the actual coefficient as shown in Fig. 15. Considering that the largest differential glucose absorption is only in this work, a deviation of would result in an error of 14%. Table 4Example of inferring unknown optical thickness.
This error does not just occur on the slope of the log transmittance, but also on its amplitude. Thus, for cases with fixed-thickness samples, the error on transmittance amplitude is similarly as severe and correction is necessary. 6.Appendix B: Approach for Fixed or Unknown Sample Effective Optical ThicknessIn general, the transmission or diffuse backscattered return signal can be expressed as: where is a loss factor that may or may not be known, and is a sum of absorptions:where , is the baseline absorption coefficient of the medium (which can be a solvent other than water), can represent single or round-trip, and for to denotes concentration of substance of interest.In controlled applications with fixed, known and , Eq. 24 can be written: and a solution similar to Eq. 8 is obtained:If is unknown, Eq. 24 can be written:If is negligible, Eq. 27 leads to Eq. 14. If is unknown, the problem cannot be solved. However, usually represents optics effects such as light collection loss, Fresnel loss, scattering loss, or chromatic dispersion, and therefore can be a slowly varying function of wavelength. For the sake of discussion, suppose that it can be parametrized to the first order:then Eq. 27 can be modified to yield:which allows inferring all three unknown parameters , , and .A crucial requirement relevant to the wavelengths used in this work is that should have low uncertainty. Let each row or column of be denoted as a vector . If the wavelength range is too narrow, various vector can become close to each other, can be nearly degenerate with large uncertainty, and the measurements would have poor confidence. A nominal criterion for the vectors to be “far” from each other is to use the quasi-Mahalanobis distance: where are the vector elements, are their uncertainties, and is a proportional factor.Applying this criterion in this case explains the reason why the 5.4- and wavelength combination produced the best results in Sec. 3.4. The quasi-Mahalanobis distance for columns with 5.4 and is 93, compared to only to 6 for any set of wavelengths in the 9.5 to range. AcknowledgmentThis work was supported in part by the DARPA/ARO under contract DAAD-190010361, the Air Force Research Laboratory under contract F29601-00-2-0058, and the State of Texas through the THECB and the Texas Center for Advanced Materials and Center for Superconductivity. We thank the late Pradeep Johnson, C. Peng, G. P. Luo, and A. Seetheraman for valuable assistance. ReferencesC. M. Peterson, Diabetes Management in the ’80s: The Role of Home Blood Glucose Monitoring and New Insulin Delivery Systems, Praeger, New York (1982). Google Scholar
M. Rohrscheib,
R. Robinson, and
R. Philip Eaton,
“Non-invasive glucose sensors and improved informatics—the future of diabetes management,”
Diabetes, Obesity Metabolism, 5
(5), 280
–281
(2003). Google Scholar
D. C. Klonoff,
“Noninvasive blood glucose monitoring,”
Diabetes Care, 20
(3), 433
–437
(1997). https://doi.org/10.2337/diacare.20.3.433 0149-5992 Google Scholar
O. S. Khalil,
“Spectroscopic and clinical aspects of noninvasive glucose measurements,”
Clin. Chem., 45 165
–177
(1999). 0009-9147 Google Scholar
G. Cote,
“Non-invasive and minimally invasive optical monitoring technologies,”
J. Nutr., 131 1596S
–1604S
(2001). 0022-3166 Google Scholar
H. M. Heise,
R. Marbach,
G. Janatsch, and
J. D. Kruse-Jarres,
“Multivariate determination of glucose in whole blood by attenuated total reflection infrared spectroscopy,”
Anal. Chem., 61 2009
–2015
(1989). https://doi.org/10.1021/ac00193a004 0003-2700 Google Scholar
D. M. Haaland,
“Multivariate calibration methods applied to quantitative FT-IR analysis,”
Practical Fourier Transform Infrared Spectroscopy, 395
–486 Eds.Academic Press, New York (1990). Google Scholar
P. Bhandare,
Y. Mendelson,
E. Stohr, and
R. A. Peura,
“Glucose determination in simulated plasma blood serum solutions by FTIR spectroscopy: investigation of spectral interferences,”
Vib. Spectrosc., 6 363
–378
(1994). https://doi.org/10.1016/0924-2031(93)E0059-B 0924-2031 Google Scholar
R. Vonach,
J. Buschmann,
R. Falkowski,
R. Schindler,
B. Lendl, and
R. Keller,
“Application of mid-infrared transmission spectrometry to the direct determination of glucose in whole blood,”
Appl. Spectrosc., 52 820
–822
(1998). https://doi.org/10.1366/0003702981944553 0003-7028 Google Scholar
Y. C. Shen,
A. G. Davies,
E. H. Linfleld,
T. S. Elsey,
P. F. Taday, and
D. D. Arnone,
“The use of Fourier-transform infrared spectroscopy for the quantitative determination of glucose concentration in whole blood,”
Phys. Med. Biol., 48
(13), 2023
–2032
(2003). https://doi.org/10.1088/0031-9155/48/13/313 0031-9155 Google Scholar
Y. J. Kim,
S. Hahn, and
G. Yoon,
“Determination of glucose in whole blood samples by mid-infrared spectroscopy,”
Appl. Opt., 42
(4), 745
–749
(2003). https://doi.org/10.1038/nj6959-745a 0003-6935 Google Scholar
W. B. Martin,
S. Mirov, and
R. Venugopalan,
“Using two discrete frequencies within the middle infrared to quantitatively determine glucose in serum,”
J. Biomed. Opt., 7
(4), 613
–617
(2002). https://doi.org/10.1117/1.1501893 1083-3668 Google Scholar
M. Rhiel,
P. Ducommun,
I. Bolzonella,
I. Marison, and
U. von Stockar,
“Real-time in situ monitoring of freely suspended and immobilized cell cultures based on mid-infrared spectroscopic measurements,”
Biotechnol. Bioeng., 77
(2), 174
–185
(2002). https://doi.org/10.1002/bit.10134 0006-3592 Google Scholar
C. D. Malchoff,
K. Shoukri,
J. I. Landau, and
J. M. Buchert,
“A novel noninvasive blood glucose monitor,”
Diabetes Care, 25 2268
–2275
(2002). https://doi.org/10.2337/diacare.25.12.2268 0149-5992 Google Scholar
M. A. Arnold,
“Non-invasive glucose monitoring,”
Curr. Opin. Biotechnol., 7
(1), 46
–49
(1996). https://doi.org/10.1016/S0958-1669(96)80093-0 0958-1669 Google Scholar
M. A. Arnold and
C. W. Small,
“Noninvasive glucose sensing,”
Anal. Chem., 77
(17), 5429
–5439
(2005). https://doi.org/10.1021/ac050429e 0003-2700 Google Scholar
J. W. Hall and
A. Pollard,
“Near-infrared spectrophotometry: a new dimension in clinical chemistry,”
Clin. Chem., 38
(9), 1623
–1631
(1992). 0009-9147 Google Scholar
M. R. Robinson,
R. P. Eaton,
D. M. Haaland,
G. W. Koepp,
E. V. Thomas,
B. R. Stallard, and
P. L. Robinson,
“Noninvasive glucose monitoring in diabetic patients: a preliminary evaluation,”
Clin. Chem., 38
(9), 1618
–1622
(1992). 0009-9147 Google Scholar
H. M. Heise,
R. Marbach,
T. Koschinsky, and
F. A. Gries,
“Noninvasive blood glucose sensors based on near-infrared spectroscopy,”
Artif. Organs, 18
(6), 439
–447
(1994). 0160-564X Google Scholar
R. Marbach,
T. H. Koschinsky,
F. A. Gries, and
H. M. Heise,
“Noninvasive blood glucose assay by near-infrared diffuse reflectance spectroscopy of the human inner lip,”
Appl. Spectrosc., 47 875
–881
(1993). https://doi.org/10.1366/0003702934415354 0003-7028 Google Scholar
V. Saptari and
K. Youcef-Toumi,
“Design of a mechanical-tunable filter spectrometer for noninvasive glucose measurement,”
Appl. Opt., 43
(13), 2080
–2088
(2004). https://doi.org/10.1364/AO.43.002680 0003-6935 Google Scholar
M. J. McShane and
G. L. Cote,
“Near-infrared spectroscopy for determination of glucose, lactate, and ammonia in cell culture media,”
Appl. Spectrosc., 52 1073
–1078
(1998). https://doi.org/10.1366/0003702981944968 0003-7028 Google Scholar
J. Burd, Google Scholar
K. M. Quan,
G. B. Christison,
H. A. MacKenzie, and
P. Hodgson,
“Glucose determination by a pulsed photoacoustic technique: an experimental study using a gelatin-based tissue phantom,”
Phys. Med. Biol., 38 1911
–1922
(1993). https://doi.org/10.1088/0031-9155/38/12/014 0031-9155 Google Scholar
Y. Shen,
Z. Lu,
S. Spiers,
H. A. MacKenzie,
H. S. Ashton,
J. Hannigan,
S. S. Freeborn, and
J. Lindberg,
“Measurement of the optical absorption coefficient of a liquid by use of a time-resolved photoacoustic technique,”
Appl. Opt., 39 4007
–4012
(2000). 0003-6935 Google Scholar
B. Guo,
Y. Wang,
C. Peng,
G. Luo, and
H. Q. Le,
“Multi-wavelength mid-infrared micro-spectral imaging using semiconductor lasers,”
Appl. Spectrosc., 57 811
–822
(2003). https://doi.org/10.1366/000370203322102906 0003-7028 Google Scholar
G. P. Luo,
C. Peng,
H. Q. Le,
S. S. Pei,
W. Y. Hwang,
B. Ishaug,
J. Um,
J. N. Baillargeon, and
C. H. Lin,
“Grating-tuned external-cavity quantum cascade semiconductor lasers,”
Appl. Phys. Lett., 78 2834
–2836
(2001). https://doi.org/10.1063/1.1371524 0003-6951 Google Scholar
G. Luo,
C. Peng,
H. Q. Le,
S. S. Pei,
H. Lee,
W. Y. Hwang, and
B. J. Zheng,
“Broadly wavelength-tunable external cavity mid-infrared quantum cascade lasers,”
IEEE J. Quantum Electron., 38 931
(2002). 0018-9197 Google Scholar
C. Peng,
G. P. Luo, and
H. Q. Le,
“Broadband, continuous, and fine-tuning properties of external-cavity TE-stabilized mid-IR quantum cascade laser,”
Appl. Opt., 42 4877
–4882
(2003). 0003-6935 Google Scholar
C. Peng,
H. L. Zhang, and
H. Q. Le,
“Mid-infrared external-cavity two-segment quantum cascade laser,”
Appl. Phys. Lett., 83 4098
–4100
(2003). https://doi.org/10.1063/1.1627486 0003-6951 Google Scholar
B. Guo,
Y. Wang,
C. Peng,
H. L. Zhang,
G. P. Luo,
H. Q. Le,
C. Gmachl,
D. L. Sivco,
M. L. Peabody, and
A. Y. Cho,
“Laser-based mid-infrared reflectance imaging of biological tissues,”
Opt. Express, 12 208
(2004). https://doi.org/10.1364/OPEX.12.000208 1094-4087 Google Scholar
|