Expand this Topic clickable element to expand a topic
Skip to content
Optica Publishing Group

Rapid, accurate, and precise concentration measurements of a methanol–water mixture using Raman spectroscopy

Open Access Open Access

Abstract

Here we design, construct, and characterize a compact Raman-spectroscopy-based sensor that measures the concentration of a methanol–water mixture with 0.5% accuracy. The sensor measures the concentration with a precision of 0.2% with a 1 second measuring time, and with longer measurement times, the precision reaches as low as 0.006%. We characterize the long-term stability of the instrument over an 11-day period of constant measurement, and confirm that systematic drifts are on the level of 0.02%. We describe methods to improve the sensor performance, providing a path towards accurate, precise, and reliable concentration measurements in harsh environments. This sensor should be adaptable to other water–alcohol mixtures, or other small-molecule liquid mixtures.

© 2018 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

1. Introduction

Raman spectroscopy is a powerful analytical technique that provides detailed information on molecular vibrations. While Raman scattering is typically a weaker effect than directly probing vibrational transitions in the infrared (IR) spectral region (~3000 to ~20,000 nm), Raman spectroscopy can be accomplished using light in the visible or near-IR region. The ability to use shorter wavelengths can be a significant advantage, because laser sources, spectrographs, and detectors that operate in the visible/NIR region are typically less expensive and more compact than comparable devices in the mid- or long-wave-IR region. Being a convenient and field-deployable technique, Raman spectroscopy has found diverse applications, such as imaging biological samples [1], detecting trace levels of explosives at a long distance [2], and identifying pharmaceutical compounds inside an intact capsule [3].

Since a Raman measurement only requires optical access to the sample, and can be completed using relatively inexpensive lasers and spectrometers, it is ideally suited for probing the precise concentrations of mixtures in a laboratory, field, or industrial setting. While Raman spectroscopy has indeed been demonstrated for sensing the concentrations of solutions [4–8], liquid mixtures [9–12], and gas mixtures [13], the accuracy is typically observed to be on the level of a few percent, the precision is often not quantified, and long-term systematic uncertainties are not investigated. Moreover, demonstrations are typically completed in a laboratory setting and with bulky, expensive equipment. Commercial Raman probes are often expensive, and their ability to measure concentrations is usually not described quantitatively. Consequently, it is not clear if a simple, inexpensive Raman spectrometer can be used to measure concentrations with high levels of accuracy, precision, repeatability, reliability, and immunity from environmental perturbations.

Here we investigate the use of Raman spectroscopy to measure the concentration of a methanol–water mixture. Being a mixture of two widely available small molecules, a methanol–water mixture serves as a straightforward benchmark of concentration-sensing ability. In addition, methanol–water mixtures are of considerable interest for several important applications. For example, Rebound Technologies is currently operating refrigeration systems that are driven by the mixing and separation of a freeze-point suppressant like methanol and water [14]. This new style of refrigeration technology offers both efficiency and operational benefits but requires accurate and inexpensive ways of sensing concentration in harsh environments. Additionally, the measurement of methanol concentration is important for the operation of direct methanol fuel cells [15].

Here we design, construct, and thoroughly characterize a low-cost, low-power-consumption, Raman-spectroscopy-based sensor to precisely measure the concentration of a methanol–water mixture. Our Raman-spectroscopy-based concentration sensor (RCS) utilizes a excitation wavelength of 532 nm, which allows the use of a low-optical-power (4 mW) laser and an efficient silicon-detector-based spectrometer (Fig. 1). By fitting the Raman spectrum of a methanol–water mixture as a linear combination of the spectra of the two constituents, we extract a concentration estimate that is largely immune to variations in laser power, optical alignment, and fluctuations in the absorption of the sample. This simple approach provides an estimate of the methanol concentration to within 3% accuracy, but with high repeatability. Moreover, by measuring several mixtures with known concentration, we prepare an empirical correction curve, which can improve the accuracy of the sensor to 0.5%. We characterize the precision of the sensor on timescales ranging from 0.02 to 600,000 seconds. We find that the sensor provides a precision of 0.2% with a 1-second averaging time, and 0.006% with a 1000-second averaging time.

 figure: Fig. 1

Fig. 1 The experimental scheme for the Raman-spectroscopy-based concentration sensor (RCS). A 532-nm laser is focused into the sample. Some photons experience Stokes Raman scattering, and are emitted in all directions. The back-scattered Raman-shifted photons (yellow) are collected by the same lens, but transmit through the longpass dichroic mirror. After passing through a long-pass filter to further reject elastically scattered 532-nm light from the sample, the Raman-shifted light is couped into a fiber which transports it into the grating-based spectrometer.

Download Full Size | PDF

2. Experiment

2.1. Spectrometer design

The RCS (Fig. 1) is constructed with all commercially available components, with an emphasis on a compact, and durable design. Our apparatus roughly follows the design of Mohr et al. [16], and the simple backreflection geometry ensures that the strong un-scattered excitation beam is directed away from the detector, relaxing requirements on the long-pass filters. This geometry relies on a dichroic mirror (Thorlabs DMLP567R) to send the laser light to the sample and direct the red-shifted light (generated through spontaneous Stokes Raman scattering) to the spectrometer. The laser light is focused into the sample using a 5-mm-focal length lens, which also collects the scattered light. The sample is contained in a cuvette tube (Thorlabs CV10Q3500) made from 1-mm thickness quartz.

The excitation laser (Thorlabs CPS532) is a compact diode-pumped–solid-state (DPSS) neodymium-doped yttrium aluminium garnet (Nd:YAG), with an internal frequency doubling crystal to provide 4.5 mW (manufacturer’s specification) of continuous-wave 532-nm light. We note that this is not a “single frequency laser”, and changes in the laser frequency and power as a function of time and temperature are expected.

The energy range of Raman shifts that can be measured by our spectrometer includes the important C–H stretch (∼3000 cm−1), O–H stretch (∼3400 cm−1), and C–O stretch (~1000 cm−1) regions. The range is limited on the low-energy side by our choice of long-pass filters, and on the high energy side by the range of the spectrometer. A 550-nm long-pass filter (Thorlabs FEL0550) is used to reject elastically (Rayleigh) scattered 532-nm light from reaching the spectrometer. This longpass filter sets a hard cut-off to the observable Raman shift at 615 cm−1. The dichroic mirror has a nominal cut-off wavelength of 567 nm, but, in fact, provides a slow decrease in transmission from 580 to 550 nm, meaning that our Raman spectrometer should exhibit a decreasing response for Raman shifts lower than 1550 cm−1.

The spectrometer (Ocean Optics USB4000) is a compact grating-based spectrometer with a silicon array-detector and measures the spectrum at 3647 pixels from 340.2 to 1036.3 nm. This allows us to measure Stokes Raman shifts up to 9150 cm−1 and (if a notch filter were used rather than the long-pass filter) anti-Stokes Raman light up to 10,600 cm−1, a range that is more than sufficient for most applications. The USB4000 is a general purpose spectrometer that does not utilize a thermo-electric cooled sensor; better performance may be obtainable by using a cooled detector. At 532 nm, we observe that the spectrometer provides a resolution of approximately 2 nm, which corresponds to 70 cm−1. For samples with many narrow peaks, a higher resolution spectrometer may allow for better separation of the Raman peaks.

The entire apparatus, including the laser, optics assembly, and spectrometer is housed in a steel box, which reduces sensitivity to ambient light, air currents, vibration, and dust. We note that, in industrial applications, the cuvette could be replaced with a glass window to provide access to a vessel or pipe containing the liquid of interest. Beam steering mirrors are held with commercial mirror mounts (Thorlabs KM100) and attached to an optical breadboard (Thorlabs MB612F).

Example spectra for water, methanol, and a mixture of water and methanol are shown in Fig. 2. For water, a broad peak is seen in the region of 3400 cm−1, which corresponds to the O–H stretching vibration. A very broad peak can also be seen across the 1000 to 2800 cm−1 region, which corresponds to a bending mode. For methanol, narrower, overlapping peaks resulting from C–H stretch are seen near 3000 cm−1, and additional peaks for C–O stretch and a CH3 bending mode are seen 1400 and 1050 cm−1 respectively. For comparison, the exact frequencies for the Raman transitions listed in the NIST WebBook [17] are shown in Fig. 2(a).

 figure: Fig. 2

Fig. 2 (a) Raman spectra for pure methanol and pure water samples. These spectra are an average of 100 spectra, each collected for 2 seconds. The various peaks in the methanol spectrum can be assigned to various vibrational transitions listed in the NIST WebBook [17]. (b) The concentration of the mixture of can be estimated by fitting the observed spectrum (green) with a linear combination (magenta) of water (blue) and methanol (red) reference spectra.

Download Full Size | PDF

2.2. Data processing

Minimal data processing was required to extract the methanol concentration from the Raman spectrum, and the complete procedure is described here. First, a 9-pixel median filter was applied to remove effects from “hot pixels” in the spectrometer CCD and to smooth the acquired spectra. The spectra were clipped to the 600 to 670-nm region, which contains the C–H stretch peak from methanol and the O–H stretch peak from water. We note that we achieved roughly comparable results when using data from the 550 to 670-nm region, but found that expanding the region could increase susceptibility to broadband photoluminescent background that was present in some methanol samples, presumably from impurities.

We fit the spectrum of the methanol-water mixture fmix(λ) as a linear combination of a pure water reference spectrum fw(λ) and a pure methanol reference spectrum fme(λ), according to

fmix(λ)=afw(λ)+bfme(λ)+cλ+d,
where a , b , c , and d are coefficients determined by the fitting procedure and λ is the wavelength. The offset (d ) and the arbitrary linear term () are not strictly necessary, but allow the fitting procedure to be resilient to changes in ambient light levels. The nonlinear fit was completed using the Levenberg–Marquardt algorithm [18] as implemented in MINPACK [19] and accessed through the scipy.optimize.curve_fit function in the Scipy Python package [20]. The “initial guess” values were a0 = 0.5, b0 = 0.5, c0 = 0, and d0 = 0. The methanol mass fraction was estimated as
wmemmemw+mmebρmeaρw+bρme,
where ρw and ρme are the densities of water and methanol. The densities must be included to convert from the Raman measurement, which measures the volume fraction, to mass-fraction, so that we can compare our measurements to the mass-fractions obtained using a scale. In principle, one could simply estimate the methanol mass fraction as b ρme, however, by dividing by the total signal in Eq. (2), we can make the concentration estimate immune to fluctuations in the total spectral intensity, which can be caused by changes in the laser intensity, alignment, or spectrometer response. Figure 2(b) shows an example spectrum of a methanol-water mixture fit as a linear combination of the water and methanol spectra. As seen in Fig. 2(b), the coefficients applied to the water and methanol spectra often do not sum to one, likely due to fluctuations in the spectral intensity resulting from fluctuations in the laser power. Taking into account the density of methanol, we can extract a weight percent of 44.4% methanol. Using the empirical correction factor discussed below, we obtain a concentration of 47.3%, which is the same value as obtained by weighing the water and methanol using a scale. While we achieve good success with this procedure for our two-component mixtures, we note that more sophisticated algorithms are available to treat multi-component mixtures [21].

2.3. Sample preparation

Samples of known mass-fraction of methanol and water were prepared using a scale, which had a precision of approximately 1 mg. The total sample volume that was prepared was approximately 1 g. We prepared the samples by adding a volume of methanol, measuring the weight, and then adding a volume of water and measuring the total weight. Consequently, in a situation where we measured 500 mg of methanol and 1000 mg of total weight, the actual mass fraction of methanol might have been 499/1001 ≈ 49.85%. Thus, we consider the error from the scale to be ±0.15%.

The samples were weighed into standard 1.5 mL test tubes and shaken vigorously before being poured into a quartz cuvette tube for analysis. We note that early trials, where the water and methanol were combined directly in the cuvette tube, displayed anomalous results due to insufficient mixing of the water and methanol. For accurate results, thorough mixing is critical.

3. Results and discussion

3.1. Accuracy

To assess the accuracy of the RCS, we prepared samples of known weight-percent methanol using a scale, as described above. We then recorded reference spectra of water, methanol, and mixtures of known concentration using a 2-second integration time and averaging 10 spectra. As shown in Fig. 3(a), the weight-percent methanol recovered from the Raman sensor is true to the real concentration to within a few percent. The agreement is best when nearly pure methanol or nearly pure water is used. But, when the concentration is close to 50%, there is a systematic offset of about 3% [Fig. 3(b)]. This accuracy level is similar to the “detection limit” for the concentration of a water-toluene mixture reported by Hashimoto et al. [22] using a sophisticated femtosecond laser system.

 figure: Fig. 3

Fig. 3 The accuracy of the Raman concentration sensor. (a) The extracted mass fraction of a methanol–water mixture extracted from the Raman spectrum versus the mass fraction measured using a scale. Each point corresponds to a separate measurement, consisting of an average of 10 spectra, each collected using a 2-second exposure time. The various colors indicate independent experiments, conducted on over three days, demonstrating the repeatability of the measurements. A systematic underestimate of the methanol concentration is observed, which becomes most pronounced for methanol concentrations near 50%. (b) The difference between the points in panel (a) and the true mass fraction [dashed line in panel (a)] highlight a systematic underestimate of the methanol concentration. The error-bars indicate the 95-percent confidence interval extracted from the fitting procedure. The empirical correction (blue line) is a simple cosine function that has been fit to the data to model the systematic error. (c) Using this empirical correction function, we can correct the observed methanol concentrations and achieve improved agreement between the observed and true methanol concentrations. We note that the error from sample preparation is estimated to be ±0.15%, which is smaller than the width of the data points.

Download Full Size | PDF

It is likely that our observed deviation from the true methanol mass fraction is a result of real changes in the mixture. Indeed, it has been shown that methanol–water [10, 12, 23, 24] and ethanol– water [24–26] mixtures form different alcohol–water complexes at different concentrations, and that these complexes display fundamentally different Raman spectra from that of the pure compound. Thus, fundamentally, the spectrum of a methanol–water mixture cannot be perfectly modeled as a linear combination of the spectra of the two components. Thus, any method that makes this assumption, including methods that use the peak heights, peak areas, or complete spectral shape (as done here) will experience errors depending on the nature of the complexes that form in the mixture. We note that the O–H stretching peaks are likely the most sensitive to the formation of different complexes in the mixture, and in many cases this issue of spectral distortion may be minimized by excluding the O–H stretch region from the fitting procedure. However, the only sharp Raman peak in water is the O–H stretch peak, and so this issue cannot be avoided in our case.

Nevertheless, while it may be difficult to extract an absolute concentration measurement from a Raman spectrum from first-principles, it’s possible to empirically correct for the observed deviation, as long as there is a one-to-one relationship between the extracted concentration and the true concentration. In other words, we can prepare mixtures of a known concentration and, as long as the observed concentration monotonically increases with the true concentration, then we can calculate the true concentration from an empirical dataset. For the methanol–water mixtures that we prepare, we observe that our method of fitting the spectra of the mixture as a linear combination does indeed provide a one-to-one mapping to the true methanol concentration. Thus, we can apply an empirical correction based on our recorded data. To obtain the correction function, we fit a cosine function through the data points in Fig. 3(b) and obtain a satisfactory fit. Using this as an empirical correction, we show that the RCS can provide an estimate of the true methanol concentration with a standard deviation of 0.52% [Fig. 3(c)].

This behavior was repeatable, and each color in Fig. 3 represents an independent trial, where new pure-water and pure-methanol reference spectra were collected and fresh mixture samples were prepared. These trials were collected on three separate days, separated by several weeks. Thus, after empirical calibration, the RCS provides repeatable estimates of the methanol concentration with high accuracy. We speculate that the accuracy obtainable with a RCS may be even better than we report, and we suspect that our estimated accuracy may be limited by our ability to repeatably measure and mix ∼1 mL samples. The use of better mixing techniques, a higher precision balance, or larger sample preparation volumes may allow even better accuracy to be realized. As shown in the next sections, the precision and stability of the RCS are sufficient to support measurements better than 0.1%, and, consequently, the empirically corrected concentration measurements should, in principle, translate this level of precision into accuracy.

3.2. Precision

Like many measurement instruments, the precision of the RCS depends on the time available for the measurement. If a longer time is allowable for the measurement, better precision can, in principle, be obtained. The only parameter that can be adjusted in our measurements is the integration time of the spectrometer, which can be equal to or smaller than the total measurement time. For example, if an application requires a concentration measurement every 1 second, then we can envision two ways of performing this measurement. First (Averaging Method A) we could simply set the exposure time of the spectrometer to 1 second and extract the concentration from the acquired spectrum. Second (Averaging Method B), we could set the exposure time to 0.1 seconds, collect 10 spectra, extract the concentration from each spectrum, and average the concentrations. We would expect Method A to provide an advantage in a situation where sensor noise was limiting the precision of the measurement, since it would allow for the highest signal-to-noise ratio for an individual spectrum. Method B might excel in a situation where other systematic errors dominate.

Figure 4 displays the spectra recorded with various exposure times and demonstrates that with very low integration times, the Raman features are obscured by detector noise. Consequently, we expect that a longer measurement time will have a double benefit – not only will more Raman-scattered photons be collected, but a longer integration time can be used, improving the signal-to-noise ratio of each collected spectrum. Indeed, the extracted concentration versus time is quite noisy when the 20 ms integration time is used and becomes much flatter when a 2000 ms integration time is used [Fig. 5(a)–5(d)].

 figure: Fig. 4

Fig. 4 The Raman spectra (green dots) show improved signal-to-noise ratio (SNR) as the exposure time is increased. Practically, the exposure time could not be increased indefinitely due to the finite well depth of the spectrometer. Though the 20 ms spectrum exhibits poor SNR, it is still possible to fit the spectrum and obtain an estimate of the methanol concentration. As shown in Eq. (1), the a and b parameters listed in the legend represent the intensity of the water and methanol spectra (respectively) extracted from the fitting algorithm. The values in parentheses are the standard deviations of these values. We note that negative values for some data points occur due to an electronic background subtraction performed automatically by the spectrometer.

Download Full Size | PDF

 figure: Fig. 5

Fig. 5 Precision of the Raman-based methanol sensor. (a–d) The methanol concentration acquired with various measurement times shows improved precision with increasing measurement time. A black dashed line is drawn at the mean concentration, and gray dashed lines are drawn at ±1% to facilitate comparison. (e) The overlapping Allan deviation of each time series (a–d) provides an estimate of the uncertainty in the measurement for a given measurement time. There is a clear trend that the uncertainty decreases with increasing spectrometer integration time, and with the overall measurement time. For the spectrometer integration longer than 200 ms, the precision roughly follows the trend of 0.2%×1/τ. (Note: 100-second measurement time (d) was acquired using an exposure time of 2000 ms and averaging of 50 spectra. Additionally, the data in (d) were corrected for a slow drift, as described in Section 3.3, and (d) presents the deviation from the fit.)

Download Full Size | PDF

In order to visualize the relationship between measurement time and measurement uncertainty, we employ the overlapping Allan deviation [27, 28], which is the IEEE standard method [29] for quantifying the instability of frequency references, such as atomic clocks. The Allan deviation is closely related to the standard deviation, but it provides more robust performance in situations that exhibit “random walk” noise. Its use here allows for a straightforward analysis of the achievable precision for a given measurement time. Indeed, the overlapping Allan deviation presented in Fig. 5(e), enables the comparison of the precision achieved by various exposure times of the spectrometer. For a given exposure time, the uncertainty versus the measurement time (τ) roughly follows the 1/τ trend would be expected for a signal that randomly varies with white (flat) frequency noise. In addition, the measurement becomes more precise as the spectrometer exposure time is increased. Thus, we conclude that the best precision is offered by setting the spectrometer integration time to be as long as the desired measurement time. In other words, Averaging Method A appears to be the best option.

The uncertainty of the measurement continues to decrease with increasing measurement time until it reaches a value of 0.008% at a 400 second measurement time. This behavior suggests that the sensor noise of the spectrometer is the limiting factor for measurement times less than about 400 seconds. At measurement times longer than ~1000 seconds, long-terms drifts begin to exceed the uncertainty of the measurement, and the Allan deviation [Fig. 5(e)] rises to the 0.02% level. For times longer than 10,000 seconds, the Allan deviation begins to decrease, indicating the lack of fluctuations on timescales longer than a few hours.

3.3. Long-term stability

In many applications, a sensor would be used to monitor the concentration of a sample for many hours, days, or years. Consequently, a sensor that exhibits minimal drift on long timescales is desirable. To test the stability of our measurement system, we measured the concentration methanol in a methanol–water mixture over 11.2 days (Fig. 6). The test was conducted in a non-laboratory environment where twice-daily temperature variations from approximately 12 to 20 C are programmed into the building heating system, with larger local temperature excursions likely. We saved complete spectra continuously over this period, and, to keep the data storage requirements reasonable, we set the spectrometer to average 50 spectra, each with an exposure time of 2 seconds, providing one spectra every 100 seconds. In principle, the concentration could be monitored using a much finer sampling interval with similar results.

 figure: Fig. 6

Fig. 6 The concentration of a methanol–water mixture measured constantly over an 11-day period. (a) The concentration shows a slow decrease over the sampling period, which likely results from the evaporation of methanol from the sample. Since the concentration due to evaporation should decay exponentially, an exponential function provides a good fit to the slow trend in the data. (b) With the slow drift subtracted, faster drifts(on the several hours timescale) can be seen and correspond to systematic error in the measurement. These errors are less than 0.1% for the entire measurement period and have a standard deviation of 0.02%. We speculate that these drifts in the measured concentration are associated with drifts in the laser frequency. Two data points, corresponding with the opening of the box lid to inspect the sample, were omitted. Otherwise, all data points are plotted, indicating the absence of any spurious readings across the entire measurement period.

Download Full Size | PDF

Over the 269.4 hour measurement time, the measured concentration of methanol exhibits an overall decrease from 44.3% to 42.8% methanol [Fig. 6(a)]. We attribute this slow decrease to an imperfect seal between the teflon cap and the quartz cuvette, which allowed for some evaporation of the mixture. Since methanol has a higher vapor pressure than water, it evaporates more quickly, leaving the remaining mixture with a lower concentration of methanol. Since the relative evaporation rate of methanol should be proportional to the concentration of methanol, the decay of methanol concentration should be described by

wme=w0et/tevap,
where w0 is the initial methanol concentration and tevap is the lifetime of the evaporation process. These parameters determined by a similar nonlinear fitting procedure to that described in Section 2.2. By subtracting the exponential function from the data, we can visualize the instrumental fluctuations on the hours-to-days timescale, without the concentration change [Fig. 6(b)]. The extracted lifetime of the evaporation process tevap is 315 days. Consequently, because we subtract an exponential with this time constant, we are not sensitive to systematic drifts on similar timescales. Nevertheless, Fig. 6b provides insight into systematic errors on the hours and days timescales.

Figure 6(b) shows that the long-term fluctuations of the concentration measurement are below 0.1% over the entire 11.2 day measurement time, are typically below 0.05%, and have a standard deviation of 0.021%. This is below the accuracy of the measurement estimate in Section 3.2 and is comparable to the precision achieved with an averaging time of a few-10s of seconds. Consequently, for many applications, the systematic drifts of the measurement will be hidden under the uncertainty of the measurement given by the integration time.

However, even higher precision could be achieved if these systematic fluctuations could be eliminated, so we speculate about the cause of the fluctuations seen in Fig. 6b. One explanation is that the drifts are real changes in the local methanol concentration probed by the focused laser. During the measurement, we noticed the formation of droplets on the walls of the cuvette (above the level of the liquid). It is likely that the heat from the laser, or thermal gradients within the RCS enclosure, drives an evaporation-condensation cycle within the cuvette, which causes the formation of small droplets on the walls of the cuvette. Since the vapor pressure of methanol is higher than that of water, these droplets would likely have a much higher methanol concentration than the rest of the solution. As the droplets fall back into the solution, the methanol concentration would increase.

However, we find that the fluctuations display a strong Fourier component near 1/12 hours, suggesting that the fluctuations are related to the twice-daily cycling of the building heating system. There are numerous ways that temperature fluctuations could affect the observed concentration, including changing the alignment of the optical system, affecting the response of the spectrometer, or through a change of the Raman spectrum as a function of temperature. However, we suspect that a change in the wavelength of the 532-nm laser is the most likely cause of the drift. Lasers that have not been specifically stabilized to be “single frequency” are well known to exhibit a strong dependence between the wavelength and ambient temperature, and can experience “mode-hops” where a small change in temperature results in a sudden change in the laser wavelength [30]. The use of a wavelength-stabilized laser, such as a distributed-feedback (DFB) laser could reduce these drifts and improve the precision of the measurement. A less expensive alternative might be to use a light-emitting diode (LED) as an illumination source [31, 32].

4. Future design improvements

Above, we have presented a proof-of-principle experiment showing that a compact, accurate, precise, and reliable concentration sensor can be built using relatively inexpensive components. Depending on the application, the sensor described here might provide sufficient performance, or may require further optimization. Here we discuss possible strategies that may be employed for optimizing the size/weight, power consumption, optical access, and robustness to background signals.

4.1. Size, weight, and cost

Currently, the size and weight of the RCS are limited by the use free-space optical components and their associated mounting hardware. The use of fiber-optic components (or miniaturized free-space components) would allow the footprint to be greatly decreased. The cost is currently dominated by the commercial spectrometer. As spectrometer units continue to decrease in price, the use of RCS sensors should become more economical. Alternatively, so-called “spectral sensors” provide detection at just a few points across the spectrum, but are very inexpensive, compact, and power efficient. It is possible that recording Raman spectral information at just a few wavelengths is sufficient to estimate the concentration of a simple mixture.

4.2. Power consumption

The specified power consumption of the spectrometer is 250 mA, and the laser is also specified to consume 250 mA. Both require 5 V of DC input. Thus, the nominal power consumption of the instrument is ∼2.5 watts, plus whatever power is required to run the electronics that acquires and processes the data from the spectrometer. In this case, a laptop computer was used to control the spectrometer, greatly increasing the power consumption of the entire system. However, since the processing requirements for recording and processing spectrometer data are quite modest, a micro-controller or a compact single-board computer (for example, a Raspberry Pi, which consumes ∼1 watt) could replace the computer and keep the total power consumption below 4 watts.

4.3. Suppression of photoluminescent background

In many applications, it’s likely that the measured mixtures will contain moderate amounts of impurities, which may provide strong photoluminescent background signals. Indeed, in our initial tests, we used methanol samples with lower purity, and saw significant photoluminuscent background signals.

One method of suppressing the background is to change the wavelength of the laser. Since the probability for Raman scattering scales as λ−4, our 532-nm laser provides significantly higher Raman signal per input photon than longer-wavelength lasers. However, the amount of photoluminescent background increases with shorter wavelengths. This may complicate the analysis, and prevent the reliable extraction of concentrations. Consequently, for practical applications, the use of a longer wavelength laser may provide the best performance. However, the use of longer wavelength lasers produces longer wavelength Raman signals, and detectors for longer wavelength light do not typically match the performance–cost ratio of silicon-based detectors.

Another method to suppress the photoluminescent background is shifted excitation Raman difference spectroscopy (SERDS) [33, 34]. In SERDS, the Raman spectrum is collected at two different laser wavelengths that are separated by only a few nm. The photoluminescent background will be almost exactly the same in the two spectra, but the Raman peaks will shift. Subtraction of the two spectra, combined with post-processing, can generate a Raman spectrum without the photoluminescent background.

5. Conclusion

Here we have demonstrated a proof-of-principle device for using Raman spectroscopy to measure the concentration of a mixture of water and methanol. The system consists of all commercial components, including a low-power 532-nm laser and a silicon-detector-based spectrometer. With an empirical correction factor, we can measure the concentration of methanol to better than 0.52% accuracy and to 0.006% precision. The instrument can operate over the entire range of methanol concentrations, from 0 to 100%. We operated the sensor for over 11 days of continuous measurements, and found that the systematic fluctuations remained below 0.1%. We expect that similar performance should be realized for other water-alcohol mixtures, and likely for any simple liquid mixture where the components have distinct Raman spectral signatures. We have outlined several ways that the design of the concentration sensor may be improved, providing an avenue for the development of field-deployable Raman-based concentration sensors.

Acknowledgments

We acknowledge helpful discussions with Kyle Schnitzenbaumer, Paul Vallett, David Carlson, Ron Johnson, and Carol Johnson. We thank Kate and Milo Goldfarbmuren for their moral support.

R.G., J.D., and L.E. are employed by Rebound Technologies, a company that manufactures refrigeration units based on liquid mixtures, including water and methanol. D.D.H works for Kapteyn-Murnane Laboratories, Inc., a company that manufactures femtosecond laser systems. The authors declare no other competing interests.

References

1. C. H. Camp Jr, Y. J. Lee, J. M. Heddleston, C. M. Hartshorn, A. R. H. Walker, J. N. Rich, J. D. Lathia, and M. T. Cicerone, “High-speed coherent Raman fingerprint imaging of biological tissues,” Nat. Photonics 8, 627 (2014). [CrossRef]  

2. M. T. Bremer and M. Dantus, “Standoff explosives trace detection and imaging by selective stimulated Raman scattering,” Appl. Phys. Lett. 103, 061119 (2013). [CrossRef]  

3. J. Kim, J. Noh, H. Chung, Y.-A. Woo, M. S. Kemper, and Y. Lee, “Direct, non-destructive quantitative measurement of an active pharmaceutical ingredient in an intact capsule formulation using Raman spectroscopy,” Anal. Chimica Acta 598, 280 (2007). [CrossRef]  

4. W. W. Rudolph and G. T. Hefter, “Quantitative analysis in alkaline aluminate solutions by Raman spectroscopy,” Anal. Methods 1, 132–138 (2009).

5. Z. Du, J. Chen, W. Ye, J. Guo, X. Zhang, and R. Zheng, “Investigation of two novel approaches for detection of sulfate ion and methane dissolved in sediment pore water using Raman spectroscopy,” Sensors 15, 12377–12388 (2015). [CrossRef]   [PubMed]  

6. Z. Li, M. J. Deen, S. Kumar, and P. R. Selvaganapathy, “Raman spectroscopy for in-line water quality monitoring instrumentation and potential,” Sensors 14, 17275–17303 (2014). [CrossRef]  

7. T. H. Kauffmann and M. D. Fontana, “Optical sensor of salt concentration: uncertainty evaluation,” Sens. Actuators B: Chem. 161, 21–27 (2012). [CrossRef]  

8. R. Pandey, S. K. Paidi, T. A. Valdez, C. Zhang, N. Spegazzini, R. R. Dasari, and I. Barman, “Noninvasive monitoring of blood glucose with Raman spectroscopy,” Accounts Chem. Res. 50, 264–272 (2017). [CrossRef]  

9. V. S. Tiwari, R. R. Kalluru, F. Y. Yueh, J. P. Singh, W. S. Cyr, and S. K. Khijwania, “Fiber optic Raman sensor to monitor the concentration ratio of nitrogen and oxygen in a cryogenic mixture,” Appl. Opt. 46, 3345–3351 (2007). [CrossRef]   [PubMed]  

10. C. A. Holden, S. S. Hunnicutt, R. Sánchez-Ponce, J. M. Craig, and S. C. Rutan, “Study of complexation in methanol/water mixtures by infrared and Raman spectroscopy and multivariate curve resolution–alternating least-squares analysis,” Appl. Spectrosc. 57, 483–490 (2003). [CrossRef]   [PubMed]  

11. F. M. Zehentbauer, E. J. Bain, and J. Kiefer, “Multiple parameter monitoring in a direct methanol fuel cell,” Meas. Sci. Technol . 23, 045602 (2012). [CrossRef]  

12. Y. Yu, W. Fan, Y. Wang, X. Zhou, J. Sun, and S. Liu, “C–H...O interaction in methanol–water solution revealed from Raman spectroscopy and theoretical calculations,” The J. Phys. Chem. B 121, 8179–8187 (2017). [CrossRef]  

13. J. Kiefer, T. Seeger, S. Steuer, S. Schorsch, M. C. Weikl, and A. Leipertz, “Design and characterization of a Raman-scattering-based sensor system for temporally resolved gas analysis and its application in a gas turbine power plant,” Meas. Sci. Technol . 19, 085408 (2008). [CrossRef]  

14. R. Sarhadian, “Freeze point suppression cooling system with embedded thermal energy storage for refrigeration application,” Tech. rep., Southern California Edison, Report cET15SCE1230 (2017).

15. S. Doerner, T. Schultz, T. Schneider, K. Sundmacher, and P. Hauptmann, “Capacitive sensor for methanol concentration measurement in direct methanol fuel cells (DMFC),” in IEEE SENSORS (2004), pp. 639–641.

16. C. Mohr, C. L. Spencer, and M. Hippler, “Inexpensive Raman spectrometer for undergraduate and graduate experiments and research,” J. Chem. Educ. 87, 326–330 (2010). [CrossRef]  

17. M. Jacox, “Vibrational and electronic energy levels of polyatomic transient molecules,” in NIST Chemistry WebBook, NIST Standard Reference Database Number 69, P. Linstrom and W. Mallard, eds. (National Institute of Standards and Technology, Gaithersburg, Maryland, 2018).

18. K. Levenberg, “A method for the solution of certain non-linear problems in least squares,” Q. J. Appl. Math. 2, 164–168 (1944). [CrossRef]  

19. J. J. Moré, D. C. Sorensen, K. E. Hillstrom, and B. S. Garbow, The MINPACK Project, in Sources and Development of Mathematical Software(Prentice-Hall, Inc., 1984).

20. “SciPy: Open source scientific tools for Python,” (2001–2018). [Online; accessed 2018-05-02].

21. S. Gawinkowski, A. Kaminska, T. Rolinski, and J. Waluk, “A new algorithm for identification of components in a mixture: application to Raman spectra of solid amino acids,” Analyst 139, 5755–5764 (2014). [CrossRef]   [PubMed]  

22. K. Hashimoto, M. Takahashi, T. Ideguchi, and K. Goda, “Broadband coherent raman spectroscopy running at 24,000 spectra per second,” Nat. Sci. Rep . 6, 21036 (2016). [CrossRef]  

23. D. K. Singh, S. K. Srivastava, S. Schlücker, R. K. Singh, and B. P. Asthana, “Self-association and hydrogen bonding of propionaldehyde in binary mixtures with water and methanol investigated by concentration-dependent polarized Raman study and DFT calculations,” J. Raman Spectrosc. 42, 851–858 (2010). [CrossRef]  

24. M. Nedic, T. N. Wassermann, R. W. Larsen, and M. A. Suhm, “A combined Raman- and infrared jet study of mixed methanol–water and ethanol–water clusters,” Phys. Chem. Chem. Phys. 13, 14050–14063 (2011). [CrossRef]  

25. S. Burikov, T. Dolenko, S. Patsaeva, Y. Starokurov, and V. Yuzhakov, “Raman and IR spectroscopy research on hydrogen bonding in water-ethanol systems,” Mol. Phys. 108, 2427–2436 (2010). [CrossRef]  

26. F. Li, Z. Men, S. Li, S. Wang, Z. Li, and C. Sun, “Study of hydrogen bonding in ethanol-water binary solutions by Raman spectroscopy,” Spectrochimica Acta Part A: Mol. Biomol. Spectrosc. 189, 621 – 624 (2018). [CrossRef]  

27. D. W. Allan, “Statistics of atomic frequency standards,” Proc. IEEE 54, 221–230 (1966). [CrossRef]  

28. D. A. Howe, D. U. Allan, and J. A. Barnes, “Properties of signal sources and measurement methods,” in Thirty Fifth Annual Frequency Control Symposium, (1981), pp. 669–716. [CrossRef]  

29. “IEEE standard definitions of physical quantities for fundamental frequency and time metrology–random instabilities,” IEEE Std 1139-2008 (Revision IEEE Std 1139-1999) pp. 1–50 (2009).

30. P. J. Herre and U. Barabas, “Mode switching of Fabry-Perot laser diodes,” IEEE J. Quantum Electron. 25, 1794–1799 (1989). [CrossRef]  

31. J. S. Greer, G. I. Petrov, and V. V. Yakovlev, “Raman spectroscopy with LED excitation source,” J. Raman Spectrosc . 44, 1058–1059 (2013). [CrossRef]  

32. M. A. Schmidt and J. Kiefer, “Polarization-resolved high-resolution Raman spectroscopy with a light-emitting diode,” J. Raman Spectrosc. 44, 1625–1627 (2013). [CrossRef]  

33. M. T. Gebrekidan, C. Knipfer, F. Stelzle, J. Popp, S. Will, and A. Braeuer, “A shifted-excitation raman difference spectroscopy (SERDS) evaluation strategy for the efficient isolation of Raman spectra from extreme fluorescence interference,” J. Raman Spectrosc. 47, 198–209 (2015). [CrossRef]  

34. E. Cordero, F. Korinth, C. Stiebing, C. Krafft, I. W. Schie, and J. Popp, “Evaluation of shifted excitation raman difference spectroscopy and comparison to computational background correction methods applied to biochemical Raman spectra,” Sensors 17, 1724 (2017). [CrossRef]  

Cited By

Optica participates in Crossref's Cited-By Linking service. Citing articles from Optica Publishing Group journals and other participating publishers are listed here.

Alert me when this article is cited.


Figures (6)

Fig. 1
Fig. 1 The experimental scheme for the Raman-spectroscopy-based concentration sensor (RCS). A 532-nm laser is focused into the sample. Some photons experience Stokes Raman scattering, and are emitted in all directions. The back-scattered Raman-shifted photons (yellow) are collected by the same lens, but transmit through the longpass dichroic mirror. After passing through a long-pass filter to further reject elastically scattered 532-nm light from the sample, the Raman-shifted light is couped into a fiber which transports it into the grating-based spectrometer.
Fig. 2
Fig. 2 (a) Raman spectra for pure methanol and pure water samples. These spectra are an average of 100 spectra, each collected for 2 seconds. The various peaks in the methanol spectrum can be assigned to various vibrational transitions listed in the NIST WebBook [17]. (b) The concentration of the mixture of can be estimated by fitting the observed spectrum (green) with a linear combination (magenta) of water (blue) and methanol (red) reference spectra.
Fig. 3
Fig. 3 The accuracy of the Raman concentration sensor. (a) The extracted mass fraction of a methanol–water mixture extracted from the Raman spectrum versus the mass fraction measured using a scale. Each point corresponds to a separate measurement, consisting of an average of 10 spectra, each collected using a 2-second exposure time. The various colors indicate independent experiments, conducted on over three days, demonstrating the repeatability of the measurements. A systematic underestimate of the methanol concentration is observed, which becomes most pronounced for methanol concentrations near 50%. (b) The difference between the points in panel (a) and the true mass fraction [dashed line in panel (a)] highlight a systematic underestimate of the methanol concentration. The error-bars indicate the 95-percent confidence interval extracted from the fitting procedure. The empirical correction (blue line) is a simple cosine function that has been fit to the data to model the systematic error. (c) Using this empirical correction function, we can correct the observed methanol concentrations and achieve improved agreement between the observed and true methanol concentrations. We note that the error from sample preparation is estimated to be ±0.15%, which is smaller than the width of the data points.
Fig. 4
Fig. 4 The Raman spectra (green dots) show improved signal-to-noise ratio (SNR) as the exposure time is increased. Practically, the exposure time could not be increased indefinitely due to the finite well depth of the spectrometer. Though the 20 ms spectrum exhibits poor SNR, it is still possible to fit the spectrum and obtain an estimate of the methanol concentration. As shown in Eq. (1), the a and b parameters listed in the legend represent the intensity of the water and methanol spectra (respectively) extracted from the fitting algorithm. The values in parentheses are the standard deviations of these values. We note that negative values for some data points occur due to an electronic background subtraction performed automatically by the spectrometer.
Fig. 5
Fig. 5 Precision of the Raman-based methanol sensor. (a–d) The methanol concentration acquired with various measurement times shows improved precision with increasing measurement time. A black dashed line is drawn at the mean concentration, and gray dashed lines are drawn at ±1% to facilitate comparison. (e) The overlapping Allan deviation of each time series (a–d) provides an estimate of the uncertainty in the measurement for a given measurement time. There is a clear trend that the uncertainty decreases with increasing spectrometer integration time, and with the overall measurement time. For the spectrometer integration longer than 200 ms, the precision roughly follows the trend of 0.2 % × 1 / τ. (Note: 100-second measurement time (d) was acquired using an exposure time of 2000 ms and averaging of 50 spectra. Additionally, the data in (d) were corrected for a slow drift, as described in Section 3.3, and (d) presents the deviation from the fit.)
Fig. 6
Fig. 6 The concentration of a methanol–water mixture measured constantly over an 11-day period. (a) The concentration shows a slow decrease over the sampling period, which likely results from the evaporation of methanol from the sample. Since the concentration due to evaporation should decay exponentially, an exponential function provides a good fit to the slow trend in the data. (b) With the slow drift subtracted, faster drifts(on the several hours timescale) can be seen and correspond to systematic error in the measurement. These errors are less than 0.1% for the entire measurement period and have a standard deviation of 0.02%. We speculate that these drifts in the measured concentration are associated with drifts in the laser frequency. Two data points, corresponding with the opening of the box lid to inspect the sample, were omitted. Otherwise, all data points are plotted, indicating the absence of any spurious readings across the entire measurement period.

Equations (3)

Equations on this page are rendered with MathJax. Learn more.

f mix ( λ ) = a f w ( λ ) + b f me ( λ ) + c λ + d ,
w me m me m w + m me b ρ me a ρ w + b ρ me ,
w me = w 0 e t / t evap ,
Select as filters


Select Topics Cancel
© Copyright 2024 | Optica Publishing Group. All rights reserved, including rights for text and data mining and training of artificial technologies or similar technologies.