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

Phase sensitive optical rotation measurement using the common-path heterodyne interferometry and a half-wave plate at a specific azimuth angle

Open Access Open Access

Abstract

We proposed a new method for small optical rotation measurement. The method is based on the use of a half-wave plate and the high-stability common-path heterodyne interferometry. When the azimuth angle is at 22.5° of a half-wave plate, the phase has a distinct change caused by the small polarization rotation of the test beam. The optical rotation can be obtained from the relationship between the phase and the azimuth angle of the wave plate. The resolution of polarization rotation measurement can achieve 1.6 × 10−5 degree/mm and the detection sensitivity on circular birefringence of glucose-water solution can be up to δ|nr − nl | = 5.6 × 10−11.

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

1. Introduction

In the early 19th century, Biot & Fresnel discovered that when linearly polarized light passed through an optically active medium, the light would maintain the initial polarization state and start to rotate. Optically active medium can be divided into left-handed and right-handed mediums [1]. The polarization would gradually rotate to the left as the thickness of the left-handed medium increases when the linearly polarized light enters a left-handed medium. On the other hand, the polarization would rotate to the right as linearly polarized light enters a right-handed medium. This kind of medium is ubiquitous in nature, especially in biology, where there are many molecules with optical activity. For example, the protein composed of amino acid is a kind of biological molecule with optical activity. Feng et al. [2] used the Mach-Zehnder interferometer to measure the optical rotation of the glucose concentration. The rotation angle was directly converted to the phase between the two test signals. Because the rotation angle was small for most solutions, the resolution was not high. The stability of the light intensity and the non-common optical path structure could reduce the system accuracy. Lin et al. [3] proposed an optical heterodyne polarimeter, using the linearly polarized light to reflect on the interface of a chiral liquid. When its incident angle was close to the critical angle, the coupling coefficient of the reflected light would be amplified which could be used to measure up to an accuracy of $2 \times {10^{ - 8}}$ for the chiral parameter and $1 \times {10^{ - 5}}$ for the average refractive index of a chiral liquid, respectively. Chou [4] presented an optical heterodyne interferometer using a Zeeman laser to measure the optical rotation in a scattered chiral medium solution whose scattering particle size was 1.072 $\mu m$ in diameter and the volume concentration was 4.6%. The use of a circularly polarized photon pair (CPPP) from the two orthogonal circularly polarized light waves was able to reduce the scattering effect and improve the sensitivity on optical rotation measurement.

Chou et al. [5] reused the Zeeman laser as the light source of the polarized light common-path optical heterodyne interferometer to measure the quartz Cornu depolarizer (QCD) optically active object. Its balanced detection circuit plus a differential amplifier could greatly reduce most of the noise from laser system and environment, and then obtain accurate and highly sensitive optical rotation measurement results. The detection sensitivity of the optical activity of a quartz crystal was measured as $8.5 \times {10^{ - 10}}$. Steinbacher et al. [6] presented a setup combining common-path optical heterodyne interferometry and an accumulative technique for femtosecond laser spectroscopy to measure small changes in optical rotatory power. It could be used for the photoreaction of chiral sulfoxides, and simultaneously recorded linear absorption spectra. The presented polarimeter had an experimentally determined angular resolution of 0.10 mdeg, a measurement time of only one second, an interaction length of only 250 µm, and a probe beam size (aperture) of only 30 µm. Kumagai et al. [7] proposed a Sagnac interferometer for the glucose concentration measurement. The angle of optical rotation was measured by detecting the phase difference between clockwise and counterclockwise circular polarized light that propagated in a sensing loop with a polarization-maintaining optical fiber. The angular resolution of optical rotation is about $5.4 \times {10^{ - 4}}$ degree. Liu et al. [8] proposed a low-coherence interferometer for the optical rotation measurement in phase-sensitive detection. They used the polarization maintaining fiber (PMF) and polarization optic system to create the left- and right-handed circular states. These two states experienced the same optical path in the reference arm and test arm of the interferometer. A Wollaston prism separated two interferences, the optical rotation can be achieved by the phase difference between these two signals.

Cao et al. [9] proposed an optical rotation measurement system based on the centroid algorithm. The beam of the semiconductor laser was converted into linearly polarized passes through the polarizer (Glan-Taylor prism). Then, the polarized light was received by a photodiode after transmitting through the sample tube and analyzer (Glan-Taylor prism) fixed on a step-motor rotating stage. With the photodiode, the light signals were transformed into electrical signals, taking advantage of processing light signals with or without sample to obtain the optical rotation angle through the step difference between two centroids. After calibrating with several standard quartz tubes and measuring glucose solution with different concentration, the relative error and precision of the system were determined to be 0.4% and 0.004$^\circ $. Although this method was simple, the accuracy was not high. All the above methods did not use water as the test specimen to obtain further confirmation. Goldberg et al. [10] proposed a curve-fitting polarimetric method based on a self-referenced optical rotation polarimeter. They used two high-contrast (1 : ${10^5}$) Glan–Thompson polarizers as the polarizer and analyzer and used an Electro-Optic modulator for phase modulation. The calculated accuracy is $8.1 \times {10^{ - 5}}{\; }$degree, and the calculated precision is $2.22 \times {10^{ - 4}}$ degree.

In this paper we propose a new method for measuring tiny polarization, which is based on the common-path heterodyne interferometry combining with a half-wave plate (HWP) for measurement. When the HWP is rotated to a specific azimuth angle at 22.5$^\circ $, the phase difference changes greatly, and can be used to measure the optically active medium. The analytic polarization rotation is $1.6 \times {10^{ - 5}}\textrm{degree}/\textrm{mm}$, and the detection sensitivity of $\mathrm{\delta }|{{\; }{\textrm{n}_\textrm{r}} - {\textrm{n}_\textrm{i}}} |$ is as high as $5.6 \times {10^{ - 11}}$. Such a high angle resolution is a result that is difficult to achieve by ordinary measuring instruments. In the case of high resolution of polarization rotation measurement, in addition to measuring the concentration of glucose, the optical rotation of different water qualities is also measured to further demonstrate the feasibility of this method.

2. Principles

2.1 Principles and experimental setup

The reason why we use the common-path heterodyne interference as the experimental light source is that the reference light (s-polarized) and the test light (p-polarized) are overlapping in the same path and cancel the phase change with each other so as not to be affected by environmental interference. This can improve the accuracy and reduce the error of the experiment. The structure is simple and can accurately resolve the phase. Not only can it be used to measure the optical rotation of optically active objects, it can also be used to measure long wavelength changes [11,12], small angle measurement [1315], and small displacement [1618]. The experimental structure is shown in Fig. 1, using the common-path heterodyne interferometry combining with a HWP and the test sample uses a vessel with an inner length of 5 cm to hold the solution. The experimental setup of this study is divided into two parts. Part 1 is a heterodyne light source, including a He-Ne laser (λ=632.8 nm), two polarization beam splitters (PBS1, PBS2), two acousto-optic modulators (AOM1, AOM2), two mirrors (M1, M2), and two apertures (Iris1, Iris2). Initially, in order to generate a heterodyne light source, we use PBS1 to divide the laser light into the p- and s- polarizations. The modulation frequencies of AOM1 and AOM2 are 80.00 and 80.01 MHz, respectively. After the s-polarized light is reflected by the mirror (M2) and then incidents on AOM2, the +1 order beam is reflected by PBS2. Since the modulated light will have the -1, 0 and +1 order lights, we use the aperture (Iris1or Iris2) to block the -1 and 0 order lights which are not used in this experiment, and we only take out the +1 order light. The p-polarized light from PBS1 passing through the AOM1 and reflected by the mirror (M1) is combined with the s-polarized light by PBS2. A heterodyne light source with a frequency difference of 10kHz can be obtained. Part 2 is the optical measurement architecture, which includes a beam splitter (BS), two analyzers ($\textrm{A}{\textrm{N}_\textrm{r}},{\; A}{\textrm{N}_\textrm{t}}$) with the same transmission axis at azimuth angle of 45$^\circ $, two photodetectors (${\textrm{D}_\textrm{r}},{\textrm{D}_\textrm{t}}$), a HWP (${\textrm{W}_{\mathrm{\lambda }/2}}$) and the test medium (Sample). The heterodyne light source is divided into the transmitted light and reflected light by BS. The reflected light passes through an analyzer ($\textrm{A}{\textrm{N}_\textrm{r}}$) and received by a photodetector $({\textrm{D}_\textrm{r}}$). The detected interference signal is used as the reference signal ${\textrm{I}_\textrm{r}}$. The transmitted light passes through the test medium (Sample) and the HWP (${\textrm{W}_{\mathrm{\lambda }/2}}$), and then passes through an analyzer ($\textrm{A}{\textrm{N}_\textrm{t}}$). The interference signal detected by the photodetector ${\textrm{D}_\textrm{t}}$ is regarded as the test signal ${\textrm{I}_\textrm{t}}$. Finally, a lock-in amplifier (the Stanford Research Systems, SR830) is used to capture the phase difference between the test signal ${\textrm{I}_\textrm{t}}$ and the reference signal ${\textrm{I}_\textrm{r}}$.

 figure: Fig. 1.

Fig. 1. The experimental structure. Laser: He-Ne Laser, λ=632.8 nm; PBS1, PBS2: Polarization Beam-Splitter; AOM1, AOM2: Acousto-Optic Modulator; M1, M2: Mirror; BS: Beam-Splitter; ${W_{\lambda /2}}$: half-wave plate (HWP); $A{N_r}\; ,\; A{N_t}$: Analyzer; ${D_r}$, ${D_t}$: Photodetector; Lock-in Amplifier: Stanford Research Systems, SR830

Download Full Size | PDF

Because each interference is from combining the partial components of the p- and s-polarizations on the transmission axis of the analyzer and two polarization lights are always on the same optical path, the signal is the common-path heterodyne interference signal. It is stable enough without concerning the environment disturbance. In addition, the test signal ${\textrm{I}_\textrm{t}}\; $and reference signal ${\textrm{I}_\textrm{r}}\; $are from the same heterodyne light source, the same air disturbance in the light source (Part 1) can cancel each other.

2.2 Formula derivation

This section mainly describes the derived formula. The optical measurement structure is shown in Fig. 2. In the process of derivation, the best measurement conditions are found in the relationship curve between the phase difference and the azimuth angle of the HWP. From Fig. 2, ${\textrm{E}_1}$ represents the optical field intensity reflected by BS and the $\textrm{Jon}{\textrm{e}^{{\prime}}}\textrm{s}$ Matrix of the reference signal can be written as

$${{\mathbf E}_{\mathbf r}} = {\mathbf A}{{\mathbf N}_{\mathbf r}}({{{45}^{\mathbf o}}} ){{\mathbf E}_1} = \frac{1}{2}\left[ {\begin{array}{cc} 1&1\\ 1&1 \end{array}} \right]\left[ {\begin{array}{c} {{\textrm{A}_\textrm{x}}{\textrm{e}^{\textrm{j}({{\mathrm{\omega }_1}\textrm{t}} )}}}\\ {{\textrm{A}_\textrm{y}}{\textrm{e}^{\textrm{j}({{\mathrm{\omega }_2}\textrm{t}} )}}} \end{array}} \right] = \frac{1}{2}\left[ {\begin{array}{c} {{\textrm{A}_\textrm{x}}{\textrm{e}^{\textrm{j}({{\mathrm{\omega }_1}\textrm{t}} )}} + {\textrm{A}_\textrm{y}}{\textrm{e}^{\textrm{j}({{\mathrm{\omega }_2}\textrm{t}} )}}}\\ {{\textrm{A}_\textrm{x}}{\textrm{e}^{\textrm{j}({{\mathrm{\omega }_1}\textrm{t}} )}} + {\textrm{A}_\textrm{y}}{\textrm{e}^{\textrm{j}({{\mathrm{\omega }_2}\textrm{t}} )}}} \end{array}} \right],$$
and
$${\; }{\textrm{I}_\textrm{r}} = {|{{\textrm{E}_\textrm{r}}} |^2} = \frac{1}{2}\textrm{A}_\textrm{x}^2 + \frac{1}{2}\textrm{A}_\textrm{y}^2 + {\textrm{A}_\textrm{x}}{\textrm{A}_\textrm{y}}\cos ({\mathrm{\omega }\textrm{t}} )$$
where ${\textrm{A}_\textrm{x}}$ and ${\textrm{A}_\textrm{y}}$ are the electric field amplitudes of ${\textrm{E}_1}$ in the x and y directions respectively, $\mathrm{\omega } = {\mathrm{\omega }_2} - {\mathrm{\omega }_1}$.

 figure: Fig. 2.

Fig. 2. Schematic diagram of experimental formula derivation.

Download Full Size | PDF

S represents the test specimen, and ${\textrm{E}_2}$ is the intensity of electric field of the transmitted light from the BS. Assuming that the phase retardation of the HWP is Γ, the $\textrm{Jon}{\textrm{e}^{{\prime}}}\textrm{s}$ Matrix of the test light can be written as

$$\begin{aligned} {{\mathbf E}_{\mathbf t}} &= {\mathbf A}{{\mathbf N}_{\mathbf t}}({{{45}^{\mathbf o}}} )\; \; {{\mathbf W}_{\frac{\mathrm{\lambda }}{2}}}({\mathbf \theta } ){\mathbf S}{{\mathbf E}_2} \\ & = \frac{1}{2}\left[ {\begin{array}{cc} 1&1\\ 1&1 \end{array}} \right]\left[ {\begin{array}{cc} {\textrm{co}{\textrm{s}^2}\mathrm{\theta }{\textrm{e}^{ - \textrm{j}\frac{\mathrm{\Gamma }}{2}}} + \textrm{si}{\textrm{n}^2}\mathrm{\theta }{\textrm{e}^{\textrm{j}\frac{\mathrm{\Gamma }}{2}}}}&{\textrm{cos}\mathrm{\theta} \textrm{sin}\mathrm{\theta }({\textrm{e}^{ - \textrm{j}\frac{\mathrm{\Gamma }}{2}}} - {\textrm{e}^{\textrm{j}\frac{\mathrm{\Gamma }}{2}}})}\\ {\textrm{cos}\mathrm{\theta} \textrm{sin}\mathrm{\theta }({\textrm{e}^{ - \textrm{j}\frac{\mathrm{\Gamma }}{2}}} - {\textrm{e}^{\textrm{j}\frac{\mathrm{\Gamma }}{2}}})}&{\textrm{si}{\textrm{n}^2}\mathrm{\theta }{\textrm{e}^{ - \textrm{j}\frac{\mathrm{\Gamma }}{2}}} + \textrm{co}{\textrm{s}^2}\mathrm{\theta }{\textrm{e}^{\textrm{j}\frac{\mathrm{\Gamma }}{2}}}} \end{array}} \right]\left[ {\begin{array}{cc} {{\textrm{e}^{\textrm{j}{\mathrm{\delta }_1}}}}&0\\ 0&{{\textrm{e}^{\textrm{j}{\mathrm{\delta }_2}}}} \end{array}} \right]\left[ {\begin{array}{c} {{\textrm{A}_\textrm{x}}{\textrm{e}^{\textrm{j}\left( {{\mathrm{\omega }_1}\textrm{t}} \right)}}}\\ {{\textrm{A}_\textrm{y}}{\textrm{e}^{\textrm{j}\left( {{\mathrm{\omega }_2}\textrm{t} + {\phi _0}} \right)}}} \end{array}} \right] \\ & = \frac{1}{2}\left[ {\begin{array}{c} {\textrm{{\rm M}}{\textrm{A}_\textrm{x}}{\textrm{e}^{\textrm{j}({{\mathrm{\omega }_1}\textrm{t} + {\mathrm{\delta }_1}} )}} + \textrm{{\rm N}}{\textrm{A}_\textrm{y}}{\textrm{e}^{\textrm{j}({{\mathrm{\omega }_2}\textrm{t} + {\phi_0} + {\mathrm{\delta }_2}} )}}}\\ {\textrm{{\rm M}}{\textrm{A}_\textrm{x}}{\textrm{e}^{\textrm{j}({{\mathrm{\omega }_1}\textrm{t} + {\mathrm{\delta }_1}} )}} + \textrm{{\rm N}}{\textrm{A}_\textrm{y}}{\textrm{e}^{\textrm{j}({{\mathrm{\omega }_2}\textrm{t} + {\phi_0} + {\mathrm{\delta }_2}} )}}} \end{array}} \right], \end{aligned}$$
where ${\phi _0}$ is the phase difference between ${\textrm{E}_2}$ and ${\textrm{E}_1}$, $\mathrm{\delta }$ is the phase retardation difference of the test specimen, $\mathrm{\delta } = {\mathrm{\delta }_2} - {\mathrm{\delta }_1}$ and $\mathrm{\theta }$ is the azimuth angle of the wave plate. Let
$$\textrm{{\rm M}} = \cos \frac{\mathrm{\Gamma }}{2} - \textrm{j}(\cos 2\mathrm{\theta } + \sin 2\mathrm{\theta })\sin \frac{\mathrm{\Gamma }}{2}$$
and
$$\textrm{{\rm N}} = \cos \frac{\mathrm{\Gamma }}{2} + \textrm{j}(\cos 2\mathrm{\theta } - \sin 2\mathrm{\theta })\sin \frac{\mathrm{\Gamma }}{2}$$
The test signal can be written as
$$ \begin{aligned} {\textrm{{\rm I}}_\textrm{t}} & = {|{{\textrm{E}_\textrm{t}}} |^2} \\ & = \frac{1}{2}\left\{ {co{s^2}\left( {\frac{\mathrm{\Gamma }}{2}} \right)({\textrm{A}_\textrm{x}^2 + \textrm{A}_\textrm{y}^2} )+ si{n^2}\left( {\frac{\mathrm{\Gamma }}{2}} \right)[{{({{\; }\cos 2\mathrm{\theta } + \textrm{sin}2\mathrm{\theta }} )}^2}\textrm{A}_\textrm{x}^2 + ({\cos 2\mathrm{\theta } + \textrm{sin}2\mathrm{\theta }{)^2}\textrm{A}_\textrm{y}^2} ]+ 2\mathrm{\gamma }\cos ({\textrm{a} - \textrm{b}} )} \right\} \end{aligned}$$
In this case,
$$\mathrm{\gamma } = \sqrt {{{\left( {co{s^2}\left( {\frac{\mathrm{\Gamma }}{2}} \right) - si{n^2}\left( {\frac{\mathrm{\Gamma }}{2}} \right)cos4\theta } \right)}^2} + {{({\sin \mathrm{\Gamma }\cos 2\mathrm{\theta }} )}^2}} {\textrm{A}_\textrm{x}}{\textrm{A}_\textrm{y}}$$
$$\textrm{a} = \mathrm{\omega }\textrm{t} - {\phi _0} + \mathrm{\delta }, $$
and
$$\textrm{b} = {\tan ^{ - 1}}(\frac{{\sin \mathrm{\Gamma }\cos 2\mathrm{\theta }}}{{co{s^2}\left( {\frac{\mathrm{\Gamma }}{2}} \right) - si{n^2}\left( {\frac{\mathrm{\Gamma }}{2}} \right)cos4\theta }}).$$

2.3 MATLAB simulation

According to the formula presented in Section 2-2, we can compare the phase difference between ${\textrm{I}_\textrm{t}}$ and ${\textrm{I}_\textrm{r}}$. The value is between 0 and 2π. For the convenience of explanation, we will convert 0 ∼ 2π to -180°∼180°. And the function of phase difference $\phi$ can be written as

$$\phi = {\phi _{\textrm{t} - }}{\phi _\textrm{r}} ={-} {\phi _0} + {\; }\mathrm{\delta } + {\tan ^{ - 1}}(\frac{\textrm{B}}{\textrm{A}}).$$
A and B can be expressed by
$$\textrm{A} = co{s^2}\left( {\frac{\mathrm{\Gamma }}{2}} \right) - si{n^2}\left( {\frac{\mathrm{\Gamma }}{2}} \right)cos4\theta $$
and
$$\textrm{B} = \sin \mathrm{\Gamma }\cos 2\mathrm{\theta }$$
respectively. Initially, we assume that there is no phase difference between ${\textrm{E}_1}$ and ${\textrm{E}_2}$, set ${\phi _0}$ to zero, and assume that the phase retardation δ of the test specimen is equal to zero. After setting each parameter, we use MATLAB software to perform the simulation, so that we can get the relationship curves of the phase difference $\phi $ versus the azimuth angle $\mathrm{\theta }$ for different phase retardations Γ of the wave plates.
  • (1) For $0^\circ \le |\mathrm{\Gamma } |\le 90^\circ $,
    $$\phi = {\tan ^{ - 1}}(\frac{\textrm{B}}{\textrm{A}})$$
  • (2) For $90^\circ \le |\mathrm{\Gamma } |\le 180^\circ $,
    • (i)
      $$\; \textrm{A} \ge 0, \phi = {\tan ^{ - 1}}(\frac{\textrm{B}}{\textrm{A}});$$
    • (ii)
      $$\textrm{A} < 0\; \textrm{and}\; \textrm{B} \ge 0, \phi = \mathrm{\pi } + {\tan ^{ - 1}}(\frac{\textrm{B}}{\textrm{A}});$$
    • (iii)
      $$\textrm{A} < 0\; \textrm{and}\; \textrm{B} < 0, \phi ={-} \mathrm{\pi } + {\tan ^{ - 1}}(\frac{\textrm{B}}{\textrm{A}}).$$

In Eqs. (13) to (16), the range of the arctangent function is from -π/2 to π/2. In order to calculate the phase difference conveniently, we define the phase retardation nearly equal to ${\pm} \mathrm{\pi }$ when the arctangent function has two variables. From the above relationship, we use MATLAB to draw the relationship between the phase difference $\phi $ and the azimuth angle θ. If the wave plate used is a nearly half-wave plate, we let the phase retardation Γ be a positive (175°∼179°), the result of the simulated phase difference is shown in Fig. 3(a) The largest sensitive angles are at $22.5^\circ ,{\; }67.5^\circ ,{\; }112.5^\circ ,{\; }157.5^\circ ,{\; \; }202.5^\circ ,{\; }247.5^\circ ,{\; }292.5^\circ ,{\; }$and $337.5^\circ $. Finally, using the result of Γ as a positive value in Fig. 3(a), when θ=22.5$^\circ $, we take out the phase $\phi $ versus the Γ and draw it as shown in Fig. 3(b). It reveals that the phase shift caused by the slight change of phase retardation is smaller than the change of the azimuth angle.

 figure: Fig. 3.

Fig. 3. (a). When Γ is positive, the relationship between the phase difference and the azimuth angle θ; (b) When $\mathrm{\Gamma }$ is a positive and the azimuth angle θ=22.5$^\circ $, the relationship between the phase difference $\phi {\; }$and $\mathrm{\Gamma }$.

Download Full Size | PDF

Taking out the relationship diagram between $\mathrm{\Gamma } = 175^\circ{\sim} 179^\circ $ and its azimuth angle during 22.5°±0.2 degrees, the phase difference as shown in Fig. 4(a), and the slope diagram (angular sensitivity) is shown in Fig. 4(b). The better linearity range is when the phase retardation is less than $177^\circ $.

For the same reason, if Γ is negative (-179°∼-175°), the relationship between the simulated phase difference $\phi $ and the azimuth angle θ is plotted as shown in Fig. 5(a). The largest sensitive angles are the same at $22.5^\circ ,{\; }67.5^\circ ,{\; }112.5^\circ ,{\; }157.5^\circ ,{\; }202.5^\circ ,{\; }247.5^\circ $, $292.5^\circ {\; }$, and $337.5^\circ $. Then, we use the result of Γ as a negative value in Fig. 5(a), when θ=22.5$^\circ $, taking out the phase $\phi $ versus Γ, and draw it as Fig. 5(b). It also reveals that the phase change caused by the retardation change is very small at this azimuth angle. It is almost negligible if compared with the phase change caused by the change of $\mathrm{\theta }$. It can be found from Fig. 3(b) and Fig. 5(b) that at the azimuth angle θ=22.5$^\circ $, the maximum phase change is only about 1.4$^\circ $.

 figure: Fig. 4.

Fig. 4. (a). When $\Gamma = -175^\circ \sim -179^\circ $, the relationship between $\theta = 22.5^\circ{\pm} 0.2^\circ $ and the phase difference $\phi$; (b). The relationship between $\Gamma = -175^\circ \sim -179^\circ $ and the slope (average angular sensitivity)

Download Full Size | PDF

 figure: Fig. 5.

Fig. 5. (a). When Γ is negative, the relationship between the phase difference $\phi $ and the azimuth angle θ. (b) When Γ is negative and the azimuth angle θ=22.5$^\circ $, the relationship between the phase difference $\phi {\; }$and Γ

Download Full Size | PDF

Taking out the relationship diagram between $\mathrm{\Gamma} = -175^\circ \sim -179^\circ$ and its azimuth angle during $22.5^\circ{\pm} 0.2^\circ \; $, the phase difference is shown in Fig. 6(a), and the slope diagram is shown in Fig. 6(b). The better linearity range is when the absolute phase retardation is less than $177^\circ $.

 figure: Fig. 6.

Fig. 6. (a). When $\Gamma = -175^\circ \sim -179^\circ $, the relationship between $\theta = 22.5^\circ{\pm} 0.2^\circ $ and the phase difference $\phi $; (b) The relationship between $\Gamma = -175^\circ \sim -179^\circ $ and the average slope (average angular sensitivity)

Download Full Size | PDF

3. Experimental results and discussion

3.1 Phase retardation measurement results of a half-wave plate

For ensuring the phase retardation of the HWP in this experiment. The experimental results without any test specimen are shown in Fig. 7(a). The solid line is the MATLAB simulation result of Γ=+176.78°. The dotted line is the result of the experiment. From Eqs. (10) ∼ (12), if $\theta = 0^\circ ,$ then $\phi = \mathrm{\Gamma }$, we can prove that the phase retardation Γ of the half-wave plate is +176.78$^\circ $. When the azimuth angle of the half wave plate is equal to 22.5$^\circ $, we assume ±0.2$^\circ $ for the azimuth angle tolerance to find the average slope of this line segment, as shown in Fig. 7(b). The slope of the phase difference is -98.2 (degree/degree) in the range of $\mathrm{\theta } = 22.3^\circ {\; }\sim {\; }22.7^\circ $.

 figure: Fig. 7.

Fig. 7. (a). The experimental results of a half-wave plate; (b) Graph with Γ=176.78$^\circ $ and azimuth angle $\theta $ from 22.25$^\circ $ to 22.7$^\circ $.

Download Full Size | PDF

This phase change caused by the azimuth rotation is much larger than that of the phase retardation variation. Even if the medium has a small phase retardation variation, the phase shift caused by the phase retardation is very small and can be ignored. The polarization rotation and azimuth rotation are relative motions. In the absence of azimuth rotation, the phenomenon of phase shift appears, which means that the phase shift is caused by the polarization rotation and the polarization rotation is the negative azimuth rotation ($\alpha ={-} \Delta \theta $). So we directly calculate the polarization rotation by using the equation $\alpha = \frac{{\Delta \phi }}{{98.2}}$.

3.2 When the azimuth angle of the half-wave plate is equal to 0$^\circ $, the effect of phase retardation of the glucose

We use 50mm-long glucose solution as the test specimen in this experiment. The concentration is from 100 mg/dL to 1000 mg/dL, and the measurement is performed with an increase of 100 mg/dL each time. When the azimuth of the HWP is equal to 0$^\circ $, the measurement result is shown in Fig. 8. The abscissa of Fig. 8 is the concentration and its unit is milligrams per deciliter (mg/dL), which is the unit of measurement used in the United States for blood glucose concentration. The ordinate is the total phase retardation including the retardations of the HWP and glucose solutions and its unit is degree.

 figure: Fig. 8.

Fig. 8. The total phase retardation measurement result graph at $\mathrm{\theta } = 0^\circ $

Download Full Size | PDF

The up and down jump of the phase difference indicates the change of the phase retardation of the solution, and there is only a vibration in the phase retardation of ±0.2$^\circ $. The average total phase retardation is $176.78^\circ .{\; }$Because its distribution is a jump up and down, there is no correlation with the concentration of glucose. Therefore, the decreasing (or increasing) phase retardation is only related to the fast and slow axes of the birefringence of the glucose solution. In addition, because the fast or slow axis is random, the total phase retardation will slightly jump up and down regardless of the concentration. From Fig. 8, we can confirm once again that the phase retardation is equal to $176.78^\circ $. Because its phase slope near at $\theta = 22.5^\circ $ is $ - 98.2^\circ $ (degree/degree) in Fig. 7(b), the phase change caused by the slight variation of phase retardation is much smaller than the phase change caused by the polarization rotation. However, that does not affect the measurement results of polarization rotation.

3.3 When the half-wave plate rotates to θ=22.5°, the specific rotation number of different concentrations of glucose and fructose

When the HWP is rotated to θ=22.5°, the temperature $T = 25\circ{C}$, the phase shift ($\Delta \phi = \phi - {\phi _0}$, where ${\phi _0} = 91.7^\circ $ is the initial phase) measurement results of high and low concentration glucose solutions are shown in Figs. 9(a) and 9(b), and the results of fructose solutions are shown in Fig. 9(c) and Fig. 9(d), respectively. When the concentration increases, the phase changes will increase, too. This phase variable is only caused by the optical activity of the solutions. The greater polarization rotation angle, the greater phase change will be measured.

 figure: Fig. 9.

Fig. 9. The phase shift ($\Delta \phi $) measurement results of the different concentrations in 50mm-long transparent container at $\mathrm{\theta } = 22.5^\circ $ and $T = 25^\circ C$. (a). High-concentration glucose solution and (b) Low-concentration glucose solution; (c) High-concentration fructose solution and (d) Low-concentration fructose solution$.$

Download Full Size | PDF

Because of the polarization rotation and the azimuth rotation are relative motions ($\mathrm{\alpha }\textrm{=} - \Delta \mathrm{\theta }$), and from Figs. 9(a), 9(b) and Fig. 3(b), the phase shift of glucose solution is positive $(\Delta \phi > 0)$. The meaning is that $\Delta \theta = ({\theta - 22.5^\circ } )< 0$, the wave plate is equivalent to clockwise rotation, so the polarization rotation ($\mathrm{\alpha } > 0)$ is counterclockwise and right-handed. In other words, from Figs. 9(c), 9(d), and Fig. 5(a), the phase shift of fructose solution is negative ($\Delta \phi < 0$) and the polarization rotation is left-handed ($\alpha < 0$).

The phase shift of Figs. 9(a)∼9(d) have been subtracted from the initial phase difference by $91.7^\circ $ (without any solution) and then are divided by value of 98.2 (degree/degree) of the phase retardation $\mathrm{\Gamma } = 176.78^\circ $, the results of optical rotation (or polarization rotation) measurements for the 50 mm long high- and low-concentration glucose solutions are shown in Figs. 10(a) and 10(b), respectively. And the optical rotation measurement results of high- and low-concentration fructose solutions are shown in Figs. 10(c) and 10(d), respectively.

 figure: Fig. 10.

Fig. 10. The polarization rotation measurement results for (a) high-concentration glucose solutions, (b) low-concentration glucose solutions, (c) high-concentration fructose solutions, (d) low-concentration fructose solutions. The solid line is the polarization rotation $\alpha $ and the dotted line is the specific rotation number $[\alpha ]$

Download Full Size | PDF

The solid line is the optical rotation $\alpha $ per millimeter and the dotted line is the specific rotation number $[\alpha ]$. Based on the specific rotation definition [19],

$$[\alpha ]= \frac{\alpha }{{l \times c}},\; \left( {unit:\frac{{degree}}{{\frac{g}{{mL}} \times dm}}} \right),$$
where $\alpha $ is the optical rotation, l is the solution length (the unit is $dm = 0.1m$), c is the concentration (the unit is gram per milliliter). The specific rotation of the glucose and fructose solutions are near to 47${\pm} 1.5$ and -83.5${\pm} 2.5$ ($deg \times mL\ast {g^{ - 1}}d{m^{ - 1}}$), respectively for the wavelength of 632.8 nm. Thus, the experimental error is about ${\pm} $3%. Although we use the average of all values as the specific rotation, the determined error is all values minus the average. We know the specific rotation within this concentration range and also the associated error. In addition, the specific rotation may not be a fixed value as it may change with the change of concentration.

In this paper, we used the average slope of phase curve to evaluate the optical rotation for convenience. In fact, the phase difference diagram $\phi \; $is a curve and not a straight line. Thus the nonlinear error exists and affects the calculations of azimuth angle $\theta $ and optical rotation $\alpha \; $. $\; $The theoretical calculation of the phase error is shown in Fig. 11. The maximum phase error is ${\pm} 0.26^\circ $ within the measurement range of $ - 0.2^\circ \le \alpha \le 0.2^\circ $. It is larger than the resolution of phase measurement. Thus, the accuracy of the polarization rotation is only dependent on the nonlinear error. Therefore, $\Delta {\alpha _{max}}/\textrm{mm} = \frac{{0.26^\circ }}{{98.2 \times 50}} = 5.3 \times {10^{ - 5}}{\; }$degree/mm. For the resolution of phase measurement, the standard deviation is $0.08^\circ \; $from our past experience [17], so the standard deviation of $\alpha $ per millimeter is ${\pm} 1.6 \times {10^{ - 5}}$ degree/mm and the detection sensitivity on circular birefringence of glucose-water solution up to $\delta |{{n_r} - {n_l}} |= 5.6 \times {10^{ - 11}}$.

 figure: Fig. 11.

Fig. 11. The phase error due to the nonlinearity in the theoretical calculation within the measurement range of $\theta $ (or $ - 0.2^\circ \le \alpha \le 0.2^\circ $).

Download Full Size | PDF

3.4 Optical rotation of different water qualities

When the HWP is rotated to θ=22.5°, $T = 25\circ{C},\; $the container length is 50 mm, the phase shift and optical rotation measurement results of different water are shown in Figs. 12(a) and 12(b), respectively. From Fig. 12(a), the phase shifts after the phase differences subtracting the initial phase in air by 91.7$^\circ $ are negative and then dividing the results by the average slope of -98.2 degree/degree of Γ=17$6.78^\circ $ is $\Delta \theta $, and $\alpha ={-} \Delta \theta .$ The optical rotations per millimeter of different water qualities are obtained as shown in Fig. 12(b). These optical rotations are left-handed and the average rotation is $ - 2.0 \times {10^{ - 4}}{\; }degree/mm$. The distilled and slightly alkaline water have larger phase shifts and optical rotations, and the tap water has the smallest optical rotation. It might be revealed that larger optical rotation has a higher water quality.

 figure: Fig. 12.

Fig. 12. (a). The phase shift from air of different water qualities at $\mathrm{\theta } = 22.5^\circ $, where, A: Distilled saltwater, B: Bamboo charcoal water, C: Slightly alkaline water, D: Tap water, E: Reverse osmosis water, F: Desalinated water; (b) The experimental results of the optical rotation of different water qualities

Download Full Size | PDF

4. Conclusions

The small optical rotation measurement based on the common-path heterodyne interferometry and the use of a half-wave plate was proposed in this paper. The measurable range of the polarization rotation is from 0$^\circ $ to ±0.$2^\circ $. Because the phase error can reach ${\pm} 0.26^\circ $, the accuracy of the optical rotation can be up to $5.3 \times {10^{ - 5}}\; \textrm{degree}/\textrm{mm}$. For the phase standard deviation is $0.08^\circ {\; },$ the deviation of $\alpha $ is ${\pm} 1.6 \times {10^{ - 5}}$ degree/mm and the optical rotations of different water qualities are negative (left-handed) and have been successfully measured. From the experimental results, the error of the specific rotation measurement is about ${\pm} 3\%.$ The feasibility of this method is demonstrated. In this paper, we emphasize that the proposed method can measure very tiny optical rotation even for water. The optical rotation may be another important parameter for water.

Funding

National Formosa University (EN2020-09170456074223).

Acknowledgment

M.H. Chiu thanks the MOST in Taiwan for supporting this work.

Disclosures

The authors declare no conflicts of interest.

References

1. L. Velluz, M. Legrand, and M. Grosjean, “Optical Circular Dichroism,” Academic, New York, Ch.1-Ch.3 (1965).

2. C. M. Feng, Y. C. Huang, J. G. Chang, M. Chang, and C. Chou, “A true phase sensitive optical heterodyne polarimeter on glucose concentration measurement,” Opt. Commun. 141(5-6), 314–321 (1997). [CrossRef]  

3. J. Y. Lin and D. C. Su, “A new method for measuring the chiral parameter and the average refractive index of a chiral liquid,” Opt. Commun. 218(4-6), 317–323 (2003). [CrossRef]  

4. C. Chou, W. C. Kuo, T. S. Hsieh, and H. K. Teng, “A phase sensitive optical rotation measurement in a scattered chiral medium using a Zeeman laser,”,” Opt. Commun. 230(4-6), 259–266 (2004). [CrossRef]  

5. C. Chou, W. C. Kuo, T. S. Hsieh, and H. K. Teng, “Optical activity measurement by use of a balanced detector optical heterodyne interferometer,” Appl. Opt. 45(16), 3733–3739 (2006). [CrossRef]  

6. A. Steinbacher, J. Buback, P. Nuernberger, and T. Brixner, “Precise and rapid detection of optical activity for accumulative femtosecond spectroscopy,” Opt. Express 20(11), 11838–11854 (2012). [CrossRef]  

7. T. Kumagai, Y. Tottori, R. Miyata, and H. Kajioka, “Glucose sensor with a Sagnac interference optical system,” Appl. Opt. 53(4), 720–726 (2014). [CrossRef]  

8. C. J. Liu, T. Li, and T. Akkin, “Low-coherence interferometry for phase-sensitive measurement of optical rotation,” Appl. Opt. 57(20), 5893–5898 (2018). [CrossRef]  

9. J. Cao, H. Jia, X. Shen, and S. Jiang, “Research of optical rotation measurement system based on centroid algorithm,” Rev. Sci. Instrum. 87(9), 093108 (2016). [CrossRef]  

10. D. Goldberg and Z. Weissman, “Accurate, nondifferential, curve-fitted, self-referenced optical rotation polarimeter,” Appl. Opt. 58(36), 9815–9822 (2019). [CrossRef]  

11. J. Y. Lee and D. C. Su, “Common-path heterodyne interferometric detection scheme for measuring wavelength shift,” Opt. Commun. 162(1-3), 7–10 (1999). [CrossRef]  

12. K. H. Chen, C. C. Hsu, and D. C. Su, “Measurement of wavelength shift by using surface plasmon resonance heterodyne interferometry,” Opt. Commun. 209(1-3), 167–172 (2002). [CrossRef]  

13. M. H. Chiu and D. C. Su, “Angle measurement using total-internal-reflection heterodyne interferometry,” Opt. Eng. 36(6), 1750–1753 (1997). [CrossRef]  

14. M. H. Chiu, S. F. Wang, and R. S. Chang, “Instrument for measuring small angles by use of multiple total internal reflections in heterodyne interferometry,” Appl. Opt. 43(29), 5438–5442 (2004). [CrossRef]  

15. F. Zhu, J. Tan, and J. Cui, “Common-path design criteria for laser datum based measurement of small angle deviations and laser autocollimation method in compliance with the criteria with high accuracy and stability,” Opt. Express 21(9), 11391–11403 (2013). [CrossRef]  

16. S. J. Liao, S. F. Wang, and M. H. Chiu, “Small-displacement sensing system based on multiple total internal reflection in heterodyne interferometry,” Appl. Opt. 48(13), 2566–2573 (2009). [CrossRef]  

17. M. H. Chiu, B. Y. Shih, C. W. Lai, L. H. Shyu, and T. H. Wu, “Small absolute distance measurement with nanometer resolution using geometrical optics principles and a SPR angular sensor,” Sens. Actuators, A 141(1), 217–223 (2008). [CrossRef]  

18. M. H. Chiu, W. C. Chen, and C. T. Tan, “Small displacement measurements based on an angular-deviation amplifier and interferometric phase detection,” Appl. Opt. 54(10), 2885–2890 (2015). [CrossRef]  

19. M. F. G. Wood, “Combined optical intensity and polarization methodology for analyte concentration determination in simulated optically clear and turbid biological media,” J. Biomed. Opt. 13(4), 044037 (2008). [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 (12)

Fig. 1.
Fig. 1. The experimental structure. Laser: He-Ne Laser, λ=632.8 nm; PBS1, PBS2: Polarization Beam-Splitter; AOM1, AOM2: Acousto-Optic Modulator; M1, M2: Mirror; BS: Beam-Splitter; ${W_{\lambda /2}}$ : half-wave plate (HWP); $A{N_r}\; ,\; A{N_t}$ : Analyzer; ${D_r}$ , ${D_t}$ : Photodetector; Lock-in Amplifier: Stanford Research Systems, SR830
Fig. 2.
Fig. 2. Schematic diagram of experimental formula derivation.
Fig. 3.
Fig. 3. (a). When Γ is positive, the relationship between the phase difference and the azimuth angle θ; (b) When $\mathrm{\Gamma }$ is a positive and the azimuth angle θ=22.5 $^\circ $ , the relationship between the phase difference $\phi {\; }$ and $\mathrm{\Gamma }$ .
Fig. 4.
Fig. 4. (a). When $\Gamma = -175^\circ \sim -179^\circ $ , the relationship between $\theta = 22.5^\circ{\pm} 0.2^\circ $ and the phase difference $\phi$ ; (b). The relationship between $\Gamma = -175^\circ \sim -179^\circ $ and the slope (average angular sensitivity)
Fig. 5.
Fig. 5. (a). When Γ is negative, the relationship between the phase difference $\phi $ and the azimuth angle θ. (b) When Γ is negative and the azimuth angle θ=22.5 $^\circ $ , the relationship between the phase difference $\phi {\; }$ and Γ
Fig. 6.
Fig. 6. (a). When $\Gamma = -175^\circ \sim -179^\circ $ , the relationship between $\theta = 22.5^\circ{\pm} 0.2^\circ $ and the phase difference $\phi $ ; (b) The relationship between $\Gamma = -175^\circ \sim -179^\circ $ and the average slope (average angular sensitivity)
Fig. 7.
Fig. 7. (a). The experimental results of a half-wave plate; (b) Graph with Γ=176.78 $^\circ $ and azimuth angle $\theta $ from 22.25 $^\circ $ to 22.7 $^\circ $ .
Fig. 8.
Fig. 8. The total phase retardation measurement result graph at $\mathrm{\theta } = 0^\circ $
Fig. 9.
Fig. 9. The phase shift ( $\Delta \phi $ ) measurement results of the different concentrations in 50mm-long transparent container at $\mathrm{\theta } = 22.5^\circ $ and $T = 25^\circ C$ . (a). High-concentration glucose solution and (b) Low-concentration glucose solution; (c) High-concentration fructose solution and (d) Low-concentration fructose solution $.$
Fig. 10.
Fig. 10. The polarization rotation measurement results for (a) high-concentration glucose solutions, (b) low-concentration glucose solutions, (c) high-concentration fructose solutions, (d) low-concentration fructose solutions. The solid line is the polarization rotation $\alpha $ and the dotted line is the specific rotation number $[\alpha ]$
Fig. 11.
Fig. 11. The phase error due to the nonlinearity in the theoretical calculation within the measurement range of $\theta $ (or $ - 0.2^\circ \le \alpha \le 0.2^\circ $ ).
Fig. 12.
Fig. 12. (a). The phase shift from air of different water qualities at $\mathrm{\theta } = 22.5^\circ $ , where, A: Distilled saltwater, B: Bamboo charcoal water, C: Slightly alkaline water, D: Tap water, E: Reverse osmosis water, F: Desalinated water; (b) The experimental results of the optical rotation of different water qualities

Equations (17)

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

E r = A N r ( 45 o ) E 1 = 1 2 [ 1 1 1 1 ] [ A x e j ( ω 1 t ) A y e j ( ω 2 t ) ] = 1 2 [ A x e j ( ω 1 t ) + A y e j ( ω 2 t ) A x e j ( ω 1 t ) + A y e j ( ω 2 t ) ] ,
I r = | E r | 2 = 1 2 A x 2 + 1 2 A y 2 + A x A y cos ( ω t )
E t = A N t ( 45 o ) W λ 2 ( θ ) S E 2 = 1 2 [ 1 1 1 1 ] [ co s 2 θ e j Γ 2 + si n 2 θ e j Γ 2 cos θ sin θ ( e j Γ 2 e j Γ 2 ) cos θ sin θ ( e j Γ 2 e j Γ 2 ) si n 2 θ e j Γ 2 + co s 2 θ e j Γ 2 ] [ e j δ 1 0 0 e j δ 2 ] [ A x e j ( ω 1 t ) A y e j ( ω 2 t + ϕ 0 ) ] = 1 2 [ {\rm M} A x e j ( ω 1 t + δ 1 ) + {\rm N} A y e j ( ω 2 t + ϕ 0 + δ 2 ) {\rm M} A x e j ( ω 1 t + δ 1 ) + {\rm N} A y e j ( ω 2 t + ϕ 0 + δ 2 ) ] ,
{\rm M} = cos Γ 2 j ( cos 2 θ + sin 2 θ ) sin Γ 2
{\rm N} = cos Γ 2 + j ( cos 2 θ sin 2 θ ) sin Γ 2
{\rm I} t = | E t | 2 = 1 2 { c o s 2 ( Γ 2 ) ( A x 2 + A y 2 ) + s i n 2 ( Γ 2 ) [ ( cos 2 θ + sin 2 θ ) 2 A x 2 + ( cos 2 θ + sin 2 θ ) 2 A y 2 ] + 2 γ cos ( a b ) }
γ = ( c o s 2 ( Γ 2 ) s i n 2 ( Γ 2 ) c o s 4 θ ) 2 + ( sin Γ cos 2 θ ) 2 A x A y
a = ω t ϕ 0 + δ ,
b = tan 1 ( sin Γ cos 2 θ c o s 2 ( Γ 2 ) s i n 2 ( Γ 2 ) c o s 4 θ ) .
ϕ = ϕ t ϕ r = ϕ 0 + δ + tan 1 ( B A ) .
A = c o s 2 ( Γ 2 ) s i n 2 ( Γ 2 ) c o s 4 θ
B = sin Γ cos 2 θ
ϕ = tan 1 ( B A )
A 0 , ϕ = tan 1 ( B A ) ;
A < 0 and B 0 , ϕ = π + tan 1 ( B A ) ;
A < 0 and B < 0 , ϕ = π + tan 1 ( B A ) .
[ α ] = α l × c , ( u n i t : d e g r e e g m L × d m ) ,
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.