## 1. Introduction

In order to determine amplitude and phase in a two-arm interferometer, a phase modulator is often required to clearly discriminate phase changes from amplitude changes. Such operation is straightforward when the phase modulation is a linear function of time. In this case, a standard Lock-in Amplifier (LIA) gives the required information with an optimal signal to noise ratio (SNR). Unfortunately, a number of phase modulators interesting in term of cost, achromaticity or integration offers non-linear responses, that are even sometimes coupled with unwanted amplitude modulation. A critical question that arises is “How can we extract phase and amplitude information in an optimal way when non-linear phase modulation is used?”

To solve this issue while keeping the benefits of high SNR, approaches have been proposed based on multiple lock-in detection at selected signal harmonics. These approaches were mainly employed in the case where the phase modulation is a sine function [1–6]. Such phase modulation is, for example, achieved using piezo-actuator, fiber stretchers, and other phase modulators where a sine excitation typically offers the best response. The multiple lock-in approach works fine but it is less direct and does not necessarily provide an optimal SNR or a straightforward implementation. Especially, if an amplitude modulation is present at the same frequency as that of the phase modulation. Alternately, the Generalized Lock-in Amplifier (G-LIA) technique was recently introduced [7] to solve this issue with a procedure similar to a single LIA operation. In this chapter, we first detail the principle of this method when operated in the simplest case where no amplitude modulation is present. Application is provided notably in the context of digital holography. Then we consider the case where there is an additional amplitude modulation in the signal field. The first case which is discussed is related to unbalanced interferometry where the phase modulation is achieved via a power modulation of the laser source. Finally, we also discuss the case of phase-sensitive near-field imaging.

## 2. Theory: introduction to G-LIA

In order to introduce the G-LIA technique, we need to provide an expression for the detected signal. We consider the simplest configuration of a 2-arm interferometer comprising a reference arm and a signal arm (cf. **Figure 1**). The system is illuminated by a monochromatic radiation. The detected signal intensity *I*(*t*) can be expressed as:

where *E*_{s} and *E*_{r} are, respectively, the amplitude of the field of interest and the reference field impinging on the detector. The phase difference between the two fields is Δ*φ*(*t*) while the factor *m* ≤ 1 in the interferometric term accounts for the interference contrast. Alternately, it is common to express the detected intensity as a function of the laser power *P*:

where *s* is proportional to the unknown signal amplitude *E*_{s}, considered as constant during a measurement. The unknown spatial phase *φ*_{s} of the signal field is also supposed to be constant during the measurement although a phase modulator can be included inside the signal arm. Any time dependence in Δ*φ*(*t*) is therefore arbitrarily considered as coming from the reference field:

From the expression of *I* given by Eq. (2), different strategies can be proposed in order to recover amplitude and phase information. In order to extract the two unknowns (*s*, *φ*_{S}) from the signal *I*(*t*), the phase of one of the two beams can be modulated in time by a frequency shifter or another phase modulator. When the time dependence is linear *φ*_{r}(*t*) ∝ *t* as in the former case, (*s*, *φ*_{S}) are precisely determined by a Lock-in Amplifier (LIA) locked at the single frequency component present in *I*(*t*). For other functions of time such as a sine waveforme (i.e. *φ*_{R} = *α* sin Ω*t*, where *a* is the modulation depth), the use of a conventional LIA is less trivial as the signal information is typically spread over a number of frequency components (nΩ/2π) having different weights. The G-LIA method was introduced to handle such cases, while keeping an extraction procedure very similar to a LIA. Hereafter, the two approaches are provided to highlight the similitudes and differences.

### 2.1. Amplitude and phase determination using a standard LIA

#### 2.1.1. Case of linear phase modulation (LIA)

Considering the general Eq. (2), the use of LIA is direct when Δφ(t) is linearly modulated by *φ*_{R} (*t*) = Ω*t* = 2*π*Δ*Ft*, where the phase modulation rate Ω can be induced in different ways. These ways include notably the use of a frequency shifter in one of the two arms (where Δ*F* is the frequency shift), a linear translation of one of the mirror (where Δ*F* is the associated Doppler shift), or a linear variation of the laser frequency1 if the two arms are unbalanced. Amongst the cited methods, heterodyne measurement based on the use of frequency shifters is often considered as more favorable as a purely linear displacement of the mirror is hardly achievable in practise without alignment or coherence issues, while unbalanced interferometry is subject to noise [8] induced by small wavelength fluctuation typically related to temperature drifts.

The two unknowns (*s*, *φ*_{S}) are simply the amplitude and phase of the sine waveform present in Eq. (1) which is modulated at the frequency shift Δ*F*. This determination is optimally carried out by a dual-output (*X*, *Y*) LIA locked at the frequency Δ*F* that is precisely provided by the modulator driver. Depending on the output, the detected signal *I* is multiplied by an in-phase or a quadrature sinusoids modulated at the same angular frequency Δ*F* and it is averaged over a time *t*_{int}:

From these two outputs, the quantities *LIA reference signals* built from the frequency shift precisely provided to the LIA. If the LIA is not locked exactly at the angular frequency Ω, the measured values of *φ*_{s} will drift in time, while the extracted signal amplitude will decay for long integration time.

#### 2.1.2. Case of a non-linear phase modulation (LIA)

Achieving *φ*_{r} (*t*) ∝ *t* is not possible for a number of phase modulators, notably because they have a finite range of phase modulation. Pseudo-heterodyne approaches were proposed long ago to circumvent this problem [9] by using a sawtooth modulation of the optical path, where the peak to peak amplitude of the sawtooth corresponds to an integer number of times 2π in term of phase. The approach is illustrated by **Figure 1**, in the case where the phase ramp is achieved by a piezo-actuated mirror in a balanced interferometer.

As can be seen, the detected intensity mimics the sinusoidal beating observed in heterodyne setups. Such approach is not widely used since errors are induced during the flyback time on the sawtooth edges, especially if the modulation is fast.

As mentioned, the use of sine modulation *φ*_{R} = *α* sin Ω*t* is regarded as much more desirable as most of the modulator can operate better and faster when they are sinusoidal excited. It was early highlighted in this context [10] that for such modulation, the Fourier spectrum of the signal has harmonic sidebands coming from the interferometric term *I*_{mod} ∝ *P* *s* cos (Δ *φ* (*t*)) in Eq. (2), as shown in **Figure 2**.

The amplitudes of these frequency components are obtained by developing the term in *I*_{mod }∝ *P*[cos (*φ*_{s}) cos (*φ*_{R}) + sin(*φ*_{s}) sin (*φ*_{R})] and using the Jacobi-Anger expansion of cos (*φ*_{R}) (even harmonics) and sin (*φ*_{R}) (odd harmonics) [11]. From this expansion, we see that a LIA locked at an angular frequency harmonics mΩ, with *m* ≠ 0, gives:

where *J*_{m} (*a*) is the *m*-th Bessel function. The amplitude and phase can then be extracted using both odd and even harmonics, for example, using

When only two harmonics *m* = (1, 2) are used, caution must be exercised in the choice of the modulation depth *a* in order to maximize the power density on the selected harmonics. The optimum value of *a* in this case is *a *= 2.19 rad (maximum of *J*_{1} (*a*)^{2} + *J*_{2}(*a*)^{2}).

### 2.2. G-LIA method

The main benefit of the G-LIA method is that all the weighted harmonics are used to retrieve phase and amplitude with an operation similar to that of a LIA. To introduce this method, we also remark that the interferometric term *I*_{mod} ∝ *P* *s* cos (Δ*φ*(*t*)) in Eq. (2) can be expressed as *I*_{mod} ∝ *P* *s* [cos (*φ*_{s}) *C*(*t*) + sin (*φ*_{s}) *S*(*t*)], with *C*(*t*) = *cos*(*φ*_{R}) and *S*(*t*) = *sin*(*φ*_{R}). From this expression, it appears that *C*(*t*) and *S*(*t*) can be used as relevant reference signals within a modified LIA having the two following outputs:

These references contain the same frequency components than the interferometric term since *I*_{mod}(*t*) is a function of *C*(*t*) and *S*(*t*). These frequency components are naturally weighted so that the contribution of stronger harmonics will be favored. **Figure 2(b)** exemplifies the case where the phase modulation function is a sine function.

We note that in the particular case where *φ*_{r}(*t*) = Ω*t*, *C*(*t*) = cos(Ω*t*) and *S*(*t*)=sin(Ω*t*) and the G-LIA operation degenerates to that of a standard LIA. More generally, by replacing the expression of *I*_{mod} in Eqs. (8) and (9), we see that for any phase modulation we have:

where *k*_{x} = < *C*^{2}(*t*) >, *C*(*t*) and *S*(*t*) are orthogonal, that is, *X*_{φR} and *Y*_{φR} outputs are:

The G-LIA outputs are then similar to that of the LIA in the linear case (cf. Eqs. (4)–(5)). The difference is the presence of the additional proportionality constants *k*_{x} and *k*_{y} which need to be evaluated by calculating the average values <*C*^{2}(*t*)> and <*S*^{2}(*t*)>, respectively. In the case of sine phase modulation *φ*_{R} = *a* sin Ω*t*, the constants have analytical expressions obtained from the integral representations of Bessel functions and simple trigonometric developments:

where *J*_{0} is the Bessel function of first kind. However, it is difficult to extract *I*_{mod} from *I*(*t*) to feed the G-LIA input as it requires to precisely remove the signal and reference field intensities from the detected intensity *I*(*t*). In general, it is not possible either to use directly *I*(*t*) to perform the G-LIA operation described above. The reason is that for phase modulation function such as sine or triangle, the useful term *I*_{mod} also contains a DC component. In consequence the references *C*(*t*) and/or *S*(*t*) also contain a DC term so that the constant, non-interferometric term *P* in Eq. (2) is also detected by the G-LIA if *I*(*t*) is used directly. To avoid this problem, specific phase modulation depths for which both *C*(*t*) and *S*(*t*) do not have a DC component can be used.2 For a sine phase modulation, this specific depth of modulation *a* corresponds to the zeros of the *J*_{0}(*a*), for example, *a *= 2.405 rad. We note that there is no prejudice in term of signal to noise ratio using this *a* since all the harmonic contents is detected by the G-LIA operation.

Alternately, a satisfactory solution is to filter the detected intensity to remove all DC component from the signal. In fact, such operation is easy to do and is often highly desirable to directly remove the ambient light contribution in normal conditions [12]. In this case where the signal is filtered, the G-LIA operation is:

where the *Ĩ* denotes a DC-filtered quantity,3 that is,

In the useful case of a sine modulation of the form *φ*_{R} = *a* sin Ω*t*, the constants *K*_{x} and *K*_{y} have the following analytical expressions:

where the negative extra term in *K*_{x} comes from the filtering of the DC component which is not zero for *C*(*t*). As shown in **Figure 3**, for certain phase modulation amplitude the two constants are identical and approximately equal to unity. This is obtained for *a* ≈ 2.814 rad, but any phase modulation can be used. A better signal to noise ratio (SNR) is naturally achieved when there is no DC component in *C*(*t*) and *S*(*t*), since this part is filtered. These cases correspond to the zeros of *J*_{0}(*a*) as previously mentioned (e.g. *a *= 2.405 rad), however the SNR is nearly optimum for a continuous range of values above *a *= 2 rad. Other analytical expressions can be given, for example, for a triangular modulation [7], however the constants estimation can be made numerically without difficulty for a variety of phase modulation functions.

## 3. Application examples

In this section we review and present several results of interferometric measurements performed with the G-LIA approach described in the previous section. Results include measurement with a point detector reported elsewhere and interferometric measurement with 2D detector in the framework of holographic measurement.

### 3.1. Measurement with a point detector

**Figure 4** shows measurement results adapted from the Ref. [7], where the G-LIA can be used with or without filtering to monitor an arbitrary displacement (here a triangle-shaped displacement).

#### 3.1.1. Displacement measurement

The setup is shown on **Figure 4(a)**. In this example, the phase modulation is a sine function *φ*_{R} = *a* sin Ω*t*, with a modulation depth approximately equals to *a *= 2.405 rad. As explained before, with this value it is not necessary to filter the detected intensity to extract amplitude and phase from the G-LIA operations *X* = <*I C*(*t*)> and *Y* = <*I S*(*t*)>. However, it is a common practise to filter the DC signal directly after the detector to get rid of the environmental light condition before acquisition in order to optimize the analog to digital conversion. In the Michelson configuration shown here, the value of *a* corresponds to a peak to valley oscillation of the reference mirror of about 38.3% of the wavelength. If no position sensor is present on the reference mirrors, several methods can be used to achieve the desired phase modulation depth *a* used in the references *C*(*t*) and *S*(*t*), notably by inspecting the signal *I*(*t*) whose shape changes continuously with increasing value of *a*. Alternately the modulation depth can be precisely adjusted to recover precisely a controlled displacement of the signal mirror without affecting the amplitude output.

#### 3.1.2. Sensing

Determining the phase rather than the amplitude is known to offer potential advantage in term of sensitivity in optical sensing systems [13]. More precisely, the phase detection coupled with surface plasmon resonance (SPR) is known to improve the measurement sensitive by one to several order of magnitude depending on the exact system geometry. Many different designs on combining interferometry or heterodyne detection on Kretschmann configuration-based SPR sensor have been done [14, 15].

**Figure 5** shows the demonstration setup used in [7] to demonstrate the applicability of the G-LIA for phase sensitive sensing application. The setup is similar to that of **Figure 4(a)**, except that an *I*_{2} gas cell was added in the reference arm. The phase modulation is still achieved with a piezo-actuated mirror and the wavelength of a laser diode emitting near 660 nm is ramped over 4pm across an absorption line of the gas. To avoid phase drift induced by the wavelength ramp, the length of the arms must be precisely balanced, since even minute wavelength fluctuation can create phase fluctuation in unbalanced interferometer. This balanced setting can be achieved by ramping the laser diode wavelength outside an absorption band, and adjusting the position the reference mirror until the phase output of the G-LIA remains constant.

As can be seen, the phase varies more abruptly at the absorption peak center. However, the benefit of measuring the phase for monitoring a gas concentration is not clear since the amplitude has similar variation on the two sides of the absorption peak which indicates a similar sensitivity than the phase if the detection is made where the slope is maximum on the amplitude.

The interest of phase sensitive detection in SPR-based measurement is more obvious. In fact, strong plasmonic resonances can be reached by carefully adjusting the opto-geometrical parameters of the plasmonic layer in order to obtain very sharp phase variation across a resonance. One possible combination of phase sensitivity SPR bio-sensor using G-LIA for phase extraction is proposed in **Figure 6**, where a cuvette is put on a plasmonic chip to convey a fluid on the surface of a plasmonic chip. A coupling prism makes it possible to satisfy the Kretschmann condition for which the reflectivity of a p-polarized incident beam reaches a minimum corresponding the excitation of the plasmon-polariton surface mode. In order to have a stable phase, immune to wavelength fluctuations, the length of the two arms are made equal. **Figure 6(b)** presents the numerically calculated complex reflectivity as a function of the incident beam angle in the case of a glass coupling prism coated by a gold layer of thickness *h* covered by water and excited by a red laser. The best coupling angle is close to *α* = 70° in the provided examples.

As can be seen, the phase variation across the resonance can be made very sharp by adjusting the metal thickness *h*. It should be noted however that this fast variation is associated with a strong attenuation of the reflected beam and a compromise between signal level and sensitivity can be made depending on the available laser power. For example, on **Figure 6(b)** the reflection is about 1.6% for *h =* 48 nm at Kretschmann angle, but it drops to 2.5‰ for *h* = 50 nm.

If we consider a reasonable phase resolution of 10^{−3} rad, simple calculation show that the case *h *= 50 nm shown on the figure leads to sensitivity slightly better than 8×10^{−7} RIU (refractive index unit). On the other hand, if the thickness is slightly inferior, the RIU sensitivity drops rapidly (about 2×10^{−6} for *h* = 48 nm). To obtain a similar sensitivity with the only amplitude signal, the noise level on the amplitude measured at the maximum slope should be smaller than 10^{−3}%, which is hardly achievable.

### 3.2. Digital holography

In digital holography, the holograms of a sample object are recorded on a 2D detector such as a Charge-coupled Device (CCD) or a Complementary Metal-Oxide-Semiconductor (CMOS) camera. Such system can notably be used as an optical profilometer, or for sensing applications [1, 6, 16]. **Figure 7(a)** presents the experimental setup of a lensless, compact, digital microscope working with the G-LIA extraction method.

In the provided example, a metallic grid of slit is imaged in amplitude and phase. The Lead Zirconate Titanate (PZT) oscillates in the reference arm at 10 Hz to generate the phase modulation function *φ*_{R} = *a* sin Ω*t* at a wavelength *λ* = 640 nm in the reference arm. During the piezo oscillation, a video is recorded at a frame rate of 120 Hz. Then, amplitude and phase of the detected field is obtained on the camera by performing a G-LIA operation on each pixel using a program code. In this 2D case, the operation was not made in real-time because of the non-negligible processing time (In the order of 1 min or less depending on the computing resource). Once the detected complex field is retrieved, the associated plane wave spectrum can be obtained by Fourier transform. Then each plane wave can be back-propagated numerically to any position before the CCD [17], typically up to the sample plane or surface.

In this example, the raw signal *I*(*X*,*Y*) is processed by the numeric, software-based, G-LIA, without filtering. For a correct operation it is then mandatory to use an amplitude modulation *a* of about 2.405 rad. **Figure 7(b)** shows the amplitude and phase of the complex field which is back-propagated up to the sample plane. As we are dealing with a rough surface, the recorded phase has a speckle-like distribution. Because the illumination direction is normal to the sample surface, the system is strongly sensitive to any out-of-plane displacement of the structure. To give an idea of the system sensitivity, the sample is slightly rotated. By subtracting the phase image before rotation to the phase image after rotation, we obtain the phase-shift associated with the out-of-plane displacements, as shown in **Figure 8**.

## 4. G-LIA in the case of amplitude and phase-modulated signal

In some important cases, the signal field is modulated both in phase and amplitude. The modulated term *I*_{mod} ∝ *P* *s* cos (Δ*φ*(*t*)) in the detected intensity can then be expressed as:

where *f*(*t*) accounts for the amplitude modulation function. We can mention two relevant examples in the context of *phase-sensitive nanoscopy* and *unbalanced interferometry*.

### 4.1. Unbalanced interferometry

In interferometers having a path unbalance, a phase modulation can be efficiently induced by a wavelength modulation of the emission wavelength. For this purpose a spectrally single-mode laser diode working at a central wavelength *λ*_{0} can be used. The wavelength modulation is typically obtained by a current modulation which is associated with a power modulation. In air, the phase modulation is related to the small wavelength modulation *δλ*(*t*) by: *l* is the length difference between the two arms.

For modulation frequency below the MHz range, the change in wavelength is considered to be primarily due to a change of temperature that increases with the current. Therefore, using a sawtooth function to create a quasi-linear phase change is usually not an excellent choice, as the thermal inertia of the system prevents the wavelength to precisely follow the driving excitation. On the other hand, a sine power modulation will typically induce the desired sine wavelength and phase modulation. In this case, the detected intensity within an unbalanced interferometer is:

with *φ*_{R}(*t*) = *a* sin (*φ*_{R}), where *a*_{λ} corresponding to the depth of wavelength modulation.

It is clear that in the case where the amplitude modulation is small (*µ*≪1 and *µ*≪*s*), the amplitude modulation disappears and the G-LIA method can be applied directly using the signal *C*(*t*) = *cos* (*φ*_{R}) and *S*(*t*) = *sin* (*φ*_{R}). The operations are identical to those given by Eqs. (16)–(21). However, while µ can be small compared to 1, the signal *s* of interest can also be very small in some experiments and setups so that *µ*≪*s* cannot be satisfied in general. To see how to handle this problem, we can express the filtered intensity

where we have normalized the detected intensity by the constant laser power factor. The brackets indicate the quantity is filtered from its DC component. We see that the main issue comes from the modulated term outside the bracket which is independent from the signal *s*. A direct solution to get rid of this unwanted term is to use references without harmonic component at Ω. This is achieved by choosing *J*_{1}(*a*) = 0 (e.g. *a* = 3.832 rad). With this choice of *a*, the G-LIA operation given by Eqs. (16)–(24) gives excellent results provided that µ is kept reasonably small in comparison with 1 (e.g. for µ = 0.1 the maximum error on the G-LIA *X* and *Y* output is about 2%).5 This condition is reasonably achieved in many cases. For example, if we consider a standard single mode laser diode such as a Vertical Cavity Surface Emitting Laser (VCSEL) with a tunability of 0.5 nm/mA and a bias current of 1 mA. We see that a current modulation of 4.4% is sufficient to achieve *a *= 3.8 rad for an unbalance of 10 mm at a wavelength of 850 nm. The corresponding power modulation µ depends on its *P*(*I*) but will be typically smaller than 0.1.

#### 4.1.1. Improved unbalanced configuration

Despite its advantages in term of cost, unbalanced interferometry is not currently widely used. The main reason is also related to the extreme sensitivity of the system to minute wavelengths changes. **Figure 9(a)** represents a compensation scheme to solve this issue. The idea is to illuminate the interferometer with a linear polarization at 45° with respect to the horizontal and vertical axis and to discriminate the two s and p polarization using polarization beam splitters. An additional signal arm is equipped with a fixed mirror in order to measure the phase fluctuation induced by any wavelength drifts in time. The light impinging on this mirror is s-polarized and is selectively detected by the photodiode PD1, using a polarization beam splitter in reflection. On the other hand, the p-polarized light impinging on the piezo-actuated mirror is reflected back onto the second photodiode (PD2). Both amplitude and phases are recorded with the above described G-LIA operation.

**Figure 9(b)** shows a controlled triangular displacement which is correctly determined despite the presence of intentional wavelength drifts. In this experiment, the wavelength of the VCSEL is driven sinusoidally at about 10 kHz to create the phase modulation. The important wavelength drifts are artificially created by adding a low frequency sine to this excitation signal. The compensation is obtained by plotting the phase of the p-polarized light minus the phase of the s-polarized light which is coming from the fixed mirror. Both signal phases are obtained by the G-LIA method with *a* approximately equal to 3.83 rad.

Such system is really interesting in term of performance since VCSELs are very affordable laser sources that can be driven at very sinusoidally at very high speed. In the described experiment the phase modulation frequency was only limited by the acquisition card used to perform the G-LIA measurement.

### 4.2. Phase-sensitive nanoscopy

A modulation of the amplitude at an angular frequency Ω_{A} can be introduced in the signal arm to discriminate the amplitude-modulated signal from other unwanted contributions able to interfere with the reference field. This is the case in near-field nanoscopy where the signal light is coming from a near-field probe in interaction with a surface, oscillating at an angular frequency Ω_{probe}. The situation is depicted in **Figure 10**.

In **Figure 10**, the near-field head is included in the signal arm of a Michelson interferometer. Alternately the near-field microscope can be used in the signal arm of a Mach-Zhender which is well adapted to the characterization of waveguiding photonic devices as in [18–20]. Here, the sample is scanned under a nano-tip which is precisely positioned in the focus spot of an objective lens. The light backscattered by the oscillating probe operating in tapping mode contains information on the local optical properties of the sample. This backscattered light can have a rich harmonic content due to its near-field interaction with a sample. The amplitude modulation function appearing in *I*_{mod} is therefore

We note that *I*_{mod} in this case refers to the interference term between the *near-field signal* from the probe and the phase-modulated *reference signal*. However other parasitic fields can be backscattered by the probe-sample system. When using a standard, cantilevered, AFM probe, the oscillation amplitude is typically larger than the probe radius. In this case it is often interesting to detect the contribution from a higher harmonics Ω_{A} = *k*Ω_{probe} as the unwanted part of the light which is modulated by the shaft of the probe (rather than the apex) only contributes to the first harmonic(s). This unwanted contribution modulated by the probe shaft can be referred as *modulated background contribution* (MBC) in contrast with the *unmodulated background contribution* (UBC) coming, for example, from the sample backscattering (nano-dusts, roughness, …) which is also unwanted. In less stringent configurations, especially when the oscillation is small compared to the tip radius (the near field varies almost linearly on the excursion of the probe), the MBC is not perceived6 while keeping *k=1*.

As shown in **Figure 10(b)**, because of the interference between the probe signal and the reference field, the signal is split into sidebands at *k*Ω_{probe} ± *p*Ω where Ω characterize the phase modulation frequency. In order to collect the information spread throughout all these sidebands, the G-LIA can use the following references:

where _{A}. But even in this favorable case where the MBC is easily excluded, the UBC can be especially detrimental and has been described by several authors. UBC however can be efficiently removed with interferometric detection. As noted elsewhere, the UBC is perceived because the unmodulated background interferes with the modulated signal from the tip on the detector creating signal peaks at the harmonics frequencies *k*Ω_{probe}. Therefore a direct solution consists in excluding the frequencies Ω_{A} = *k*Ω_{probe} from the references *C*(*t*) and *S*(*t*). In the case of a sine modulation, the condition *a* = 2.405 rad (*J*_{0}(*a*) = 0) is sufficient to exclude the unwanted frequencies component by removing the possible DC contribution in cos(*φ*_{R}*t*).

First examples of phase-sensitive near-field imaging based on G-LIA can be found in Ref. [7]. In **Figure 11**, a simple demonstration experiment is made by using a bare tuning oscillating fork to modulate part of the signal at an angular frequency Ω_{A}. The trace of the experimentally detected signal intensity exhibits clearly a slow modulation due the phase modulation (*φ*_{R} = *a* sin Ω*t* with a frequency of 1 kHz) as well as a sine amplitude oscillation at the oscillation frequency of the tuning fork Ω_{A} (at about 32 kHz). An additional coverglass can be added to check the system immunity to UBC. In this example, the references are those given by Eqs. (25)–(26) with *a *= 2.405 rad to exclude the UBC, the two G-LIA outputs provide:

where *s*_{1} corresponds to the amplitude of signal field modulated at Ω_{A}. For a sine phase modulation, an analytical expression can be derived for the proportionality constants *k*_{x} and *k*_{y}.These expressions can be found in the summary table given in Appendix A.

**Figure 11(b)** shows the phase determined with this method when a triangular phase modulation having a peak to peak phase modulation depth of about 2.0 rad is induced by the signal mirror. The signal phase is precisely retrieved.

The value of ψ includes the mechanical phase shift existing between the driving signal and the actual motion of the fork. In fact, in a near-field experiment, this shift can vary from one position to another on the sample depending on the material in interaction with the probe. Depending on the system, the value of ψ in

If ψ is unknown, the G-LIA can be applied twice to solve this issue with two quadrature amplitude modulation functions. In other word, we can calculate for outputs signals

The four outputs (*X*, *Y*, *X*′, *Y*′) can be evaluated numerically or analytically, they provide, respectively, *φ*_{R}(*t*).

## 5. Conclusion

We have detailed the principle of the G-LIA method, first in the case of pure phase modulations then in the case where the amplitude of the signal is also modulated. For pure phase modulations, the interest of the approach was illustrated in different contexts: position monitoring, sensing, and digital holography. In these experiments, the non-linear phase modulation was achieved by mirrors mounted on sinusoidally driven piezo-actuators. In this case, the main advantage of the G-LIA is to extract amplitude and phase information directly from all the harmonic contents created by the phase modulation function. While the examples only considered sine phase modulation functions which is often the most desirable one, the G-LIA also provides a unified treatment to handle arbitrary phase modulation function.

We have also detailed the case where an amplitude modulation can be present. This is notably the case in unbalanced interferometry where a non-negligible amplitude modulation can be perceived at the same frequency than the phase modulation. Experimentally, we considered the case of unbalanced interferometers where a fast sine phase modulation is provided by a current-driven single mode laser diode. A simple yet efficient setup was described to neutralize the impact of wavelength fluctuation on the system. Such approach offers the opportunity to develop simple and cost-efficient system without sacrificing precision. Finally, we discussed the case where the signal of interest is modulated in amplitude at a frequency different from that of the phase modulation. This case was detailed in the context of phase sensitive SNOM, where the low available signal requires to exploit all the available sidebands induced by the phase modulation. Notably, the condition to cancel the effect of the unmodulated background light was presented and attention was paid to the impact of the mechanical phase of the oscillating probe.

**Appendix A**

The table (**Figure 12**) provides a summary of case handled by the G-LIA method.