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

Optical refractometry using lensless holography and autofocusing

Open Access Open Access

Abstract

Conventional optical refractometry methods are often limited by a narrow measurement range, complex hardware, or relatively high cost. Here, we present a novel refractometry method to measure the bulk refractive index (RI) of materials (including solids and liquids) using lensless holographic on-chip imaging and autofocusing, which is simple, cost-effective, and has a large RI measurement range. As a proof of concept, two compact prototypes were built to measure the RIs of solid materials and liquids, respectively, and they were tested by measuring the RIs of a ZnSe plate and a microscopy immersion oil. Experimental results show that our devices have an average accuracy of ~3 × 10−4 RI unit (RIU) with an estimated precision of ~3 × 10−3 RIU for solids; and an average accuracy of ~1 × 10−4 RIU with an estimated precision of ~3 × 10−4 RIU for liquids. We believe that this cost-effective and portable RI measurement platform holds promise to be used in laboratory and industrial settings.

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

1. Introduction

Refractive index (RI) is a fundamental physical property of materials. The measurement of RI, i.e., refractometry, has broad applications in chemical sensing [1–3], biology [4–8], the food industry [9–11], optical fabrication [12–15], among others. In the past decades, several methods for bulk RI measurement have been developed based primarily on geometric optics [16–20] or wave optics [21–28]. On the basis of geometric optics, the Snell’s law and total internal reflection (TIR) are commonly used principles. The minimum deviation method [16], in the light of Snell’s law, is widely used for RI measurement of solid transparent materials such as glasses, which provides an accuracy of ~10−5 refractive index unit (RIU) [17]. However, as a disadvantage of this method, one must prepare prism-shaped samples, making the routine measurements cumbersome. As an alternative, RI can also be determined by measuring the critical angle of TIR that occurs at the interface between the sample and a high-RI prism [18,19]. The Abbé refractometer [20], as one of the most famous refractometry instruments employing the TIR principle, has been widely used to measure the RI of solids with flat surfaces and liquids. Although the TIR-based methods provide a very good accuracy (e.g., ~10−3-10−5 RIU) for transparent liquids and solids, their measurement range is typically limited by the RI of the reference prism [20], which poses a practical limitation.

On the basis of wave optics, interferometric methods were also developed to detect a relative phase difference with high sensitivity, which can be used to measure the RI of the sample [21–25]. One advantage of interferometric methods is that the RI measurement range is not theoretically limited, which is far wider than TIR-based methods. As an example, a standard Michelson interferometer (MI) can be used to measure the RI of transparent solid plates [24], where a transparent sample plate is rotated in the sample arm continuously to change the optical path difference, and fringe patterns at different sample orientations are measured. By analyzing the fringe patterns, the RI can be estimated with an accuracy of ~10−3 RIU given that the thickness of the sample is known. However, since the Michelson interferometer is sensitive to air flow in the light path and the vibration of each component in the setup, the accuracy of the measurement is easily influenced by the environment. The Fabry-Perot (FP) interferometer [25] is another method used to measure the RI of a sample, where the fringe pattern is determined by the phase difference between the directly transmitted light wave and the multiply reflected waves between the reflectors. Thus, the FP arrangement is more robust to environmental perturbations than MI, and the accuracy can be as good as ~10−4 RIU. However, interferometers are in general limited by their relatively high cost, complex optical setup and sensitivity to alignment. Diffraction-based methods to measure RI have been demonstrated with simplified configurations [26–28] achieving an accuracy of e.g., ~10−4 RIU, but portable devices have not been reported. Besides these techniques, various other refractometry methods were proposed, including ellipsometers [29], optofluidics-based methods [30], and surface plasmon resonance (SPR) sensors [31]. They have several advantages, such as a low sample volume requirement and high accuracy, but they entail complex configurations and require precise and expensive optical components.

Here, we present a novel method to measure the RI of solids and liquids using holographic imaging [32–38] and autofocusing [39–44]. The measurement principle is based on the relation between RI and optical diffraction: in a holographic imaging system, a change in the RI of the medium after the light hits an object will affect the effective distance of optical diffraction, thereby changing the focus distance of the object’s hologram. Therefore, the RI can be indirectly measured by automatically measuring the focus distance of a two-dimensional (2D) object using a holographic autofocusing algorithm. Because this method does not involve TIR, the theoretical measurement range of this method is not limited, and therefore it can be used to measure both low- and high-RI samples. Based on this principle, we built two portable refractometer prototypes to measure RIs of solids and liquids, respectively. A lensless holographic on-chip imaging configuration was adopted due to its simplicity, cost-effectiveness and robustness to alignment errors. A laser diode provides a partially-coherent illumination, and a complementary metal-oxide-semiconductor (CMOS) image sensor array is used to record the diffraction patterns of a 2D object (e.g., a USAF-1951 target). The sample to be measured is inserted between this 2D object and the image sensor plane, and the auto-focus distances of the 2D object before and after the sample insertion are used to calculate the RI of the sample.

We demonstrated the accuracy as well as the wide measurement range of our technique by measuring the RIs of a ZnSe glass sample and a microscopy immersion oil sample. The average accuracy was quantified to be ~3 × 10−4 RIU for solids and ~1 × 10−4 RIU for liquids; and the precision was estimated as ~3 × 10−3 RIU for solids and ~3 × 10−4 RIU for liquids. The cost-effectiveness, wide RI measurement range and the applicability of this measurement approach to different samples (solids and liquids) can make this technique a versatile tool to measure RIs in laboratory as well as industrial settings.

2. Measurement principle

2.1 Scalar wave propagation under Fresnel approximation

As shown in Fig. 1, for a given complex optical field on the plane z = 0, U(x, y; 0), denote the complex optical field after it propagates by a distance of z as U(x, y; z). Let A(fx, fy; z) be the angular spectrum of U(x, y; z), with a Fourier transform relation

A(fx,fy;z)=U(x,y;z)exp[j2π(fxx+fyy)]dxdy.
A(fx, fy; z) and A(fx, fy; 0) are related by [45]:
A(fx,fy;z)=A(fx,fy;0)H(fx,fy;z).
The transfer function, H (fx, fy; z), is given by
H(fx,fy;z)={exp[j2πnλ0z1(λ0fxn)2(λ0fyn)2],fx2+fy2fcut-off0,others,
where λ0 is the wavelength in vacuum, n is the refractive index of the medium (assumed to be uniform across the sample), and fcut-off is the cut-off frequency for free-space propagation,
fcut-off=nλ0.
Under Fresnel approximation, where |λ0fx/n|≪1 and |λ0fy/n|≪1, H can be approximated by:
H(fx,fy;z)exp(j2πnzλ0)exp[jπλ0zn(fx2+fy2)].
Here, in Eq. (5), the first part of H (fx, fy; z) is a constant phase factor that is independent of fx and fy, which cannot be sensed by a typical CMOS or charge-coupled device (CCD) based image sensor. Thus, it is ignored in our discussion here and only the second part is kept, given by:
H˜(fx,fy;z)=exp[jπλ0zn(fx2+fy2)].
It is evident from Eq. (6) that the effect of wave propagation in a medium with RI = n over a distance z is equivalent to propagation in vacuum (where RI = 1) over a distance of z/n. In other words, having a larger RI reduces the amount of effective diffraction. Therefore, we define here a quantity named as “effective diffraction distance” (EDD), denoted by lE. For wave propagation from z = z1 to z = z2 with multiple layers of materials of different RIs in between, lE is given by:

 figure: Fig. 1

Fig. 1 Schematic of scalar wave propagation.

Download Full Size | PDF

lE=z1z21n(z)dz.

In practice, for a typical lensless holographic imaging system without pixel super-resolution [32,34,35,38], the additional frequency cut-off due to the finite spacing of the image sensor’s pixels ensures that Fresnel approximation is satisfied. For example, assuming a pixel pitch of Δx = Δy = 1.67 μm, according to the Nyquist theorem, |fx, y| ≤ 1/(2Δx) = 0.299 μm−1. Further assuming λ0 = 655 nm and n = 1, |λ0 fx, y /n| ≤ 0.299 × 0.655 = 0.194 ≪ 1, satisfying the Fresnel approximation condition.

2.2 Holographic autofocusing

Holographic autofocusing algorithms are used to automatically determine the wave propagation distance (i.e., the “focus distance”) of a hologram. Various autofocusing algorithms have been previously demonstrated [29–44] and compared [40,42]. In this work, we adopt a recently proposed edge sparsity-based autofocusing criterion named “Gini index of the gradient” (GoG), which is shown to be robust and accurate for a wide range of specimens [39,40]. GoG of a complex refocused (“back-propagated”) optical field using the angular spectrum method detailed in Section 2.1 is given by

GoG(U)=Gini(|U|),
whereis the 2D gradient operator; |⋅| is the complex modulus operator; U is the back-propagated complex optical field to a given distance z; and the Gini index of a 2D real-valued non-negative image is defined as:
Gini(C)=12k=1Na[k]sum(C)(Nk+0.5N),
where N is the number of pixels in the image; a[k] is the k-th sorted entry of the 2D image C in the ascending order, k = 1, 2, …, N; and sum(C) is the sum of all pixels in the image C. The autofocusing process using GoG can be expressed as a maximization problem,
zAF=argmaxz{GoG[(G;z)]},
where zAF is the calculated focus distance of the hologram intensity G; ℘(G; z) denotes the back-propagation of G by a distance of z using the angular spectrum method. When the coherent illumination light source is viewed as a point source, the wavefront curvature should be corrected. For this, before the autofocusing step, G should be multiplied by a phase factor exp(jφ) [46], i.e.,

φ=πλ0R[(xxc)2+(yyc)2],
where (xc, yc) define the coordinates of the intersection point of the optical axis with the image sensor, and R is the distance from the light source to the image sensor. The search for zAFis conducted by a global search algorithm based on the golden-section search method [40]. If the RI is assumed to be 1 during autofocusing (that is, the RI related to the back-propagation kernel ℘(G; z) is set as 1), the calculated zAF is an estimation of lE introduced in Eq. (7).

2.3 Optical refractometry using lensless holography and autofocusing

A typical lensless holographic imaging setup is depicted in Fig. 2(a), where a partially coherent quasi-planar light wave illuminates a high-resolution planar target (e.g., a positive USAF-1951 target), and its diffraction pattern is recorded on a CMOS/CCD image sensor. The EDD, Eq. (7), between the target and the image sensor will change if a sample plate with a thickness of d and a RI of nx is inserted between the sensor and target planes. Let lE,0 be the EDD without inserting the sample, and lE,1 with the sample. They are related to each other by:

lE,1=lE,0dna+dnx,
where na is the RI of the medium (typically air) without inserting the sample. Because lE,0 and lE,1 can be estimated by the holographic autofocusing algorithm, and assuming the thickness d is known accurately, nx can be obtained by:
nx=dlE,1lE,0+dna.
Note the thickness of a solid sample plate (d) can be measured mechanically, e.g., using a micrometer or a microscope. The thickness of a liquid sample can be more difficult to obtain directly. For this reason, we developed a three-step measurement procedure for liquid samples, detailed in Sec. 3.2.

 figure: Fig. 2

Fig. 2 Working principle of optical refractometry using lensless in-line holography and holographic autofocusing.

Download Full Size | PDF

3. Results and discussion

3.1 Portable holographic device for RI measurement of solid transparent plates

Based on the working principle detailed in Fig. 2 and Section 2.3, we custom-fabricated a holographic device for the RI measurement of solid samples, shown in Fig. 3. We use a laser diode (LD) (LD201710XHXL, QiLin optics, Nantong, China) with a peak wavelength of 655 nm as the illumination source, and a 10-megapixel CMOS image sensor (1.67 μm pixel size, 3840 × 2748 pixels; DMM 27UJ003-ML, ImagingSource Corporation, Germany) for hologram recording. A positive USAF-1951 target (Thorlabs, Inc.) with a 25.4 mm diameter is adopted as the “target object” used for autofocusing. Among various patterns of the USAF-1951 target, the high-resolution part (i.e., the central part) is more important as it contains richer fine features, useful for autofocusing. The housing of the device is made from 3D-printed plastic (printing resolution 0.1 mm; service provided by Zhantong 3D, Shenzhen, China), and steel posts are used to connect the illumination part and the imaging part. In order to demonstrate that the RI measurement range of our presented technique is very wide, we chose a ZnSe plate (GL2017102145R, QiLin optics, Nantong, China) as the sample with n = 2.5778 at λ = 655 nm (data obtained from Zemax software [47]), which falls outside the measurement range of a typical commercial Abbé refractometer (typically 1.3 < n < 1.8). The thickness d of the ZnSe plate is measured by a micrometer to be 1697.0 μm.

 figure: Fig. 3

Fig. 3 A custom-fabricated device to measure the RI of solids. RI measurement of a ZnSe plate is demonstrated here. (a) Photograph of our device. (b) Schematic side view of the inset in (a).

Download Full Size | PDF

Following the procedure detailed in Section 2.3, a hologram of the USAF-1951 target is captured before and after inserting the ZnSe sample plate, i.e., Fig. 4(a) and (d), respectively, and autofocusing is performed for each case to estimate the EDDs: lE,0 = 2811.4 μm and lE,1 = 1772.1 μm, respectively (Fig. 4(b) and (e)). An area with 1024 × 1024 pixels at the center of the CMOS image sensor is chosen for autofocusing. The wave curvature from the LD source is corrected by multiplying the hologram with a phase compensation term given by Eq. (11) [46]. The RI of air is taken as na = 1.0003 [48] in Eq. (13). Based on the measurements reported in Fig. 4 (i.e., lE,0 and lE,1) and the known quantities (i.e., d and na), the unknown RI of the solid material, nx, can be calculated as 2.5802. Including this one, we repeated the RI measurement of the same ZnSe plate for a total of 10 times (by repeating all the necessary steps), and obtained nx = 2.5781 ± 0.0016, with the individual RI measurement results as 2.5802, 2.5771, 2.5799, 2.5755, 2.5779, 2.5767, 2.5795, 2.5775, 2.5798, and 2.5767 RIU. According to Zemax software, the RI of ZnSe is n = 2.5778 at λ = 655 nm [47], which confirms that the mean RI accuracy of our measurement is ~3 × 10−4 RIU.

 figure: Fig. 4

Fig. 4 Experimental RI measurement of a ZnSe plate sample. (a-c) The acquired hologram and its autofocusing result before inserting the ZnSe plate between the USAF-1951 target and the CMOS image sensor. (d-f) The hologram and its autofocusing result after a ZnSe plate is inserted between the USAF-1951 target and the CMOS image sensor. (a) Full field-of-view (FOV) hologram and a zoom-in view of the high-resolution area of the USAF-1951 target (inset), without the ZnSe plate. (b) GoG autofocusing curve as a function of z. The discrete points in the curve represent the evaluated z-distances by the search algorithm (Section 2.2). (c) The hologram back-propagated to the z-distance determined by the autofocusing algorithm. (d) The hologram recorded after the ZnSe plate is inserted. Interference pattern is visible due to slightly nonparallel front and back surfaces of the ZnSe plate. (e, f) Autofocusing result and the back-propagation of the hologram in (d).

Download Full Size | PDF

We also provide here a theoretical analysis on our measurement error to understand the effect of each parameter. Our measurement accuracy is mainly limited by the accuracies of the micrometer, the autofocusing algorithm [40], and the RI of air [48]. The estimated uncertainties of these parameters are listed in Table 1.

Tables Icon

Table 1. Estimated uncertainties in the measured parameters (solid sample).

According to the Eq. (13), we can estimate the standard error of nx,

σ(nx)=(nxdδd)2+(nxlE,0δlE,0)2+(nxlE,1δlE,1)2+(nxnaδna)2,
where δd, δlE,0, δlE,1 and δna are uncertainties of d, lE,0, lE,1 and na, respectively. The partial derivatives in Eq. (14) are given by

nxd=lE,1lE,0(lE,1lE,0+dna)2=na2(lE,1lE,0)(nalE,1nalE,0+d)2,
nxlE,0=d(lE,1lE,0+dna)2=na2d(nalE,1nalE,0+d)2,
nxlE,1=d(lE,1lE,0+dna)2=na2d(nalE,1nalE,0+d)2,
nxna=d2na2(lE,1lE,0+dna)2=d2(nalE,1nalE,0+d)2,

After inserting the known and measured values into Eqs. (15)-(18) (i.e., na = 1.0003, d = 1697.0 μm, lE,0 = 2811.4 μm and lE,1 = 1772.1 μm), we obtain |nxdδd|=2.4×103, |nxlE,0δlE,0|=3.9×103, |nxlE,1δlE,1|=3.9×103, and |nxnaδna|=6.7×106. Based on these, we can estimate the total standard error of nx as: σ(nx)=6.1×103.This standard error is approximately 4-fold larger than our measurement standard deviation (0.0016 RIU); this might be partially due to that fact that the above analysis assumes that lE,0 and lE,1 are independent from each other. In reality, the autofocusing algorithm can output a z-distance that is consistently offset by a certain amount from the true z-distance, such that δlE,0 and δlE,1 can partially cancel each other when subtracted in Eq. (13). This is confirmed by the relatively small standard deviation (0.0016 RIU) of 10 individual experiments, which would otherwise be significantly larger if the autofocusing algorithm’s error (~1 μm) is entirely random. As a result, we take the experimentally measured 0.0016 RIU as the random error related to uncertainties in δlE,0, δlE,1 and na. Another major source of error comes from the micrometer’s inaccuracy. Here, we used a micrometer with an accuracy of 1 μm, which, by itself, results in a δnx of 0.0024 RIU. Because we used the same d throughout the 10 measurements, this error source is not reflected in the standard deviation of our experiments, i.e., 0.0016 RIU. Considering the error introduced by the micrometer is 0.0024 RIU and the measurement random error is 0.0016 RIU, the total random error (i.e., precision) of the current device for measuring solids can be estimated to be 0.0029 RIU, i.e., approximately 3 × 10−3 RIU.

The additional error that is related to the tilting of the test sample is not considered in these calculations, but it can be avoided or mitigated through the design and precise fabrication of the device. In our design, we placed the solid sample plate into a slot inside the sample holder, and placed the device in an upright (i.e., vertical) position such that the sample lies flat in the sample holder. This ensures that the sample is parallel to the image sensor and the target, mitigating errors due to sample tilting. The RI error related to thickness inaccuracy can also be partially mitigated by increasing the thickness of the sample. Moreover, because most of the errors are random errors, we believe it is possible to reduce the error significantly by averaging results from multiple experiments. In summary, we estimate that our current prototype achieves a mean accuracy of ~3 × 10−4 RIU with an estimated precision of ~3 × 10−3 RIU for the RI measurement of solids.

3.2 Portable holographic device for RI measurement of liquids

The photograph and design schematic of a portable refractometer for measuring the RIs of liquid samples is presented in Fig. 5, which is very similar to the one shown in Fig. 3. The illumination is similarly provided by a LD. A custom-made liquid chamber, which can be inserted in/out, is placed between the USAF-1951 target and the CMOS image sensor. The front surface (close to the USAF-1951 target) and the back surface (close to the CMOS image sensor) of the liquid chamber are both transparent glass cover slips (18 × 18 × 0.1 mm). The inner thickness of the liquid chamber is roughly 3 mm. Owing to the design of the liquid chamber, the device in Fig. 5 is in a horizontal position during the measurements.

 figure: Fig. 5

Fig. 5 Photograph and schematic of a portable device for measuring the RI of liquid samples based on lensless holographic on-chip imaging and autofocusing. (a) Photograph of our portable device. (b) 3D design schematic. (c) Zoom-in top-view of the device excluding the laser diode used for illumination.

Download Full Size | PDF

We tested this prototype using a standard immersion oil (Type 300, Cargille Labs, US) under room temperature (T = 25°C). Because the exact thickness of the liquid chamber is unknown, a three-step measurement procedure was developed, which also involves calibration of the chamber itself:

  • 1) When the chamber is empty (filled with air, na = 1.0003), a hologram H0 is recorded by the CMOS image sensor (left image in Fig. 6(a)). The EDD lE,0 is calculated from H0 by autofocusing. In this experiment, lE,0 was measured as 4539.1 μm (right curve in Fig. 6(a)).
  • 2) Next, the chamber is filled with deionized (DI) water (RI: nw = 1.3310 [49]). A
 figure: Fig. 6

Fig. 6 Experimental results for measuring the RI of liquid sample. (a-c) Hologram (left) and autofocus curve (right) when the chamber is with (a) air; (b) DI water; and (c) immersion oil to be tested.

Download Full Size | PDF

hologram H1 is recorded (left pattern in Fig. 6(b)), and the EDD lE,1 is calculated. In this experiment, lE,1 was measured as 3831.1 μm (right curve in Fig. 6(b)).

  • 3) Lastly, the chamber is filled with the liquid sample to be tested (an immersion oil is tested here). A hologram H2 is recorded (left pattern in Fig. 6(c)), and the EDD lE,2 is calculated. In this experiment, lE,2 was measured as 3575.1 μm (right curve in Fig. 6(c)).

The RI of the unknown liquid sample that is calculated by this method is given by (see the Appendix for derivation):

nx=nanwlE,0nanwlE,1nwlE,2nalE,2+nalE,0nwlE,1,
where na and nw are RIs of air and DI water, respectively. According to Eq. (19), nx of this experiment (Fig. 6) is calculated to be 1.5117 (λ = 655 nm). Based on 10 independent measurements, we experimentally obtained nx = 1.5115 ± 0.0003 (with the individual RI measurement results as 1.5117, 1.5114, 1.5114, 1.5117, 1.5108, 1.5117, 1.5115, 1.5112, 1.5116, and 1.5119 RIU), while the RI data provided by Cargille Labs is 1.5116 [50]. Based on these measurements, the mean error of our device (i.e., accuracy) was estimated to be ~1 × 10−4 RIU, with an experimental standard deviation of ~3 × 10−4 RIU.

Similar to the error analysis reported earlier for solid RI measurements, we conducted an error analysis to understand the effect of each parameter. Table 2 lists the estimated uncertainties of the measured parameters. According to Eq. (19), the estimated standard error of the measured liquid RI (nx) is given by

σ(nx)=(nxlE,0δlE,0)2+(nxlE,1δlE,1)2+(nxlE,2δlE,2)2+(nxnaδna)2+(nxnwδnw)2,
where δlE,0, δlE,1, δlE,2, δna and δnw are uncertainties of lE,0, lE,1, lE,2, na and nw, respectively. The partial derivatives in Eq. (20) are given by:

Tables Icon

Table 2. Estimated uncertainties in measured parameters (liquid sample).

nxlE,0=nanw2lE,2na2nwlE,2nanw2lE,1+na2nwlE,1(nwlE,2nalE,2+nalE,0nwlE,1)2,
nxlE,1=na2nwlE,2na2nwlE,0+nanw2lE,0nanw2lE,2(nwlE,2nalE,2+nalE,0nwlE,1)2,
nxlE,2=na2nwlE,0na2nwlE,1nanw2lE,0+nanw2lE,1(nwlE,2nalE,2+nalE,0nwlE,1)2,
nxna=nw2lE,0lE,2nw2lE,1lE,2nw2lE,0lE,1+nw2lE,12(nwlE,2nalE,2+nalE,0nwlE,1)2,
nxnw=na2lE,02+na2lE,1lE,2na2lE,0lE,2na2lE,0lE,1(nwlE,2nalE,2+nalE,0nwlE,1)2.

After inserting the known and measured values into Eq. (21)-(25) (i.e., na = 1.0003, nw = 1.3310, lE,0 = 4539.1 μm, lE,1 = 3831.1 μm and lE,2 = 3575.1 μm), we obtain |nxlE,0δlE,0|=2.9×104, |nxlE,1δlE,1|=1.1×103, |nxlE,2δlE,2|=8.02×104, |nxnaδna|=8.3×107 and |nxnwδnw|=1.76×104. Thus, the estimated standard error of nx is σ(nx)=1.4×103, which is ~5 times larger than the standard deviation that is experimentally measured (i.e., ~3 × 10−4 RIU). However, analogous to the analysis in Sec. 3.1, we also assumed here that δlE,0, δlE,1 and δlE,2 are independent. In fact, δlE,0, δlE,1 and δlE,2 may cancel each other significantly as indicated by the low standard deviation of our RI measurements (~3 × 10−4 RIU), and the other errors due to δna and δnw are only on the order of ~10−4. Therefore, self-calibration of the sample thickness using reference liquids (e.g., DI water) is advantageous compared to direct thickness measurement of the sample using e.g., a micrometer, which improves our precision in liquid measurements compared to solid materials. In summary, we estimate that our current prototype achieves a mean accuracy of ~1 × 10−4 RIU and a precision of ~3 × 10−4 RIU for RI measurements of liquids.

The main components in our devices are a LD (0.05 USD), a positive USAF-1951 target (~200 USD) and a CMOS image sensor (~400 USD). For industrial applications, the positive USAF-1951 target can be substituted by other high-resolution planar patterns that may even be printed by a high-resolution printer, which will significantly reduce the cost of our instrument. When custom-designing a target pattern, one should note that the pattern should have various fine spatial features and sharp edges for accurate autofocusing [39,40], but the feature size must be larger than the resolution of the imaging system. Furthermore, various low-cost monochromatic CMOS/CCD image sensors can be used for image acquisition, which might altogether bring down the total cost of the parts of our platform to <100 USD under large volume manufacturing.

4. Conclusion

In this manuscript, we present a novel bulk refractometry principle based on lensless holographic on-chip imaging and autofocusing, which has the advantages of cost-effectiveness, simplicity, and a wide RI measurement range. As a proof of concept, we designed two compact and portable prototypes to measure the RIs of solids and liquids, respectively, and demonstrated RI measurement of a solid ZnSe plate and a microscopy immersion oil using these prototypes. We show that our devices can achieve an average RI measurement accuracy of ~3 × 10−4 RIU for solid samples, and ~1 × 10−4 RIU for liquid samples. This presented refractometry principle is promising to be widely applicable in laboratory and industrial settings.

Appendix Derivation of Eq. (19)

Let us denote the inner thickness of the liquid chamber as t; denote the EDD outside the liquid chamber (after light hits the 2D object) as lE,S, which should be constant. Thus, lE,0, lE,1 and lE,2 in Sec. 3.2 can be expressed as:

lE,0=lE,S+tna,
lE,1=lE,S+tnw,
lE,2=lE,S+tnx.

According to Eqs. (26) and (27), the inner thickness of the chamber for liquid, t, is

t=(lE,0lE,1)nanwnwna.

By replacing t in Eq. (28) with Eq. (29) and moving nx to the left and the other items to the right, Eq. (19) can be derived.

Funding

National Science Foundation; Howard Hughes Medical Institute (HHMI); National Natural Science Foundation of China (61575175, 61427818); National Key Research and Development Program of China (2016YFB1001502); China Scholarship Council (No. 201606320051).

Disclosures

The authors declare that there are no conflicts of interest related to this article.

References

1. C. R. Parkerson, “Determination of Saturation Temperatures of Inorganic Salt Solutions. Refractive Index Measurements,” Anal. Chem. 23(4), 610–613 (1951). [CrossRef]  

2. C. Y. Tan and Y. X. Huang, “Dependence of Refractive Index on Concentration and Temperature in Electrolyte Solution, Polar Solution, Nonpolar Solution, and Protein Solution,” J. Chem. Eng. Data 60(10), 2827–2833 (2015). [CrossRef]  

3. M. Yahya and M. Z. Saghir, “Prediction and Experimental Measurement of Refractive Index in Ternary Hydrocarbon Mixtures,” J. Chem. Eng. Data 60(8), 2329–2342 (2015). [CrossRef]  

4. W. Lukosz, “Integrated optical chemical and direct biochemical sensors,” Sens. Actuators B Chem. 29(1-3), 37–50 (1995). [CrossRef]  

5. R. Barer and D. A. T. Dick, “Interferometry and refractometry of snail amoebocytes,” Exp. Cell Res. 16(2), 285–291 (1959). [CrossRef]   [PubMed]  

6. T. R. Corner, T. C. Beaman, J. T. Greenamyre, and P. Gerhardt, “Spectrophotometric determination of refractive index increment for bacterial cells,” J. Microbiol. Methods 11(3-4), 255–260 (1990). [CrossRef]  

7. N. Lue, G. Popescu, T. Ikeda, R. R. Dasari, K. Badizadegan, and M. S. Feld, “Live cell refractometry using microfluidic devices,” Opt. Lett. 31(18), 2759–2761 (2006). [CrossRef]   [PubMed]  

8. R. Barer, “Refractometry and Interferometry of Living Cells,” J. Opt. Soc. Am. 47(6), 545–556 (1957). [CrossRef]   [PubMed]  

9. F. Rogerson and C. Symington, “A method for the estimation of alcohol in fortified wines using hydrometer Baume and refractometer Brix,” Am. J. Enol. Vitic. 57, 486–490 (2006).

10. M. L. Dongare, P. B. Buchade, and A. D. Shaligram, “Refractive index based optical Brix measurement technique with equilateral angle prism for sugar and Allied Industries,” Optik (Stuttg.) 126(20), 2383–2385 (2015). [CrossRef]  

11. A. F. Rawle, “Refractive Index Measurements,” in Encyclopedia of Spectroscopy and Spectrometry (Third Edition), J. C. Lindon, G. E. Tranter, and D. W. Koppenaal, eds. (Academic Press), pp. 936–945.

12. S. Mitachi and T. Miyashita, “Refractive-index dispersion for BaF2-GdF3-ZrF4-AlF3 glasses,” Appl. Opt. 22(16), 2419–2425 (1983). [CrossRef]   [PubMed]  

13. D. Malacara, Optical Shopping Test, 3rd ed. (John Wiley & Sons, Inc.).

14. D. Malacara and Z. Malacara, Handbook of Optical Design, 2nd ed. (Marcel Dekker, Inc.).

15. A. J. Werner, “Methods in High Precision Refractometry of Optical Glasses,” Appl. Opt. 7(5), 837–843 (1968). [CrossRef]   [PubMed]  

16. M. C. Simon and P. A. Larocca, “Minimum deviation for uniaxial prisms,” Appl. Opt. 34(4), 709–715 (1995). [CrossRef]   [PubMed]  

17. H. Hattori, H. Yamanaka, H. Kurniawan, S. Yokoi, and K. Kagawa, “Using minimum deviation of a secondary rainbow and its application to water analysis in a high-precision, refractive-index comparator for liquids,” Appl. Opt. 36(22), 5552–5556 (1997). [CrossRef]   [PubMed]  

18. H. Li and S. Xie, “Measurement method of the refractive index of biotissue by total internal reflection,” Appl. Opt. 35(10), 1793–1795 (1996). [CrossRef]   [PubMed]  

19. G. H. Meeten, “Refractive index errors in the critical-angle and the Brewster-angle methods applied to absorbing and heterogeneous materials,” Meas. Sci. Technol. 8(7), 728–733 (1997). [CrossRef]  

20. https://www.amazon.com/Azzota-AR-1-Abbé-Monocular-Refractometer/.

21. M. S. Shumate, “Interferometric Measurement of Large Indices of Refraction,” Appl. Opt. 5(2), 327–331 (1966). [CrossRef]   [PubMed]  

22. G. Coppola, P. Ferraro, M. Iodice, and S. De Nicola, “Method for measuring the refractive index and the thickness of transparent plates with a lateral-shear, wavelength-scanning interferometer,” Appl. Opt. 42(19), 3882–3887 (2003). [CrossRef]   [PubMed]  

23. L. Chen, X. Guo, and J. Hao, “Lens refractive index measurement based on fiber point-diffraction longitudinal interferometry,” Opt. Express 21(19), 22389–22399 (2013). [CrossRef]   [PubMed]  

24. G. D. Gillen and S. Guha, “Refractive-index measurements of zinc germanium diphosphide at 300 and 77 K by use of a modified Michelson interferometer,” Appl. Opt. 43(10), 2054–2058 (2004). [CrossRef]   [PubMed]  

25. G. D. Gillen and S. Guha, “Use of Michelson and Fabry-Perot interferometry for independent determination of the refractive index and physical thickness of wafers,” Appl. Opt. 44(3), 344–347 (2005). [CrossRef]   [PubMed]  

26. S. H. Lu, S. P. Pan, T. S. Liu, and C. F. Kao, “Liquid refractometer based on immersion diffractometry,” Opt. Express 15(15), 9470–9475 (2007). [CrossRef]   [PubMed]  

27. M. T. Tavassoly and A. Saber, “Optical refractometry based on Fresnel diffraction from a phase wedge,” Opt. Lett. 35(21), 3679–3681 (2010). [CrossRef]   [PubMed]  

28. M. T. Tavassoly, R. R. Naraghi, A. Nahal, and K. Hassani, “High precision refractometry based on Fresnel diffraction from phase plates,” Opt. Lett. 37(9), 1493–1495 (2012). [CrossRef]   [PubMed]  

29. M. Stchakovsky, Y. Battie, and A. E. Naciri, “An original method to determine complex refractive index of liquids by spectroscopic ellipsometry and illustrated applications,” Appl. Surf. Sci. 421, 802–806 (2017). [CrossRef]  

30. C. Prokop, N. Irmler, B. Laegel, S. Wolff, A. Mitchell, and C. Karnutsch, “Optofluidic refractive index sensor based on air-suspended SU-8 grating couplers,” Sens. Actuators A Phys. 263, 439–444 (2017). [CrossRef]  

31. L. Ji, X. Sun, G. He, Y. Liu, X. Wang, Y. Yi, C. Chen, F. Wang, and D. Zhang, “Surface plasmon resonance refractive index sensor based on ultraviolet bleached polymer waveguide,” Sens. Actuators B Chem. 244, 373–379 (2017). [CrossRef]  

32. W. Luo, A. Greenbaum, Y. Zhang, and A. Ozcan, “Synthetic aperture-based on-chip microscopy,” Light Sci. Appl. 4(3), e261 (2015). [CrossRef]  

33. C. Zuo, J. Sun, J. Zhang, Y. Hu, and Q. Chen, “Lensless phase microscopy and diffraction tomography with multi-angle and multi-wavelength illuminations using a LED matrix,” Opt. Express 23(11), 14314–14328 (2015). [CrossRef]   [PubMed]  

34. W. Luo, Y. Zhang, A. Feizi, Z. Göröcs, and A. Ozcan, “Pixel super-resolution using wavelength scanning,” Light Sci. Appl. 5(4), e16060 (2016). [CrossRef]   [PubMed]  

35. Y. Zhang, S. Y. C. Lee, Y. Zhang, D. Furst, J. Fitzgerald, and A. Ozcan, “Wide-field imaging of birefringent synovial fluid crystals using lens-free polarized microscopy for gout diagnosis,” Sci. Rep. 6(1), 28793 (2016). [CrossRef]   [PubMed]  

36. Y. Zhang, Y. Wu, Y. Zhang, and A. Ozcan, “Color calibration and fusion of lens-free and mobile-phone microscopy images for high-resolution and accurate color reproduction,” Sci. Rep. 6(1), 27811 (2016). [CrossRef]   [PubMed]  

37. M. Rostykus and C. Moser, “Compact lensless off-axis transmission digital holographic microscope,” Opt. Express 25(14), 16652–16659 (2017). [CrossRef]   [PubMed]  

38. A. Greenbaum and A. Ozcan, “Maskless imaging of dense samples using pixel super-resolution based multi-height lensfree on-chip microscopy,” Opt. Express 20(3), 3129–3143 (2012). [CrossRef]   [PubMed]  

39. M. Tamamitsu, Y. Zhang, H. Wang, Y. Wu, and A. Ozcan, “Comparison of Gini index and Tamura coefficient for holographic autofocusing based on the edge sparsity of the complex optical wavefront,” arXiv:1708.08055 [physics.optics] (2017).

40. Y. Zhang, H. Wang, Y. Wu, M. Tamamitsu, and A. Ozcan, “Edge sparsity criterion for robust holographic autofocusing,” Opt. Lett. 42(19), 3824–3827 (2017). [CrossRef]   [PubMed]  

41. Z. Ren, Z. Xu, and E. Y. Lam, “Learning-based nonparametric autofocusing for digital holography,” Optica 5(4), 337–344 (2018). [CrossRef]  

42. E. S. R. Fonseca, P. T. Fiadeiro, M. Pereira, and A. Pinheiro, “Comparative analysis of autofocus functions in digital in-line phase-shifting holography,” Appl. Opt. 55(27), 7663–7674 (2016). [CrossRef]   [PubMed]  

43. Z. Ren, N. Chen, and E. Y. Lam, “Automatic focusing for multisectional objects in digital holography using the structure tensor,” Opt. Lett. 42(9), 1720–1723 (2017). [CrossRef]   [PubMed]  

44. P. Gao, B. Yao, J. Min, R. Guo, B. Ma, J. Zheng, M. Lei, S. Yan, D. Dan, and T. Ye, “Autofocusing of digital holographic microscopy based on off-axis illuminations,” Opt. Lett. 37(17), 3630–3632 (2012). [CrossRef]   [PubMed]  

45. J.W.Goodman, Introduction to Fourier Optics (Roberts).

46. A. Khmaladze, M. Kim, and C.-M. Lo, “Phase imaging of cells by simultaneous dual-wavelength reflection digital holography,” Opt. Express 16(15), 10900–10911 (2008). [CrossRef]   [PubMed]  

47. https://www.zemax.com/.

48. P. E. Ciddor, “Refractive index of air: new equations for the visible and near infrared,” Appl. Opt. 35(9), 1566–1573 (1996). [CrossRef]   [PubMed]  

49. G. M. Hale and M. R. Querry, “Optical Constants of Water in the 200-nm to 200-microm Wavelength Region,” Appl. Opt. 12(3), 555–563 (1973). [CrossRef]   [PubMed]  

50. http://www.cargille.com/immeroilspecs.shtml.

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 Schematic of scalar wave propagation.
Fig. 2
Fig. 2 Working principle of optical refractometry using lensless in-line holography and holographic autofocusing.
Fig. 3
Fig. 3 A custom-fabricated device to measure the RI of solids. RI measurement of a ZnSe plate is demonstrated here. (a) Photograph of our device. (b) Schematic side view of the inset in (a).
Fig. 4
Fig. 4 Experimental RI measurement of a ZnSe plate sample. (a-c) The acquired hologram and its autofocusing result before inserting the ZnSe plate between the USAF-1951 target and the CMOS image sensor. (d-f) The hologram and its autofocusing result after a ZnSe plate is inserted between the USAF-1951 target and the CMOS image sensor. (a) Full field-of-view (FOV) hologram and a zoom-in view of the high-resolution area of the USAF-1951 target (inset), without the ZnSe plate. (b) GoG autofocusing curve as a function of z. The discrete points in the curve represent the evaluated z-distances by the search algorithm (Section 2.2). (c) The hologram back-propagated to the z-distance determined by the autofocusing algorithm. (d) The hologram recorded after the ZnSe plate is inserted. Interference pattern is visible due to slightly nonparallel front and back surfaces of the ZnSe plate. (e, f) Autofocusing result and the back-propagation of the hologram in (d).
Fig. 5
Fig. 5 Photograph and schematic of a portable device for measuring the RI of liquid samples based on lensless holographic on-chip imaging and autofocusing. (a) Photograph of our portable device. (b) 3D design schematic. (c) Zoom-in top-view of the device excluding the laser diode used for illumination.
Fig. 6
Fig. 6 Experimental results for measuring the RI of liquid sample. (a-c) Hologram (left) and autofocus curve (right) when the chamber is with (a) air; (b) DI water; and (c) immersion oil to be tested.

Tables (2)

Tables Icon

Table 1 Estimated uncertainties in the measured parameters (solid sample).

Tables Icon

Table 2 Estimated uncertainties in measured parameters (liquid sample).

Equations (29)

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

A ( f x , f y ; z ) = U ( x , y ; z ) exp [ j 2 π ( f x x + f y y ) ] d x d y .
A ( f x , f y ; z ) = A ( f x , f y ; 0 ) H ( f x , f y ; z ) .
H ( f x , f y ; z ) = { exp [ j 2 π n λ 0 z 1 ( λ 0 f x n ) 2 ( λ 0 f y n ) 2 ] , f x 2 + f y 2 f cut-off 0 , others ,
f cut-off = n λ 0 .
H ( f x , f y ; z ) exp ( j 2 π n z λ 0 ) exp [ j π λ 0 z n ( f x 2 + f y 2 ) ] .
H ˜ ( f x , f y ; z ) = exp [ j π λ 0 z n ( f x 2 + f y 2 ) ] .
l E = z 1 z 2 1 n ( z ) d z .
GoG ( U ) = Gini ( | U | ) ,
Gini ( C ) = 1 2 k = 1 N a [ k ] sum ( C ) ( N k + 0.5 N ) ,
z A F = arg max z { G o G [ ( G ; z ) ] } ,
φ = π λ 0 R [ ( x x c ) 2 + ( y y c ) 2 ] ,
l E, 1 = l E, 0 d n a + d n x ,
n x = d l E, 1 l E, 0 + d n a .
σ ( n x ) = ( n x d δ d ) 2 + ( n x l E,0 δ l E,0 ) 2 + ( n x l E,1 δ l E,1 ) 2 + ( n x n a δ n a ) 2 ,
n x d = l E,1 l E,0 ( l E,1 l E,0 + d n a ) 2 = n a 2 ( l E,1 l E,0 ) ( n a l E,1 n a l E,0 + d ) 2 ,
n x l E,0 = d ( l E,1 l E,0 + d n a ) 2 = n a 2 d ( n a l E,1 n a l E,0 + d ) 2 ,
n x l E,1 = d ( l E,1 l E,0 + d n a ) 2 = n a 2 d ( n a l E,1 n a l E,0 + d ) 2 ,
n x n a = d 2 n a 2 ( l E,1 l E,0 + d n a ) 2 = d 2 ( n a l E,1 n a l E,0 + d ) 2 ,
n x = n a n w l E,0 n a n w l E,1 n w l E,2 n a l E,2 + n a l E,0 n w l E,1 ,
σ ( n x ) = ( n x l E,0 δ l E,0 ) 2 + ( n x l E,1 δ l E,1 ) 2 + ( n x l E,2 δ l E,2 ) 2 + ( n x n a δ n a ) 2 + ( n x n w δ n w ) 2 ,
n x l E, 0 = n a n w 2 l E,2 n a 2 n w l E,2 n a n w 2 l E,1 + n a 2 n w l E,1 ( n w l E,2 n a l E,2 + n a l E,0 n w l E,1 ) 2 ,
n x l E , 1 = n a 2 n w l E,2 n a 2 n w l E,0 + n a n w 2 l E,0 n a n w 2 l E,2 ( n w l E,2 n a l E,2 + n a l E,0 n w l E,1 ) 2 ,
n x l E, 2 = n a 2 n w l E,0 n a 2 n w l E,1 n a n w 2 l E,0 + n a n w 2 l E,1 ( n w l E,2 n a l E,2 + n a l E,0 n w l E,1 ) 2 ,
n x n a = n w 2 l E,0 l E,2 n w 2 l E,1 l E,2 n w 2 l E,0 l E,1 + n w 2 l E,1 2 ( n w l E,2 n a l E,2 + n a l E,0 n w l E,1 ) 2 ,
n x n w = n a 2 l E,0 2 + n a 2 l E,1 l E,2 n a 2 l E,0 l E,2 n a 2 l E,0 l E,1 ( n w l E,2 n a l E,2 + n a l E,0 n w l E,1 ) 2 .
l E,0 = l E,S + t n a ,
l E,1 = l E,S + t n w ,
l E,2 = l E,S + t n x .
t = ( l E,0 l E,1 ) n a n w n w n a .
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.