1. Introduction
The refractive index is a key optical parameter for transparent and translucent materials. For liquids, it is closely related to parameters such as the wavelength of the light, the composition, the concentration, and the atmosphere temperature. The empirical expressions of the refractive index in terms of the concentration and temperature were derived for different liquids [
1]. To measure the refractive indices of liquids, a variety of non-destructive, sensitive refractometers [
2,
3,
4], reflectors [
5,
6], and other highly precise optical systems [
7,
8,
9] have been proposed and designed.
There have been many reviews and research studies on utilizing optical interferometers to measure the indices of the refraction of liquids with high accuracy [
10,
11,
12,
13,
14]. The Fabry–Perot interferometer is a good candidate due to its simplicity and versatility [
15,
16,
17]. In addition, a single-element interferometer to measure the refractive index was designated by rotating a rectangular optical glass cell containing the sample liquid and air simultaneously and by calculating the interference fringe shift number from an interferogram. The measurement accuracy for the refractive index was up to 10
−4 [
14]. A dual scanning interferometric technique was applied to measure the temperature and concentration dependence of the refractive indices of liquids at concentrations smaller than 0.03% and at a high precision of the order of 1
for a 20 mm optical path [
18]. High-precision refractive index measurements of liquids were also carried out for wavelengths from 0.5 to 2 μm using Rayleigh interferometry [
19].
Interferometers such as the Mach–Zehnder interferometer and the Michelson interferometer also play a prominent role in the field of precise measurement [
20]. The fiber-based Mach–Zehnder interferometric sensors were applied to measure the refractive index and temperature of transparent liquids from the wavelength shift of the interference spectrum [
21,
22,
23,
24]. A very common strategy to obtain the refractive index of a sample is to apply fringe count techniques to measure the fringe shift due to the introduction of the sample in one arm of an interferometer. For example, a hybrid optical fiber interferometer, which consists of a Mach–Zehnder interferometer and a Michelson interferometer, was fabricated to achieve accurate measurement of the refractive indices of different liquids [
25]. The standard deviation of the refractive index was about 10
−5 RIU. The refractive index of water was 1.33157 at 25 °C when the wavelength was 632 nm, very close to 1.33158 in ref. [
26]. The refractive index for NaCl solutions with concentrations smaller than 0.2% was measured. Using orbital angular momentum beams in a twisted Michelson interferometer opens the possibility for measurements of refractive index changes down to 10
−6 RIU. The refractive index was measured from the rotation of interference patterns at very small concentrations [
27]. The Michelson interferometer was applied to measure the refractive indices of liquids from measuring the fringe shifts, and the relationships among the solution concentration, temperature, and refractive index were established for different solutions [
28]. In addition, high-resolution refractive index measurement was achieved by using differential interferometry, based on a white light Mach–Zehnder configuration [
29] and by employing ultra-compact double-slot hybrid plasmonic waveguide as an active sensing arm in the interferometer [
30].
A few inspired research works have reported interpreting the indices of the refraction of samples from the interference patterns of the free-space Mach–Zehnder interferometer. An early experiment was conducted to measure the local indices of the refraction of solutions from the relationship between the relative optical path length difference and the interference fringe bending [
31]. The refractive index measurement error was about
RIU. Recently, some researchers have reported employing a white light Mach–Zehnder interferometer to generate a series of interferograms to determine the refractive indices of materials [
32]. A research work published in the same year has reported measuring the indices of the refraction of scintillation crystals from the width shift of interference patterns based on the Mach–Zehnder interferometer [
33]. Their method is also suitable for measuring the indices of the refraction of liquids in a rectangular cuvette. A change in the direction of the interference stripes was witnessed in the interference patterns for different crystals. In this work, to understand the underlying physics in fringe spacing and direction changes, we propose a novel method to resolve the indices of the refraction of transparent and translucent liquids from the interference patterns of the Mach–Zehnder interferometer.
To verify the proposed method, the NaCl solution is adopted as one of the test examples to measure the refractive indices of NaCl solutions with various concentrations. For the NaCl solutions, large-scale experiments have been carried out to measure the refractive index at different wavelengths and temperatures. A brief review of the relevant works is given in the following sentences. The refractive indices of NaCl solutions was measured using the Abbe refractometer [
34]. The fit curve for the refractive indices of NaCl solutions was
(R
2 = 0.9984) at 588 nm and at 28 °C, with
c representing the concentration of NaCl solutions. The Abbe refractometer was also used to measure the refractive indices of NaCl solutions of different concentrations and temperatures at the wavelength of 589.3 nm [
1]. The empirical expression for the refractive index was
with
c representing the concentration and
T representing the absolute temperature. Refs. [
35,
36] reported the refractive indices of NaCl solutions with concentrations ranging from 0 to 16.6667% in the spectral range of 0.3–1.0 μm at 20 °C and with concentrations from 0 to 4% at 22 °C and at the He-Ne laser wavelength, i.e., 632.8 nm, respectively. In ref. [
28], the relationships between the solution concentration, temperature, and refractive indices for NaCl solutions were established at 632.8 nm, with
, where
T represents the temperature in Celsius. The refractive indices of NaCl solutions was measured for concentrations from 0 to 20% at 632.8 nm [
37]. The fit curve for the refractive index was
(
= 0.9993). The refractive indices of NaCl solutions was measured by using the image correlation coefficient under the He-Ne laser wavelength. The data obtained agreed well with that measured using the Abbe refractors [
38]. In addition, the refractive indices of NaCl solutions was measured for concentrations from 0 to 1.8% at the wavelength of 850 nm [
39], for concentrations from 5.5 to 33.2% at 543 nm and at 20 °C [
40], for concentrations smaller than 0.03% [
18], for concentrations up to 30% at 1550 nm [
6], for concentrations from 0 to 25% at 1550 nm [
41], for various concentrations at the wavelengths of 589 and 636.94 nm at 25 °C [
14], and for very small concentrations at 532 nm [
27].
This article is organized as follows. In
Section 2, we illustrate in detail the derivation process of building up the analytic relations between the index of refraction and the fringe spacings, spatial frequencies, and directions of the interference patterns. In
Section 3, the experimental results are presented and analyzed to verify the proposed method.
Section 4 presents the conclusions of this study.
2. Theory and Methods
The schematic of the Mach–Zehnder interferometer measurement system is presented in
Figure 1. A light beam from a He-Ne laser passing through a linear polarizer P is expanded and collimated using two lens L1 and L2. In this case, the light beam can be approximated by a monochromatic plane wave. Then, it is split into two paths by the beam splitter BS1. One light beam path first goes to mirror M2 and passes through the cuvette and then propagates to the beam splitter BS2. The other light beam path first goes to mirror M1 and is then reflected to the beam splitter BS2 to meet the first light beam. The photodetector CCD is placed after BS2 to capture the interference patterns generated due to the optical path difference between the two paths in the interferometer.
When the sensor surface of the CCD is placed perpendicular to the optical axis, i.e., the
z axis in
Figure 1, the optical path difference between the two paths is a linear function of position coordinates (
x,
y). When the cuvette is removed from the stage, the interference pattern on the CCD is simply from interference between two plane waves in the two paths. This optical path difference is
This optical path difference also contains values from the beam splitters and the mirrors. Since they are highly uniform planar optical elements, the quadratic and higher-order polynomials in the optical path can be safely neglected. When the empty rectangular cuvette is placed on the stage, the interference pattern on the CCD will change due to an additional optical path length introduced in this setup. Considering the good uniformity of the front and back quartz glass plates, the extra optical path difference induced by the empty cuvette is approximated by
Thus, the total optical path length difference between the two paths in the Mach–Zehnder interferometer becomes
a and
b are the structure coefficients that are responsible for the fringe spacings and directions of the interference pattern for the empty cuvette. When the cuvette is filled with liquids, as the front and back plates of the glass cuvette are not perfectly parallel to each other, there is always a very small angle between them due to machining errors, forming a liquid wedge. In fact, we take advantage of this property of the commercial rectangular quartz cuvette. The optical path difference from the liquid wedge is approximately by
Here,
and
are independent of the liquids and are the structure coefficients for the liquid wedge in the cuvette.
n is the index of refraction of the liquid to be measured, and
is the refractive index of air since the system is placed in free space. Therefore, when the cuvette contains liquids, the total optical path difference between the two paths is
By applying the constructive or destructive interference conditions (
or
with
k = 0, ±1, 0, ±2,… the number of wavelength between the two paths), we can obtain the fringe spacings between adjacent bright or dark fringes and the fringe directions of the interference pattern. For the empty cuvette, the fringe spacings in the
x and
y directions are
in which
and
are the spatial frequencies of the interference fringes along the two directions. Thus, when the wavelength of the laser is known, after measuring the fringing spacings or spatial frequencies along the two perpendicular directions, the absolute value of structure coefficients
a and
b in Equation (
3) for the optical path difference of the empty cuvette can be obtained. To determine the sign of these two coefficients, the direction of the interference fringes is required. If the angle of the fringes to the positive direction of the horizontal
x axis is
, the fringe direction is described by the tangent of this angle
Note that the sign of the tangent in Equation (
7) indicates the direction of the interference fringes. When either
a or
b is initially set positive or negative, once the direction of the fringes is known, the sign of the other coefficient is definite. The fringe spacings
and
, the spatial frequencies
and
, and the tangent of the fringe direction
for the interference fringes when the cuvette is filled with liquids can be similarly obtained from Equation (
5), i.e.,
To obtain the two unknown structure coefficients
and
, we can fill in the cuvette with a liquid of a known refractive index
. For example, the index of refraction of water is
= 1.3325 at the room temperature. After measuring the fringe spacings
and
(or the fringe spatial frequencies of
and
) of the interference fringes for the cuvette with water,
and
are obtained
The fringe direction satisfies
and is used to determine the sign in Equation (
9). After determining all of the four structure coefficients
a,
b,
, and
, we can use this setup to measure the refractive index for a liquid from its fringe spacings
and
(or from fringe spatial frequencies
and
) and fringe direction
. The analytic expressions for the liquid refractive indices are
with
. It is clear that the refractive index of a liquid is proportional to the spatial frequencies of the interference patterns. In principle, the above three equations will give us identical results for the index of refraction for a single liquid. However, due to the external noise, the simplification of the plane wave model and the limited resolution of the detectors, there is slight difference between these three values.
To show the influence of the liquid index of refraction on the fringe spacings (or spatial frequencies) and direction of the interference fringes, the simulation results for an illustrating example are presented in
Figure 2 with
a = 0.0002,
b = −0.0002,
= 0.0005,
= 0.0005,
= 1,
= 1.3325, and
n = 1.34, 1.36, and 1.38, respectively. The height and width of the each figure is 3600 × 3600 pixels. The pixel size is 2.4 × 2.4 μm.
In this example, the interference fringes of the empty cuvette in
Figure 2a and the fringes of the water wedge in
Figure 2b are perpendicular to each other. The combination of them for the cuvette filled with water and with the liquids of refractive indices of 1.34, 1.36, and 1.38 are shown in
Figure 2c–f, respectively. The fringe spacings and the corresponding spatial frequencies in both the
x and
y directions and the fringe direction are changing when the cuvette is filled with liquids of different refractive indices. When the direction of the empty cuvette interference fringes are not parallel to that of the liquid wedge, both the change in the fringe spacings and fringe direction can be observed when the index of refraction of the liquid varies. Sometimes, the naked eyes fail to detect this variation, but the highly precise photodetectors with a much larger resolving power can detect it.
It is worth mentioning that the above derivation process is also valid for other interferometers. For example, when we consider the Michelson interferometer, the only difference is that the optical path difference between the two paths is twice of that for the Mach–Zehnder interferometer since the light beam passes through the cuvette twice in the sensing path. When the number of fringes of the empty cuvette and of the liquid wedge are comparable, the changes in the fringe direction are the most obvious. In the case of the two sets fringes for the empty cuvette and the liquid wedge are parallel, only the fringe spacings can be used to determine the indices of refraction for the liquids since the direction of the fringes are fixed. In the future, we will consider a more general case where the plane wave approximation of the optical path fails. In that case, we cannot expect very straight fringes with identical fringe widths and directions in the view field. By introducing a more complex optical path difference model, the possibility to simultaneously measure the local refractive indices of liquids is open.
3. Experimental Results and Discussion
To verify our proposed method of measuring the refractive indices of liquids, we conduct the relevant experiments to measure the indices of refraction of NaCl solutions by using the system in
Figure 1. The atmosphere temperature is 23 ± 3 °C The wavelength of the He-Ne laser is
= 632.8 nm. The experimental system is shown in
Figure 3. For the beam expansion and collimation system, the focal length of the concave lens is −9.8 mm and that of the convex lens is 100 mm. Thus, the magnification of this beam expander is about 100/9.8 ≈ 10 when the distance between the two lenses is around 99.2 mm. The dimension of the CCD is 5488 × 3672 pixels. Each pixel size is 2.4 × 2.4 μm. The quartz glass cuvette is rectangular in shape, and its optical path is around 10 mm. The distance between the front and back plates of the cuvette is also around 10 mm. To make the liquid in the cuvette stable in the small cuvette and to reduce the measurement error, the waiting time for each measurement after changing the liquid is at least 2 h.
To obtain the indices of refraction of the liquids from Equation (
11), the four structure coefficients should be determined in advance. Among them,
a and
b are determined by the fringe spacings, spatial frequencies and direction of the interference pattern for the empty cuvette in
Figure 4a from Equations (
6) and (
7). The interference pattern is closely related to the cuvette position on the stage and the two mirrors. To extract the fringe spacing from the interference pattern, make sure a sufficient number of fringes appearing in the interference area. This is realized by adjusting the two mirrors. In addition, the number of the fringes should be controlled in a certain range; otherwise, the CCD cannot detect the variation in the fringe spacings from the interference patterns of liquids with very close concentrations. In this case, the sensitivity and accuracy of measuring the refractive index with this experimental system is very limited. The spatial frequencies of the interference fringes are obtained from the spectrum by taking the 2D Fourier transform of the corresponding interference patterns. Compared with the fast Fourier transform algorithm, the former Fourier transform method is applied to increase the precision of spatial frequency with the price of an increase in the computation time. After taking the reciprocal of the spatial frequencies, one can obtain the fringe spacings along the
x and
y directions. The advantage of this technique is noise insensitivity. In
Figure 4a, the fringe spacings in
x (horizontal) and
y (vertical) direction are
mm and
mm, and the tangent of the fringe direction is
. Substituting them into Equations (
6) and (
7), the structure coefficients for the empty cuvette are
and
.
To obtain the structure coefficients of the liquid wedge, the interference pattern for the cuvette with a liquid of a known refractive index is required. In
Figure 4b, the interference pattern for the cuvette filled with water is presented. The index of refraction of water was accurately measured at different wavelengths and temperatures [
4,
42,
43]. We take
1.3315 at the considered wavelength and at room temperature in this work. The fringe spacings are
mm and
mm, receptively. The tangent of the fringe direction is
. Inputting them into Equation (
8), we have the two structure coefficients for the liquid wedge
and
. Since
a and
are both positive, the fringe spacing in the
x direction decreases as the index of refraction of the liquid increases.
b and
have different signs, and the value of
becomes smaller than that of
b. Thus, the fringe spacing in the
y direction increases as the refractive index of the liquid increases. Meanwhile, the tangent of the fringe direction is increasing with an increase in the refractive index, i.e., the fringes become more vertical for liquids with larger indices of refraction. Such trends are obvious in
Figure 4.
After gathering all four structure coefficients and filling the cuvette with the liquids to be measured, the CCD was used to capture the corresponding interference patterns. The fringe spacings
and
, spatial frequencies
and
, and direction
are first extracted from the interference patterns; then, Equation (
11) was applied to obtain the indices of refraction of the liquids. The interference patterns for
,
,
,
, and
NaCl solutions are presented in
Figure 5a–e. Here, the concentration refers to the weight fraction between NaCl and the NaCl solution. The refractive indices of these NaCl solutions together with the corresponding fringe spacings
and
, spatial frequencies
and
, and direction
are presented in
Table 1.
and
are the refractive indices obtained from the fringe spacings and spatial frequencies.
means the refractive index from the fringe direction. The second last row gives the average of the above three indices of refraction. In the last row, the data taken from ref. [
38] are presented for comparison.
In
Table 1, there is small variation between the three refractive indices for a NaCl solution. As we pointed out before, the error mainly comes from the external noise, the instability of the laser source, the plane wave approximation, the plane geometry assumption for the planar optical elements, and the limited resolving power of the CCD. Similar to ref. [
33], the measurement accuracy of our method reaches around 10
−3.
To test the stability and repeatability of the measured refractive indices of NaCl solutions, we conduct a new experiment to measure the refractive indices of NaCl solutions by using the same cuvette. The atmosphere temperature is 27 ± 3 °C. The resultant interference patterns for the empty cuvette, cuvette filled with water, and cuvette filled with NaCl solutions with concentrations of 1, 5, 10, 15, 20, and 25% are shown in
Figure 6. The interference pattern for the empty cuvette is very different from that in
Figure 4, resulting in the interference patterns for water and NaCl solutions differing dramatically with those in
Figure 5. From the fringe spacings or spatial frequencies of
Figure 6a,b, the four structure coefficients can be obtained according to Equations (
6)–(
9). The fringe spacings and spatial frequencies of the interference pattern for the empty cuvette and water are shown in
Table 2. They are used to calculate the structure coefficients for the empty cuvette and liquid wedge. The relevant data are also given in
Table 2.
In
Table 2, it is clear that the four structure coefficients are different from that of the previous experiments. The slight difference between the structure coefficients of the liquid wedge indicates a slight change in the position of the cuvette for the new experiment. The fringe spacings, spatial frequencies, and directions of the interference patterns for NaCl solutions with different concentrations are given in
Table 3. The refractive indices for these NaCl solutions are obtained by applying Equation (
11). They are also listed in this table.
To make the relation between the refractive indices of NaCl solutions and the concentration more explicit and to better compare the data obtained by the two experiments, we draw the concentration and refractive index data points in
Figure 7. Besides the linear fitting curves for the refractive index, the least squares method are plotted for the two experiments, and the function forms are
and
with
c the concentration of the NaCl solution. The refractive indices of NaCl solutions obtained by the two individual experiments are very close to each other, confirming the stability and repeatability of the proposed method. To compare our results with that in ref. [
38], the data taken from this reference are also plotted in
Figure 7, with the fit curve for the refractive index reaching
. There is decent agreement between our experiments and the experiment in ref. [
38]. The slight difference between the data of the experiments mainly comes from the difference in the atmosphere temperatures when the experiments were carried out. In principle, the refractive indices of both water and NaCl solutions at higher temperatures is smaller than the refractive indices at lower temperatures. Thus, it is apparent that the actuary of the proposed method is slightly influenced by the accuracy of the refractive index of water. The influence of the atmosphere temperature on the refractive indices of NaCl solutions is reflected in the deviation between the fitting curve and the measured data in
Figure 7. When the concentration of the liquids is of concern, the relative trend of the refractive index matters, and the influence of the slight variation in the water index of refraction is reduced.
Here, we consider the sensitivity and resolution of our method. It is closely related to the structure of the experiment system, including the position of the cuvette and the positions of the two mirrors. The refractive indices of liquids are proportional to the spatial frequencies of the interference fringes according to Equation (
11). This is confirmed in
Figure 8, where the spatial frequencies in the
x and
y directions as a function of the refractive index (the data are taken from
Table 1 and
Table 3) are plotted for the two experiments. The sensitivities are defined as the absolute values of the slopes of these curves. They are 1752/RIU in the
x direction and 5454/RIU in the
y direction for the first experiment, and 2694/RIU in the
x direction and 5589/RIU in the
y direction for the second experiment, respectively. Thus, the sensitivity of the latter experiment is improved. The resolution of the refractive index is defined as the ratio between the standard deviation of the spatial frequency and the sensitivity. After fitting the four curves in
Figure 8, we have the root mean square error of 1.851 and 0.921 for
and
for the first experiment and 0.9408 and 0.6731 for
and
for the second experiment. Therefore, the resolutions for these four measurements are 0.0010, 0.00016, 0.00034, and 0.0012 RIU, respectively. The resolutions for these two experiments are not very high. The proposed method offers a new perspective to measure the refractive indices of liquids from the fringe spacings, spatial frequencies, and directions of interference patterns. The advantage of this method is that the sensitivity and resolution can be adjusted and controlled. In addition, the degree of freedom on adjusting the two mirrors in the interferometer ensures a large measurement range of the refractive indices of liquids.
Finally, to show the capacity of applying our method to measure the refractive indices of translucent liquids, we conducted an experiment on measuring the refractive index of Coca Cola by using the experimental system in
Figure 3 at the atmosphere temperature of 28 ± 2 °C. The interference patterns for the empty cuvette, the cuvette filled with water, and the cuvette filled with Coca Cola are presented in
Figure 9. To apply Equation (
11) to obtain the refractive index of Coca Cola, we shall first calculate the structure coefficients of the empty cuvette and liquid wedge from the fringe spacings and spatial frequencies of
Figure 9a,b. The corresponding results are given in
Table 4. Since we adopted the same cuvette as the above two experiments, the structure coefficients for the liquid wedge in this experiment were very close to the results in
Table 2. The slight difference between them comes from the slight difference in the positions of the cuvette between the two experiments. The fringe spacings, spatial frequencies, and directions of the interference pattern for Coca Cola in
Figure 9c are
mm,
mm,
= 1348.4,
, and
, respectively. Thus, inputting the relevant data into Equation (
11), we obtain the refractive indices of Coca Cola, i.e.,
,
, and
. The average of these three refractive indices is
, very close to 1.3472 in ref. [
44], where the refractive index of Coca Cola was measured by the minimum deviation method at the wavelength of 650 nm and at room temperature, and close to 1.3487 and 1.3489 in ref. [
45], where the refractive index of Coca Cola was measured using the Abbe refractometer and the designed apparatus at the wavelength of 532 nm and at the lab temperature of 20 ± 0.5 °C. It is the difference in the measurement wavelength and temperature that leads to the difference between our results and the results in the above-mentioned studies. Therefore, our method is applicable in measuring the real parts of the refractive indices of translucent liquids. It is beyond the capacity of our method to obtain the imaginary part of the refractive index of a translucent liquid to reflect the absorption of the liquid from the spacings, spatial frequencies, and directions of the interference fringes. In addition, for liquids with high turbidity, the applicability of our method is greatly reduced since no obvious interference fringes are generated due to the large scattering and absorption of light by the liquids.