Next Article in Journal
What Is Significant in Modern Augmented Reality: A Systematic Analysis of Existing Reviews
Next Article in Special Issue
Deep Learning Network for Speckle De-Noising in Severe Conditions
Previous Article in Journal
Digital Hebrew Paleography: Script Types and Modes
Previous Article in Special Issue
Phase Retardation Analysis in a Rotated Plane-Parallel Plate for Phase-Shifting Digital Holography
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Generation of Ince–Gaussian Beams Using Azocarbazole Polymer CGH

1
Doctor’s Program of Materials Chemistry, Graduate School of Science and Technology, Kyoto Institute of Technology, Kyoto 606-8585, Japan
2
Bachelor’s Program of Materials Chemistry, Graduate School of Science and Technology, Kyoto Institute of Technology, Kyoto 606-8585, Japan
3
Materials Innovation Laboratory, Kyoto Institute of Technology, Kyoto 606-8585, Japan
4
Faculty of Materials Science and Engineering, Kyoto Institute of Technology, Kyoto 606-8585, Japan
*
Authors to whom correspondence should be addressed.
Submission received: 30 March 2022 / Revised: 16 May 2022 / Accepted: 17 May 2022 / Published: 21 May 2022
(This article belongs to the Special Issue Digital Holography: Development and Application)

Abstract

:
Ince–Gaussian beams, defined as a solution to a wave equation in elliptical coordinates, have shown great advantages in applications such as optical communication, optical trapping and optical computation. However, to ingress these applications, a compact and scalable method for generating these beams is required. Here, we present a simple method that satisfies the above requirement, and is capable of generating arbitrary Ince–Gaussian beams and their superposed states through a computer-generated hologram of size 1 mm2, fabricated on an azocarbazole polymer film. Other structural beams that can be derived from the Ince–Gaussian beam were also successfully generated by changing the elliptical parameters of the Ince–Gaussian beam. The orthogonality relations between different Ince–Gaussian modes were investigated in order to verify applicability in an optical communication regime. The complete python source code for computing the Ince–Gaussian beams and their holograms are also provided.

1. Introduction

The solution to the paraxial wave equation in elliptical coordinates can be expressed in terms of an Ince polynomial. The electromagnetic wave generated in accordance with this solution is known as an Ince–Gaussian (IG) beam. These beams form a complete orthogonal basis in Hilbert space and can be converted into other structural beams such as a Laguerre–Gaussian (LG) beam (in cylindrical coordinates) and Hermite–Gaussian (HG) beam (in cartesian coordinates), by varying their elliptical parameters [1]. Hence, IG beams are considered as one of the most generalized beams. IG beams have advantages over LG and HG beams due to their higher spatial degree of freedom and these beams can also be used as a fundamental mode which can generate other complex vector optical fields [2]. IG beams find application in various fields, such as quantum information [3], optical communication [4], optical storage [5], biological medicine [6], optical trapping [7] and non-linear optics [8]. In order to harness the benefits of IG beams, a simple and efficient method to generate these beams is important. A variety of different tools and methods for the generation of IG beams have been reported. An IG beam was first observed by Schwarz et al. [9] and Bandres et al. [10] in a stable resonator. Later, Otsuka et al. [11] reported generation by breaking the symmetry of a cavity in a solid state laser under a tightly focused pump beam. Ren et al. [12] reported the arbitrary generation of Ince–Gaussian beams using a digital micromirror device (DMD), which is an array of 1024 × 768 micromirrors with a switching frequency of 5.2 kHz. Aguirre et al. [13] used kinoform phase elements for the generation of IG and HG beams. The generation of a helical Ince–Gaussian beam was reported by Bentley et al. [14], who used a complex amplitude and phase mask encoded on a liquid crystal display. Wang et al. [15] demonstrated Ince–Gaussian-beam array generation by using a computer-generated hologram (CGH) and spatial light modulator (SLM). A second harmonic Ince–Gaussian beam generation was reported by Wang et al. [16] by employing a binary non-linear CGH. Ohtomo et al. [17] demonstrated a method to generate vortex, HG and IG beams, and also the conversion of one structural beam to another using an astigmatism mode converter (AMC).
All the existing methods, tools and techniques for generating Ince–Gaussian and other similar structural beams require a bulky system assembled on an optical bench or a complex fabrication method (e.g., metamaterials) [18]. We have recently shown that compact vortex beam generators can be realized by using azocarbazole CGH [19], which is simple and compact. In this paper, we further extend the method to incorporate the more general IG beams and their generation. The method combines a digital hologram printing (DHP) technique, CGH calculation and photoisomeric polymer fabrication, to generate arbitrary IG beams and their superposition states. We printed a CGH of an IG beam on an azocarbazole polymer film of size 1 mm 2 , which is simple and small compared to the existing popular devices such as SLM. The printed hologram can be used to reconstruct an IG beam and its superposition states without any complex optical setup. In summary, the following advantages in the generation of IG beams can be expected when comparing to the conventional methods that use SLM/DMD:
  • The two methods reported for the generation of arbitrary IG beams are LCoS-SLM and DMD, which have a pixel pitch of 3 microns. This results in a small diffraction angle and they are always required to have a 4f Fourier-filtering setup to separate out the zero-order beam from the diffracted beam. The digitally printed hologram on azopolymers has a pixel pitch of 0.8 microns and does not require any optics in the reconstruction setup. This provides a significant reduction in system size footprint compared to SLM/DMD.
  • The printed holograms are just 1 mm in size and can be used in an integrated optical device (after printing), which cannot be achieved using SLM/DMD.
  • Digital hologram printing allows to print large-sized holograms (25 cm) at a pixel pitch of 0.8 microns, thereby offering a huge space–bandwidth product. This can be beneficial in multiplexing 1000’s of beam modes, which is a crucial requirement in information processing and communication applications. To do the same using SLM/DMD requires spatially tiling SLMs which significantly increases the cost and size footprint.
It is noted that azocarbazole polymer also differs from the other popular holographic polymer called a photorefractive polymer. A photorefractive polymer is a real-time read/write material, where both the writing beams and reading beams must be present all the time, which increases system complexity. Since the purpose of the research reported in this paper is to develop a compact system without any optics in the readout process, photorefractive polymers were not considered. Moreover, the photoinduced birefringence exhibited by the azo polymers assists in obtaining a high diffraction efficiency of 30%, which is indeed helpful. Complications with respect to polarization states can be avoided, if the light beams present during recording have a reasonably pure polarization state.
In order to generate the CGH of an IG beam, it is first required to compute the complex amplitudes (of the beam cross-section) corresponding to each mode number at a particular propagation distance. Earlier reports used proprietary software such as MATLAB [20] and Mathematica [21] to realize these computations. Here, we implement the same using Python3 [22], which is open source and freely available. We also provide the Python source, in Supplementary Materials, that can reproduce all the results presented in the paper. Intensity and phase profiles corresponding to different IG beams and their derivatives (LG and HG) were successfully generated and verified experimentally. The superposition states of different IG beams that correspond to different elliptical parameters have also been generated and their orthogonality relations were investigated.

2. Theory and Methods

2.1. Theory

2.1.1. Ince–Gaussian Beam

In this section, we will discuss the theoretical foundation of an IG beam and some of its characteristics. The mathematical origin of this beam can be traced back to the Helmholtz equation (also known as the wave equation), which is defined as follows,
2 U + k 2 U = 0
where k is a wave number. Now, for the beam propagating in the z direction, the solution of Equation (1) can be described as
U ( x , y , z ) = U 0 ( x , y , z ) e i k z
Combining Equations (1) and (2) under slow-varying approximation of wave function ( U 0 ), the Helmholtz equation changes to
2 U 0 x 2 + 2 U 0 y 2 + 2 i k U 0 z = 0
The lowest order solution to Equation (3) yields a Gaussian beam defined as
U ( r ) = A W 0 W ( z ) e r 2 W 2 ( z ) e i ( k z + k r 2 2 R ( z ) ψ ( z ) )
where A is normalization constant, r = x 2 + y 2 , W(z) is the spot size of the beam W ( z ) = W 0 1 + Z Z R 2 , such as Z R is the Rayleigh range, and W 0 is beam width at z = 0, R(z) is radius of curvature R ( z ) = Z + Z R 2 Z and ψ ( z ) is the Gouy phase ( ψ ( z ) = arctan Z Z R ) . The higher order solution to Equation (1) in elliptical coordinate system under paraxial approximation is given by
I G p , m o ( r , ϵ ) = S W 0 W ( z ) S p m ( i ζ , ϵ ) S p m ( η , ϵ ) e r 2 W 2 ( z ) e i ( k z + k r 2 2 R ( z ) ( p + 1 ) ψ ( z ) )
I G p , m e ( r , ϵ ) = C W 0 W ( z ) C p m ( i ζ , ϵ ) C p m ( η , ϵ ) e r 2 W 2 ( z ) e i ( k z + k r 2 2 R ( z ) ( p + 1 ) ψ ( z ) )
where S and C are the normalization constants, indices o and e denote the odd and even Ince–Gaussian polynomial and S p , m and C p , m are odd and even Ince polynomials. These polynomials are defined with order p and degree m, such that p and m must have the same parity ( ( 1 ) p m = 1 ) and obey the condition 0 m p for odd Ince polynomials and 1 m p for even Ince polynomials. In above equations, ζ and η denote the radial and angular elliptical coordinates, respectively, which are defined as x = fcosh ζ cos η and y = fsinh ζ sin η , such that ζ [ 0 , ) , η [ 0 , 2 π ) and ϵ is elliptical parameter defined as ϵ = 2 f 2 W 2 , where f is the semifocal separation. The beam thus generated and corresponding to the Ince–Gaussian polynomial is known as the Ince–Gaussian beam. The intensity and phase profile of even and odd Ince–Gaussian beams are shown in Figure 1.
The normalization constant of IG beam can be obtained using the orthogonality relation given as
( I G p , m σ ( ζ , η ) ) × ( I G p , m σ ( ζ , η ) ) d s = δ σ σ δ p p δ m m
The linear combination of odd and even Ince–Gaussian beams form a helical Ince–Gaussian beam, which is defined as
I G p , m h ( ζ , η ) = I G p , m e ( ζ , η ) + i I G p , m o ( ζ , η )
Some of the helical Ince–Gaussian beams can form an elliptical vortex followed by elliptical rings. The total number of rings of these kinds of helical Ince–Gaussian beams, ( I G p , m h ( ζ , η ) ), are equal to 1 + p m 2 . The intensity and phase profile of these particular kinds of helical Ince–Gaussian beams are shown in Figure 2.
The IG beam can be transformed into other structural beams, such as HG and LG, by varying its elliptical parameter ( ϵ ). When ϵ = 0, the I G p , m σ beam changes to LG ( L G n , l σ ) beam, such as l = m and n = p m 2 and for ϵ = the I G p , m σ beam changes to HG ( H G n x , n y σ ) beam, such as n x = m 1 and n y = p m + 1 for σ = odd and n x = m and n y = p m for σ = even. The transformations of IG beam to LG and HG beam for different elliptical parameter values are shown in Figure 3. All these structural beams form a complete set in Hilbert space, so any one beam (e.g., LG, HG, IG) can be written in terms of other beams, as shown in Equation (9).
I G p , m σ ( ζ , η ) = n , l A n , l ( L G l , m σ ( r , ϕ ) ) = a , b B a , b ( H G n x , n y ( x , y ) )
where A n , l and B a , b are the coefficients that can be obtained by calculating the overlap integral, which is given as
A n , l = L G l , m σ ( r , ϕ ) | I G p , m σ ( ζ , η )
B a , b = H G n x , n y ( x , y ) | I G p , m σ ( ζ , η )
such that the value of coefficients A n , l and B a , b will always be unique for a given IG beam and will follow the condition given in Equation (12).
n , l A n , l 2 = a , b B a , b 2 = 1
Hence, any one beam that forms a complete set in a given coordinate system can be transformed into another complete beam defined in another coordinate system by using the transformation matrix [10], whose elements can be calculated by using the overlap integral defined in Equations (10) and (11).

2.1.2. Azocarbazole Polymer

Azocarbazole is a photoswitchable chemical compound which exhibits high efficiency in optical data storage. The cis-trans photoisomerization of azocarbazole polymer makes it a suitable polymeric material for hologram recording [23,24]. When this sample is exposed to a light pattern with alternate bright and dark regions, the molecules in brighter regions undergo isomerization and relax in a particular orientation once the light is turned off. The direction of the molecular orientation depends on the polarization state of the input beam, resulting in photoinduced birefringence. The molecules in the darker region are unaffected and stay randomly oriented. Therefore, a polarized beam experiences a difference in refractive index between the bright (exposed) and dark (unexposed) areas when passing through the sample. This difference in refractive index corresponds to the recorded hologram pattern. The recorded pattern remains intact until the polymer is heated or strongly radiated. The polarization of input (writing) beams make this azocarbazole polymer anisotropic and, hence, optically birefringent. More details on the performance of the material for hologram recording and reconstruction can be found in Kinashi et al. [24].
Figure 3. Transformation of Ince–Gaussian beam (center) into Hermite–Gaussian beam (left) and Laguerre–Gaussian beam (right) at different elliptical parameter ( ϵ ) values. (A) Conversion of I G 3 , 1 e for ϵ = 2 into H G 1 , 2 for ϵ = and L G 1 , 1 for ϵ = 0 . (B) Conversion of I G 4 , 2 o for ϵ = 2 into H G 1 , 3 for ϵ = and L G 2 , 1 for ϵ = 0 . (C) Conversion of I G 6 , 6 h for ϵ = 2 into H G 6 , 0 for ϵ = and L G 0 , 6 for ϵ = 0 .
Figure 3. Transformation of Ince–Gaussian beam (center) into Hermite–Gaussian beam (left) and Laguerre–Gaussian beam (right) at different elliptical parameter ( ϵ ) values. (A) Conversion of I G 3 , 1 e for ϵ = 2 into H G 1 , 2 for ϵ = and L G 1 , 1 for ϵ = 0 . (B) Conversion of I G 4 , 2 o for ϵ = 2 into H G 1 , 3 for ϵ = and L G 2 , 1 for ϵ = 0 . (C) Conversion of I G 6 , 6 h for ϵ = 2 into H G 6 , 0 for ϵ = and L G 0 , 6 for ϵ = 0 .
Jimaging 08 00144 g003

2.2. Methods

2.2.1. Calculation of Computer-Generated Hologram(CGH)

The computer-generated hologram of IG beam is calculated using the scripting language Python. In order to compute CGH, we need to generate the phase profile of IG beam of interest. For calculating phase profile, we first define a meshgrid of size n × n cm in cartesian coordinates which is divided into N × N number of pixels. We note that a rectangular aperture is usually used to simplify the plots, though Ince–Gaussian beams are elliptically symmetrical and correspond to an elliptical aperture. Using the values of rectangular aperture defined in cartesian coordinates, an elliptical coordinate ( ζ , η ) is defined for a fixed value of elliptical parameter ( ϵ ) [25]. After defining the elliptical coordinate, we divide the simulation into two parts: (i) odd Ince beam and (ii) even Ince beam. Each part is further divided into two sections: (i) odd indices (e.g., I G 3 , 1 σ ) and (ii) even indices (e.g., I G 4 , 2 σ ). In each section we generate the Ince polynomial by computing the coefficients in the summation and normalization constant [10]. After generating the Ince polynomial, we set various parameters such as beam waist, propagation distance, wavelength of beam, etc. Using these parameters, odd and even IG beams are generated, as defined by Equations (5) and (6). Helical Ince–Gaussian beam can easily be computed by adding even and odd IG beams. In order to compute the hologram of different IG beams, we can rewrite Equations (5), (6) and (8) in simplified form, as
I G ( ζ , η ) = A ( x , y ) e i ϕ ( x , y )
where A ( x , y ) is the amplitude and ϕ ( x , y ) is the phase of the IG beam in cartesian coordinates. A linear phase grating ( e i 2 π x d ) is multiplied with IG beam which serves as a carrier, where d is the period of grating. The hologram is computed by adding the phases of the IG beam and carrier beam, assuming its amplitude to be uniform across the beam cross section. Therefore, what we are generating is known as a ‘phase hologram’ and should not be confused with ‘phase-only’ hologram that is generated by iterative algorithms such as Gerchberg—Saxton (GS). The holograms of different IG beams are simulated for a meshgrid dimension of 1 cm × 1 cm with pixels values of 1200 × 1200. The elliptical parameter ( ϵ ) is 2.0 and propagation distance is z = 0 m. The wavelength of the beam is 532 nm and beam waist is 3 mm; hence, the calculated value of f(z) at z = 0 m is given as f(0) = 0.003 m. The hologram is generated for a carrier frequency of 10 pixels per cycle and is shown in Figure 4. These computer-generated holograms will be printed on azocarbazole polymers using digital hologram printing.
Commented Python3 source files for calculating phase, intensity and hologram of the different Ince–Gaussian beams and their superposition are provided as a Supplementary Materials. They include two files: (i) I G b e a m . p y —the core module, which includes the necessary function to generate IG beams (this file needs no modification), and (ii) e x a m p l e . p y —the file where user sets the parameters and then executes it to obtain the results. The provided ‘example.py’ is pre-set to generate amplitude and phase patterns shown in Figure 1d and the corresponding hologram pattern shown in Figure 4c, when executed without any modification. To change to a different pattern or to generate the superposition states, the user may follow the comments inlined in the source code.

2.2.2. Sample Synthesis

The azocarbazole polymer film is prepared using a four-step process. In the first step, poly(CACzE-MMA), which is composed of 3-((4-cynophenyl)azo)-9H-carbazole-9-ethanol (CACzE) and methyl metacrylate(MMA), is mixed with CACzE and Diphenyl Phthalate(DPP) in ratio of 9:3:1 by weight. In the second step, we dissolve the above prepared sample into Tetrahydrofuran(THF) and stir for 48 h. The third step involves drying the sample at 70 C for another 48 h. In the last step, the sample is melted at 180–185 C and pressed between glass slides, where the thickness is controlled by placing polyimide spacer between the slides. The method of synthesis of azocarbazole polymer film is summarized in Figure 5. The synthesized azocarbazole polymer film is shown in Figure 6a. The thickness of the sample ranges from 30 μ m to 35 μ m .The sample shown in Figure 6a has thickness of 32 μ m and the total thickness of sample is 1.87 mm (including thickness of glass slides).

2.2.3. Digital Hologram Printing (DHP)

Digital hologram printing (DHP) is a technique used to print digitally calculated holograms (such as CGH) on a material medium using an optical setup (e.g., fringe printers) [26]. A simple schematic of the hologram printing setup which is used to print the CGH computed in Section 2.2.1 on the azocarbazole polymer film synthesized in Section 2.2.2 is shown in Figure 7. In this setup, the CGH is displayed on the spatial light modulator (SLM), where the physical size of hologram is 1 cm × 1 cm. The SLM being employed is the LC-R 1080 reflective LCOS SLM from HOLOEYE Photonics AG (HOLOEYE Photonics AG, Berlin, Germany), which has a resolution of 1920 × 1200 pixels, a pixel pitch of 8.1 μ m and a frame rate of 60 Hz. A collimated beam coming out from the laser is split by a polarizing beam splitter (PBS) and the reflected beam illuminates the SLM. The phase hologram (phase modulation) pattern displayed on the SLM (phase only) is converted to amplitude modulation when the reflected beam from SLM passes through the PBS. A 10× demagnifying lens arrangement transfers the intensity pattern passing through the PBS onto the sample as refractive index modulation through a process called photoinduced birefringence. The sample is mounted on a high-precision 3-axis motion stage for accurate position of the sample w.r.t the illuminating beam (for more details about the printing setup, one may refer to Jackin et al. [19]). The laser power at the sample plane is measured to be 32 mW/mm 2 and the recording is performed with a 5 s exposure. Multiple holograms are printed on different areas on a sample of size 3 cm × 3 cm and a photograph of the printed sample is shown in Figure 6b. The printed area of one hologram on the sample is 1 mm 2 . The diameter of the whole sample is measured to be 2 cm and we can print upto 36 holograms (avoiding the overlap between two holograms) on one sample. The sample shown in Figure 6b has a diameter of 1.87 cm and 9 holograms were printed on it. Figure 6c shows the phase contrast optical microscopic image of one of the printed holograms on azocarbazole polymer film. The refractive index changes are clearly visible as bright and dark areas in the sample.

3. Results

The printed hologram can be reconstructed using the simple setup shown in Figure 8. A collimated beam from a laser of wavelength 640 nm was passed through the printed area of the sample and the intensity pattern of the exiting radiation was observed using a CCD. The reconstructed intensity pattern for the different odd IG beams ( I G p , m o ) and even IG beams ( I G p , m e ) are shown in Figure 9a–j. The experimentally observed intensity patterns were compared with the simulated intensity patterns and were found to be in good agreement.
The ring-type helical Ince–Gaussian beams ( I G p , m h ) were also generated using the same setup. The experimental verification of generated intensity profiles for different helical Ince–Gaussian beams are shown in Figure 10.
The diffraction efficiency ( η ( % ) ) for the generated IG beams was calculated using Equation (14), given as:
η ( % ) = I + 1 I t × 100
where I + 1 is the intensity of +1 order of the diffracted beam and I t is the total transmitted intensity, which can be measured by passing a collimated beam through the unprinted sample. The calculated value of diffraction efficiency is 30%. This is in close agreement with the theoretical maximum efficiency of the Raman–Nath regime, which is 33.9%. Hence, we can conclude that the efficiency obtained from our sample is very high. The retention time of the printed hologram on the sample is nearly 50 days, which makes the azocarbazole polymer suitable for writing once and later reading multiple times (permanent storage). However, with the application of heat on purpose, the written pattern can be erased, and the sample be reused for another recording (unlike photopolymers, which are non-reusable).
IG beams were transformed into other structural beams by changing their elliptical parameter ( ϵ ) values. Initially, an even IG beam ( I G 2 , 2 e ) and an odd IG beam ( I G 3 , 3 o ) were generated for ϵ = 2 and then, by changing the value of ϵ to 10 6 , the I G 2 , 2 e beam changes to H G 2 , 0 e and the I G 3 , 3 o beam changes to H G 2 , 1 e . For ϵ = 0, the I G 2 , 2 e beam changes to L G 0 , 2 e and I G 3 , 3 o changes to L G 0 , 3 e . The corresponding generated beams are shown in Figure 11.
The superposition states of IG beams were also generated, which is important in the domain of optical communication. We generated the superposition of an IG beam for ϵ = 2 and ϵ = 0 and they are shown in Figure 12. Figure 12a shows the superposition of I G 3 , 1 e and I G 3 , 3 e , while Figure 12b shows the superposition of I G 4 , 4 o and I G 4 , 2 o , for ϵ = 2 . The superposition of I G 9 , 9 o and I G 9 , 9 e and the superposition of I G 9 , 3 o and I G 9 , 3 e for ϵ = 0 is shown in Figure 12c,d, respectively. Figure 12c,d corresponds to the generation of a vortex and a vortex with circular rings, respectively. Information can be encoded into these superposition states and can be transmitted to the receiver who can use the orthogonality relation between these beams to decode the information. On the other hand, plane electromagnetic beams do not possess orthogonality and hence are difficult to separate once they are superposed (unless angular, polarization and wavelength multiplexing, etc., are employed). Moreover, certain modes of IG beams also possess orbital angular momentum (vortex modes), which can be used to encode information directly in terms of topological charges.
The orthonormality constant was calculated by evaluating Equation (7) numerically in Python. The process involved in the calculation of an orthonormality constant consists of four steps. In the first step, the dot product between two simulated IG beams was calculated; then, in the second step, the total number of pixels occupied by this dot product in the region defined in the meshgrid was calculated. In the third step, the area occupied by 1 pixel was calculated by using formula ( L N 1 ) 2 , such that L is the length of the meshgrid and N is the total number of pixels. Finally, in the last step, the orthonormality constant was calculated by multiplying the total number of pixels occupied by the dot product with the area of one pixel.
The investigated orthonormality relation for different configurations is shown in Figure 13. Figure 13a,b shows the orthogonality relation between even and odd indices. The variation in orthogonality relation with different parameters such as beam waist ( W o ), meshgrid dimension (n × n), pixels numbers (N × N) and elliptical parameter ( ϵ ) was analysed and the results thus obtained are shown in Figure 13c–f. The logarithmic variation in orthogonality constant with meshgrid sizes for modes I G 5 , 3 e and I G 3 , 1 e for ϵ = 2 can be seen from Figure 13c for a given beam waist of 3 mm and pixel number of 512 × 512. We found that for a meshgrid size greater than 1.3 cm × 1.3 cm, the orthogonality relation is mostly satisfied. Figure 13d shows the variation in orthogonality constant between modes I G 5 , 3 e and I G 3 , 1 e for ϵ = 2 with the pixel number for a meshgrid size of 1 cm × 1 cm and a beam waist of 3 mm. We find that the orthogonality relation in this case deteriorates when we increase the pixel number. Figure 13e manifests the variation in the normalization constant (factor of 1) between modes I G 4 , 2 o and I G 4 , 2 o for ϵ = 2 with a varying beam waist for a meshgrid size of 1 cm × 1 cm and pixel number of 1000 × 1000. Here, we find that, for beam waist size less than 4 mm, the normalization relation is well satisfied and after that, the normalization relation starts violating the ideal case. Figure 13f conveys the variation in normalization constant with elliptical parameter between modes I G 4 , 4 o and I G 4 , 4 o for a meshgrid size of 1 cm × 1 cm, pixels number 1000 × 1000 and beam waist 3 mm. We observe that the normalization constant is well-defined at all values of elliptical parameters, which indicates the validation of the orthonormality relation even when the IG beam changes to an HG or LG beam. Hence, from all these observations, we can conclude that a perfect combination of meshgrid dimension, pixel numbers and beam waist is important for generating IG beams and their superposed state to access their application in optical communication.

4. Conclusions

A combination of digital hologram printing and an azocarbazole polymer was used to overcome the high cost, large footprint and non-scalablity of existing IG beam generators. Intensity, phase and the corresponding CGH of different IG beams were computed using the open source Python programming language. The simulated and experimentally generated intensity profiles of IG beams were found to be in good agreement, which determines that the beams were generated successfully. The conversion of an IG beam into other structural beams, such as LG and HG, is also demonstrated by simulation and was verified experimentally. In order to access the applicability of IG beams in optical communication, the superposed states of different IG and LG beams were generated and their orthogonality relations were examined. The reported compact and scalable IG-beam-generation technique could open up new possibilities in structural-beam science. Further improvement to this technique can be made by increasing the printing area of a hologram on polymer film, which will possibly increase the diffraction efficiency and assist larger mode generation. Other structural beams defined in different coordinates, such as Bessel–Gaussian beams [27], Mathieu–Gauss beams [28], Airy–Gauss beams [29], hypergeometric Gaussian beams [30], etc., can also be generated using the proposed technique.

Supplementary Materials

Author Contributions

Conceptualization, B.J.J. and S.K.S.; Methodology, S.K.S. and H.H.; Validation, S.K.S., B.J.J. and K.K.; Investigation, B.J.J. and K.K.; Software, Visualization and formal analysis, S.K.S.; writing—original draft preparation, S.K.S.; writing—review and editing, B.J.J., K.K., N.T. and W.S.; Supervision, B.J.J., K.K., N.T. and W.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

S.K.S. acknowledges the Ministry of Education, Sports and Culture of Japan (MEXT) scholarship provided by the Japanese Government for his studies in Kyoto Institute of Technology.

Conflicts of Interest

The authors declare no conflict of interest and certify that they have no affiliations with or involvement in any organization or entity with any financial interest or non-financial interest in the subject matter or materials discussed in this manuscript.

Abbreviations

The following abbreviations are used in this manuscript:
IGInce–Gaussian
HGHermite Gaussian
LGLaguerre Gaussian
CGHComputer generated hologram
DHPDigital Hologram printing
CCDCharge coupled device
SLMSpatial light modulator
MMAMethyl metacrylate
DPPDiphenyl pthalate
THFTetrahydrofuran
PBSPolarizing beam splitter

References

  1. Bandres, M.A.; Gutiérrez-Vega, J.C. Ince–Gaussian beams. Opt. Lett. 2004, 29, 144–146. [Google Scholar] [CrossRef] [PubMed]
  2. Jiaqing, L.; Jiming, W.; Chongjun, H.; Tong, W.; Yuangang, L.; Haixia, M.; Youwen, L. Research on the Generation of Ince-Gaussian Vector Optical Field. Acta Opt. Sin. 2019, 39, 0826001. [Google Scholar]
  3. Baghdasaryan, B.; Fritzsche, S. Enhanced entanglement from Ince-Gaussian pump beams in spontaneous parametric down-conversion. Phys. Rev. A 2020, 102, 052412. [Google Scholar] [CrossRef]
  4. Sakpal, S.; Milione, G.; Li, M.J.; Nouri, M.; Shahoei, H.; LaFave, T.; Ashrafi, S.; MacFarlane, D. Stability of Ince–Gaussian beams in elliptical core few-mode fibers. Opt. Lett. 2018, 43, 2656–2659. [Google Scholar] [CrossRef] [PubMed]
  5. Yu, Y.; Chen, Y.; Wang, C.; Wang, J.; Sun, Z.; Cao, M.; Gao, H.; Li, F. Optical storage of Ince–Gaussian modes in warm atomic vapor. Opt. Lett. 2021, 46, 1021–1024. [Google Scholar] [CrossRef] [PubMed]
  6. Gather, M.; Yun, S. Single-cell biological lasers. Nat. Photonics 2011, 5, 406. [Google Scholar] [CrossRef]
  7. Yang, Y.; Ren, Y.; Chen, M.; Arita, Y.; Rosales-Guzmán, C. Optical trapping with structured light: A review. Adv. Photonics 2021, 3, 034001. [Google Scholar] [CrossRef]
  8. Pires, D.; Rocha, J.; Jesus-Silva, A.; Fonseca, E. Optical mode conversion through nonlinear two-wave mixing. Phys. Rev. A 2019, 100, 043819. [Google Scholar] [CrossRef]
  9. Schwarz, U.T.; Bandres, M.A.; Gutiérrez-Vega, J.C. Observation of Ince–Gaussian modes in stable resonators. Opt. Lett. 2004, 29, 1870–1872. [Google Scholar] [CrossRef] [Green Version]
  10. Bandres, M.A.; Gutiérrez-Vega, J.C. Ince–Gaussian modes of the paraxial wave equation and stable resonators. JOSA A 2004, 21, 873–880. [Google Scholar] [CrossRef] [Green Version]
  11. Otsuka, K.; Nemoto, K.; Kamikariya, K.; Miyasaka, Y.; Chu, S.C. Linearly polarized single-frequency oscillations of laser-diode-pumped microchip ceramic Nd: YAG lasers with forced Ince–Gaussian mode operations. Jpn. J. Appl. Phys. 2007, 46, 5865. [Google Scholar] [CrossRef]
  12. Ren, Y.X.; Fang, Z.X.; Gong, L.; Huang, K.; Chen, Y.; Lu, R.D. Dynamic generation of Ince-Gaussian modes with a digital micromirror device. J. Appl. Phys. 2015, 117, 133106. [Google Scholar] [CrossRef]
  13. Aguirre-Olivas, D.; Mellado-Villaseñor, G.; Sánchez-de-la Llave, D.; Arrizón, V. Efficient generation of Hermite–Gauss and Ince–Gauss beams through kinoform phase elements. Appl. Opt. 2015, 54, 8444–8452. [Google Scholar] [CrossRef] [PubMed]
  14. Bentley, J.B.; Davis, J.A.; Bandres, M.A.; Gutiérrez-Vega, J.C. Generation of helical Ince-Gaussian beams with a liquid-crystal display. Opt. Lett. 2006, 31, 649–651. [Google Scholar] [CrossRef] [PubMed]
  15. Wang, D.; Li, Y.; Hu, X.; Jin, L. Tailoring Ince-Gaussian arrays based on a computer-generated hologram. Appl. Opt. 2021, 60, 7172–7176. [Google Scholar] [CrossRef]
  16. Wang, M.Y.; Tang, J.; Wang, H.J.; Ming, Y.; Zhang, Y.; Cui, G.X.; Lu, Y.Q. Generation of second-harmonic Ince-Gaussian beams. Appl. Phys. Lett. 2018, 113, 081105. [Google Scholar] [CrossRef]
  17. Ohtomo, T.; Chu, S.C.; Otsuka, K. Generation of vortex beams from lasers with controlled Hermite-and Ince-Gaussian modes. Opt. Express 2008, 16, 5082–5094. [Google Scholar] [CrossRef]
  18. Jin, J.; Pu, M.; Wang, Y.; Li, X.; Ma, X.; Luo, J.; Zhao, Z.; Gao, P.; Luo, X. Multi-Channel Vortex Beam Generation by Simultaneous Amplitude and Phase Modulation with Two-Dimensional Metamaterial. Adv. Mater. Technol. 2017, 2, 1600201. [Google Scholar] [CrossRef]
  19. Jackin, B.J.; Shirai, M.; Haginaka, H.; Kinashi, K.; Tsutsumi, N.; Sakai, W. Compact and Scalable Large Vortex Array Generation Using Azocarbazole Polymer and Digital Hologram Printing Technique. Nanoscale Res. Lett. 2022, 17, 44. [Google Scholar] [CrossRef]
  20. Bandres, M.A. Ince Gaussian Beam. 2022. Available online: https://www.mathworks.com/matlabcentral/fileexchange/46222-ince-gaussian-beam (accessed on 20 February 2022).
  21. Zeleny, E. Ince Polynomial. 2014. Available online: http://demonstrations.wolfram.com/IncePolynomials/ (accessed on 20 February 2022).
  22. Van Rossum, G.; Drake, F.L. Python 3 Reference Manual; CreateSpace: Scotts Valley, CA, USA, 2009. [Google Scholar]
  23. Kinashi, K.; Fukami, T.; Yabuhara, Y.; Motoishi, S.; Sakai, W.; Kawamoto, M.; Sassa, T.; Tsutsumi, N. Molecular design of azo-carbazole monolithic dyes for updatable full-color holograms. NPG Asia Mater. 2016, 8, e311. [Google Scholar] [CrossRef] [Green Version]
  24. Kinashi, K.; Nakanishi, I.; Sakai, W.; Tsutsumi, N. Material Design of Azo-Carbazole Copolymers for Preservation Stability with Rewritable Holographic Stereograms. Macromol. Chem. Phys. 2019, 220, 1800456. [Google Scholar] [CrossRef]
  25. Sun, C. Explicit equations to transform from Cartesian to Elliptic coordinates. Math. Model. Appl. 2017, 2, 43–46. [Google Scholar] [CrossRef] [Green Version]
  26. Yoshikawa, H.; Tachinami, M. Development of direct fringe printer for computer-generated holograms. In Proceedings of the Practical Holography XIX: Materials and Applications, San Jose, CA, USA, 25–27 January 2005; International Society for Optics and Photonics: Bellingham, WA, USA, 2005; Volume 5742, pp. 259–266. [Google Scholar]
  27. Chu, X.; Sun, Q.; Wang, J.; Lü, P.; Xie, W.; Xu, X. Generating a Bessel-Gaussian beam for the application in optical engineering. Sci. Rep. 2015, 5, 18665. [Google Scholar] [CrossRef] [PubMed]
  28. Ez-Zariy, L.; Nebdi, H.; Chafiq, A.; Belafhal, A. Mathieu-Gauss beam and its propagation. Phys. Chem. News 2012, 64, 42–52. [Google Scholar]
  29. Bandres, M.A.; Gutiérrez-Vega, J.C. Airy-Gauss beams and their transformation by paraxial optical systems. Opt. Express 2007, 15, 16719–16728. [Google Scholar] [CrossRef] [Green Version]
  30. Karimi, E.; Zito, G.; Piccirillo, B.; Marrucci, L.; Santamato, E. Hypergeometric-gaussian modes. Opt. Lett. 2007, 32, 3053–3055. [Google Scholar] [CrossRef] [Green Version]
Figure 1. (A) Intensity (first row) and phase (second row) profile of even Ince−Gaussian beams (a) I G 1 , 1 e , (b) I G 3 , 1 e , (c) I G 5 , 5 e , (d) I G 8 , 6 e , (e) I G 12 , 2 e for ϵ = 2 . (B) Intensity (first row) and phase (second row) profile of odd Ince−Gaussian beams (f) I G 2 , 2 o , (g) I G 3 , 1 o , (h) I G 5 , 5 o , (i) I G 10 , 2 o , (j) I G 12 , 10 o for ϵ = 2 .
Figure 1. (A) Intensity (first row) and phase (second row) profile of even Ince−Gaussian beams (a) I G 1 , 1 e , (b) I G 3 , 1 e , (c) I G 5 , 5 e , (d) I G 8 , 6 e , (e) I G 12 , 2 e for ϵ = 2 . (B) Intensity (first row) and phase (second row) profile of odd Ince−Gaussian beams (f) I G 2 , 2 o , (g) I G 3 , 1 o , (h) I G 5 , 5 o , (i) I G 10 , 2 o , (j) I G 12 , 10 o for ϵ = 2 .
Jimaging 08 00144 g001
Figure 2. Intensity (first row) and phase profile (second row) profile of helical Ince−Gaussian beams (a) I G 8 , 8 h , (b) I G 10 , 8 h , (c) I G 12 , 8 h , (d) I G 14 , 8 h for ϵ = 2 .
Figure 2. Intensity (first row) and phase profile (second row) profile of helical Ince−Gaussian beams (a) I G 8 , 8 h , (b) I G 10 , 8 h , (c) I G 12 , 8 h , (d) I G 14 , 8 h for ϵ = 2 .
Jimaging 08 00144 g002
Figure 4. Computer-generated hologram for different Ince–Gaussian beams. (a) I G 5 , 3 o , (b) I G 3 , 1 e , (c) I G 8 , 6 e , (d) I G 6 , 4 o , (e) I G 12 , 8 h for ϵ = 2 .
Figure 4. Computer-generated hologram for different Ince–Gaussian beams. (a) I G 5 , 3 o , (b) I G 3 , 1 e , (c) I G 8 , 6 e , (d) I G 6 , 4 o , (e) I G 12 , 8 h for ϵ = 2 .
Jimaging 08 00144 g004
Figure 5. Process involved in the synthesis of azocarbazole polymer film.
Figure 5. Process involved in the synthesis of azocarbazole polymer film.
Jimaging 08 00144 g005
Figure 6. (a) Photograph of the synthesized azocarbazole polymer film sample. (b) Photograph of the sample with 9 holograms printed on film. (c) Phase contrast image of the printed hologram (scalebar is 50 μ m in length).
Figure 6. (a) Photograph of the synthesized azocarbazole polymer film sample. (b) Photograph of the sample with 9 holograms printed on film. (c) Phase contrast image of the printed hologram (scalebar is 50 μ m in length).
Jimaging 08 00144 g006
Figure 7. Schematic of optical setup for digital hologram printing on azocarbazole polymer film.
Figure 7. Schematic of optical setup for digital hologram printing on azocarbazole polymer film.
Jimaging 08 00144 g007
Figure 8. Optical setup for reconstruction of printed hologram on azocarbazole polymer film.
Figure 8. Optical setup for reconstruction of printed hologram on azocarbazole polymer film.
Jimaging 08 00144 g008
Figure 9. Intensity patterns of simulated and experimentally generated Ince–Gaussian beams (a) I G 6 , 6 e , (b) I G 6 , 4 o , (c) I G 5 , 3 o , (d) I G 4 , 2 o , (e) I G 3 , 3 o , (f) I G 2 , 2 o , (g) I G 5 , 5 e , (h) I G 1 , 1 e , (i) I G 2 , 2 e , (j) I G 5 , 5 o for ϵ = 2 .
Figure 9. Intensity patterns of simulated and experimentally generated Ince–Gaussian beams (a) I G 6 , 6 e , (b) I G 6 , 4 o , (c) I G 5 , 3 o , (d) I G 4 , 2 o , (e) I G 3 , 3 o , (f) I G 2 , 2 o , (g) I G 5 , 5 e , (h) I G 1 , 1 e , (i) I G 2 , 2 e , (j) I G 5 , 5 o for ϵ = 2 .
Jimaging 08 00144 g009
Figure 10. Intensity pattern of simulated and experimentally generated helical Ince–Gaussian beams (a) I G 8 , 8 h , (b) I G 10 , 8 h , (c) I G 12 , 8 h , (d) I G 14 , 8 h for ϵ = 2 .
Figure 10. Intensity pattern of simulated and experimentally generated helical Ince–Gaussian beams (a) I G 8 , 8 h , (b) I G 10 , 8 h , (c) I G 12 , 8 h , (d) I G 14 , 8 h for ϵ = 2 .
Jimaging 08 00144 g010
Figure 11. (A) Conversion of I G 2 , 2 e beam (second column) into H G 2 , 0 e beam (first column) and L G 0 , 2 e beam (third column) for different ϵ values. (B) Conversion of I G 3 , 3 o beam (second column) into H G 2 , 1 o beam (first column) and L G 0 , 3 o beam (third column) for different ϵ values.
Figure 11. (A) Conversion of I G 2 , 2 e beam (second column) into H G 2 , 0 e beam (first column) and L G 0 , 2 e beam (third column) for different ϵ values. (B) Conversion of I G 3 , 3 o beam (second column) into H G 2 , 1 o beam (first column) and L G 0 , 3 o beam (third column) for different ϵ values.
Jimaging 08 00144 g011
Figure 12. Holograms (first row) of superposition of different IG beams. Simulation (second row) and experimental (Third Row) intensity profile of superposition of IG beams (a) I G 3 , 3 e + I G 3 , 1 e for ϵ = 2 , (b) I G 4 , 4 o + I G 4 , 2 o for ϵ = 2 , (c) I G 9 , 9 o + I G 9 , 9 e for ϵ = 0 , (d) I G 9 , 3 o + I G 9 , 3 e for ϵ = 0 .
Figure 12. Holograms (first row) of superposition of different IG beams. Simulation (second row) and experimental (Third Row) intensity profile of superposition of IG beams (a) I G 3 , 3 e + I G 3 , 1 e for ϵ = 2 , (b) I G 4 , 4 o + I G 4 , 2 o for ϵ = 2 , (c) I G 9 , 9 o + I G 9 , 9 e for ϵ = 0 , (d) I G 9 , 3 o + I G 9 , 3 e for ϵ = 0 .
Jimaging 08 00144 g012
Figure 13. (a) Orthogonality relation between Ince–Gaussian beam of even indices for ϵ = 2. (b) Orthogonality relation between Ince–Gaussian beam of odd indices for ϵ = 2. (c) Variation ins log(orthogonality constant) for beams I G 5 , 3 e and I G 3 , 1 e with meshgrid size (L) for ϵ = 2. (d) Variation in orthogonality constant for beams I G 5 , 3 e and I G 3 , 1 e with pixels number (N) for ϵ = 2. (e) Variation in normalization constant for beams I G 4 , 2 o and I G 4 , 2 o with beam waist ( W 0 ) for ϵ = 2. (f) Variation in normalization constant for beams I G 4 , 4 o and I G 4 , 4 o with elliptical parameter ( ϵ ).
Figure 13. (a) Orthogonality relation between Ince–Gaussian beam of even indices for ϵ = 2. (b) Orthogonality relation between Ince–Gaussian beam of odd indices for ϵ = 2. (c) Variation ins log(orthogonality constant) for beams I G 5 , 3 e and I G 3 , 1 e with meshgrid size (L) for ϵ = 2. (d) Variation in orthogonality constant for beams I G 5 , 3 e and I G 3 , 1 e with pixels number (N) for ϵ = 2. (e) Variation in normalization constant for beams I G 4 , 2 o and I G 4 , 2 o with beam waist ( W 0 ) for ϵ = 2. (f) Variation in normalization constant for beams I G 4 , 4 o and I G 4 , 4 o with elliptical parameter ( ϵ ).
Jimaging 08 00144 g013
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Singh, S.K.; Haginaka, H.; Jackin, B.J.; Kinashi, K.; Tsutsumi, N.; Sakai, W. Generation of Ince–Gaussian Beams Using Azocarbazole Polymer CGH. J. Imaging 2022, 8, 144. https://0-doi-org.brum.beds.ac.uk/10.3390/jimaging8050144

AMA Style

Singh SK, Haginaka H, Jackin BJ, Kinashi K, Tsutsumi N, Sakai W. Generation of Ince–Gaussian Beams Using Azocarbazole Polymer CGH. Journal of Imaging. 2022; 8(5):144. https://0-doi-org.brum.beds.ac.uk/10.3390/jimaging8050144

Chicago/Turabian Style

Singh, Sumit Kumar, Honoka Haginaka, Boaz Jessie Jackin, Kenji Kinashi, Naoto Tsutsumi, and Wataru Sakai. 2022. "Generation of Ince–Gaussian Beams Using Azocarbazole Polymer CGH" Journal of Imaging 8, no. 5: 144. https://0-doi-org.brum.beds.ac.uk/10.3390/jimaging8050144

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