Next Article in Journal
Formation of a 3D Particle Array Actuated by Ultrasonic Traveling Waves in a Regular Polygon Resonator
Next Article in Special Issue
Propagation Losses Estimation in a Cationic-Network-Based Hydrogel Waveguide
Previous Article in Journal
Research on Intelligent Robot Point Cloud Grasping in Internet of Things
Previous Article in Special Issue
Validation of a Hyperspectral Imaging System for Color Measurement of In-Vivo Dental Structures
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Contribution to the Improvement of the Correlation Filter Method for Modal Analysis with a Spatial Light Modulator

by
David Benedicto
1,
María Victoria Collados
1,*,
Juan C. Martín
1,
Jesús Atencia
1,
Omel Mendoza-Yero
2 and
Juan A. Vallés
1
1
Departamento de Física Aplicada, Instituto de Investigación en Ingeniería de Aragón (I3A), Facultad de Ciencias, Universidad de Zaragoza, Pedro Cerbuna 12, 50009 Zaragoza, Spain
2
Institut de Noves Tecnologies de la Imatge (INIT), Universitat Jaume I, 12080 Castelló, Spain
*
Author to whom correspondence should be addressed.
Micromachines 2022, 13(11), 2004; https://doi.org/10.3390/mi13112004
Submission received: 14 October 2022 / Revised: 11 November 2022 / Accepted: 16 November 2022 / Published: 17 November 2022

Abstract

:
Modal decomposition of light is essential to study its propagation properties in waveguides and photonic devices. Modal analysis can be carried out by implementing a computer-generated hologram acting as a match filter in a spatial light modulator. In this work, a series of aspects to be taken into account in order to get the most out of this method are presented, aiming to provide useful operational procedures. First of all, a method for filter size adjustment based on the standard fiber LP-mode symmetry is presented. The influence of the mode normalization in the complex amplitude encoding-inherent noise is then investigated. Finally, a robust method to measure the phase difference between modes is proposed. These procedures are tested by wavefront reconstruction in a conventional few-mode fiber.

1. Introduction

Characterization of optical fields by means of modal decomposition (MD) is key for the analysis, design, and optimization of multimode waveguide-based devices [1,2,3]. Interest in these techniques is increasing in parallel with the growing research on applications based on multimode structures such as high-power large-mode-area fiber lasers [4,5] and space division multiplexing techniques, aimed at widening the channel capacity of optical communication systems [6,7,8]. Modal decomposition facilitates the study of modal competition [9,10], modal instabilities limiting the maximum power emission in fiber lasers [11,12], nonlinear optical effects where the interaction of many transversal modes is involved [13,14,15], and the design of optical devices based on optical power transfer [16,17] and selective mode excitation [18,19]. This analysis has proven to be a key tool to study the transmission properties of waveguides regarding its modal behavior due to the physical insight that can be obtained [20].
In general, MD methods can be classified as numerical or experimental. On one hand, several numerical MD algorithms based on measured optical intensity distributions have been reported [21]. Iterative methods such as the GS algorithm [22], the SPGD algorithm [23], and the genetic algorithm [24] show a high accuracy but are highly sensitive to the initial values due to the local minima problem and, in some cases, the necessary iterative process can be computationally intensive [25]. The newly emerged neural network methods [26,27,28] outperform the iterative methods in decomposition speed without having the local minima problems and showing a highly accurate performance, but require high-performance computers, a large amount of memory, and a long time for training the neural network. Recently, a non-iterative method based on a sophisticated analytical model has been published [29], achieving a very fast MD performance and without any optimization or training process, as long as the signal-to-noise ratio of the acquired images remains particularly low.
On the other hand, MD can be directly performed through experiments, such as by the use of ring-resonators [30], the evaluation of multimode interference commonly known as S2 [31] (spatially and spectrally resolved imaging) and C2 [32] (frequency domain cross-correlated imaging), low-coherence interferometry [33], the application of fiber Bragg gratings [34] or diffraction gratings [35], and the correlation filter method (CFM) [36]. The latter allows us to retrieve the full optical field information by measuring only a few modal amplitudes and phases. Compared to numerical methods, this procedure requires a more complex experimental setup, but it provides a fast and accurate MD without the need for developing any software for the analysis of the intensity distribution images.
In the CFM, the waveguide output beam illuminates a computer-generated hologram (CGH), acting as a match filter, which performs the decomposition of the output field into the waveguide propagation modes by means of a specific previously designed transmittance function. Hence, a priori knowledge of the field distribution of the waveguide-under-study set of modes is required. Therefore, the designed CGH would be limited to analyzing only the corresponding type of waveguide. This limitation can be overcome by implementing the CGH into a spatial light modulator (SLM), allowing the user to adapt the decomposition mode set so that any waveguide can be investigated [37].
Thanks to the recently emerged liquid-crystal SLM technology, the experimental procedure has been significantly simplified, and the CFM has been used for mode analysis [36,37,38], mode-resolved bend loss analysis [39], mode-resolved gain analysis [10], wavefront reconstruction [40], and mode excitation [20]. However, there are only a few works describing the details of this method [41], which could be related both to the novelty of the technology and to the complexity of the method compared to other MD approaches. Taking this into account, we provide in this work a series of operational procedures, not detailed to date in any paper as far as we know, in order to carry out this method correctly and to get the most out of it. This way, part of its complexity becomes reduced.
Correlation filters require light amplitude and phase modulation. Nevertheless, commercial SLMs can only induce a phase shift or an amplitude modulation. Both phase and amplitude modulation can be achieved with a phase-only SLM by conveniently grouping the SLMs pixels into what have been called macropixels [42,43,44,45,46], thus allowing the correlation filter method implementation. However, this technique gives rise to an output field presenting a noise term [45], which not only depends on the filter but also on the input-field distribution. In this work, we study the effect of the mentioned noise term, inherent to the use of macropixels when an SLM is employed in an MD setup.
While implementing the CFM, we have studied some factors that can condition and worsen its performance (e.g., system alignment, CGH adjustment, a priori mode computation, laser instabilities), for which we offer a set of techniques that can be useful in order to lessen its influence. On the one hand, we propose a CGH-scale-adjustment technique based on the waveguide mode symmetry in order to easily adjust the match-filter size. On the other hand, a different way of measuring the phase difference between modes is presented, reducing the system instabilities effects and improving the MD performance.
The manuscript is organized as follows. For the sake of completeness, Section 2 is devoted to revise the MD technique based on the CFM implemented in a phase-only SLM. In Section 3, the experimental setup is explained. Section 4 contains the mentioned noise-term analysis, inherent to the necessary encoding technique, while Section 5 and Section 6 include the proposed scale adjustment and phase-retrieval techniques, respectively. Finally, in Section 7 we test the experimental performance of the modal-analysis procedure by a wavefront reconstruction, which leads to the conclusions in Section 8.

2. The Correlation Filter Method for Modal Decomposition

2.1. Basics of the Correlation Filter Method

The notation to be employed is established in Figure 1, which shows a basic scheme of the setup necessary for modal decomposition using the CFM. Lenses L1 and L2, positioned following a 4f configuration and act as a beam expansor, where the fiber output end is placed at the focus of lens L1; the distance between lenses is the sum of their focal distances and the CGH is located at the focal plane of lens L2. Through this beam expansor, the fiber output electric field U ( ε , μ ) is imaged at the CGH plane, U ( x , y ) . T ( x , y ) is the CGH transmittance implemented in the SLM, W 0 ( x , y ) is the electric field after the SLM, and W f ( x , y ) is the field at the CCD camera position, which is placed at the focal plane of lens L3.
The multimode waveguide output field to be analyzed constitutes the input CGH field
U ( x , y ) = n = 1 N c n ϕ n ( x , y ) .
where x and y are orthogonal coordinates in the CGH plane, N is the total number of propagation modes allowed by the waveguide under study [47], ϕ n ( x , y ) is the nth-mode normalized distribution of the transversal electric field, and c n is its complex expansion coefficient, which can be expressed in terms of a modal weight, c n , and a relative phase term, φ n ,
c n = c n e i φ n ,
where i is the imaginary unit and where all mode relative phases are expressed with respect to the 0-mode phase ( φ 0 = 0 ) .
The electric fields W 0 ( x , y ) and U ( x , y ) are related by:
W 0 ( x , y ) = U ( x , y ) T ( x , y ) .
The Fourier transform of the output distribution, W f ( u , v ) , is obtained at the L3 lens focal plane. Specifically, W f and W 0 are related by means of [48]:
W f ( u , v ) = 1 i λ f [ W 0 ( x , y ) ] ,
where denotes Fourier transform, λ is the input-light wavelength, f is the lens focal distance, and the coordinates at the focal plane ( u , v ) are related to the Fourier-transform frequency-space coordinates f x and f y by means of u = λ ff x , v = λ ff y . In order to measure the amplitude of mode p, a filter with the following transmittance should be employed [36]:
T p ( x , y ) = ϕ p ( x , y ) ,
where the asterisk denotes the complex conjugate. After some operations and taking into account orthogonality between mode field distributions [36]:
W f ( u = 0 , v = 0 ) = 1 i λ f n = 1 N c n δ pn   = 1 i λ f c p ,
where δ pn is the Kronecker delta. If a CCD detector is placed on axis in the focal plane of L3, the use of a filter, as given by Equation (5), allows one to obtain information regarding the p-mode amplitude c p at the specific focal plane coordinates ( u = 0 ,   v = 0 ) , the intensity at this point being proportional to the modal weight: I ( 0 , 0 ) =   W f ( 0 , 0 ) 2 c p 2 .
A simultaneous procedure based on angular multiplexing can be carried out. It consists of employing a filter whose transmittance is a superposition of different Tp filters, associating to each p-mode a different wave vector ( k x , p , k y , p ) , so that the total transmittance function is:
T ( x , y ) = p = 1 N ϕ p ( x , y ) · e i ( xk x , p + yk y , p ) .
By a convenient choice of these wave vectors, information regarding the different mode modal weights appears at separated-enough points in the focal plane, u p = λ fk x , p and v p = λ fk y , p being the coordinates of each one of them.
For some applications, it is enough to obtain the modal weights, c p 2 . However, sometimes it is required to determine the phase difference between modes. According to [36], in order to measure it, it is necessary to use two transmittance functions, defined as
T p cos ( x , y ) = 1 2 [ ϕ 0 ( x , y ) + ϕ p ( x , y ) ] · e i ( xk x + yk y )
and
T p sin ( x , y ) = 1 2 [ ϕ 0 ( x , y ) + i ϕ p ( x , y ) ] · e i ( xk x + yk y ) ,
so that the p-mode relative phase, φ p , can be determined with respect to the 0-mode one. In these cases, the intensity at the CCD plane is proportional to the cosine and sine, respectively, of the phase difference [36]. This allows one to obtain an unambiguous solution for the phase differences.

2.2. Double-Phase Method for Complex Amplitude Encoding

When implementing the CGH in a phase SLM, all the transmittance functions defined in the previous subsection require the encoding of a complex amplitude in an only-phase device. Among the different proposed procedures for pixel grouping, the double-phase CGHs have been chosen [45] due to their simple implementation and maximum resolution, as in this technique macropixels are composed of only two pixels. This method is based on the principle that any complex value inside the unit circle can be resolved into the sum of two complex values with unit modules.
Consider a CGH implemented in a commercial SLM that only admits phase modulation of its pixels, and a combination of them in macropixels that may offer an approximate complex modulation. Specifically, as explained in [45], a CGH filter with an approximate transmittance to the theoretical one can be achieved by means of grouping pairs of pixels into a macropixel. Considering that each couple of consecutive pixels constitute the macropixel with indexes ( m ,   n ) , if the desired complex transmittance of the macropixel ( m ,   n ) is T mn = T mn e i α mn , the following respective phase shifts are assigned to the two consecutive pixels [45]:
α mn ( 1 ) = α mn     Δ mn   and   α mn ( 2 ) = α mn + Δ mn ,
where T mn     1 , 0     α mn     2 π and
Δ mn = cos 1 ( T mn ) .
Denoting by Q the transmittance function of the phase-only device where the double-phase CGH are implemented, its Fourier transform can be expressed as:
Q = Q S + Q N ,
in such a way that [45]
Q S = [ T ( x , y ) ] .
Contribution Q S , referred to as signal, is the only one that should be expected if the double-phase solution provided an exact T ( x , y ) transmittance. Nevertheless, another contribution to Q also appears, due to the double-phase implementation: Q N , referred to as noise.

3. Modal Analysis Setup

The required experimental setup is depicted in Figure 2. An SMF28 optical fiber supporting six propagation LP-modes is illuminated by a He-Ne laser ( λ = 632 . 8   nm ) . The first section of optical fiber is repeatedly curved in order to eliminate modes other than the fundamental one, LP 01 . We use this fiber to illuminate a second optical fiber with the same characteristics. This way, when both fibers are placed closely one after another over the same longitudinal axis, the LP 01 mode is the only one that propagates along the second section of optical fiber. It can be used to carry out the system alignment. By sideways displacing the first fiber with respect to the second one, different transversal distributions are generated, due to different combinations of excited modes, which is useful for testing the MD performance. Figure 2A,B show two examples: the first one with both fiber sections totally aligned, used to calibrate the setup, and the second one at an arbitrary position, employed to test the proposed technique’s performance.
In order to increase the effective resolution, the fiber output beam size has been magnified at the SLM display using a microscope objective (DIN × 40) and a lens (f1 = 50 cm) following a 4f configuration (theoretical magnification 125). The 4f-lens combination images the fiber output field distribution through the beam splitter (BS) both in the CCD1 camera and in the phase-only SLM operating in reflection mode (PLUTO VIS Phase Only Spatial Light Modulator from Holoeye [49]). As it only operates in the horizontal polarization, we select this state of polarization from the input beam with a linear polarizer (LP). In order to obtain a complete MD, the same procedure should be performed for the vertical polarization. In order to do that, the horizontal LP should be replaced by a vertical one, followed by a half-wave plate. The reflected field from the SLM, where the CGH is implemented, is Fourier transformed by lens L2 (f2 = 15 cm), and this signal is detected by the CCD2 camera, placed at the L2 lens focal plane, where the modal weights are obtained by measuring the intensity in the specified coordinates according to Equation (6). In order to avoid interferences, The CCD2 and the SLM are covered when measuring intensity profiles with the CCD1 camera. The SLM was previously calibrated following the procedure described in [50].

4. Correlation-Filter Size adjustment

The MD performance is based on the overlap in the SLM between the incident light beam and the implemented CGH. Inaccuracies in the position or size of the encoded CGH with respect to the incident light can result in the detection of erroneous modes or in the incorrect determination of the modal weights. Notwithstanding its importance, these issues have not been studied in depth, except for in a recently published tutorial [41]. Their proposed approach to transversally align the CGH position consists of displacing the filter position in such a way that, for a mismatched overlap, the on-axis null appears centered regarding the input light beam. If the first and second optical fiber sections are aligned so as to achieve single LP 01 mode propagation in the second one, there should be zero on-axis intensity when implementing any other mode into the match filter. We can benefit from LP 11 e and LP 11 o mode symmetry in order to center the CGH vertically and horizontally, respectively. The simulated transversal intensity distributions at the CCD2 plane, when the LP 11 e and LP 11 e modes are implemented in the SLM, are shown in Figure 3a,b, respectively, for different cases of transversal centering. As one can see, when the filters are centered with respect to the incident field, there is zero intensity in the optical axis point (marked with a red cross). As the filters decenter, the intensity at the measurement position increases.
We propose that this same idea can be used to adjust the filter size. Despite it being straightforward to calculate the theoretical 4f-magnification, this value can be slightly modified in the experimental implementation, mainly because it is not easy to place the optical fiber output end at the exact focal length distance from the microscope principal plane. In Ref. [51], a method for filter-size adjustment is proposed, for which it is required to obtain some parameters, such as the beam quality factor and the second moment. Here, we present a different approach, based on the LP 02   mode symmetry. If we implement the LP 02 mode into the CGH and illuminate it with the LP 01 mode, only when both their sizes are matched will there be zero intensity at the optical axis, as shown in the simulations from Figure 4.
We can use this behavior to find the correct scale just by minimizing the optical axis intensity as a function of the theoretical magnification value. This is shown experimentally in Figure 5. Overlaying the graph, three measured transversal intensity distributions are shown for a smaller, adjusted, and larger size (from left to right) of the filter with respect to the incident light. The measured intensity distributions follow a similar behavior than the simulated ones in Figure 4.
According to the minimum intensity, the magnification value is 141, which has a 13% relative error with respect to the theoretical one. This error can be explained due to the inherent uncertainties of the optical system positioning. This method can be used as long as the set of modes has the required symmetry, which is common in cylindrical-like structures.

5. Double Phase Method Noise Term Analysis

Consider a double-phase CGH with transmittance Q ( x , y ) (Equation (18)), implemented in order to emulate an ideal CGH filter with complex transmittance T ( x , y ) (Equation (13)). According to Equations (3), (4), and (12), this double-phase filter yields a field distribution in the L2 back focal plane, such as:
W f = 1 i λ f [ U ( Q S + Q N ) ] ,
where U ˜ = [ U ( x , y ) ] and the ( ) sign represents convolution product. We consider two different contributions to W f , called here signal and noise terms, respectively:
W f , S = 1 i λ f [ U Q S ] ,
W f , N = 1 i λ f [ U Q N ] ,
Because Q ˜ S = [ T ( x , y ) ] , W f , S represents the field distribution that would be obtained if an ideal CGH filter with complex transmittance T ( x , y ) was employed. Nevertheless, a superimposed noise term W f , N is also present. As stated in Equation (16), this term contains the convolution of two factors. One of them, Q N , is intrinsic to the filter, while the other, U , is directly related to the input electric field distribution, which is obviously unpredictable as it is the object of analysis. Therefore, the magnitude and distribution of the noise term is expected to be very different depending on each working condition. For this reason, its impact cannot be analyzed by means of a general mathematical treatment.
A case in which the noise factor strongly affects the MD is shown in Figure 6. Suppose an incident beam with all six LP-modes equally present and a CGH in which all the match filters (one for each mode) are simultaneously multiplexed with a different grating to spatially separate the signals at the Fourier plane. Specifically, the vertexes of a regular hexagon, as shown in Figure 6a–c, show the simulated intensity at the L2 back focal plane, separated into the signal term W f , S (Equation (15)) and the noise term W f , N (Equation (16)), respectively.
By comparing the modal weight percentages with and without noise terms, relative errors up to 50% are obtained, as shown in Table 1, thus preventing the correct performance of the MD. The modal weight in percentage and the relative error are computed by c l 2 ( % ) = 100 c l 2 / n = 1 N c n 2 and ε r = 100 ( c Exp 2 c Teo 2 ) / c Teo 2 , respectively, where c Exp 2 and c Teo 2 are the experimental and theoretical modal weights, the last one being 16.67% as all modes are equally present.
In view of the findings, it is essential to reduce the noise term influence on the MD performance. One way to achieve this is to realize the MD sequentially rather than simultaneously. However, it is not possible to obtain all modal weights in real time with this procedure, which can be an important disadvantage depending on the system dynamics. Another approach is to normalize the transmittance function (Equation (7)). Theoretically, any transmittance value between zero and one is valid when implementing the double phase method. Nevertheless, in order to use the whole range of the SLM and reduce the noise impact, the maximum value of the transmittance function amplitude should be as close as possible to one. This affects parameter Δ mn from Equation (11), which is necessary for the implementation of the double phase method. The normalization has a direct impact on the signal to noise ratio, as shown in Figure 7, where we have used the same configuration than the one in Figure 6 and summed the intensities at the six information positions for both the signal (Equation (15)) and the noise (Equation (16)) terms.
As shown in Figure 7, by adjusting the maximum value of the transmittance amplitude to one, we improve the signal to noise ratio, while if we keep the set of modes orthonormality (situation marked with the red arrows), the noise intensity cannot be neglected with respect to the signal, thus affecting the MD. When all the modes are measured using the same transmittance function (simultaneous MD case), the renormalization appears as a constant factor equal for all measured modal weights. However, the renormalization factors must be considered when measuring different modal weights with different transmittance functions in order not to break the mode orthonormality.

6. Robustness Improvement in Relative Phase Measurements

In theory, the phase difference between modes can be unambiguously determined from two measurements by using Equations (8) and (9), each one having two possible solutions. Nevertheless, the mode relative phase measurements can be easily influenced by system instabilities (i.e., laser power fluctuations, micro-positioners looseness, etc.) or by misalignments due to position uncertainties (i.e., lenses, fiber, CCD, SLM). As a consequence, there may be an error in the arccosine and arcsine result that could prevent correct phase retrieval. For such cases, we propose a different transmittance function in order to measure the relative phases and also to study the phase determination uncertainty.
Suppose we implement in the SLM the following transmittance function:
T ( x , y ) = 1 2 [ ϕ 0 ( x , y ) + ϕ l ( x , y ) · e i θ ] ,
as a function of a variable phase θ . The particular cases θ = 0 and θ = π / 2 recover Equations (8) and (9), respectively. The intensity at the CCD2 measurement position as a function of θ is
W f   ( u = λ fk x , v = λ fk y , θ ) 2 c 0 2 + c l 2 + 2 c 0 c l · cos ( φ l + θ ) .
In order to improve the robustness of the relative phase determination, we perform a series of intensity measurements as a function of θ and fit them to the function in Equation (18).
As an example, Figure 8 shows the so mentioned measurements of the phase difference between the LP 21 e and LP 01   modes, when the input distribution is the one shown in Figure 2b, together with its best fit by least squares. It is clear to see that, despite the fact that the experimental dots follow a similar behavior than the cosine function, there are some significant differences with respect to the fitting function. Both instabilities (i.e., laser fluctuations, fiber disturbances) and experimental inaccuracies (i.e., devices alignments) may be responsible for the lack of agreement between the measurements and the theoretical function.
Superimposed to the graph in Figure 8, four points have been highlighted. First, two points, θ = 0 and θ = π / 2 , are the particular cases from Equations (8) and (9). In this case, the points show a good agreement with the fitting curve. In fact, following the two-measurements procedure, a phase difference of 0.23π rad is obtained, which is almost the same value as the one determined through the fitting approach: 0.24π rad. Based on this result, one could think that the fitting approach does not improve the phase determination. However, the relative phase obtained through the two-measurements procedure may depend heavily on the selected points, as one can see in Figure 8. Suppose we choose the particular cases θ = π and θ = 3 π / 2 , whose experimental values are far from the fitting curve, and compute the phase difference. In this case, the obtained angles do not match each other. The obtained phases that fall in the same quadrant are 0 and 0.39π rad, respectively. It is clear to see that, in general, the function fit allows us to obtain a relative phase less influenced by possible instabilities and uncertainties. In the next section, we observe its impact by performing a wavefront reconstruction.

7. Wavefront Reconstruction

In order to test the proposed size adjustment technique and phase retrieval procedure, together with the MD method itself, we have reconstructed the transversal intensity distribution shown in Figure 2b. Figure 9 shows the comparison between the experimental distribution (a) and the numerically reconstructed ones (b–d). The reconstructed distributions have been obtained by measuring all six LP-mode weights and relative phases. However, in the reconstructions shown in Figure 9b,c, the proposed CGH size adjustment and phase retrieval techniques have not been used, respectively. The reconstruction shown in Figure 9b has been obtained by employing the proposed phase retrieval fitting approach but not the CGH size adjustment. Instead, its size has been computed theoretically with the 4f magnification. On the contrary, for the reconstruction shown in Figure 9c, the CGH size adjustment technique has been used, but the relative phases have been determined through the two-measurements method. Finally, the distribution shown in Figure 9d has been computed using both of the proposed techniques. As one can see, a good agreement is shown between Figure 9a,d distributions, highlighting the importance of both the size adjustment and the phase retrieval. The residual intensity patterns, computed where m and r mean measured and reconstructed, are respectively provided, together with each reconstructed image.
By comparing the transversal intensity distributions one can see that the CGH size adjustment is quite critical and that, although a look alike distribution is obtained without determining the relative phases with our proposed procedure, a better reconstruction is obtained when employed. To quantitatively compare the intensity distributions from Figure 9, the root-mean-square error (RMSE) has been computed for each of the reconstructed images with respect to the experimentally registered one, all four of them previously normalized to the unit, obtaining the RMSE values 2.4 × 10−2, 1.8 × 10−2, and 1.7 × 10−3 for the Figure 9b–d intensity distributions, respectively. Distribution from Figure 9d gives a RMSE one order of magnitude lower than the ones from Figure 9b,c, showing the good performance of the reconstruction when using the proposed procedures. Table 2 summarizes the modal weights and relative phases obtained by using the proposed techniques, with which distribution from Figure 9d is computed.

8. Conclusions

We have presented a series of practical procedures to perform an accurate modal decomposition of light, and we tested them satisfactorily by wavefront reconstruction in a few-mode fiber, allowing us to verify both the good performance of the modal analysis procedure and the proposed techniques. In order to reduce the SNR when performing a simultaneous MD and exploit the SLM range, the set of modes orthonormality needs to be broken, and thus a scale factor taken into account. It has been shown that it is possible to benefit from the LP-mode symmetry, not only to transversally center the CGH, but also to adjust its size with respect to the incident beam. This procedure is not restricted to the LP-modes, and it can be used as long as the set of modes provides the necessary symmetry, which is usually fulfilled in cylindrical-like waveguides. When setup instabilities make the correct phase difference determination through the classical approach difficult, a possibility to improve phase retrieval robustness consists of implementing a transmittance as a function of an angle parameter and adjust the set of measurements. The good performance of these techniques has been successfully tested by wavefront reconstruction.

Author Contributions

Conceptualization and investigation, D.B., M.V.C., J.C.M., J.A., and J.A.V.; software development, D.B. and O.M.-Y.; experimental performance, D.B. and M.V.C.; data analysis, D.B., M.V.C., and J.C.M.; writing the manuscript, D.B., M.V.C., and J.C.M.; reviewing the manuscript, J.A. and J.A.V.; project administration and funding acquisition, M.V.C. and J.A.V. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported in part by the Spanish Ministry of Science and Innovation under Grant PID2019−108598GB-I00/AEI/10.13039/501100011033 and in part by the Department of Industry and Innovation (Government of Aragon) through the research group (Cofinanciando con Feder 2014–2020: Construyendo Europa desde Aragón) under Grant E44-20R. This work has also been funded by project PID2019-110927RB-I00, financed by MCIN/AEI/10.13039/501100011033.

Data Availability Statement

Data available on request.

Acknowledgments

The authors are also very grateful to Gladys Mínguez Vega for her help in the understanding and calibration of the spatial light modulator.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bachmann, M.; Besse, P.A.; Melchior, H. General self-imaging properties in N× N multimode interference couplers including phase relations. Appl. Opt. 1994, 33, 3905–3911. [Google Scholar] [CrossRef] [PubMed]
  2. Smith, C.M.; Venkataraman, N.; Gallagher, M.T.; Müller, D.; West, J.A.; Borrelli, N.F.; Allan, D.C.; Koch, K.W. Low-loss hollow-core silica/air photonic bandgap fibre. Nature 2003, 424, 657–659. [Google Scholar] [CrossRef] [PubMed]
  3. Schmidt, O.A.; Schulze, C.; Flamm, D.; Brüning, R.; Kaiser, T.; Schröter, S.; Duparré, M. Real-time determination of laser beam quality by modal decomposition. Opt. Express 2011, 19, 6741–6748. [Google Scholar] [CrossRef] [PubMed]
  4. Ji, J.; Raghuraman, S.; Huang, X.; Zang, J.; Ho, D.; Zhou, Y.; Benudiz, Y.; Ami, U.B.; Ishaaya, A.A.; Yoo, S. 115 W fiber laser with an all solid-structure and a large-mode-area multicore fiber. Opt. Lett. 2018, 43, 3369–3372. [Google Scholar] [CrossRef]
  5. Li, L.; Schülzgen, A.; Chen, S.; Temyanko, V.L.; Moloney, J.V.; Peyghambarian, N. Phase locking and in-phase supermode selection in monolithic multicore fiber lasers. Opt. Lett. 2006, 31, 2577–2579. [Google Scholar] [CrossRef] [PubMed]
  6. Li, G.; Bai, N.; Zhao, N.; Xia, C. Space-division multiplexing: The next frontier in optical communication. Adv. Opt. Photonics 2014, 6, 413–487. [Google Scholar] [CrossRef] [Green Version]
  7. Saridis, G.M.; Alexandropoulos, D.; Zervas, G.; Simeonidou, D. Survey and evaluation of space division multiplexing: From technologies to optical networks. IEEE Commun. Surv. Tutor. 2015, 17, 2136–2156. [Google Scholar] [CrossRef] [Green Version]
  8. Trichili, A.; Park, K.H.; Zghal, M.; Ooi, B.S.; Alouini, M.S. Communicating using spatial mode multiplexing: Potentials, challenges, and perspectives. IEEE Commun. Surv. Tutor. 2019, 21, 3175–3203. [Google Scholar] [CrossRef] [Green Version]
  9. Napartovich, A.P.; Vysotsky, D.V. Theory of spatial mode competition in a fiber amplifier. Phys. Rev. A 2007, 76, 063801. [Google Scholar] [CrossRef] [Green Version]
  10. Jollivet, C.; Mafi, A.; Flamm, D.; Duparré, M.; Schuster, K.; Grimm, S.; Schülzgen, A. Mode-resolved gain analysis and lasing in multi-supermode multi-core fiber laser. Opt. Express 2014, 22, 30377–30386. [Google Scholar] [CrossRef]
  11. Eidam, T.; Wirth, C.; Jauregui, C.; Stutzki, F.; Jansen, F.; Otto, H.J.; Schmidt, O.; Schreiber, T.; Limpert, J.; Tünnermann, A. Experimental observations of the threshold-like onset of mode instabilities in high power fiber amplifiers. Opt. Express 2011, 19, 13218–13224. [Google Scholar] [CrossRef] [PubMed]
  12. Stutzki, F.; Otto, H.J.; Jansen, F.; Gaida, C.; Jauregui, C.; Limpert, J.; Tünnermann, A. High-speed modal decomposition of mode instabilities in high-power fiber lasers. Opt. Lett. 2011, 36, 4572–4574. [Google Scholar] [CrossRef] [PubMed]
  13. Wright, L.G.; Sidorenko, P.; Pourbeyram, H.; Ziegler, Z.M.; Isichenko, A.; Malomed, B.A.; Menyuk, C.R.; Christodoulides, D.N.; Wise, F.W. Mechanisms of spatiotemporal mode-locking. Nat. Phys. 2020, 16, 565–570. [Google Scholar] [CrossRef] [Green Version]
  14. Wright, L.G.; Christodoulides, D.N.; Wise, F.W. Spatiotemporal mode-locking in multimode fiber lasers. Science 2017, 358, 94–97. [Google Scholar] [CrossRef] [Green Version]
  15. Mohammadzahery, Z.; Jandaghi, M.; Aghayari, E.; Nabavi, H. Observation of spatial nonlinear self-cleaning in a few-mode step-index fiber for special distributions of initial excited modes. Sci. Rep. 2021, 11, 1–7. [Google Scholar] [CrossRef]
  16. Eckhardt, R.; Ulrich, R. Mode-beating spectroscopy in a few-mode optical guide. Appl. Phys. Lett. 1993, 63, 284–286. [Google Scholar] [CrossRef]
  17. Vallés, J.A.; Benedicto, D. Optimized active multicore fiber bending sensor. Opt. Mater. 2019, 87, 53–57. [Google Scholar] [CrossRef] [Green Version]
  18. Andrews, D.L. Structured Light and Its Applications: An Introduction to Phase-Structured Beams and Nanoscale Optical Forces; Academic Press: Cambridge, MA, USA, 2011. [Google Scholar]
  19. Forbes, A. Laser Beam Propagation: Generation and Propagation of Customized Light; CRC Press: Boca Raton, FL, USA, 2014. [Google Scholar]
  20. Flamm, D.; Schulze, C.; Naidoo, D.; Schroter, S.; Forbes, A.; Duparré, M. All-Digital Holographic Tool for Mode Excitation and Analysis in Optical Fibers. J. Light. Technol. 2013, 31, 1023–1032. [Google Scholar] [CrossRef]
  21. Shapira, O.; Abouraddy, A.F.; Joannopoulos, J.D.; Fink, Y. Complete modal decomposition for optical waveguides. Phys. Rev. Lett. 2005, 94, 143902. [Google Scholar] [CrossRef]
  22. Gerchberg, R.W. A practical algorithm for the determination of plane from image and diffraction pictures. Optik 1972, 35, 237–246. [Google Scholar]
  23. Lü, H.; Zhou, P.; Wang, X.; Jiang, Z. Fast and accurate modal decomposition of multimode fiber based on stochastic parallel gradient descent algorithm. Appl. Opt. 2013, 52, 2905–2908. [Google Scholar] [CrossRef] [PubMed]
  24. Li, L.; Leng, J.; Zhou, P.; Chen, J. Multimode fiber modal decomposition based on hybrid genetic global optimization algorithm. Opt. Express 2017, 25, 19680–19690. [Google Scholar] [CrossRef] [PubMed]
  25. Brüning, R.; Gelszinnis, P.; Schulze, C.; Flamm, D.; Duparré, M. Comparative analysis of numerical methods for the mode analysis of laser beams. Appl. Opt. 2013, 52, 7769–7777. [Google Scholar] [CrossRef]
  26. An, Y.; Huang, L.; Li, J.; Leng, J.; Yang, L.; Zhou, P. Learning to decompose the modes in few-mode fibers with deep convolutional neural network. Opt. Express 2019, 27, 10127–10137. [Google Scholar] [CrossRef] [PubMed]
  27. Rothe, S.; Zhang, Q.; Koukourakis, N.; Czarske, J. Intensity-only mode decomposition on multimode fibers using a densely connected convolutional network. J. Light. Technol. 2021, 39, 1672–1679. [Google Scholar] [CrossRef]
  28. Zhang, Q.; Rothe, S.; Koukourakis, N.; Czarske, J. Learning the matrix of few-mode fibers for high-fidelity spatial mode transmission. APL Photonics 2022, 7, 066104. [Google Scholar] [CrossRef]
  29. Manuylovich, E.S.; Dvoyrin, V.V.; Turitsyn, S.K. Fast mode decomposition in few-mode fibers. Nat. Commun. 2022, 11, 1–9. [Google Scholar] [CrossRef]
  30. Andermahr, N.; Theeg, T.; Fallnich, C. Novel approach for polarization-sensitive measurements of transverse modes in few-mode optical fibers. Appl. Phys. B 2008, 91, 353–357. [Google Scholar] [CrossRef]
  31. Nicholson, J.W.; Yablon, A.D.; Ramachandran, S.; Ghalmi, S. Spatially and spectrally resolved imaging of modal content in large-mode-area fibers. Opt. Express 2008, 16, 7233–7243. [Google Scholar] [CrossRef]
  32. Schimpf, D.N.; Barankov, R.A.; Ramachandran, S. Cross-correlated (C2) imaging of fiber and waveguide modes. Opt. Express 2011, 19, 13008–13019. [Google Scholar] [CrossRef]
  33. Ma, Y.Z.; Sych, Y.; Onishchukov, G.; Ramachandran, S.; Peschel, U.; Schmauss, B.; Leuchs, G. Fiber-modes and fiber-anisotropy characterization using low-coherence interferometry. Appl. Phys. B 2009, 96, 345–353. [Google Scholar] [CrossRef]
  34. Daniel, J.M.O.; Chan, J.S.P.; Kim, J.W.; Sahu, J.K.; Ibsen, M.; Clarkson, W.A. Novel technique for mode selection in a multimode fiber laser. Opt. Express 2011, 19, 12434–12439. [Google Scholar] [CrossRef] [PubMed]
  35. Stutzki, F.; Jauregui, C.; Voigtlaender, C.; Thomas, J.; Nolte, S.; Limpert, J.; Tuennermann, A. Real-time monitoring of the modal content of monolithic large-mode-area fiber lasers. In Proceedings of the 2010 Conference on Optical Fiber Communication (OFC/NFOEC), collocated National Fiber Optic Engineers Conference, San Diego, CA, USA, 21–25 March 2010. [Google Scholar]
  36. Kaiser, T.; Flamm, D.; Schröter, S.; Duparré, M. Complete modal decomposition for optical fibers using CGH-based correlation filters. Opt. Express 2009, 17, 9347–9356. [Google Scholar] [CrossRef] [PubMed]
  37. Flamm, D.; Naidoo, D.; Schulze, C.; Forbes, A.; Duparré, M. Mode analysis with a spatial light modulator as a correlation filter. Opt. Lett. 2012, 37, 2478–2480. [Google Scholar] [CrossRef]
  38. Gervaziev, M.D.; Zhdanov, I.; Kharenko, D.S.; Gonta, V.A.; Volosi, V.M.; Podivilov, E.V.; Babin, S.A.; Wabnitz, S. Mode decomposition of multimode optical fiber beams by phase-only spatial light modulator. Laser Phys. Lett. 2020, 18, 15101. [Google Scholar] [CrossRef]
  39. Schulze, C.; Lorentz, A.; Flamm, D.; Hartung, A.; Schröter, S.; Bartelt, H.; Duparré, M. Mode resolved bend loss in few-mode optical fibers. Opt. Express 2013, 21, 3170–3181. [Google Scholar] [CrossRef]
  40. Schulze, C.; Naidoo, D.; Flamm, D.; Schmidt, O.A.; Forbes, A.; Duparré, M. Wavefront reconstruction by modal decomposition. Opt. Express 2012, 20, 19714–19725. [Google Scholar] [CrossRef] [Green Version]
  41. Pinnell, J.; Nape, I.; Sephton, B.; Cox, M.A.; Rodríguez-Fajardo, V.; Forbes, A. Modal analysis of structured light with spatial light modulators: A practical tutorial. JOSA A 2020, 37, C146–C160. [Google Scholar] [CrossRef]
  42. Lee, W.H. Sampled Fourier Transform Hologram Generated by Computer. Appl. Opt. 1970, 9, 639–643. [Google Scholar] [CrossRef]
  43. Burckhardt, C.B. A simplification of Lee’s method of generating holograms by computer. Appl. Opt. 1970, 9, 1949. [Google Scholar] [CrossRef]
  44. Hsueh, C.K.; Sawchuk, A.A. Computer-generated double-phase holograms. Appl. Opt. 1978, 17, 3874–3883. [Google Scholar] [CrossRef] [PubMed]
  45. Arrizón, D.; Sánchez-de-la-Llave, D. Double-phase holograms implemented with phase-only spatial light modulators: Performance evaluation and improvement. Appl. Opt. 2002, 41, 3436–3447. [Google Scholar] [CrossRef] [PubMed]
  46. Mendoza-Yero, O.; Mínguez-Vega, G.; Lancis, J. Encoding complex fields by using a phase-only optical element. Opt. Let. 2014, 39, 1740–1743. [Google Scholar] [CrossRef] [Green Version]
  47. Gloge, D. Weakly guiding fibers. Appl. Opt. 1971, 10, 2252–2258. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Goodman, J.W. Introduction to Fourier Optics; McGraw-Hill: New York, USA, 1968. [Google Scholar]
  49. Holoeye, PLUTO. Available online: https://holoeye.com/spatial-light-modulators/slm-pluto-phase-only/ (accessed on 19 September 2022).
  50. Reichelt, S. Spatially resolved phase-response calibration of liquid-crystal-based spatial light modulators. Appl. Opt. 2013, 52, 2610–2618. [Google Scholar] [CrossRef] [PubMed]
  51. Schulze, C.; Ngcobo, S.; Duparré, M.; Forbes, A. Modal decomposition without a priori scale information. Opt. Express 2012, 20, 27866–27873. [Google Scholar] [CrossRef]
Figure 1. Scheme of the modal decomposition setup. L1, L2, and L3 are lenses. SLM is the spatial light modulator where the CGH is implemented. U ( ε , μ ) is the waveguide output electric field, U ( x , y ) is the electric field at the CGH plane after passing through the beam expansor in a 4f configuration, T ( x , y ) is the SLM transmittance, W 0 ( x , y ) is the electric field after the SLM, and W f ( x , y ) is the electric field at the CCD camera plane.
Figure 1. Scheme of the modal decomposition setup. L1, L2, and L3 are lenses. SLM is the spatial light modulator where the CGH is implemented. U ( ε , μ ) is the waveguide output electric field, U ( x , y ) is the electric field at the CGH plane after passing through the beam expansor in a 4f configuration, T ( x , y ) is the SLM transmittance, W 0 ( x , y ) is the electric field after the SLM, and W f ( x , y ) is the electric field at the CCD camera plane.
Micromachines 13 02004 g001
Figure 2. Scheme of the modal analysis setup. FMF: few mode fiber, MP: micro-positioners, MO: microscope objective, L1 and L2: lenses, LP: linear polarizer, BS: beam splitter, CCD: cameras, SLM: spatial light modulator; A and B are the measured transversal intensity distributions in order to align the setup and test the proposed techniques, respectively.
Figure 2. Scheme of the modal analysis setup. FMF: few mode fiber, MP: micro-positioners, MO: microscope objective, L1 and L2: lenses, LP: linear polarizer, BS: beam splitter, CCD: cameras, SLM: spatial light modulator; A and B are the measured transversal intensity distributions in order to align the setup and test the proposed techniques, respectively.
Micromachines 13 02004 g002
Figure 3. Simulated transversal intensity distributions at the CCD2 plane, when the CGH is illuminated with the LP 01   mode and the (a) LP 11 e and (b) LP 11 e modes are implemented in the SLM, for different transversal displacements of the CGH. The optical axis point is marked with a red cross.
Figure 3. Simulated transversal intensity distributions at the CCD2 plane, when the CGH is illuminated with the LP 01   mode and the (a) LP 11 e and (b) LP 11 e modes are implemented in the SLM, for different transversal displacements of the CGH. The optical axis point is marked with a red cross.
Micromachines 13 02004 g003
Figure 4. Simulated transversal intensity distributions at the L2 back focal plane when the CGH is illuminated with the LP 01   mode and the LP 02   mode is implemented in the SLM for different magnification values of the CGH.
Figure 4. Simulated transversal intensity distributions at the L2 back focal plane when the CGH is illuminated with the LP 01   mode and the LP 02   mode is implemented in the SLM for different magnification values of the CGH.
Micromachines 13 02004 g004
Figure 5. Measured intensity at the CCD2 camera optical axis point as a function of the LP 02   mode magnification implemented in the SLM when illuminated with the LP 01   mode. Three measured transversal intensity distributions are superimposed for three specified points.
Figure 5. Measured intensity at the CCD2 camera optical axis point as a function of the LP 02   mode magnification implemented in the SLM when illuminated with the LP 01   mode. Three measured transversal intensity distributions are superimposed for three specified points.
Micromachines 13 02004 g005
Figure 6. (a) Simultaneous MD spatial configuration scheme; (b,c) normalized transversal intensity distributions W f , S and W f , N , respectively, when all six LP-modes are equally present in the incident light.
Figure 6. (a) Simultaneous MD spatial configuration scheme; (b,c) normalized transversal intensity distributions W f , S and W f , N , respectively, when all six LP-modes are equally present in the incident light.
Micromachines 13 02004 g006
Figure 7. Signal and noise intensity evolution as a function of the maximum value in the transmittance function (Equation (7)) for the configuration case shown in Figure 6.
Figure 7. Signal and noise intensity evolution as a function of the maximum value in the transmittance function (Equation (7)) for the configuration case shown in Figure 6.
Micromachines 13 02004 g007
Figure 8. Relative phase measurements of the LP 21 e mode together with its cosine least-squares best fit, based on Equation (18), for the transversal intensity distribution shown in Figure 2b.
Figure 8. Relative phase measurements of the LP 21 e mode together with its cosine least-squares best fit, based on Equation (18), for the transversal intensity distribution shown in Figure 2b.
Micromachines 13 02004 g008
Figure 9. Arbitrary measured transversal intensity distribution (a) together with their numerical reconstructions (b) when the size is not readjusted through the proposed technique in Section 4, (c) when the phases are determined without using the fitting procedure presented in Section 6, and (d) when both proposed techniques are employed. The residual intensity patterns are shown below each reconstructed image.
Figure 9. Arbitrary measured transversal intensity distribution (a) together with their numerical reconstructions (b) when the size is not readjusted through the proposed technique in Section 4, (c) when the phases are determined without using the fitting procedure presented in Section 6, and (d) when both proposed techniques are employed. The residual intensity patterns are shown below each reconstructed image.
Micromachines 13 02004 g009
Table 1. Simulated modal weights (obtained with the W f function) and relative errors (obtained by comparing them with the ones determined through the W f , S function) following the simultaneous configuration of Figure 6 when all six LP-modes are equally present in the incident light.
Table 1. Simulated modal weights (obtained with the W f function) and relative errors (obtained by comparing them with the ones determined through the W f , S function) following the simultaneous configuration of Figure 6 when all six LP-modes are equally present in the incident light.
ModeWeight (%)Relative Error (%)
LP 01   23.440
LP 11 e 14.1−15
LP 11 o 7.8−53
LP 21 e 17.23
LP 21 o 17.23
LP 02   20.322
Table 2. Measured modal weights and relative phases to that of the LP01 mode from the output intensity distribution shown in Figure 9a.
Table 2. Measured modal weights and relative phases to that of the LP01 mode from the output intensity distribution shown in Figure 9a.
ModeWeight (%)Phase (π rad)
LP 01   330.0
LP 11 e 251.0
LP 11 o 101.9
LP 11 e 140.2
LP 11 o 30.9
LP 02   160.9
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Benedicto, D.; Collados, M.V.; Martín, J.C.; Atencia, J.; Mendoza-Yero, O.; Vallés, J.A. Contribution to the Improvement of the Correlation Filter Method for Modal Analysis with a Spatial Light Modulator. Micromachines 2022, 13, 2004. https://doi.org/10.3390/mi13112004

AMA Style

Benedicto D, Collados MV, Martín JC, Atencia J, Mendoza-Yero O, Vallés JA. Contribution to the Improvement of the Correlation Filter Method for Modal Analysis with a Spatial Light Modulator. Micromachines. 2022; 13(11):2004. https://doi.org/10.3390/mi13112004

Chicago/Turabian Style

Benedicto, David, María Victoria Collados, Juan C. Martín, Jesús Atencia, Omel Mendoza-Yero, and Juan A. Vallés. 2022. "Contribution to the Improvement of the Correlation Filter Method for Modal Analysis with a Spatial Light Modulator" Micromachines 13, no. 11: 2004. https://doi.org/10.3390/mi13112004

APA Style

Benedicto, D., Collados, M. V., Martín, J. C., Atencia, J., Mendoza-Yero, O., & Vallés, J. A. (2022). Contribution to the Improvement of the Correlation Filter Method for Modal Analysis with a Spatial Light Modulator. Micromachines, 13(11), 2004. https://doi.org/10.3390/mi13112004

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop