## 1.

## Introduction

The coherence-controlled holographic microscope (CCHM)^{1}^{,}^{2} is an innovative system particularly designed for quantitative phase imaging (QPI) and measurement of live cell dynamics.^{3} An achromatic off-axis interferometer based on the diffraction grating is used in CCHM. The spatial and temporal coherence of illumination can be widely varied, and in this way, imaging properties of the microscope are substantially modifiable. The coherence gating effect induced by low coherence makes possible formation of optical sections of the sample^{4} (in reflection mode) or imaging through turbid media^{5}^{,}^{6} (in transmission mode). Low coherence also improves lateral resolution and the imaging in general.^{7} For exploiting these effects, CCHM works with a broad polychromatic light source, which makes all the difference to other off-axis holographic microscopes usually equipped with a laser light source.^{8}^{,}^{9} This is because the high-coherence light source leads to the formation of unwanted artifacts in QPI as a consequence of coherence noise, random interferences, and diffraction of light. Low-coherence illumination, however, requires precise alignment and high stability of the system to be maintained during long-lasting time-lapse QPI studies of activity and reaction of living cells. And just these advantages stand for the contribution of CCHM to cell biology research. For alignment of highly sensitive interferometric systems, a secondary light source is often used. The interferometer state detection is carried out by an auxiliary detector^{10} or by a system of detectors.^{11}12.^{–}^{13} These methods are unsuitable for CCHM, because they would cause photo-toxic effect on living cells,^{10}^{,}^{11}^{,}^{13} and because the alignment of the many additional optical elements would be very complicated.^{12} For this reason, an original method based on the measurement of the modulus $\overline{w}$ of the reconstructed holographic signal was elaborated for assessing the instrument state and guiding the optimization. The basic realignment three-dimensional algorithm (BReTA) method is suitable for CCHM because it does not require additional optical components and light sources. Moreover, it is amenable to full automation because the value of the reconstructed holographic signal is available from the online image processing and the operation can be robotized. Finally, the verification of the BReTA method applicability is presented.

## 2.

## Methods and Materials

## 2.1.

### Optical Setup and Image Processing

The optical assembly of CCHM is shown in Fig. 1. A halogen lamp is used as a light source; its light is guided by an optical fiber to the plane light source. This plane is imaged by the relay lens RL into the rear focal plane of condensers ${\mathrm{C}}_{1}$ and ${\mathrm{C}}_{2}$ so that the sample SP and the reference object RO are Köhler illuminated. The light beam is split by beam splitter ${\mathrm{BS}}_{1}$ and mirrors ${\mathrm{M}}_{1}$, ${\mathrm{M}}_{4}$, and ${\mathrm{M}}_{5}$ to the optical paths of object and reference arm. Changeable aperture stop AS is placed into the beam for setting spatial coherence of the light and bandpass filter F for setting its temporal coherence. The specimen and the reference object are imaged from the object plane to the output plane OP by the objective lens ${\mathrm{O}}_{1}$ and ${\mathrm{O}}_{2}$, tube lens ${\mathrm{TL}}_{1}$ and ${\mathrm{TL}}_{2}$, and output lens ${\mathrm{OL}}_{1}$ and ${\mathrm{OL}}_{2}$, respectively. The hologram is captured by the CCD camera D in the output plane OP. Diffraction grating DG and mirror ${\mathrm{M}}_{2}$ are placed in the intermediate image planes behind the tube lenses ${\mathrm{TL}}_{1}$ and ${\mathrm{TL}}_{2}$ to form an interference pattern with the same spatial frequency of fringes (by the same carrier frequency ${f}_{c}$) for all wavelengths of light and for all points of the extended source LS in OP.^{2} In principle, the system provides off-axis holography by the interference of the zeroth and first diffraction orders. Achromaticity and spatial invariance of the interference pattern follow from the fact that the diffraction grating is imaged onto output plane OP by ${\mathrm{OL}}_{2}$.^{14}

The microscope is capable of imaging at various depths in a sample. Two approaches are available depending on the coherence of illumination.^{15}^{,}^{16} Full numerical refocusing is possible with highly coherent illumination, while with partially coherent illumination, the axial range of the refocusing gets reduced.^{15}^{,}^{16} On the other hand, low spatial coherence leads to the formation of an optical section in a broader sense,^{15} the thickness of which corresponds to the depth of field of the objective.^{17} Then, refocusing outside the section is only possible optically. This imaging mode is suitable for observation of adhering cells on a coverslip.

The hologram captured by the CCD is processed by the computer in the following way, see Fig. 2. First, the Fourier transform of the hologram is computed. Then, the spatial-frequency spectrum of the object image is obtained by the windowing operation on the sideband of the spatial-frequency spectrum of the hologram (around the carrier frequency ${f}_{c}$). The zeroth spatial frequency is shifted to the center of the window and the spectrum is multiplied by the Hanning weight function. Finally, the image complex amplitude ${w}_{D}(i,j)$ is computed by the inverse Fourier transform, which provides the image phase $\mathrm{Arg}[{w}_{D}(i,j)]$ and the image amplitude $|{w}_{D}(i,j)|$. A quantitative phase image is obtained after applying an unwrapping procedure and background compensation.^{18}^{,}^{19}

## 2.2.

### Theory

Spatially incoherent illumination requires correct lateral alignment of the images in both arms of the microscope. Temporally incoherent light demands adjusting the arms to the same length. In order to prove that the modulus of the reconstructed holographic signal is of significant value for the microscope alignment, we had to describe its dependence on the optical path length difference, $\mathrm{\Delta}L$, of the arms and on the mutual displacement of the images in the output plane OP that is described by the corresponding displacement vector ${\mathit{q}}_{f}=({x}_{f},{y}_{f})$ in the object plane. The theoretical reconstructed holographic signal can be described by the following formula [see Eq. (3.20) in Ref. 15 for zero defocus]:

## (1)

$$w({\mathit{q}}_{t})={\iint}_{-\infty}^{\infty}{T}_{t}({\mathit{Q}}_{t}){H}_{t}({\mathit{Q}}_{t})\mathrm{exp}(2\pi \mathrm{i}{\mathit{Q}}_{t}\xb7{\mathit{q}}_{t}){\mathrm{d}}^{2}{\mathit{Q}}_{t},$$## (2)

$$w({\mathit{q}}_{t})={K}^{-2}{H}_{t}(\mathrm{0,0})={K}^{-2}{\iint}_{-\infty}^{\infty}{S}_{t}^{*}({\mathit{K}}_{t};K){P}_{t}({\mathit{K}}_{t};K){\mathrm{d}}^{2}{\mathit{K}}_{t},$$^{20}with the support of the radius 1, ${S}_{t}({\mathit{K}}_{t};K)={i}_{t}({\mathit{K}}_{t})\mathrm{circ}[({K}_{t}n)/(K{\mathrm{NA}}_{s})]$ is the 2-D effective pupil function of the illumination, where ${\mathrm{NA}}_{s}$ is the lowest of the numerical apertures of the condenser lenses and the objective lens ${\mathrm{O}}_{2}$ in the reference arm, the function ${i}_{t}({\mathit{K}}_{t})$ describes the distribution of the light intensity in the plane of the light source LS for the Köhler illumination,

^{15}${\mathit{K}}_{t}$ is the transverse part of $\mathit{K}$, ${K}_{t}=|{\mathit{K}}_{t}|$, and $|\mathit{K}|=K=1/\lambda $, $\lambda $ is the wavelength and $n$ is the refraction index in the object space of the objective lenses ${\mathrm{O}}_{1}$ and ${\mathrm{O}}_{2}$. Because ${\mathrm{NA}}_{s}\le \mathrm{NA}$, the function ${P}_{t}({\mathit{K}}_{t};K)$ can be removed from the integrand in Eq. (2). If the image fields of both arms are mutually shifted by the nonzero displacement vector ${\mathit{q}}_{f}=({x}_{f},{y}_{f})$, the function $\mathrm{exp}(-2\pi \mathrm{i}{\mathit{K}}_{t}\xb7{\mathit{q}}_{f})$ must be added to Eq. (2) [see Eq. (2.25) in Ref. 15]. Then,

## (3)

$$w({\mathit{q}}_{t};{\mathit{q}}_{f})={K}^{-2}{\iint}_{-\infty}^{\infty}{i}_{t}({\mathit{K}}_{t})\mathrm{circ}\left(\frac{{K}_{t}n}{K{\mathrm{NA}}_{s}}\right)\mathrm{exp}(-2\pi \mathrm{i}{\mathit{K}}_{t}\xb7{\mathit{q}}_{f}){\mathrm{d}}^{2}{\mathit{K}}_{t}.$$## (4)

$$w({\mathit{q}}_{t};{\mathit{q}}_{f})=2\pi {\int}_{0}^{{\mathrm{NA}}_{s}/n}{i}_{T}(\kappa ){J}_{0}(2\pi K\mathit{\kappa}\xb7{\mathit{q}}_{f})\kappa \mathrm{d}\kappa ,$$Equation (4) describes the reconstructed holographic signal for a spatially incoherent source characterized by the intensity distribution ${i}_{T}(\kappa )$. To also take into account the limited temporal coherence, suppose that the spectral properties of the source are described by the function ${i}_{K}(K)$, and the difference of the optical length of the arms in the image space is $\mathrm{\Delta}L$. The phase shift of the beam inclined by $\psi $ is then $\mathrm{\Delta}{\mathrm{\Phi}}_{\psi}=2\pi K\mathrm{\Delta}L\text{\hspace{0.17em}}\mathrm{cos}\text{\hspace{0.17em}}\psi $ for the refractive index $n\approx 1$. Assuming that the sine condition holds for the microscope objectives, we get $\mathrm{sin}\text{\hspace{0.17em}}\psi ={\mathrm{NA}}_{s}/M\le \mathrm{NA}/M$ for the beam of the maximum inclination $\psi $ in the image space, where $M$ is the magnification between the object plane and the output plane OP. The maximum ratio $\mathrm{NA}/M$ for the objectives used in CCHM has been found for NIKON CFI S Fluor $20\times /0.75$, where $\psi =0.038\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{rad}$. Supposing $\lambda =650\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{nm}$ and maximum $\mathrm{\Delta}L=100\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mu \mathrm{m}$, the difference of the phase shifts for the central beam and the maximum-inclination beam is $\mathrm{\Delta}{\mathrm{\Phi}}_{0}-\mathrm{\Delta}{\mathrm{\Phi}}_{\psi}=2\pi K\mathrm{\Delta}L(1-\mathrm{cos}\text{\hspace{0.17em}}\psi )=0.21\pi \text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{rad}$. Hence, we can approximate the phase shift by $\mathrm{\Delta}{\mathrm{\Phi}}_{0}$ for any inclination of the beam. After completing Eq. (4) by the spectral function ${i}_{K}(K)$ of the source and by the complex exponential depending on $\mathrm{\Delta}{\mathrm{\Phi}}_{0}$ and by the integration over $K$, we get

## (5)

$$w({\mathit{q}}_{t};{\mathit{q}}_{f},\mathrm{\Delta}L)=2\pi {\int}_{0}^{\infty}{i}_{K}(K)\mathrm{exp}(-2\pi \mathrm{i}K\mathrm{\Delta}L){\int}_{0}^{{\mathrm{NA}}_{s}}{i}_{T}(\kappa ){J}_{0}(2\pi K\mathit{\kappa}\xb7{\mathit{q}}_{f})\kappa \mathrm{d}\kappa \mathrm{d}K.$$## 2.3.

### Experiment

For verification of the BReTA methodology, the modulus $\overline{w}$ of the theoretically reconstructed holographic signal described by Eq. (5) was compared with experimental data. No specimen was inserted on the object plane. The measurement was performed with a spatially broad and spectrally narrowband light source. The source spectral function ${i}_{K}(K)$ was given by the manufacturer data (Thorlabs) of the interference filter FB650-10 ($\lambda =650\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{nm}$, $\mathrm{FWHM}=10\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{nm}$); its radial intensity distribution ${i}_{T}(\kappa )$ was approximated by a Gaussian distribution^{21} with the proportional reciprocal standard deviation $\alpha =2.5$ that was found by fitting the theoretical curve to the experimental data. To eliminate the noise, the modulus of the measured reconstructed holographic signal was averaged over the whole image field as follows:

## (6)

$${\overline{w}}_{D}=\sum _{i=0}^{{N}_{I}-1}\sum _{j=0}^{{N}_{J}-1}\frac{|{w}_{D}(i,j)|}{{N}_{I}{N}_{J}},$$Figure 3(a) shows the comparison of dependences of reconstructed holographic signal moduli $\overline{w}$ and ${\overline{w}}_{D}$ on $|{\mathit{q}}_{f}|$. Measured curves slightly differ mutually and from the theoretical curve. It is due to the optical aberrations of the output lens ${\mathrm{OL}}_{2}$, which is located behind the DG in the reference arm. Light diffracted to the first order passes through the border of aperture ${\mathrm{OL}}_{2}$ (in comparison to the axial propagation in the object arm). Moreover, shifting the objective lens ${\mathrm{O}}_{2}$ laterally leads to transversal displacement of the diffracted beam in the aperture ${\mathrm{OL}}_{2}$. In contrast, the measured values of modulus of the reconstructed holographic signal on Fig. 3(b) perfectly fit the theoretical values. Measurement was performed for the same light source and for $|{\mathit{q}}_{f}|=0\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mu \mathrm{m}$ by changing $\mathrm{\Delta}L$.

The theoretical reconstructed holographic signal modulus $\overline{w}$ independence between $|{\mathit{q}}_{f}|$ and $\mathrm{\Delta}L$ is shown in Fig. 3(c). Because of its symmetry, it is displayed only in the first quadrant. Its global maximum ${\overline{w}}_{\mathrm{max}}$ is apparent in the origin. Local maximum ${\overline{w}}_{l,i}$ along the axis $\mathrm{\Delta}L$ is the consequence of the form of the source spectral function ${i}_{K}(K)$. Behavior in the direction of the axis $|{\mathit{q}}_{f}|$ is smooth with no side lobes as a result of the Gaussian form of the intensity ${i}_{T}(\kappa )$. The modulus $\overline{w}$ has the maximum values on the axis $\mathrm{\Delta}L$ for constant $|{\mathit{q}}_{f}|$ and on the axis $|{\mathit{q}}_{f}|$ for constant $\mathrm{\Delta}L$, see also experimental data in Fig. 4.

To demonstrate that the measured reconstructed holographic signal modulus, ${\overline{w}}_{D}$, has a strong maximum at $|{\mathit{q}}_{f}|=0\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mu \mathrm{m}$ and $\mathrm{\Delta}L=0\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mu \mathrm{m}$ even with broadband light source, a measurement was carried out similar to the previous case, but with the narrowband filter removed. Figure 5 compares measurements containing the maximum value ${\overline{w}}_{D,\mathrm{max}}$ for the case of spectrally narrowband and broadband light source. Broadband illumination leads to a broader peak around ${\overline{w}}_{D,\mathrm{max}}$ [Fig. 5(b)] in comparison with spectrally narrowband illumination [Fig. 5(a)]. Extension of the peak is caused by the superposition of peaks related to different wavelengths, which do not overlap ideally due to chromatic aberration of the optical imaging system. Therefore, the measured reconstructed holographic signal modulus ${\overline{w}}_{D}$ is a significant enough value for BReTA method in the case of spectrally narrowband light, as well as broadband light.

## 2.4.

### Basic Realignment Three-Dimensional Algorithm Method

## 2.4.1.

#### Setup adjustment to basic realignment three-dimensional algorithm method

Several essential elements of the microscope setup had to be motorized (see Fig. 1) for experimental implementation of the BReTA method. The XY linear stage of the objective ${\mathrm{O}}_{2}$ was motorized to vary ${\mathit{q}}_{f}$; the holder of the mirror ${\mathrm{M}}_{2}$ was mounted on a linear stage for longitudinal movement $\mathrm{\Delta}L$. A sighting pattern was placed in the plane conjugated with the image plane, and computer-controlled shutters ${\mathrm{S}}_{1}$ and ${\mathrm{S}}_{2}$ were placed in the reference and object arms. The algorithm of the BReTA method described in this part is based on the measurement of the reconstructed holographic signal modulus ${\overline{w}}_{D}$ (hereafter the signal ${\overline{w}}_{D}$), see Fig. 6.

## 2.4.2.

#### Initial alignment

This procedure is used for correction of large misalignment typical for the microscope startup. First, the proper shifts of the mirror ${\mathrm{M}}_{2}$ and objective lens ${\mathrm{O}}_{2}$ are found with the aim to approach the area of the global maximum of ${\overline{w}}_{D}$. The final value of the signal must be higher than the threshold value ${\overline{w}}_{D}^{\prime}=k{\overline{w}}_{D,{l}_{1}}$, where ${\overline{w}}_{D,{l}_{1}}$ is the value of the first side maximum and $k$ is the safety factor. In the beginning of the process, the sighting pattern P is inserted into the field plane and imaged to the object plane (the plane of the specimen SP and of the reference object RO) and finally to the output plane OP. The centers of images are mutually shifted by

## (7)

$${\mathit{q}}_{f}^{\prime}=({x}_{{\mathrm{P}}_{\mathrm{r}}}^{\prime}-{x}_{{\mathrm{P}}_{\mathrm{o}}}^{\prime},{y}_{{\mathrm{P}}_{\mathrm{r}}}^{\prime}-{y}_{{\mathrm{P}}_{\mathrm{o}}}^{\prime}).$$## (8)

$${\mathit{q}}_{f}=\frac{{\mathit{q}}_{f}^{\prime}}{M}=\frac{({x}_{{\mathrm{P}}_{\mathrm{r}}}^{\prime}-{x}_{{\mathrm{P}}_{\mathrm{o}}}^{\prime},{y}_{{\mathrm{P}}_{\mathrm{r}}}^{\prime}-{y}_{{\mathrm{P}}_{\mathrm{o}}}^{\prime})}{M},$$After this process, the value of the signal ${\overline{w}}_{D}$ is tested. If ${\overline{w}}_{D}$ is greater than ${\overline{w}}_{D}^{\prime}$, the process ends. If not, $\mathrm{\Delta}L$ is changed in the direction of increasing signal value ${\overline{w}}_{D}$ by moving the mirror ${\mathrm{M}}_{2}$, until the condition is met.

## 2.4.3.

#### Advanced alignment procedure

In this procedure, the elements ${\mathrm{O}}_{2}$ and ${\mathrm{M}}_{2}$ are moved sequentially to search for the maximum value of signal ${\overline{w}}_{D,\mathrm{max}}$. Because the area of the global maximum is smooth, it is possible to choose any sequence of independent processes for search of the value ${\overline{w}}_{D,\mathrm{max}}$, which are described in this section.

The first process deals with searching the greatest value of signal ${\overline{w}}_{D}$ by changing the microscope objective ${\mathrm{O}}_{2}$ position while maintaining $\mathrm{\Delta}L$ constant. It can be performed in many ways. For example, by a 2-D scanning of objective lens ${\mathrm{O}}_{2}$ around its current position; the resulting scans can be seen in Fig. 7(b). This is a very easy and robust way. Another faster possibility is to use a heuristic algorithm with appropriate termination condition. Algorithm of this process is illustrated in Fig. 7(a). It determines the value of signal ${\overline{w}}_{D,1}$ in the initial position of ${\mathrm{O}}_{2}$. Its position ${\mathit{q}}_{f}$ is then changed by a small defined step in any direction, and the obtained value of signal ${\overline{w}}_{D,2}$ is compared with ${\overline{w}}_{D,1}$. If ${\overline{w}}_{D,2}<{\overline{w}}_{D,1}$ is true, algorithm executes the next step right toward the original step direction; otherwise, the step is performed in the original direction. If the exit condition is not true, the value of signal ${\overline{w}}_{D,3}$ is obtained and again compared with the previous value of the signal ${\overline{w}}_{D,2}$. The termination condition can be selected according to the situation, whichever occurs first. Optionally, the combination of the conditions may be required. The termination can be based on multiple detection of the greatest signal ${\overline{w}}_{D}$ in the same position, when the algorithm is repeatedly going through the same coordinates, or on the calculation of the variational coefficient from the last several values of signal ${\overline{w}}_{D}$. The variation coefficient must be close to one.

The second process finds the value of the signal ${\overline{w}}_{D,\mathrm{max}}$ by moving the mirror ${\mathrm{M}}_{2}$ in the direction of the increased values of the signal ${\overline{w}}_{D}$ [see Fig. 7(b)]. It seeks such a position of mirror ${\mathrm{M}}_{2}$ that corresponds to $\mathrm{\Delta}L=0\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mu \mathrm{m}$.

Both processes together lead to the alignment defined by $\mathrm{\Delta}L=0\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mu \mathrm{m}$ and $|{\mathit{q}}_{f}|=0\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mu \mathrm{m}$.

## 2.4.4.

#### Testing of both procedures

For testing the procedures, Nikon objective lenses $10\times /0.30$ were used. Broadband light from the source was filtered by the interference filter ($\lambda =650\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{nm}$, $\mathrm{FWHM}=10\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mathrm{nm}$). Holograms were captured by camera XIMEA MR4021MC-BH. The choice of the minimum alignment steps in the lateral (${\mathit{q}}_{f}$) and axial ($\mathrm{\Delta}L$) directions depends on the targeted signal level, according to the graph in Fig. 3(c). As it will be shown in Sec. 2.4.5, the sufficient signal level resulting in a good-quality QPI is $\sim 90\%$. The corresponding acceptable misalignment is $0.5\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mu \mathrm{m}$ in ${\mathit{q}}_{f}$ and $10\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mu \mathrm{m}$ in $\mathrm{\Delta}L$ [see Fig. 3(c)]. Hence, for the initial alignment, we chose the step $10\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mu \mathrm{m}$ of $\mathrm{\Delta}L$, while ${\mathit{q}}_{f}$ is set approximately to zero by overlapping the sighting pattern images formed in reference and object arm. The robust scanning method runs laterally within $5\times 5$ positions with the steps $0.5\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mu \mathrm{m}$ of ${x}_{f}$ and ${y}_{f}$, while the axial step of $\mathrm{\Delta}L$ could be fined down to $1\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mu \mathrm{m}$. As the movement of the mirror ${\mathrm{M}}_{2}$ is linear and unidirectional, this correction did not noticeably prolong the axial scanning, in contrast to the situation in lateral directions.

The testing method consists of repeated misalignment of the microscope and subsequent activation of both alignment procedures. The microscope was randomly misaligned by the shift of the objective lens ${\mathrm{O}}_{2}$ and by moving the mirror ${\mathrm{M}}_{2}$ by a distance $100\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mu \mathrm{m}$. Twenty independent measurements were performed. The average time required to approach the aligned state of the microscope was 136 s. Figure 8 shows the results of automated alignment in comparison with the manually obtained (i.e., the most precise) result. The relative values of signal ${\overline{w}}_{D,s}=({\overline{w}}_{D,\mathrm{max}}/{\overline{w}}_{D,\mathrm{max},m})\xb7100\%$, where ${\overline{w}}_{D,\mathrm{max},m}$ is the value of the signal obtained by manual alignment, are on the $x$ axis. The $y$ axis is the count axis.

It is obvious that tested procedures of the BReTA method always resulted without fail with a high success rate.

## 2.4.5.

#### Alignment procedure for a long-term experiment

During the long-term experiments, the value of signal ${\overline{w}}_{D}$ is influenced by temperature changes or vibrations. Therefore, we extended the BReTA method for long-term maintenance of the maximum value of the signal ${\overline{w}}_{D}$. The procedure periodically activates both the transverse stages for the movement of the microscope lens ${\mathrm{O}}_{2}$ and the axial stage to move the mirror ${\mathrm{M}}_{2}$. They are moved symmetrically around the current positions of active elements ${\mathrm{O}}_{2}$ and ${\mathrm{M}}_{2}$ by given steps. The procedure compares the original value of signal ${\overline{w}}_{D,\mathrm{max}}$ and values ${\overline{w}}_{D}$ in adjacent positions of an active element. The active element is always moved to the position with the maximum value of signal ${\overline{w}}_{D}$. In Fig. 9(a), the time dependence of the signal values ${\overline{w}}_{D}$ during the run of the procedure is shown. The procedure changes the positions of the microscope objective ${\mathrm{O}}_{2}$ with the step 100 nm of ${x}_{f}$ and ${y}_{f}$ and the optical path difference $\mathrm{\Delta}L$ of arms (the position of the mirror ${\mathrm{M}}_{2}$) with the step 200 nm. As the procedure is running during observations, the minimum repeatable steps of the linear stages have been chosen to avoid visible changes of the image recorded. Figure 9(b) shows the time dependence of the signal ${\overline{w}}_{D}$ without the alignment procedure running. The signal is inconsistent and too low most of the time. The decrease of the signal makes the phase noise higher, thereby impairing the QPI resolution. We define the QPI resolution as $r=2{\sigma}_{\varphi}$, where ${\sigma}_{\varphi}$ is the standard deviation of the background phase values measured in the window drawn in Fig. 9(c). The decrease of the signal to 90% has no measurable effect on the resolution; for its lower values, the resolution is elevated measurably [see Fig. 9(b)]. The increase of the noise in the phase image of real object is demonstrated in Figs. 9(d) and 9(e).

The microscope was subjected to temperature changes of environment with fluctuations of tenths of centigrade. It is obvious that the procedure successfully maintained the aligned state of the microscope.

## 3.

## Conclusions

We have developed a unique set of procedures constituting the original BReTA. The BReTA method allows for automated alignment of CCHM based on maximizing the value of modulus $\overline{w}$ of the measured holographic signal. For exerting control over this parameter, some alignment elements of the original CCHM setup had to be motorized. The method consists in the initial alignment of the microscope in order to find the required minimal interference signal. The holographic signal is then optimized by searching the best alignment corresponding to the maximum signal. This maximum can be subsequently maintained by small changes in the alignment during a long-term experiment. All procedures were programmed in LabView and C++, and they are being used in the multimodal holographic microscopes produced by TESCAN ORSAY HOLDING a.s. The automatic BReTA method offers easy alignment of the microscope for ordinary users. However, it is essential for managing long-term QPI observations/measurements of live cells activity, which is the primary assignment of holographic microscopy. The BReTA method described in this article is patent pending.^{22}

## Acknowledgments

The authors thank their colleagues from the Experimental Biophotonics Group (CEITEC) for helpful discussions, especially Jana Collakova, Pavel Vesely, and Vera Kollarova. This work was supported by CEITEC-Central European Institute of Technology (CZ.1.05/1.1.00/02.0068) from the European Regional Development Fund.

## References

P. Kolman and R. Chmelik, “Coherence-controlled holographic microscope,” Opt. Express 18, 21990–22003 (2010).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.18.021990Google Scholar

T. Slaby et al., “Off-axis setup taking full advantage of incoherent illumination in coherence-controlled holographic microscope,” Opt. Express 21, 14747–14762 (2013).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.21.014747Google Scholar

M. Mir et al., “Progress in optics,” in Quantitative Phase Imaging, and E. Wolf, Ed., pp. 133–217, Elsevier, Amsterdam (2012).Google Scholar

R. Chmelik and Z. Harna, “Parallel-mode confocal microscope,” Opt. Eng. 38(10), 1635–1639 (1999).http://dx.doi.org/10.1117/1.602217Google Scholar

M. Lostak et al., “Diffuse light imaging with a coherence controlled holographic microscope,” Proc. SPIE 7746, 77461N (2010).PSISDG0277-786Xhttp://dx.doi.org/10.1117/12.882198Google Scholar

M. Lostak et al., “Coherence-controlled holographic microscopy in diffuse media,” Opt. Express 22, 4180–4195 (2014).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.22.004180Google Scholar

L. Lovicar, J. Komrska and R. Chmelík, “Quantitative-phase-contrast imaging of a two-level surface described as a 2D linear filtering process,” Opt. Express 18, 20585–20594 (2010).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.18.020585Google Scholar

N. T. Shaked et al., “Quantitative analysis of biological cells using digital holographic microscopy,”in Holography, Research and Technologies, and J. Rosen, Ed., pp. 219–236, InTech, Rijeka (2011).Google Scholar

F. Charrie et al., “Cell refractive index tomography by digital holographic microscopy,” Opt. Lett. 31, 178–180 (2006).OPLEDP0146-9592http://dx.doi.org/10.1364/OL.31.000178Google Scholar

H. Iwai et al., “Quantitative phase imaging using actively stabilized phase-shifting low-coherence interferometry,” Opt. Lett. 29, 2399–2401 (2004).OPLEDP0146-9592http://dx.doi.org/10.1364/OL.29.002399Google Scholar

D. Kalamatianos et al., “Active alignment for two-beam interferometers,” Rev. Sci. Instrum. 77(1), 013103 (2006).RSINAK0034-6748http://dx.doi.org/10.1063/1.2150823Google Scholar

R. Schwiesow, “Stabilization mechanism for optical interferometer,” US Patent 4, 444, 501 (1984).Google Scholar

V. V. Krishnamachari et al., “An active interferometer-stabilization scheme with linear phase control,” Opt. Express 14, 5210–5215 (2006).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.14.005210Google Scholar

E. N. Leith and J. Upatnikes, “Holography with achromatic-fringe systems,” J. Opt. Soc. Am. 57, 975–979 (1967).JOSAAH0030-3941http://dx.doi.org/10.1364/JOSA.57.000975Google Scholar

R. Chmelik et al., “Progress in optics,” in The Role of Coherence in Image Formation in Holographic Microscopy, and E. Wolf, Ed., pp. 133–217, Elsevier, Amsterdam (2014).Google Scholar

M. Tyc et al., “Numerical refocusing in digital holographic microscopy with extended-sources illumination,” Opt. Express 21, 28258–28271 (2013).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.21.028258Google Scholar

J. Collakova et al., “Coherence-controlled holographic microscopy enabled recognition of necrosis as the mechanism of cancer cells death after exposure to cytopathic turbid emulsion,” J. Biomed. Opt. 20(11), 111213 (2015).JBOPFO1083-3668http://dx.doi.org/10.1117/1.JBO.20.11.111213Google Scholar

R. M. Goldstein, H. A. Zebken and C. L. Werner, “Satellite radar interferometry: two-dimensional phase unwrapping,” Radio Sci. 23(4), 713–720 (1988).RASCAD0048-6604http://dx.doi.org/10.1029/RS023i004p00713Google Scholar

D. C. Ghiglia and M. D. Pritt, Two-Dimensional Phase Unwrapping: Theory, Algorithms and Software, Wiley-Interscience, New York (1998).Google Scholar

J. W. Goodman, Introduction to Fourier Optics, 3rd ed., Roberts & Company, Greenwood Village (2005).Google Scholar

F. Harris, “On the use of windows for harmonic analysis with the discrete Fourier transform,” Proc. IEEE 66, 51–83 (1978).IEEPAD0018-9219http://dx.doi.org/10.1109/PROC.1978.10837Google Scholar

R. Chmelik et al., “Interferometric system with variable optics for incoherent source of radiation,” Czech Republic Patent Application No. PV 2014-538 (2014).Google Scholar

## Biography

**Zbynek Dostal** is a PhD student at the Brno University of Technology. He received his BS degree in mechanical engineering in 2007 and MS degree in optics and precise mechanics in 2009 from the Brno University of Technology. He is a coauthor of the coherence-controlled holographic microscope patent for which he received the Werner von Siemens Excellence Award 2013. His current research interests include holographic microscopy, automation in microscopy, and optical and mechanical design.

**Tomas Slaby** is head of light microscopy R&D in TESCAN Brno company. He received his MS and PhD degrees in holographic microscopy from the Brno University of Technology in 2008 and 2015. He is involved in the development of a new generation of coherence-controlled holographic microscope.

**Lukas Kvasnica** is a PhD student at the Brno University of Technology. He received his MS degree in optics and precise mechanics in 2008 from the Brno University of Technology. Since 2012, he has been working as a C++ programmer in TESCAN Brno company, where he is involved in the development of the coherence-controlled holographic microscope.

**Martin Lostak** is a physicist in TESCAN Brno company. He received his MS and PhD degrees in holographic microscopy from the Brno University of Technology in 2008 and 2015. He is involved in developing the new generation of coherence-controlled holographic microscope.

**Aneta Krizova** is a PhD student at the Brno University of Technology. She received her BS and MS degrees in physical engineering from the Brno University of Technology in 2010 and 2012. Her current research interests include light microscopy, especially holographic microscopy and its applications.

**Radim Chmelik** is a professor of applied physics at the Brno University of Technology. He received his MS degree in solid-state physics from the Masaryk University in Brno in 1989 and his PhD degree in physical and materials engineering from the Brno University of Technology in 1997. He is the author of more than 40 journal papers. His current research interests include wave optics, imaging theory, advanced and three-dimensional light microscopy, and holographic microscopy.