• Photonics Research
  • Vol. 7, Issue 8, 875 (2019)
Yingfei Pang1、2、3, Axiu Cao1, Jiazhou Wang1, Hui Pang1, Wei Yan1, Xiangdong Wu2, Lifang Shi1、4、*, and Qiling Deng1、5、*
Author Affiliations
  • 1Institute of Optics and Electronics, Chinese Academy of Sciences, Chengdu 610209, China
  • 2School of Mechanical Engineering, Southwest Jiaotong University, Chengdu 610031, China
  • 3University of Chinese Academy of Sciences, Beijing 100049, China
  • 4e-mail: shilifang@ioe.ac.cn
  • 5e-mail: dengqiling@ioe.ac.cn
  • show less
    DOI: 10.1364/PRJ.7.000875 Cite this Article Set citation alerts
    Yingfei Pang, Axiu Cao, Jiazhou Wang, Hui Pang, Wei Yan, Xiangdong Wu, Lifang Shi, Qiling Deng, "Design and experimental verification of a monolithic complete-light modulator based on birefringent materials," Photonics Res. 7, 875 (2019) Copy Citation Text show less

    Abstract

    This paper presents a method to design a monolithic complete-light modulator (MCLM) that fully controls the amplitude, phase, and polarization of incident light. The MCLM is made of birefringent materials that provide different refractive indices to orthogonal eigen-polarizations, the ordinary o and extraordinary e states. We propose an optimization method to calculate the two relief depth distributions for the two eigen-polarizations. Also, a merging algorithm is proposed to combine the two relief depth distributions into one. The corresponding simulations were carried out in this work and the desired light distribution, including information on amplitude, phase, and four polarization states, was obtained when a laser beam passed through a 16-depth-level micro-structure whose feature size is 8 μm. The structure was fabricated by common photolithography. An experimental optical system was also set up to test the optical effects and performances of the MCLM. The experimental performance of the MCLM agrees with the simulation results, which verifies the validity of the algorithms we propose in this paper.

    1. INTRODUCTION

    A light wave contains multi-dimensional parameters: amplitude, phase, and polarization. Optical modulation technologies can be used to modulate these parameters to obtain desired optical effects, which have numerous practical applications, such as particle acceleration [15], laser fabrication [6,7], beam shaping [814], and holographic projection [1524]. At present, most light field control techniques realize only partial control of the parameters of the incident light field. However, many experts and scholars have found that there are potential applications for a complete light control technique that can simultaneously control the amplitude, phase, and polarization of incident beams. For example, this technique can be used to structure the profiles of vector fields for special purposes, such as enabling a light beam to exert force and torque on illuminated objects [25,26]. The most attractive aspect of this technique is that it can be used to enhance the longitudinal component field [27] to obtain an optical needle field [28] or focus field [29], making it useful for particle capture and particle manipulation [3032].

    Davis et al. [33] designed polarization-multiplexed diffractive optical elements (DOEs) by using two cascaded liquid-crystal displays (LCDs). The first LCD encoded the two multiplexed phase-only DOEs. The second LCD acted as a pixelated polarization rotator to change the polarization state for each of these two DOEs. This approach was able to fully control all the parameters of the spatial distributions of an optical field. However, its optical system is difficult to align precisely. Plewicki et al. [34,35] proposed a method that can independently manipulate over the major axis orientation and the axis ratio of the polarization ellipse by integrating a 4f shaper setup in both arms of a Mach–Zehnder interferometer and rotating the polarization by 90° in one of the arms before overlaying the beams. Moreno et al. [36] used a single transmissive spatial light modulator (SLM) to generate complex light fields. They divided an SLM into two regions, using the dual modulation of two transmissive subsystems to completely control the incident beam. Han et al. [37] proposed a cascaded modulation method that divided two reflective SLMs into four regions to form four reflective subsystems. This method generated various optical fields containing amplitude, phase, and polarization modulations. Chen et al. [38] and Liu et al. [39] designed different methods to generate two orthogonal complex fields. These two complex amplitude distributions generated desired light fields by coaxial superposition. The former used a phase-only SLM based on a 4f system, while the latter divided a phase-type SLM into two regions, which respectively control the ordinary (o) ray and extraordinary (e) ray of the incident light. All of the above methods were based on SLMs to achieve full control of all parameters of incident light. It has to be acknowledged that these methods have great flexibility in light field modulation, but have the common disadvantages of complex systems and relatively high cost. Wang et al. [40] obtained a high-purity optical needle field by fabricating slot antennas on a metasurface structure. Each slot antenna can be regarded as a linear polarizer to fully control the light field. The system is greatly simplified, but in order to obtain the metasurface structure, the manufacturing cost is extremely high.

    Our group previously designed a polarization multiplexing phase slice [41], which realized polarization multiplexing imaging on a single yttrium vanadate substrate (YVO4). In this manuscript, we report further developments of the modulation technique of a monolithic complete-light modulator (MCLM) and the design method for an MCLM that can yield a complex light field containing amplitude, phase, and polarization modulation. All parameters of the incident beam can be modulated without any other optical element except one MCLM. The system is simple and stable with little external interference and relatively low fabrication cost. This manuscript is organized with six sections. After this introduction, Section 2 demonstrates the principle of how to design the MCLM. Section 3 describes the design process about calculating the relief distribution of the MCLM. Section 4 shows the simulation result. In Section 5, we show the experimental results, which largely coincide with the simulation results. Section 6 is the summary.

    2. PRINCIPLE

    The design process of an MCLM includes the following three major steps. First, the desired light field E(x,y) containing amplitude, phase, and polarization information can be decomposed into two orthogonal components of o-ray and e-ray light fields Eo(x,y) and Ee(x,y) according to the Jones calculus. Second, a modified Gerchberg–Saxton (GS) algorithm [42] with both amplitude and phase limits is profiled to obtain two pure phase distributions ϕo(x,y) and ϕe(x,y), which can be used to generate the phases of Eo(x,y) and Ee(x,y). Third, by adopting point-by-point calculation, we can obtain a single relief depth distribution h(x,y) of birefringent material, which can obtain two phase delay distributions ϕo(x,y) and ϕe(x,y) of o- and e-ray light propagating in the birefringent material. The aim is to find an appropriate depth distribution h(x,y) of birefringent material that minimizes the error between ϕi(x,y) and ϕi(x,y) (i=o,e). After searching for a suitable h(x,y), the relief depth distribution of the MCLM can be obtained by quantizing h(x,y) to 16 depth levels.

    When the incident light containing equal o-ray and e-ray light (i.e., linearly polarized light with a polarization orientation of 45° relative to the o-ray and e-ray axes of a birefringent material) passes through the MCLM, the o-ray and e-ray can simultaneously be modulated to generate Eo(x,y) and Ee(x,y). The vector superposition of the two components generates the desired light field distribution E(x,y).

    3. DESIGN PROCESS

    In the design process, the following parts should be considered carefully to obtain the depth distribution h(x,y), which includes how to decompose the desired optical field into components of two orthogonally polarized complex light fields of proper amplitudes, how to retrieve two pure phase distributions of the two light fields by the modified GS algorithm, and how to obtain the relief depth distribution of birefringent material that can satisfy two phase distributions. The above three parts will be described in detail below.

    A. Decomposition of the Desired Light Field

    According to the Jones calculus, the desired light field containing complex information can be expressed as E(x,y)=A(x,y)exp[iφ(x,y)]{cos[θ(x,y)]sin[θ(x,y)]exp[iδ(x,y)].Equation (1) can be divided into two equations: {Eo(x,y)=A(x,y)cos[θ(x,y)]exp[iφ(x,y)]Ee(x,y)=A(x,y)sin[θ(x,y)]exp[iφ(x,y)]exp[iδ(x,y)].E(x,y) in Eq. (1) stands for the desired light field, which can be regarded as the vector superposition of two orthogonal components of complex light fields Eo(x,y) and Ee(x,y) in Eq. (2). φ(x,y) and δ(x,y) are the common phase and phase difference of Eo(x,y) and Ee(x,y), respectively. δ(x,y) determines the polarization state distribution of desired field, and θ(x,y) represents the polarization orientation.

    An example of a desired light field of four regions containing amplitude, phase, and polarization is shown in Fig. 1.

    Example of a light field containing four regions as divided by the purple dotted lines. Distributions of the light field parameters are: (a) amplitude A(x,y), (b) polarization θ(x,y), (c) common phase φ(x,y), (d) phase differience δ(x,y), (e) amplitude Ao(x,y) of Eo(x,y), (f) amplitude Ae(x,y) of Ee(x,y), (g) phase φo(x,y) of Eo(x,y), and (h) phase φe(x,y) of Ee(x,y).

    Figure 1.Example of a light field containing four regions as divided by the purple dotted lines. Distributions of the light field parameters are: (a) amplitude A(x,y), (b) polarization θ(x,y), (c) common phase φ(x,y), (d) phase differience δ(x,y), (e) amplitude Ao(x,y) of Eo(x,y), (f) amplitude Ae(x,y) of Ee(x,y), (g) phase φo(x,y) of Eo(x,y), and (h) phase φe(x,y) of Ee(x,y).

    Figure 1(a) shows A(x,y), which contains patterns of different amplitude distributions of A, B, C, and D in different regions. Figure 1(b) shows that the polarization state distributions in four regions are linearly polarized, and the polarization orientation θ(x,y) is indicated by four blue two-way arrows oriented at 0°, 45°, 90°, and 135°, respectively. Figure 1(c) represents the common phase φ(x,y), and Fig. 1(d) shows the phase difference δ(x,y).

    In this example, the polarization state is linearly polarized, with θ(x,y)(0,π/2) in Fig. 1(b), the phase difference δ(x,y)=2kπ(k=0,1,2,) in Fig. 1(d); if θ(x,y)(π/2,π) in Fig. 1(b), the phase difference δ(x,y)=(2k+1)π(k=0,1,2,) in Fig. 1(d); if θ(x,y)=0,π/2 or π in Fig. 1(b), the phase difference δ(x,y)=kπ(k=0,1,2,) in Fig. 1(d). Substituting A(x,y), φ(x,y), δ(x,y), and θ(x,y) into Eq. (2), the two components Eo(x,y) and Ee(x,y) can be obtained as {Eo(x,y)=Ao(x,y)exp[iφo(x,y)]Ee(x,y)=Ae(x,y)exp[iφe(x,y)].By comparing Eqs. (2) and (3), we can obtain the amplitude and phase distributions of Eo(x,y) and Ee(x,y), as shown in Figs. 1(e)1(h).

    Figures 1(e) and 1(f) represent the amplitude distributions Ao(x,y) and Ae(x,y) of Eo(x,y) and Ee(x,y), respectively. Figure 1(e) shows that the amplitude of regions A, B, C, and D is equal to 1, 2/2, 0, and 2/2, respectively. Figure 1(f) shows that the amplitude of regions A, B, C, and D is equal to 0, 2/2, 1, and 2/2, respectively.

    Figures 1(g) and 1(h) represent the phase distributions φo(x,y) and φe(x,y) of Eo(x,y) and Ee(x,y), respectively. As can be seen from Eq. (2), the phase φo(x,y) of Eo(x,y) is equal to the common phase φ(x,y). Meanwhile, the phase φe(x,y) of Ee(x,y) is equal to φ(x,y)+δ(x,y).

    B. Calculation of Phase Distributions

    The common GS algorithm and most of the modified GS algorithm [43,44] can be used to calculate the desired light field only with amplitude distribution, while the phase distribution cannot be considered. However, in this paper, since Eo(x,y) and Ee(x,y) contain determinate amplitude distributions and phase distributions, it is necessary to optimize the calculation by using the modified GS algorithm for both amplitude and phase limitations. In the iterative process, the output plane is divided into the signal domain and the noise domain, as shown in Fig. 2. To avoid the influence of the center zero-order on the experimental results, the signal domain is placed at a position offset from the center of the output plane.

    Flow chart of the modified GS algorithm with both amplitude and phase limits.

    Figure 2.Flow chart of the modified GS algorithm with both amplitude and phase limits.

    Taking Ee(x,y) as an example, the specific steps can be described as follows. The initial phase distribution of design phase plane is random value from 0 to 2π, denoted as ϕe(0)(x,y).

    The process of modified GS algorithm is shown in Fig. 2. First, the amplitude of the design phase plane is replaced by 1 to obtain complex amplitude U=exp[iϕe(k1)(x,y)], k=1,2,3,. Then U propagates from the design phase plane to the output plane, which generates a failed complex amplitude distribution E1(x,y). Ee(x,y) is the modified Fresnel-domain constraint function to correct the failed results, which can be expressed as Ee(x,y)={rAe(x,y)exp[iφe(x,y)],(x,y)SE1(x,y),(x,y)N,where Ae(x,y) and φe(x,y) are the amplitude distribution and phase distribution of Ee(x,y), respectively. r is the weight factor. S and N refer to the signal domain and the noise domain, respectively. Second, when the modified function propagates back to the design phase plane, the complex amplitude distribution U(x,y) will be generated. To obtain a phase-only distribution, we retain only the phase distribution ϕe(k1)(x,y), while the amplitude |U(x,y)| is eliminated. Finally, ϕe(k1)(x,y) becomes the phase distribution at the beginning of the next iteration of the design phase plane.

    During the iteration, we need to determine whether the failed result E1(x,y) meets the error requirement, which can be represented by the root mean square error, as shown in Eq. (5). If so, the iteration is terminated and the pure phase distribution ϕe(k1)(x,y) [i.e., ϕe(x,y)] becomes the output of the iteration; otherwise, the iteration continues. Meanwhile, ϕo(x,y) can also be obtained by using the modified GS algorithm. RMSE=[|E1(x,y)||Ed(x,y)|]2[|Ed(x,y)|]2,(x,y)S,where RMSE stands for the root mean square error of amplitude between the failed result E1 and the desired result Ed.

    Adopting the proposed algorithm, ϕo(x,y) and ϕe(x,y) can be calculated. Then, ϕo(x,y) and ϕe(x,y) will be used to search the relief depth distribution of the MCLM.

    C. Search of Relief Depth Distribution

    According to the refractive indices no and ne of the o ray and e ray propagating in birefringent material, the relief depth distributions ho(x,y) and he(x,y) can be calculated by ϕ(x,y)=2π(n1)λh(x,y),where n is the refractive index of a light wave whose wavelength is λ, and h(x,y) is the relief depth of the DOE.

    However, ho(x,y) and he(x,y) have entirely different distributions. To simultaneously modulate o and e light by a single MCLM to obtain the two orthogonal components of the desired light field, it is necessary to optimize the two different relief depth distributions to a unique relief depth distribution h(x,y).

    The calculation process is as follows. Suppose the array size of the MCLM is M×M, and the relief depth of the pixel (x,y) on the MCLM is h(x,y), where x,y[1,M]. Utilizing Eq. (6), the phase delays of h(x,y) for the o and e light can be calculated as ϕo(x,y) and ϕe(x,y). Therefore, the phase errors between ϕi(x,y) and ϕi(x,y) (i=o,e) can be calculated as Δϕo and Δϕe. However, as we all know, due to the periodicity of the phase, it is difficult to calculate the actual phase error. Therefore, Δf(k) (k=1,2) is used here to characterize Δϕi: {Δf(1)(x,y)=modλ[h(x,y)×(no1)]ho(x,y)×(no1)Δf(2)(x,y)=modλ[h(x,y)×(ne1)]he(x,y)×(ne1),where moda(b) represents the remainder of b divided by a. We regard the square sum of Δf(1)(x,y) and Δf(2)(x,y) as the equivalent judgment equation of the total phase error: ΔF(x,y)=k=12[Δf(k)(x,y)]2.Obviously, as the h(x,y) changes, the errors also change. When ΔF(x,y) is the smallest, a suitable h(x,y) can be found, which is the relief depth of the MCLM. In this paper, equal-interval search is adopted to calculate h(x,y). Taking the e-ray as a reference, h(x,y)=p(x,y)×λne1.In Eq. (9), p(x,y){0,1l,2l,,Nl} (N is an integer), where l is also an integer used to determine search accuracy. Each search interval of h(x,y) is 1l×λne1. Note that the larger l is, the higher the search accuracy is, but the longer the calculation time will be. In general, we take l=100. Substituting p=0,1l,2l,,Nl into Eq. (9), we can obtain N+1 values of h(x,y), which are 0,1l×λne1××Nl×λne1. Substituting the N+1 values into Eq. (7), and then utilizing Eq. (8), we can find a suitable h(x,y) that minimizes ΔF(x,y).

    Each pixel is optimized in the same manner, and the final relief depth distribution h can be obtained. Theoretically, if the value range of h(x,y) is not limited (N), it is certain that an appropriate h(x,y) can be found to make ΔF(x,y)0. However, it is necessary to consider the actual processing situation. In reality, if the depth is too shallow, it is difficult to find a suitable h(x,y), while if the depth is too deep, it is not conducive to machining. Meanwhile, if the depth falls outside a certain depth range, the quantized result will become worse. Therefore, when determining the value range of h(x,y) or the value of N, it is necessary to comprehensively consider the above situations.

    In this paper, to simplify the calculation, the search range of depth is chosen as an integer multiple of λne1, which provides a phase modulation of integer multiple of 2π for the e ray [see Eq. (6)]. Since the maximum value of p(x,y) in Eq. (9) is Nl, for satisfying the above condition, it requires N to be an integral multiple of l. First, let N=l=100. The optimal depth distribution is calculated by the equal-interval search algorithm, which is quantized to a 16-level structure. Two orthogonal complex fields will be calculated by diffraction calculation for the o and e rays, which will be substituted into Eq. (5) for calculating the two root mean square errors between the two calculated complex fields and the desired orthogonal complex fields, respectively. Second, in the same manner, multiple optimized depth distributions respectively corresponding to N=200,300,,1000 will be searched, and multigroup root mean square errors can be calculated. Finally, according to the criterion of minimizing the sum of root mean square errors, we determine that the 16-level relief depth distribution corresponding to N=400 is the most appropriate.

    4. SIMULATION

    To verify the correctness of the above algorithm, we took a desired light field, with its parameters as shown in Figs. 3(a)3(d), to carry out the simulation, where Fig. 3(a) is the intensity distribution I(x,y) of the desired light field that was used to characterize the amplitude distribution A(x,y). It is composed of four characters from the movie “Journey to the West.” Figure 3(b) represents the polarization state. The polarization state of each character is linear polarization, with θ(x,y)=0°, 45°, 90°, and 135°, as indicated by four blue arrows. Figure 3(c) shows that the common phase φ(x,y) of the first two characters was equal to 0.3π and the rest are equal to 0. Figure 3(d) shows that the phase difference δ(x,y) of the last two characters is equal to π and the rest are equal to 0.

    Parameters of the light field in the four regions: (a) intensity I(x,y), (b) linear polarization θ=0°, 45°, 90° and 135°, (c) common phase φ(x,y), and (d) phase differience δ(x,y).

    Figure 3.Parameters of the light field in the four regions: (a) intensity I(x,y), (b) linear polarization θ=0°, 45°, 90° and 135°, (c) common phase φ(x,y), and (d) phase differience δ(x,y).

    The wavelength of incident light was 632.8 nm. The birefringent material we chose was YVO4, which had a large birefringence with refractive indices no and ne of 1.9929 and 2.2154, respectively. The array size used in the simulation is 1024 by 1024, with a pixel size of 8  μm×8  μm. The imaging distance was 0.7 m. We have compiled the modified GS algorithm and relief depth search algorithm, and simulated the light field modulation effect of the MCLM by using the software MATLAB. The relief depth distribution of MCLM which has 16 levels is shown in Fig. 4. The deepest relief depth of the MCLM is 1.95 μm. The depth difference between adjacent levels is 0.13 μm.

    Relief depth distribution of the MCLM that has 16 levels. The depth difference between adjacent levels is 0.13 μm.

    Figure 4.Relief depth distribution of the MCLM that has 16 levels. The depth difference between adjacent levels is 0.13 μm.

    The designed MCLM is examined through simulation using a polarization analyzer for incident light with equal amounts of the two orthogonal eigen-polarizations, and the resultant light field distributions are shown in Figs. 5(a)5(e).

    Simulation results: (a) without polarizer; with polarizer oriented at (b) 0°, (c) 45°, (d) 90°, and (e) 135°.

    Figure 5.Simulation results: (a) without polarizer; with polarizer oriented at (b) 0°, (c) 45°, (d) 90°, and (e) 135°.

    Figures 5(a)5(e) represent the simulation result with 16 levels of quantization. Without a polarizer, all four characters appear as shown in Fig. 5(a). Moreover, the simulated intensity distribution is consistent with the original desired intensity distribution, which proves that the method can achieve the control of the amplitude of the incident light. Then, we simulated the case where a linear polarizer analyzer was added to detect the polarization state of the reconstructed light field. Figures 5(b)5(e) simulated that when the linear polarizer analyzer was oriented at 0°, 45°, 90°, and 135°, the four characters appeared or disappeared in turn. For example, when the linear polarizer analyzer is oriented at 0°, the third character completely disappears [Fig. 5(b)]. This confirmed that the area where the third character is located has linearly polarized light with its polarization orientation at 90°. Therefore, it could be seen from Figs. 5(b)5(e) that the light in the areas of each of the four characters is all linearly polarized with their polarization orientations at 0°, 45°, 90°, and 135°.

    As can be seen from the simulation results, the modulation of amplitude and polarization by a single MCLM can be realized. Meanwhile, according to the description of Section 3.A, the determined polarization state indicates that the phase difference between the o ray and the e ray is determined. Therefore, the modulation of phase can also be achieved by the MCLM. Furthermore, this method can also be utilized to completely control a more complex light field.

    5. PREPARATION AND EXPERIMENTAL RESULTS OF MCLM

    To verify the above design method and the simulation results, the designed MCLM was fabricated and the optical effects were tested by experiments. Here, a photolithography technique was used to carry out the fabrication. As the designed MCLM was quanticized to 16 levels of depth, multiple exposure and etching processes were carried out. The fabrication process of the 16-level structure is illustrated in Fig. 6.

    Fabrication process of the MCLM: (a) the first etching with etching depth of 130 nm, (b) the second etching with etching depth of 260 nm, (c) the third etching with etching depth of 520 nm, and (d) the fourth etching with etching depth of 1024 nm.

    Figure 6.Fabrication process of the MCLM: (a) the first etching with etching depth of 130 nm, (b) the second etching with etching depth of 260 nm, (c) the third etching with etching depth of 520 nm, and (d) the fourth etching with etching depth of 1024 nm.

    Four masks used in the process are shown in Fig. 7. The masks were used to carry out the exposure one by one to form the 16-level structure. Also, etching processes were carried out four times to transfer the structure into the YVO4 birefringent crystal material. The etching depths in the four etching processes were 130, 260, 520, and 1040 nm. We ranked the four masks according to the etching depth from shallow to deep. The area within the blue dashed line represents the graphic area of the masks, as shown in Fig. 7(a). As can be seen from the structure observed by microscope of the four masks, the graphic area was divided into two parts, which were black and white. During the exposure process, the black and white parts determined whether light could pass through the mask.

    Four masks used in the MCLM fabrication process: (a) mask one, (b) mask two, (c) mask three, and (d) mask four.

    Figure 7.Four masks used in the MCLM fabrication process: (a) mask one, (b) mask two, (c) mask three, and (d) mask four.

    The photoresist material AZ9260 was used to carry out the exposure. The photoresist was spin-coated on the YVO4 crystal substrate with rotation speed of 5000 rad/min. Process parameters of coating time, prebake temperature, prebake time, and resist thickness were 30 s, 100°C, 10 min, and 3 μm, respectively. An exposure system was employed to perform the exposure process, with a illumination light source that is a Hg lamp with a central wavelength of 365 nm. The total exposure time was 20 s. Process parameters of development, after-bake temperature, and after-bake time were 30 s, 120 °C, and 60 min, respectively. Ion beam etching was carried out to transfer the structure into the substrate. The etching times in four etching processes were 15, 30, 60, and 120 min. The resulting average depths of the four etchings were 141, 240, 539, and 1050 nm, respectively. The practical and theoretical depth errors were +11, 20, +19, and +24  nm, respectively. The pictures of the obtained elements are shown in Fig. 8. Figure 8(a) shows the photograph of the fabricated MCLM and Fig. 8(b) shows the photo taken by microscope, from which we can see the structure detail.

    (a) Fabricated MCLM. (b) Structure of the MCLM observed under microscope.

    Figure 8.(a) Fabricated MCLM. (b) Structure of the MCLM observed under microscope.

    The experimental setup and light beam path were built to test the MCLM, as shown in Fig. 9. The laser beam first propagated through a polarizer oriented at 45° in order to obtain incident linearly polarized light containing equal-amplitude o-ray and e-ray light. Then, the obtained linear polarized light passed through MCLM to produce the desired light field, of which polarization state was detected by a linear polarizer analyzer. The distance from the CCD to the location of the MCLM was 0.7 m, which was the same as in the simulation process.

    Experimental setup and light beam path for testing the MCLM performance.

    Figure 9.Experimental setup and light beam path for testing the MCLM performance.

    The desired intensity pattern was received by a CCD (Fig. 10). Figure 10(a) shows the reconstructed light field, of which the intensity distribution was consistent with the desired light field. When there was no linear polarizer analyzer, four characters appeared in the light field. By using the polarizer analyzer to detect the polarization state of the desired light field, we have obtained the patterns shown in Figs. 10(b)10(e). When the linear polarizer analyzer was oriented at 0°, 45°, 90°, and 135°, the four characters disappeared one after another. As shown in Fig. 10(b), the third character completely disappeared when the polarizer analyzer was oriented at 0°. This indicated that the area where the third person was located was 90° linearly polarized. Detecting by the linear polarizer analyzer, it was easy to see that the polarized orientations of the four characters are 0°, 45°, 90°, and 135°.

    Experimental results. The desired intensity pattern is received by a CCD. (a) No polarizer analyzer. In (b)–(e), the polarizer analyzer was oriented at 0°, 45°, 90°, and 135°, respectively.

    Figure 10.Experimental results. The desired intensity pattern is received by a CCD. (a) No polarizer analyzer. In (b)–(e), the polarizer analyzer was oriented at 0°, 45°, 90°, and 135°, respectively.

    Meanwhile, it could be seen from the experimental results that the light field generated by the MCLM was basically consistent with the theoretical simulation results shown in Figs. 5(a)5(e), with only a few deviations caused by residual manufacturing errors in the MCLM. It could be considered that the experimental results were in great agreement with the simulation results. Therefore, by utilizing only one MCLM, we have successfully controlled the phase, amplitude, and polarization of incident light.

    6. CONCLUSION

    In this paper, we have proposed a new type of optical modulation element called MCLM. The simulation and experimental results have shown that the MCLM successfully realizes complete control of incident light amplitude, phase, and polarization to obtain a tunable complex light field. The method proposed in this paper can achieve any optical response controlled by only one single MCLM component. The system is simple and stable with little external interference. Moreover, it is expected to be used in the areas of optical fabrication and light trapping.

    References

    [1] J. R. Fontana, R. H. Pantell. A high-energy, laser accelerator for electrons using the inverse Cherenkov effect. J. Appl. Phys., 54, 4285-4288(1983).

    [2] R. D. Romea, W. D. Kimura. Modeling of inverse Čerenkov laser acceleration with axicon laser-beam focusing. Phys. Rev. D, 42, 1807-1818(1990).

    [3] A. Achterberg, Y. A. Gallant, J. G. Kirk, A. W. Guthmann. Particle acceleration by ultra relativistic shocks: theory and simulations. Mon. Not. R. Astron. Soc., 328, 393-408(2001).

    [4] C. Varin, M. Piché. Acceleration of ultra-relativistic electrons using high-intensity TM01 laser beams. Appl. Phys. B, 74, s83-s88(2002).

    [5] F. A. Aharonian, A. G. Akhperjanian, K. M. Aye, A. R. Bazer-Bachi, M. Beilicke, W. Benbow, D. Berge, P. Berghaus, K. Bernlöhr, O. Bolz, C. Boisson, C. Borgmeier, F. Breitling, A. M. Brown, J. Bussons Gordo, P. M. Chadwick, V. R. Chitnis, L. M. Chounet, R. Cornils, L. Costamante, B. Degrange, A. Djannati-Ataï, L. O’C. Drury, T. Ergin, P. Espigat, F. Feinstein, P. Fleury, G. Fontaine, S. Funk, Y. A. Gallant, B. Giebels, S. Gillessen, P. Goret, J. Guy, C. Hadjichristidis, M. Hauser, G. Heinzelmann, G. Henri, G. Hermann, J. A. Hinton, W. Hofmann, M. Holleran, D. Horns, O. C. de Jager, I. Jung, B. Khélifi, A. Konopelko, I. J. Latham, R. Le Gallou, M. Lemoine, A. Lemière, N. Leroy, T. Lohse, A. Marcowith, C. Masterson, T. J. L. McComb, M. de Naurois, S. J. Nolan, A. Noutsos, K. J. Orford, J. L. Osborne, M. Ouchri, M. Panter, G. Pelletier, S. Pita, M. Pohl, G. Pühlhofer, M. Punch, B. C. Raubenheimer, M. Raue, J. Raux, S. M. Rayner, I. Redondo, A. Reimer, O. Reimer, J. Ripken, M. Rivoal, L. Rob, L. Rolland, G. Rowell, V. Sahakian, L. Saugé, S. Schlenker, R. Schlickeiser, C. Schuster, U. Schwanke, M. Siewert, H. Sol, R. Steenkamp, C. Stegmann, J.-P. Tavernet, C. G. Théoret, M. Tluczykont, D. J. van der Walt, G. Vasileiadis, P. Vincent, B. Visser, H. J. Völk, S. J. Wagner. High-energy particle acceleration in the shell of a supernova remnant. Nature, 432, 75-77(2004).

    [6] C. Hnatovsky, V. Shvedov, W. Krolikowski, A. Rode. Revealing local field structure of focused ultrashort pulses. Phys. Rev. Lett., 106, 123901(2011).

    [7] K. Lou, S. X. Qian, Z. C. Ren, C. Tu, Y. Li, H. T. Wang. Femtosecond laser processing by using patterned vector optical fields. Sci. Rep., 3, 2281(2013).

    [8] A. Chakraborty, B. Das, G. Sanyal. Beam shaping using nonlinear phase distribution in a uniformly spaced array. IEEE. T. Antennas Propag., 30, 1031-1034(1982).

    [9] J. S. Liu, M. R. Taghizadeh. Iterative algorithm for the design of diffractive phase elements for laser beam shaping. Opt. Lett., 27, 1463-1465(2002).

    [10] Z. Sun, H. K. Kim. Refractive transmission of light and beam shaping with metallic nano-optic lenses. Appl. Phys. Lett., 85, 642-644(2004).

    [11] D. L. Shealy, J. A. Hoffnagle. Beam shaping profiles and propagation. Proc. SPIE, 5876, 58760D(2005).

    [12] J. S. Liu, A. J. Caley, M. R. Taghizadeh. Diffractive optical elements for beam shaping of monochromatic spatially incoherent light. Appl. Opt., 45, 8440-8447(2006).

    [13] H. Wang, L. Shi, B. Lukyanchuk, C. Sheppard, C. T. Chong. Creation of a needle of longitudinally polarized light in vacuum using binary optics. Nat. Photonics, 2, 501-505(2008).

    [14] H. Ma, P. Zhou, X. Wang, Y. Ma, F. Xi, X. Xu, Z. Liu. Near-diffraction-limited annular flattop beam shaping with dual phase only liquid crystal spatial light modulators. Opt. Express, 18, 8251-8260(2010).

    [15] J. Wang, L. Liu, A. Cao, H. Pang, C. Xu, Q. Mu, J. Chen, L. Shi, Q. Deng. Generation of color images by utilizing a single composite diffractive optical element. Micromachines, 9, 508(2018).

    [16] T. Ito, K. Okano. Color electroholography by three colored reference lights simultaneously incident upon one hologram panel. Opt. Express, 12, 4320-4325(2004).

    [17] T. Shimobaba, T. Takahashi, N. Masuda, T. Ito. Numerical study of color holographic projection using space-division method. Opt. Express, 19, 10287-10292(2011).

    [18] M. Makowski, M. Sypek, A. Kolodziejczyk. Colorful reconstructions from a thin multi-plane phase hologram. Opt. Express, 16, 11618-11623(2008).

    [19] M. Makowski, M. Sypek, I. Ducin, A. Fajst, A. Siemion, J. Suszek, A. Kolodziejczyk. Experimental evaluation of a full-color compact lensless holographic display. Opt. Express, 17, 20840-20846(2009).

    [20] M. Makowski, M. Sypek, A. Kolodziejczyk, G. Mikuła. Three-plane phase-only computer hologram generated with iterative Fresnel algorithm. Opt. Eng., 44, 125805(2005).

    [21] R. G. Dorsch, A. W. Lohmann, S. Sinzinger. Fresnel ping-pong algorithm for two-plane computer-generated hologram display. Appl. Opt., 33, 869-875(1994).

    [22] M. Sypek. Light propagation in the Fresnel region. New numerical approach. Opt. Commun., 116, 43-48(1995).

    [23] M. Makowski. Iterative design of multiplane holograms: experiments and applications. Opt. Eng., 46, 045802(2007).

    [24] M. Sypek, C. Prokopowicz, M. Górecki. Image multiplying and high-frequency oscillations effects in the Fresnel region light propagation simulation. Opt. Eng., 42, 3158-3164(2003).

    [25] D. B. Ruffner, D. G. Grier. Optical forces and torques in nonuniform beams of light. Phys. Rev. Lett., 108, 173602(2012).

    [26] M. I. Marqués. Beam configuration proposal to verify that scattering forces come from the orbital part of the Poynting vector. Opt. Lett., 39, 5122-5125(2014).

    [27] C. J. R. Sheppard, A. Choudhury. Annular pupils, radial polarization and superresolution. Appl. Opt., 43, 4322-4327(2004).

    [28] S. Zhou, S. Wang, J. Chen, G. Rui, Q. Zhan. Creation of radially polarized optical fields with multiple controllable parameters using a vectorial optical field generator. Photon. Res., 4, B35-B39(2016).

    [29] J. J. Hao, Z. L. Yu, Z. Z. Chen, H. Chen, J. P. Ding. Shaping of focal field with controllable amplitude, phase, and polarization. Chin. Opt. Lett., 12, 090501(2014).

    [30] G. Rui, J. Chen, X. Wang, B. Gu, Y. Cui, Q. Zhan. Synthesis of focused beam with controllable arbitrary homogeneous polarization using engineered vectorial optical fields. Opt. Express, 24, 23667-23676(2016).

    [31] O. M. Maragò, P. H. Jones, P. G. Gucciardi, G. Volpe, A. C. Ferrari. Optical trapping and manipulation of nanostructures. Nat. Nanotechnol., 8, 807-819(2013).

    [32] Q. Zhan. Trapping metallic Rayleigh particles with radial polarization. Opt. Express, 12, 3377-3382(2004).

    [33] J. A. Davis, G. H. Evans, I. Moreno. Polarization-multiplexed diffractive optical elements with liquid-crystal displays. Appl. Opt., 44, 4049-4052(2005).

    [34] M. Plewicki, S. M. Weber, F. Weise, A. Lindinger. Independent control over the amplitude, phase, and polarization of femtosecond pulses. Appl. Phys. B, 86, 259-263(2007).

    [35] M. Plewicki, F. Weise, S. M. Weber, A. Lindinger. Phase, amplitude, and polarization shaping with a pulse shaper in a Mach-Zehnder interferometer. Appl. Opt., 45, 8354-8359(2006).

    [36] I. Moreno, J. A. Davis, T. M. Hernandez, D. M. Cottrell, D. Sand. Complete polarization control of light from a liquid crystal spatial light modulator. Opt. Express, 20, 364-376(2012).

    [37] W. Han, Y. Yang, W. Cheng, Q. Zhan. Vectorial optical field generator for the creation of arbitrarily complex fields. Opt. Express, 21, 20692-20706(2013).

    [38] Z. Chen, T. Zeng, B. Qian, J. Ding. Complete shaping of optical vector beams. Opt. Express, 23, 17701-17710(2015).

    [39] S. Liu, S. Qi, Y. Zhang, P. Li, D. Wu, L. Han, J. Zhao. Highly efficient generation of arbitrary vector beams with tunable polarization, phase, and amplitude. Photon. Res., 6, 228-233(2018).

    [40] S. Wang, D. C. Abeysinghe, Q. Zhan. Generation of vectorial optical fields with slot-antenna-based metasurface. Opt. Lett., 40, 4711-4714(2015).

    [41] J. Wang, A. Cao, H. Pang, M. Zhang, G. Wang, J. Chen, L. Shi, Q. Deng, S. Song. Vector optical field generation based on birefringent phase plate. Opt. Express., 25, 12531-12539(2017).

    [42] S. Tao, W. Wu. Beam shaping of complex amplitude with separate constraints on the output beam. Opt. Express, 23, 1052-1062(2015).

    [43] R. W. Gerchberg, W. O. Saxton. A practical algorithm for the determination of phase from image and diffraction plane pictures. Optik, 35, 1-6(1972).

    [44] J. S. Liu, A. J. Caley, M. R. Taghizadeh. Symmetrical iterative Fourier-transform algorithm using both phase and amplitude freedoms. Opt. Commun., 267, 347-355(2006).

    Yingfei Pang, Axiu Cao, Jiazhou Wang, Hui Pang, Wei Yan, Xiangdong Wu, Lifang Shi, Qiling Deng, "Design and experimental verification of a monolithic complete-light modulator based on birefringent materials," Photonics Res. 7, 875 (2019)
    Download Citation