Detection of Interference Phase by Digital Computation of Quadrature Signals in Homodyne Laser Interferometry

We have proposed an approach to the interference phase extraction in the homodyne laser interferometry. The method employs a series of computational steps to reconstruct the signals for quadrature detection from an interference signal from a non-polarising interferometer sampled by a simple photodetector. The complexity trade-off is the use of laser beam with frequency modulation capability. It is analytically derived and its validity and performance is experimentally verified. The method has proven to be a feasible alternative for the traditional homodyne detection since it performs with comparable accuracy, especially where the optical setup complexity is principal issue and the modulation of laser beam is not a heavy burden (e.g., in multi-axis sensor or laser diode based systems).


Introduction
The laser interferometry represents the most precise technique for measurement of geometrical quantities today. Among currently known methods, it provides the best resolution and accuracy and ultimate (theoretically an infinite) dynamic range. The basic resolution of the laser interferometer corresponds to the wavelength of the laser source, but for many applications such an accuracy is insufficient. The resolution is further improved by various methods, commonly referred to as fringe subdivision or detection techniques [1,2].
Generally speaking, the detection methods usually expand the basic concept of interferometric measurement by adding extra information to the interference signal that helps to extract the phase information with an increased precision [3,4]. A typical example is the homodyne detection that handles the polarisation of the laser beam in the optical hardware so that a compound wave is produced instead of the interference signal. The phase information is then extracted using additional electro-optic system that produces a pair of interference signals in quadrature (see details in the next section). The phase is then extracted from these signals with significantly increased precision.
In this paper, we present a novel detection technique, related to a traditional homodyne detection, that relies on a digital computation instead of rather complex optical processing. It employs a frequency modulated laser source on the input of a non-polarising laser interferometer and a simple photodetector that observes the output. Then a combination of synchronous demodulation, linear transforms, mixing, phase shifting and trigonometric transforms is used to reconstruct the same pair of quadrature signals. The detection is further improved by a scale linearization techniques that mitigate the measurement error introduced when the method is applied in a real measurement system.
The presented approach brings several advantages, especially reduced demands on the precise construction, reliability and cost of laser interferometry based measurement systems. Our aim was not to present a novel technique that would allow for better precision, accuracy or bandwidth-especially the bandwidth is limited by the frequency modulation bandwidth of the laser source (we deal with DC−100 Hz in our experimentation). Instead, we proposed a method that achieves a performance comparable to that of the current methods, i.e., decent resolution (sub-nanometre).
The rest of this paper is organised as follows. Section 2 summarises the principles of the homodyne detection and presents an analytic description of the proposed method. Section 3 describes the experiments we have carried out to validate the performance of our method in comparison with a homodyne detection, while Section 4 presents the results and holds the discussion. Finally, Section 5 summarises and concludes the paper.

Homodyne Detection with Polarising Optics and Quadrature Detection
In a typical Michelson setup with homodyne detection, the interferometer setup contains a polarising beam splitter to produce the measurement and reference beams with a different polarisation plane. As a result, there is a complex light beam compound from the waves with mutually perpendicular polarisation on the output of the interferometer rather than a simple interference. That compound beam is further processed by an optoelectronic detection system, usually referred to as a quadrature detection unit (schematically shown in Figure 1), that produces a pair of interference signals I x and I y that are in quadrature, i.e., mutually phase shifted by π/2. Figure 1. Polarising laser interferometer with quadrature detection unit for homodyne detection: the interferometer employs the polarising divider, the reference corner cube reflector CC ref and measurement reflector CC meas . The output beam is compound from two beams with a different polarisation. The detection unit splits the beam into two detection branches on the non-polarising divider NP. In the first branch, the polarising beam splitter PBS1 adjusts the polarisation to a single plane so that the interference occurs and is observed by photodetectors D1, D1' (mutually in an opposite phase). In the second branch, one of the wave components is retarded by the retarder plate RP. The photodetectors D2, D2' observe the interference that occurs due to rotation of the polarisation plane on the polarising beam splitter PBS2. The subtractor S1 produces the signal I x as a difference between signals D1 and D1', the subtractor S2 produces I y from D2, D2' in the same way. The immediate value of the interference phase φ(t) is extracted by a quadrature detection, i.e., the interference phase is calculated using an inverse function from immediate values of the X-and Y-axis signals I x and I y .
One of our objectives was to propose a method that would be simpler in terms of optical hardware complexity. This complexity is represented by, e.g., the polarising beam splitters, the corner reflectors, semi-transparent mirrors and the demands for optical system adjustment. This is why using exclusively non-polarising interferometer with planar reflectors and a simple photodetector can be considered interesting direction towards cost-effective measurement systems, despite the necessity of laser beam modulation. Several extensions of the basic scheme of the homodyne detection are currently used in modern measurement systems, e.g., the use of planar mirrors [5] or advanced retroreflectors (MI-series interferometer, SIOS Messtechnik) instead of cube corner reflectors or the multichannel photodetector capable of producing the quadrature signals (Renishaw).
The next section presents the analytical description of the interference phenomenon and the signals produced by the quadrature detection unit. Note that the laser beam is considered as a set of monochromatic plane waves in vacuum since an ideal source of coherent light is assumed.

Quadrature Signals
The electromagnetic wave generated by a laser source can be described at a given point by a vector of electric field E that oscillates in the plane perpendicular to the direction of propagation. If the direction of oscillations is fixed we call the wave linearly polarised and we can describe it by the magnitude of electric field E. For definiteness, assume that the laser source is placed at the origin of the Cartesian coordinate system and the wave propagates along the x axis.
The electric field at the source is given by the real part of the complex function where E 0 is amplitude and φ(t) is a general phase. By definition, the phase is related to the instantaneous angular frequency ω(t) by Hence the phase can be calculated from the known angular frequency as The propagation of the wave from the source to the point x = l takes the time τ = l/c, where c is a velocity of propagation. Therefore at the point x = l the electric field is where φ(t − τ ) is a phase shifted by time interval τ . The intensity observed by a detector placed at the point x = l is proportional to the squared absolute value of electric field and is given by where the asterisk denotes the complex conjugation and q is a proportionality constant. Now assume there are two waves propagating along a single axis (practically representing two waves that travel along two arms of the interferometer). At the output of the interferometer the two waves are recombined, the interference occurs and the compound wave, usually denoted as an interference signal, is observed. Let l 1 be the distance travelled from the source to the detector by the wave that propagates through the reference arm and l 2 the distance through the measurement arm. Then the electric field of the compound wave at the detector is where τ 1 = l 1 /c and τ 2 = l 2 /c. The intensity at the point of observation can be calculated using Equation (6). Denoting I 1 = qE 2 01 and I 2 = qE 2 02 and using the Euler's identity for complex numbers we obtain the observed intensity Changing the variables t − τ 1 → t and denoting τ 2 − τ 1 = τ , we rewrite the observed intensity as We see that the observed intensity depends only on amplitudes of two interfering waves and their phase difference. The latter is caused by the path difference l = l 2 − l 1 between the paths travelled by two interfering waves.
When the laser source works on single frequency f , the angular frequency of the wave is constant and can be calculated as where λ is wavelength. From Equation (4) we get the phase of the first wave where φ 0 is the initial phase. Similarly, the shifted phase of the second wave is The intensity I x observed by the first branch of the detection unit can be calculated by means of Equation (9). We get I x (τ, t) = I 1 + I 2 + 2 I 1 I 2 cos(ωτ ) which is a well-known result that allows one to observe interference fringes when the time delay τ = (l 2 − l 1 )/c is not constant. This can be easily achieved by changing the length l 2 of the measurement arm. The intensity I y observed by the second branch of the detection unit (the Y-signal) that involves the retarder plate is shifted from I x by the phase factor π/2 and is given by On the basis of analytic model of the quadrature signals, the next subsection presents the approach to computational reconstruction of quadrature signals from a signal observed by simple detector.

Reconstruction of Quadrature Signals
As mentioned above, we employ a frequency modulated laser beam instead of a monochromatic one. Similar approach has been already employed, e.g., in phase-shifting and phase-modulation based interferometry [6][7][8], frequency modulated continuous wave interferometry [9] or for absolute distance measurement [4,10], nonetheless the follow-up evaluation is principally different from our approach.
Particularly, we apply the sinusoidal frequency modulation. First, we denote the mean value of angular frequency of the wave ω s and the amplitude of the sinusoidal modulation ω m . Then the angular frequency of the first wave is where Ω is modulation angular frequency. The phase of the first wave can be obtained by integration of the angular frequency. We get Similarly, the phase of the second wave shifted by the time delay τ is The intensity I x observed by the first branch of the detection unit can be calculated using Equation (9). For clarity, we first calculate the phase difference φ(t) − φ(t − τ ). Using the standard relations for trigonometric functions we obtain This expression can be further simplified on the basis of the following considerations. Due to the typical dimensions of commonly used interferometers and a large value of the speed of light, the time delay τ is negligible compared to the time t. Moreover, the term Ωτ 1. Hence we can make the approximations The phase difference is then Inserting this expression into Equation (9) we finally get the intensity I x observed by the first branch of the detection unit I x (τ, t) = I 1 + I 2 + 2 The intensity I y observed by the second branch of the detection unit differs from I x only by global phase factor π/2, therefore we have We see that the observed intensity is variable even if the time delay τ is constant. This is obviously caused by the modulation of the angular frequency of the wave emitted by the source.
The intensities I x and I y observed by particular branches of detection unit given by Equations (21) and (22) were derived for the general case of sinusoidal modulation of the angular frequency. In these calculations the depth of modulation ω m and time delay τ were not subjected to any restrictions. However, we show that under certain assumptions for these two variables, the Equations (21) and (22) can be substantially simplified.
Since the absolute value I 1 + I 2 and the factor 2 √ I 1 I 2 are common for both I x and I y , we omit them for better readability of the following text. Mathematically, this corresponds to the normalisation of functions I x and I y by means of relations Then we start with normalised intensities Using standard formulae for trigonometric functions we can rewrite them tõ Since the maximum absolute value of the term sin(Ωt) equals to unity, assuming that ω m τ 1 we can make the following approximations and for observed intensities we get formulaẽ Both observed intensitiesĨ x andĨ y are now expressed as a sum of two terms. We explain the meaning of particular terms for intensityĨ x . The discussion forĨ y can be done in a similar way.
The first term ofĨ x expressed as cos(ω s τ ) depends only on the mean value of angular frequency ω s and time delay τ . Hence it is identical to the corresponding term calculated for interference of two monochromatic waves, see Equation (13). The second term ofĨ x is the product of amplitude ω m τ sin(ω s τ ) and oscillating function sin(Ωt), which corresponds to the modulation of angular frequency. Since we have assumed that ω m τ 1, the amplitude of the second term ofĨ x is much smaller than the amplitude of the first term.
For further considerations we make the following assumption. When the mirror in measurement arm is moving we demand that the function sin(Ωt) oscillates much faster than cos(ω s τ ). Mathematically, this corresponds to the relation d(Ωt) dt For illustration, considering the motion of the mirror in measurement arm with a constant speed v, the last condition is expressed as Ω ω s v/c. The mean values Ĩ x and Ĩ y of observed intensities can be calculated using the relations where T = 2π/Ω is a period of frequency modulation. The stated assumption Equation (28) means that the value of functions sin(ω s τ ) and cos(ω s τ ) can be considered constant in time interval from 0 to T . Moreover, since the integral of sine and cosine function over the period is zero, we get the mean values We see that under the stated assumptions, the mean values of observed intensities when the frequency is modulated are identical to the intensities observed when the frequency is constant, see Equations (13) and (14). Note that the suitable depth of modulation is crucial for the reconstruction of the quadrature signals.
Now we would like to show the procedure that allows us to calculate the mean value of intensity Ĩ y from the observed intensityĨ x . We start with intensityĨ x at the time t given bỹ and its value at neighbouring time t − ∆t, which is Calculation of their difference ∆Ĩ x gives us Next, multiplying the difference ∆Ĩ x by cos(Ωt) (note that it is a phase-shifted copy of modulation signal, originally sin(Ωt)) and calculating the mean value we get the derived signal Finally, from comparison with Equation (30) the derived signal I d is We see that the derived signal I d is proportional to the mean value of intensity Ĩ y . Hence, we are able to reconstruct the mean intensity Ĩ y from intensityĨ x .
The last part of the novel detection method, presented in the final part of this section, addresses the issues with the practical implementation.

Scale Linearity Compensation
We have identified several sources of inaccuracies related to a practical implementation of the proposed method: (i) the measured intensities I x and I y are generally burdened with a measurement error, e.g., noise or the quantisation error (i.e., the error caused by the discretization of analog measure), (ii) the currently available laser sources are liable to a residual amplitude modulation (RAM) [11], i.e., the modulation of their optical frequency has an unwanted influence on the intensity and spatial characteristics, (iii) the computation over discrete representation of the analog signals causes additional inaccuracy, e.g., numeric calculation of derivation, round-off error.
Note that the list of influences is far from being complete, nonetheless these were the most significant ones that we met and found important to be addressed.
Currently, the problems with scale linearity (i) are well explored [12,13] and there are several methods to deal with them [14][15][16].
The RAM effect (ii) is superimposed onto the effect of the laser frequency modulation. As a consequence, there occurs an additional phase shift between I x and I y besides the desired π/2 when our detection method is used. This results in an unwanted offset in the detected interference phase that depends on the difference of path-lengths and the modulation depth. The effect is additionally amplified due to the discrete computation (iii), especially due to the calculation of the numeric difference over a discrete set of sampled points.
The phase shift is typically corrected by the scale linearization techniques, nonetheless our preliminary experimentation indicated that a compensation of the phase shift between the quadrature signals before the elliptic-fitting based linearization technique leads to better overall performance of our method. On the basis of this experience, we have employed both techniques: the phase shift removal is described in the following section and for the scale linearization technique we used an existing method [14]-the latter method is not described in this article.

Removal of Unwanted Phase Shift between Quadrature Signals
The principle of our phase shift removal method is the following: after calculation of I y from the observed I x we detect their mutual phase shift δ a (caused mainly by RAM), then we transform I y in an appropriate manner to correct this phase shift and after that we use the linearization method.
The following description presents the method for a general case, i.e., we assume unequal amplitudes of the quadrature signals, even though we have operated with unit amplitudes of I x and I y before.
We start again with intensities I x and I y in the form where we assume that their mean value is zero and K x , K y are real amplitudes of I x , I y . Note that the amplitudes are known from the observed time course of I x and I y in practical applications. Next, we introduce a quantity K xy that corresponds to the amplitude of the sum I x + I y , formally written as The value K xy is also known from observed data. We illustrate the functions I x and I y given by Equation (36) in a phasor diagram using two phasors K x and K y of lengths K x and K y , respectively (see Figure 2(a)). These phasors together form a constant angle π/2 − δ a resulting from Equation (36), but they both rotate around the origin as the phase φ changes. In such a diagram the values I x and I y are simply given by orthogonal projection of the phasors K x and K y to the x-axis, respectively. Figure 2. Detection of the unwanted phase shift δ a : the known instantaneous values I x , I y are projected to phasors K x , K y that form the angles φ (equal to the immediate interference phase) and φ = φ − π/2 + δ a with axis x (a); a combination of K x and K y forms phasor K xy (b). Since there exists such φ such that K xy is parallel to axis x, the real angle α = π/2 + δ a can be calculated using law of cosines.
The geometrical meaning of quantity K xy is the following: since it was defined as the amplitude of sum I x + I y , it represents a length of phasor K xy = K x + K y in the phasor diagram, see Figure 2(b). To derive the expression between the lengths K x , K y , K xy and phase difference δ a , we use the law of cosines. Since the angle α shown in Figure 2(b) equals to δ a + π/2, we get It should be emphasised that the function K xy (δ a ) is monotonic in the range −π/2 ≤ δ a ≤ π/2 for arbitrary values of K x and K y , which can be easily verified by calculating the derivative of K xy with respect to δ a . This allows us to solve Equation (38) with respect to δ a . We obtain Since the value K xy and amplitudes K x , K y can be easily determined from measured data, we are able to calculate the phase shift δ a and eliminate the unwanted phase shift between the quadrature signals. For this purpose we now need to transform intensities I x , I y given by Equation (36) to another pair I x , I y given by Here I x is chosen to be identical to I x , for I y we need to derive the transformation relation expressed in terms of I x , I y and δ a . Therefore we rewrite Equation (36) as Solving the second equation with respect to I y we obtain the transformation relations The transformed intensities I x and I y are now deprived from the additional phase shift δ a caused by amplitude instability of the laser source. As a next step, the signals are scaled to unit amplitude, the scale linearization is applied and then the phase is extracted.
Note that the algorithm is principally based on a single data-point determined from the measured data (maximum value of I x , I y ,...), so that some outlying data points might cause a substantial error in the calculations. Since the K x and K y (and δ a as well) can be considered constant, the method can rely on some kind of statistical information (and actually has within our experiments).

Experimental Verification
We have carried out several experiments to verify the performance of the proposed technique. We have assembled an optical setup that allowed us to make the displacement measurements at the nanometre scale referenced to an interferometer with the homodyne detection technique, i.e., quadrature detection unit.
The reference was physically provided by a calibrated quadrature detection unit in combination with the support logic. The reference delivers the reference phase information with measurement uncertainty of N/1, 024 (where N denotes single interference fringe), i.e., 0.35 • [17] (including the linearity compensation).

Experimental Setup
The experimental setup, schematically shown in Figure 3, comprised a Michelson interferometer with polarising optics and two detection branches: the reference branch with a quadrature detection unit and the testing branch that performed the proposed technique for interference phase detection.
A solid-state laser (frequency doubled Nd:YAG at 532 nm with linearly polarised beam; PROMETHEUS by Innolight, Germany) was used. The output beam was collimated to particular diameter (3 mm) that optimally fits the detection unit and the photodetector. A half-wave plate (λ/2) was used to adjust the polarisation plane. The measuring reflector was mounted on a piezo actuator and on a motorised stage. This combination allowed for the measurement on an operational range of 600 mm. Note that the polarising variant of the beam splitter is used, as it is a necessary condition for the reference interferometer (the detection apparatus, not to be confused with the reference arm of the interferometer). The proposed method does not require preserving different polarisation in individual arms.
A non-polarising splitter (NP) is used to form two detection branches: the testing branch that involves our novel detection method with an ordinary photodetector (PD) and the reference one that uses the traditional detection unit (denoted QDU). . Experimental setup scheme: the interferometer employs the polarising divider, the reference corner cube CC ref and measurement corner cube CC meas . The laser input is attenuated by a neutral density absorption filters ND, collimated by an optical expander EXP and its polarisation plane is rotated by the half-wave plate λ/2. The interferometer output is split into two detection branches by the non-polarising divider NP. In the first branch, the interference occurs on the polariser POL and the intensity is observed by a photodetector PD. In the second branch, the signal from the interferometer is processed by the reference quadrature detection unit QDU. The acquisition and control hardware drives the modulation of the laser source optical frequency and evaluates the signals from both the photodetector and the quadrature detection unit. The bandwidth of both detection methods is limited by the combination of several factors: response of the photodetector, signal amplifiers, analog-to-digital conversion sample rate and processing capability. In our particular case, the available bandwidth was 50 kHz for the QDU and 100 Hz for the novel method (i.e., 100 phase readings by QDU and 200 phase readings by tested method per second). The optical line-up in the testing branch is extended by a polariser plate that alters the orthogonal polarisations into single plane so that testing branch observes the interference signal instead of compound wave (that is processed by the reference branch). The bandwidth of our novel method is relatively low, nonetheless the sense was to demonstrate the principle and compare the linearity with a reference method.
The signal generation (for laser source modulation) and acquisition was done by a dedicated DAQ device (U2531 by Agilent). The generation part produced synchronised sine and cosine wave-the sine controlled the modulation, the cosine was synchronously sampled (acquired) together with the output of the QDU and the PD output and used for the quadrature signal reconstruction. The reason to generate the modulation signal and its phase shifted copy and then sample the copy again was to achieve precise synchronisation between the generated outputs and signal inputs.

Design of Experiments
The measurement arm hardware was successively placed into different positions (measurement points) {0, 1, 2.5, 5, 7.5, 10, 15, 25, 37.5, 50, 62.5, 75, 100, 150, 200, 250, 300, 400, 500, 600} mm to cover the operating range ∆L ∈ [0, 600] mm. At each point, the measurement arm was continuously displaced (with approximately constant velocity using the piezo actuator) over the trajectory of ≈0.6 µm (more than two complete fringes) and the detection outputs were acquired. The selection of points was suggested by preliminary experiments that indicated that the correct distance determination is more problematic at the beginning of the operating range. The zero position (∆L = 0) actually denotes a closest achievable position of the measurement arm, given by the mechanical construction of interferometer components and corresponds actually to the path length difference of approximately 2 mm. The laser source has been modulated by a free running sine wave (at 1 kHz). The depth of modulation was empirically estimated and manually configured for each point-the optical frequency was tuned approximately by ±20...50 MHz.
The outputs from the quadrature detection unit I x , I y , the photodetector signal I p and modulation signal M in were synchronously acquired (at sample rate of 100 kHz). One phase reading of the novel method took a time period that corresponds to five complete modulation periods (5 ms). Three measurement cycles were recorded at each measurement point. The evaluation was performed afterwards over the entire cycles at once.
We have inspected the difference between the interference phase information extracted by the reference and by the testing method, referred to as the phase determination error (PDE). In the reference branch, the X-and Y-axis signals from the quadrature detection unit were recorded and the scale linearization was applied. The photodetector output I p was used for reconstruction (standing for the I x when related to the previous analytic description) of the quadrature signals using the presented method.

Results
A total of 60 measurement cycles in 20 positions was completed, where each cycle contains 890 phase readings. A sample evaluation of a single cycle is shown in Figure 4.
There is an apparently good coincidence between the quadrature signals from reference and testing detection branch (compare Figure 4(a) to Figure 4(b)) (the slight fluctuations are caused by changes in refractive index of air, but since the optical paths are shared, the influence is subtracted within the comparison). The phase extraction error plot (see Figure 4(c)) indicates that the PDE keeps within ±1 • range (for the particular sample).
The linearity of our method is summarised in Figure 5. There is an observable periodic error, and the mean PDE fluctuates within ±0.15 • . The Figure 6 shows the overall dependence of the mean PDE on the path length difference when the results from all three cycles measured at the same distance point were combined so that the repeatability aspect is incorporated. The overall mean PDE stays below ±0.5 • and the standard deviation stays below ±0.75 • . The jitter of the mean error over the operational distance can be attributed to the optical adjustment. The combination caused a slightly greater scatter of the phase determination error.
Several related issues were omitted within the experimentation, e.g., the refraction index of air influence [18], the uncertainty of the laser source stability or a precise temperature stabilisation. In our case, when two detector paths were compared, these influences were suppressed by the differential measurement, since the optical paths and signal acquisition within the detection chain were common for both methods. The results show that the novel method has a comparable accuracy to the reference method. Exploring the limits of our proposed method is still a subject to a future research and a comparison with a more accurate method is required to state the absolute precision and uncertainty. Figure 4. Evaluation of single experimental cycle: the X-and Y-axis signals from quadrature detection unit, i.e., the Ĩ x (τ, t) and Ĩ y (τ, t) (refer to Equation (30)), with scale linearity compensation applied (a); X-and Y-axis signals reconstructed by proposed method, i.e., Ĩ x (τ, t) and I d (Equation (35)

Discussion
The results indicate that the proposed method is a suitable alternative to the traditional homodyne detection, especially in situation when the accuracy is not the most critical aim but the construction simplicity and scalability is. There are several issues that need to be considered to judge the novel approach: its dynamic range, the trade-off between optical complexity and the requirement for laser modulation, the response bandwidth limitations and the uncertainty introduced by the computations.
As the homodyne detection has theoretically infinite dynamic range, it allows for the measurement with the same precision even on long distances (several kilometres). The practical limitation is given only by the construction issues, laser stability and influences of the environment. The novel method is principally limited by the available depth of laser frequency modulation [19], but using an approach for an adaptive control of the modulation depth the method would be probably applicable on a comparable operational range.
There are also higher demands posed on the laser source-there is the need for the laser beam frequency modulation. The modulation bandwidth limits the response of the measurement systems since the duration of at least one complete modulation period (five in case of our experiments) is processed to produce a single phase reading. Consequently, velocity of measurement is significantly slower than the optical processing. The demand for laser with a fast modulation capability can be fulfilled, to a certain extent, by cost-effective semiconductor lasers [20,21].
On the other hand, the method requires less complex optical setup to achieve comparable results-there is no need for polarising beam splitters and semi-transparent mirrors within the optical lineup. As a consequence, the method also scales well, e.g., when several measurement axes would be required within a measurement system. In such a case, the added complexity per each added axis is minimal in comparison to the traditional means-it has been already demonstrated in a practical application for surface diagnostics [22] and multi-axis geometric measurements [23]. An extra axis requires only a small amount of additional optical components and poses no substantial additional requirement on the laser source (assuming there is enough power available, but it is usually not an issue, e.g., for green Nd:YAG lasers).
In our experiments, the modulation bandwidth was limited by available high-voltage amplifiers to 200 phase readings per seconds (corresponds to 1 kHz modulation). With the particular laser source, a 100× higher response would be achievable; further enhancement would be possible using and acousto-optic modulator.
In the basic principle, similar computational methods are also sensitive to the noise issues. In case of our novel method, the issue is well mitigated by a massive averaging and suitable signal filtering without increased demands on system performance as the operations are essential parts of the phase extraction method.
There were also several assumptions we stated (in sake of simplicity) in the physical description of the novel approach. We assume the time delay τ -the difference in wave propagation along unequal optical paths in the arms of the interferometer-is negligible (see the paragraph above Equation (19)). In the borderline case, the path length difference l = 1, 200 mm, so the time delay τ ≈ 4 × 10 −9 and the term Ωτ ≈ 4 × 10 −6 . We also assume that the product of τ and the modulation depth ω m τ 1 (Equation (26)). In the case when l = 1, 200 mm and ω m = 20 MHz, the product ω m τ = 0.08, which corresponds to the limit of validity of our assumption ω m τ 1. Finally, we assume that the phase change (due to displacement) is slower than the modulation of the laser, i.e., Ω ω s v/c (Equation (28)). The displacement velocity v was v ≈ 2.5 × 10 −7 m/s, the Ω = 1 kHz, thus for the mean frequency ω s ≈ 563.25 THz the ratio is 1 × 10 3 7.5 × 10 −6 .

Conclusions
We have proposed and verified a novel detection technique that represents a novel approach to the homodyne detection in laser interferometry.
A detailed analytic description of the principles has been presented and the validity and performance were evaluated by experimental means.
The experimental results indicate that the novel approach performs comparably to the traditional detection technique. Since it requires simpler optical hardware setup-e.g., exclusively a non-polarising interferometer, single photodetector, but frequency modulated laser beam and computational power-it can be considered a suitable alternative to the techniques that rely on additional optical processing.