Abstract
In the paper, we explain an automated LabVIEW-controlled setup that enables interferometric measurement of refractive indices in crystalline materials using a laser light source. The setup combines a mechanical system, a microcomputer-controlled gearless drive, a Michelson interferometer, an optical detector, a data acquisition system, and a LabVIEW virtual instrument for an accurate nondestructive determination of the refractive index of given plane-parallel samples.
We explain the concept, implementation, and hardware/software peculiarities of the developed system. Test experiments on different crystals yielded results that are in good agreement with available reference data. The range of potential applications of the proposed setup extends from fundamental optical research to biophotonics instrumentation, where efficient delivery of light is of crucial importance and reliable automated probing tools are needed for optical component characterization.
Introduction
The refractive index, n, is a fundamental characteristic that describes the interaction of electromagnetic waves with a material. The refractive index is a function of wavelength, and the development of devices that are designed to work in the specific range of frequencies requires accurate and nondestructive measurements of n.
Different methods for the measurement of n for transparent materials have been proposed,1–4 but experimental techniques are still being improved to increase accuracy, speed up measurement processes, and extend the spectral range.5–8 In the visible range, the most widely used methods for precise materials characterization are the ellipsometry method,9,10 the method of immersion interferometry, and the method of minimum deflection (method of the prism). 11 Each of these methods has its advantages and disadvantages as well as a corresponding field of application. Most optical samples or elements for technological application are fabricated as rectangular plates, often plane parallel, with a certain thickness in the micrometer to millimeter range. The most suitable method for the refractive index measurement of such plates in the visible spectrum range is the interferometric-turning method,3,12 which also can be extended to the millimeter-wave range. 13 Recently, heightened interest in the investigation of plane-parallel nanoporous membranes in the visible range has been driven by their possible applications in biomedicine14,15 and as active elements for electro-optics and photonics. 16 Thus, new approaches for the characterization of such materials need to be developed and tested.
It is required to have at least 10−3 accuracy of refractive index determination for most optical media to achieve competitive results in new materials characterization.17–19 To reach this level of precision, one needs to eliminate most of the factors possibly influencing the experimental results. The automation and precise computerized control of the measuring procedure is an important step toward high accuracy and reliability of the experiment. Furthermore, an automated system powered with appropriate software can provide the capabilities of real-time control and data processing, significantly reducing characterization time.
Optical setup automation can be done in different ways using specific software and a variety of external hardware to collect and analyze measured data. National Instruments LabVIEW (Austin, TX) is one of the most advanced solutions to date, as it combines powerful tools for instrument control and data acquisition with a considerable set of functions for signal processing and sophisticated calculations.20,21
The present paper discusses the integration of LabVIEW-powered automation techniques, an interferometric setup, and a gearless motor drive solution for the accurate measurement of the refractive index of crystalline materials using a laser light source.
Materials and Methods
Optical Setup
The optical setup, including the necessary automation elements for measuring the refractive index of crystalline materials, is shown schematically in Figure 1 . The sample to be measured is placed in the test arm of the Michelson interferometer.22–24 Light from the laser, 633 nm, is split into two beams by the semitransparent dividing prism. After splitting, one of the beams reflects from the mirror and the other passes through the polarizer; the sample then reflects from the other mirror and goes through the same sample once again. Both beams interfere in the semitransparent dividing prism to form the interference pattern that is focused on the Thorlabs PDA100A2 (Maisons-Laffitte, France) photodetector using the lens. The photodetector is equipped with an adjustable gain amplifier, and the area of the pattern that can projected onto it is about 10 × 4 mm. The pattern is a traveling series of high- and low-intensity regions, so the adjustable slit is utilized to detect only the region corresponding to a maximum or minimum intensity of light. As the sample rotates, the output signal of the photodetector changes in a sinusoidal fashion and is recorded into PC memory via an analog-to-digital converter (ADC)/digital input/output (I/O) module. This module also serves the purpose of controlling the rotation mechanism.

Experimental setup allowing automated measurement of the refractive index in the visible spectrum range. 1, laser; 2, semitransparent dividing prism; 3 and 4, mirrors; 5, polarizer; 6, lens; 7, photodetector; 8, rotation mechanism; 9, ADC and digital I/O module (USB-6009); 10, personal computer; 11, sample.
Figure 2 illustrates the concept of the interference pattern change registration (to measure the interference shift) when rotating the sample by a given angle. The width of the slit on the photodetector, t, is chosen so that only one maximum or minimum of the interference pattern can be detected while the sample is rotated. When the maximum of the interference pattern is in the middle of the slit ( Fig. 2a ), the maximum of the signal intensity corresponding to the maximum value is observed on the photodetector. However, if one continues to rotate the sample, the maximum of the interference pattern will be replaced by the minimum, and at the photodetector, we will observe the minimal signal intensity ( Fig. 2b ).

Change of the photodetector signal intensity: (
Sample Positioning Mechanism
A crucial element of the described setup for refractive index measurements is the precise sample positioning mechanism. We use a gearless drive (GD) based on a synchronous machine with permanent magnets (SMPM) to ensure precise control of the rotation angle. This type of drive is used in robotic manipulators, guidance systems, and so forth, 25 as it provides the advantages of better mechanical rigidity, higher driving power, and increased reliability due to eliminated components such as couplings and bearings.25,26
Figure 3 shows the functional diagram of the system (left) and the actual design of the drive with its control board. Switched stator windings of the motor (M) are realized using a bridge transistor inverter (I; MOSFET type transistor IRF3205) and three-phase sinusoidal pulse width modulation (3xPWM) that is synchronized with the rotor angular position using the absolute encoder (E) from Kübler (Optic 5876; Kiev, Ukraine).

Left: A functional diagram of GD. Right: The actual design of the GD with the control board. M, motor; I, bridge transistor inverter; DC-DC, DC-DC converter; E, encoder.
The GD system control module utilizes an inexpensive 8-bit microcontroller (Atmega128-16AU). The microcontroller is programmed to run the algorithm of control and pulse formation to drive the transistors of the inverter and DC-DC converter. The firmware contains sine-wave lookup tables, which allow creating up to 500,000 steps for one full turn. The absolute encoder indicates a rotor position and speed (rpm). The encoder’s output signal is 12-bit Gray code (4096 bits per one full turn). The subroutine that performs the conversion of encoder signals from Gray code to straight binary code is included in the firmware. The binary code is then transmitted to the personal computer using the interface described below.
Control Interface
The overall control of the experimental setup, the adjustment of various measurements parameters, and data readout are realized using a personal computer via a National Instruments multifunction data acquisition unit (USB-6009) connected to the computer’s universal serial bus. Setting the direction of the sample rotation is performed by sending a bit to the input port of the microcontroller through one of the USB-6009 digital output lines, and the final sample position is monitored by reading encoder bits via USB-6009 digital inputs.
In total, seven digital lines were used for the bidirectional communication between a PC and the measurement setup. In this work, we developed software that defines the number of steps required to rotate the sample from the starting position to the end position. The number of steps depends on the angle of rotation. Also, it is possible to set the speed of motor rotation in real time while keeping the smooth, stable, and desired direction of rotation of the motor (clockwise/counterclockwise). To keep track of motor positioning, the encoder’s feedback data can be displayed onscreen and stored on the computer hard drive.
The signal from the photodetector is acquired through the analog channel of the USB-6009 and digitized by a 14-bit ADC.
Data Processing
Refractive index measurement experiment data processing tools were developed using National Instruments LabVIEW. The virtual instrument (LabVIEW program) allows one to monitor experimental runs, check the rotation angle of the sample, track changes of the signal at the photodetector, and perform digital signal filtering and other calculations in real time.
Results
Automated Refractive Index Measurements
The refractive index is calculated based on the interferometric pattern shift with respect to the “zero” position, in which the laser beam is strictly perpendicular to the sample surface. Figure 4 schematically illustrates the changes of the photodetector signal as the sample rotates from the initial position going through the zero angular position along the way.

Schematic representation of (a) changes in the photodetector signal when (b) approaching the zero position in the process of clockwise rotation, (c) the angular zero position is reached, and (d) rotating clockwise from the zero position.
In general, zero position determination is a nontrivial task. The procedure applied in this work is to find an assumed zero position based on the frequency change in the recorded signal. This assumed position is refined by calculating a refractive index based on the region between a certain angle (5°–7° from the possible zero position that is taken from the recorded data of the signal change) and this position. Ideally, if the zero position is correctly chosen, indices calculated for various angles should be the same (or deviate within the error of the method). That said, the algorithm of the measurement procedure ( Fig. 4 ) is given by the flowchart in Figure 5 .

A flowchart explaining the process of automated refraction index measurements.
The value of the refractive index of the crystalline sample can be calculated using the following formula: 27
where ϕ defines an angle in degrees by which the sample was rotated from the initial position, λ is the wavelength in nanometers of the laser used in the experiment, d is the thickness in millimeters of the sample, K denotes the interference fringe shift, and nair is the refractive index of the environment. Equation 1 is applicable to calculate the refractive index of crystalline materials only for parallel-plate samples in the visible range. In the case of nonparallel samples, the methodology described in Andrushchak et al. 28 can be used.
Interference Pattern Shift Analysis
Refractive index determination using this interferometric technique requires reliable data on the interference pattern shift. A typical data set corresponding to interference fringes obtained when rotating a 21.37 mm thick plane-parallel cut from an aluminum oxide (Al2O3) crystal is shown in Figure 6 (upper part).

The typical shift of the interference fringes recorded by a photodetector when rotating the plane-parallel sample from an initial position for a small angle (about 2°). The upper part shows the raw (unfiltered) signal, and the bottom part shows a filtered signal.
When we rotate the plane-parallel sample in one arm of the interferometer, it is possible to observe the oscillations of signal intensity on the photodetector, and according to previous research,29,30 the interference fringe shift (the shift that occurs when the sample placed in the measurement channel rotates perpendicular to the laser beam) value can be found by counting the exact number of intensity peaks (interference maxima). Whereas intensity signal fluctuations might occur while running the experiment (as is the case in the upper part of Fig. 6 ), for accurate calculations, especially on the larger angles of rotation, where these peaks are changing very rapidly, filtering needs to be applied before further data analysis. The lower part of Figure 6 shows the result of applying the second-order Butterworth infinite impulse response digital bandpass filter to the raw measured data set. Filter specifications were adjusted to cut off both high-frequency noise and the DC component so that the filtered signal contains data centered around the zero line and without random irrelevant fluctuations. Circles indicate zero-crossing points found by the developed LabVIEW virtual instrument. Each pair of zero-crossing points corresponds to a peak counting toward a total number of interference maxima. The example in Figure 5 shows four peaks within less than a 2° range of angle. One has to note that such a narrow range of angles is useful only for demonstration purposes and a much wider range is used in actual measurements.
Results
Crystalline samples of Al2O3 and BK7 crown glass were selected to test the automated setup. Samples were in the form of parallel plates with different thicknesses to see how signal processing procedures should be adjusted depending on the optical path of the beam. In both cases, the range of considered angles was set to 10°.
Based on the measured data and using eq 1, the refraction index of Al2O3 at room temperature was found to be 1.774. This value is in fair agreement with the data of Malitson et al. 31 (1.765 at 644 nm and 24 °C).
The second test was performed on a 10.93 mm thick plate of standard BK7 crown glass. The calculated refractive index was 1.526. The refractive index for crown glass is typically between 1.52 and 1.62. Taking into account the absolute accuracy of the investigated sample thickness determination as δd = 1 µm and the uncertainty of the He-Ne LHN-302 laser wavelength as 3·10−8 µm for λ = 0.6328 µm, 32 then for the real accuracy of interference fringe shift determination δk = 1 and angle of sample rotation δϕ = 1.4·10−3 the accuracy of the refractive index determination using the interferometer-turning method is 0.005 for the sapphire (Al2O3) sample and 0.003 for the BK7 sample. The results are summarized in Table 1 .
Comparison of Measured Data for the Test Samples Using Our Method and Data from References.
Aiming to increase accuracy, we have also calculated a refractive index for a wider angle of sample rotation starting from the initial position. For example, consideration of a 20° range of interferometric data leads to results that are closer to the reference data. The index change, in this case, is 0.002, which corresponds to a decrease of error to approximately 0.1%. Thus, a more extensive range of rotation angles slightly increases the accuracy but requires longer measurement times and a more significant number of points to be processed, making the entire experiment more time-consuming. Still, LabVIEW signal acquisition and processing capabilities provide the necessary resources not only to control the experimental runs but also to process data on the fly.
Conclusion
We developed a fully automated experimental setup to accurately measure the refractive indices of parallel optical plates. The setup uses software developed in National Instrument LabVIEW. The hardware includes a high-precision GD based on an SMPM to achieve a minimum rotation angle of 0.003°.
All-in-one software tools developed in National Instruments LabVIEW ensure complete experimental control, data collecting, and processing from setting the rotation speed to calculating the refractive index.
The setup and software were tested in the measurement of the refractive index of aluminum oxide and crown glass. Comparison of obtained values with reference data verifies that the automated method produces reliable results.
The proposed technique is suitable for express laboratory or industrial use as the time needed for a single measurement is within 5–7 min.
The hardware discussed in this paper can be adapted to work with portable lasers, and the software can be ported to an embedded system, therefore making the entire experimental system more compact. The proposed technique and respective setup are intended to be most useful to researchers in the fields of photonics and biophotonics who deal with the optical characterization of plane-parallel samples for their applications.
Footnotes
Declaration of Conflicting Interests
The authors declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The authors disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This result of the investigation is part of a project that received funding from the European Union’s Horizon 2020 research and innovation program under the Marie Sklodowska-Curie grant agreement no. 778156. This work was also supported by the Ministry of Education and Science of Ukraine in the frame of the “SubTera” project (registration no. 0119U100609).
