Abstract
Greek ladder is a technique for approximating by rational numbers where n is a positive integer and C is a positive real number. For the classical Greek ladder, the value is. Based on the continued fraction theory and algebraic equation, the classical Greek ladder in a special case can be reduced to the generalized Fibonacci sequence. By means of proper switching and binary, ternary or quaternary phase modulation, here we have successfully designed the various kinds of nano-photonic devices to produce three-dimensional array foci whose focusing properties satisfy the above mathematical characteristics. With this technology, the diffraction-limited array foci are freely designed or distributed under the requirement at the desired multiple focal planes.
© 2015 Optical Society of America
1. Introduction
Three-dimensional array foci, as the name implies, mean that two-dimensional array diffraction-limited foci are distributed as required at multi-focal planes. Focusing and imaging of light has many applications in ophthalmology, biological cell, virtual display technique, high-resolution microscopy, spectroscopy, and lithography. Unfortunately, a conventional refractive lens cannot focus soft x-rays and EUV as a consequence of the strong absorption of solid materials in the x-ray and EUV spectral regions. Although a traditional Fresnel zone plate (FZP) can be used for this kind of focusing, it has inherent limitations [1–5]. Some aperiodic zone plates, generated by fractal Cantor set, have been proposed to overcome these limitations [6–8]. In 2001, a new concept – photon sieve (PS) was proposed [9], which is a FZP with the transparent zones replaced by a great number of completely separated pinholes to overcome the disadvantages of traditional zone plate. Then many kinds of PS were designed and studied in detail [10–19].
Actually there is only one focus at a single focal plane no matter for FZP or PS. In 2013, an aperiodic Fibonacci sequence was introduced into the planar diffractive optical elements (DOEs) and generated two foci with a definite relationship along the optical axis [20–24]. The other interesting array illuminator is the Dammann grating [25–28]. We find that the former functionality is achieved by quasi-periodic Fibonacci structure, but the latter Dammann grating is a periodic structure. In our previous work, we proposed the generalized Fibonacci PS, and produced two adjustable foci on axis according to the mathematical character of generalized Fibonacci sequences [24]. This technology can produce multiple foci on the optical axis.
In this paper, we propose a new family of mathematical sequence named by Greek ladder [29–31], and introduce the multi-dimensional quasi-periodic Greek ladder to nano-photonic device, furthermore design the corresponding device to generate two, three and four foci on axis, respectively. By analogy, more than four-focal device can be designed. Meanwhile the focal distances can be freely adjusted. Then under some special conditions, the two-dimensional quasi-periodic classical Greek ladder can be reduced to one-dimensional quasi-periodic generalized Fibonacci sequences. By means of proper encoding and binary, ternary or quaternary phase modulation, the Greek-ladder devices do generate the diffraction-limited array foci at different multi-focal planes. Essentially, three-dimensional array focusing can be achieved because the periodic and symmetric device is replaced by the aperiodic and asymmetric structure which leads to the advent of quasi-steady intermediate field between the near field and the far field.
2. Mathematical description of Greek ladders
For traditional FZP, the optical path difference (OPD) between two adjacent zones is equal to λ/2, where λ is the incident wavelength. If the OPD λ/2 is replaced with Kλ/2 where K is a positive real number and denotes the optical path difference scaling factor (OPDSF), we can get the following OPD in the case of a monochromatic plane wave incidence
Where f0 is the expected focal length, wavelength λ; rm is the radius of mth zone; K denotes the OPDSF.When a light source is incident on the optical device, the diffracted field can be given by Rayleigh-Sommerfeld diffraction integral
Where R = [(x-ζ)2 + (y-η)2 + z2]1/2, E0(ζ,η,0) denotes the incident light, t(ζ,η) is the corresponding transmission function of each transparent region.Now let us introduce the Greek ladder which has the initial seeds
Where αj is a real number.And has the following linear recursion relations
Where n is a positive integer. rn,j(j = 1,2,…N) denotes the nth Greek-ladder vector.We can easily get the math limit of Eq. (3) by the following definition
In order to determine the initial seeds of Greek ladder, we introduce the best asymptotical nth convergent by use of the finite continued fraction. In number theory, any one of irrational number can be approximated by the continued fraction. So for the classical Greek ladder, we can get the best asymptotic fraction by the method proposed by Luogeng Hua [32].
Obviously, the nth convergent approximates the value of. In a particular case r1,j = 1, j = 1, 2, the Greek ladder can approximate the simple continued fractionFrom the above analysis, we can come to a conclusion that a continued fraction of this type is eventually periodic of period 2. There is a possibility that the Greek-ladder device can be divided into two incommensurable periods in the case of proper encoding.
3. Longitudinal focusing of Greek-ladder devices
Figure 1 shows the schematic of a Greek-ladder nano-photonic device, and Table 1 gives the first 5th Greek-ladder terms. The diameter of each pinhole is equal to the width of the zone where the pinhole is located. The 1st and the 3rd quadrants have 94 opaque and 103 transparent zones respectively, while the 2nd and the 4th quadrants have 101 opaque and 111 transparent zones respectively. The total transmission function is the linear superposition of each transparent area. The calculation parameters are as follow: the incident wavelength equals to 632.8nm, K = 0.5, the expected focal length f0 = 3.5cm, C = 1.25, the math limit , the initial seeds (α1,α2,α3) = (1,2,3) and the encoded seeds (α1,α2,α3) = (1,01,010). The two focal distances are f1 = 3.567cm and f2 = 3.439cm, respectively. Obviously, the ratio of the two focal distances equals to f1/f2≈1.0372 approximating.
When r1,j = αj, j = 1,2, the above expression is reduced to the traditional Greek ladder or Theon’s ladder which (circa 140 A.D.) had described a remarkably simple way to calculate rational approximations to. Actually due to the concatenation operation of the encoded seeds, the former term is the subset of the latter one, so we can choose any column as the encoded switching. The simulation parameters are similar to the calculation in Fig. 1. For this condition, Greek-ladder zone plate (ZP) with square root of 1.25 has 85 transparent and 95 opaque zones. The initial seeds are (α1,α2) = (1,1). The result is illustrated in Table 2 and Fig. 2. The two focal distances are f1 = 3.707cm and f2 = 3.315cm, respectively. The ratio of the two focal distances equals to f1/f2≈1.1182 approximating. The other results of Greek-ladder devices with square root of 2, 3 and 4 are also shown in Table 2 and Fig. 2. To clarify, the maximum intensities of axial foci for four cases are nearly equal.
The number of axial foci will increase with the increment of the OPDSF, and the emerging foci satisfy the following Eq. (7). Taking the Greek-ladder device with square root of 2 into account, the results are indicated in Table 3 and Fig. 3. Although the values of K are different, the maximum intensity on axis keeps the same value under the condition of the plane wave incidence. The similar conclusions are also applied to other cases mentioned above with respect to the different OPDSF.
Where [x] denotes the nearest integers less than or equal to x, mod(X,2) represents the modulus after division 2, and m denotes the mth axial focus.Equation (7) represents that there are another two techniques to generate multiple equidistant foci with equal intensity on axis. One is to scale proportionally focal length as required by OPDSF. The other is to alter the expected focal length while keeping the same OPDSF. Taking trifocal distances into account, the Greek-ladder device with square root of 2 has two focal distances which equals to f4 = 2.987cm and f2 = 4.225cm, respectively. The mean of the two focal distances is 3.606cm. Compared with Greek-ladder device with K = 0.5 in the 1st and the 3rd quadrants, the value of K equals to 3.606/3.5*0.5 = 0.5151 in the 2nd and the 4th quadrants. As for four-focal distances, the Greek-ladder device with square root of 2 is also arranged in the 1st and the 3rd quadrants. If we expect that the third focal length f3 is 3.606cm, then the fourth focal length f1 satisfies the relation f1-f2 = f3-f4. The value of f1 is 4.844cm. On the basis of Eq. (7), the expected focal length is 4.1343cm. In the 2nd and the 4th quadrants, we can obtain the Greek-ladder device with the square root of C = (f1/f3)2 = 1.8045. If so, the four-focal nano-photonic device has been successfully designed. The results of the two examples are all given in Table 4 and Fig. 4, and the axial foci are located at 2.984cm, 3.635cm, 4.221cm and 4.783cm, respectively.
4. Transverse focusing at multi-focal planes of Greek-ladder devices
Based on the concept of Fourier optics, the different phase stands for the different angular spectrum in the spatial domain. In this case, when the amplitude modulation is replaced by the phase modulation in the right way, the single focal spot may be split into multiple spots. Figure 5 shows the schematic of the Greek-ladder device with three-dimensional array foci. The value of zero still denotes the opaque pinhole, but the transparent pinhole is completely replaced with binary phase. The red pinhole denotes the zeros phase, while the blue pinhole denotes the π phase. For Cartesian Greek-ladder device, the transparent pinholes is set to zeros phase if the absolute value of the center of pinhole is less than the quarter full width of device; otherwise, the transparent pinholes is set to π phase. For polar Greek-ladder device, the transparent pinholes is set to zeros phase if the polar radius of pinhole is less than the 0.574 times the radius of device; otherwise, the transparent pinholes is set to π phase. The intensity distributions corresponding to the model are illuminated in Figs. 6(a)-6(d). The Greek-ladder device does generate 2*2 array foci at two focal planes. The locations of the two focal planes are at 2.987cm and 4.225cm, respectively. For polar Greek-ladder device, the two focal planes are at the same location.
In order to get the anisotropic focal spots, the area of zeros phase along the x-axis must be different from that along the y-axis. Taking Cartesian Greek-ladder device into account, the transparent pinholes is set to zeros phase when the absolute value of the center of pinhole is less than the quarter full width of device along the x-axis and 0.35 or 0.40 times full width of device along the y-axis. The array foci are shown in Fig. 7. As the increment of the area of zeros phase along the y-axis, the two focal spots gradually blend together along the vertical direction.
For direct comparison to (2*2)*2 array foci under the same model in Cartesian coordinates, the transparent pinholes is set to zeros phase if the absolute value of the center of pinhole is less than the 0.325 times full width of device; otherwise, the transparent pinholes is set to π phase. In addition, the packed density is set to the value of 0.9 in the area of zeros phases. Figures 8(a) and 8(b) shows the intensity contour at two focal planes. The two focal planes are located at 2.987cm and 4.225cm, respectively. Again for the model of Fig. 5(b), the transparent pinholes is set to zeros phase if the polar radius of pinhole is less than 0.26 times and larger than 0.66 times the radius of device; otherwise, the transparent pinholes is set to π phase. Three-ring spots are obtained at two focal planes located at 2.987cm and 4.225cm, respectively [see Figs. 8(c) and 8(d)].
5. Three-dimensional array foci from Greek ladder to generalized Fibonacci
Finally let us consider a special form of classical Greek ladder,
The characteristic equation of Eq. (8) is a quadratic equation. That isFrom Eqs. (8) and (9), we know .Based on the algebraic equation, we can easily get another equivalent recursion relation named generalized Fibonacci sequence
In particular, all the methods in the analysis process of Greek-ladder devices are all applied to the generalized Fibonacci structures.In Fig. 7, the anisotropic focal spots can be achieved by adjusting the area of zeros phases in Cartesian coordinates. Here we propose another method to get the anisotropic focal spots by multiple Fibonacci sequences distributed in different regions in polar coordinates. One Fibonacci device with (m,k) = (1,1) has 92 transparent and 149 opaque zones in the 1st and the 3rd quadrants, the other Fibonacci device with (m,k) = (−2,0.3) has 77 transparent and 164 opaque zones in the 2nd and the 4th quadrants. The initial seeds are all (2,3). In the 1st and 3rd quadrants, the transparent pinholes is set to zeros phase if the absolute value of the center of pinhole is less than the 0.307 times full width of device; otherwise, the transparent pinholes is set to π phase. The same principle also applies to the quadrants (2,4) except that (0, π) phases are replaced with the corresponding (0.5π,1.5π). The four focal planes are located at 2.567cm, 2.830cm, 4.584cm and 5.494cm, respectively. Figures 9(a)-9(d) shows the intensity contour at four focal planes, and there are four foci at each focal plane. When the transparent pinholes in quadrants (2,4) are all set to the value of 0.5π phase and the quadrants (1,3) remain the same, at this time we successfully obtain different number of foci at each focal planes as shown Figs. 9(e)-9(h). It should be pointed out that when the 0.5π phase in quadrants (2,4) is replaced with 1.5π phase, the spots at each focal plane will be integrally rotated by 90 degree, but the diffraction pattern remain the same.
Taking the diffraction efficiency into account, the amplitude-only DOE is worse than the phase-only DOE. For binary optical element (BOE) with phase modulation, if the opaque areas in the same partition are replaced with the opposite phase, the diffraction efficiency will be improved greatly. In addition, it should be noted that the phase modulation in visible wavelengths cannot be fit for EUV and soft x-rays in that their wavelengths are too short. Fortunately, the optical theory tells us that the phase equals to the product of the wave number and its own OPD. Obviously, the OPD of two pinholes, which are located at the two adjacent Fresnel zone, is equal to half wavelength. So we can move the pinholes in odd Fresnel zone to that in even Fresnel zone in the partition as required, and vice versa. In this way, phase modulation for EUV and soft x-rays can be realized by amplitude modulation.
Furthermore, the characteristic equation of Fibonacci sequence is a quadratic equation; the math limit is proportional to the square roots. While for Greek ladder, the math limit is proportional to the nth roots. It means that the Greek ladder converges at a faster rate than Fibonacci sequences under the same condition. At this point the device with a very short distance between two focal planes can be easily constructed by Greek ladder. The submicro-to-micro sized distance between two focal planes is more suitable for biological cell, microscopy, lithography, etc.
6. Conclusion
We have proposed a kind of new technology to generate array foci at the multi-focal planes. Different from the mono-focal Fibonacci diffractive lenses with the ratio of two focal distances approaching the golden mean, the Greek ladder is introduced into the diffractive optics for the first time, and successfully produces the multiple foci on axis with the ratio of them approaching the nth root of C, where C is a positive real number. Meanwhile we can be free to adjust the focal locations as required by OPDSF. A special form of the Greek ladders can be reduced to the generalized Fibonacci sequences. With these technologies and the multi-nary phase coding, we can be free to design various kinds of multi-focal nano-photonic devices with any number and any location. In visible wavelengths, this technology has many applications in ophthalmology, virtual display technique, array imaging, laser micromachining and spectroscopy. In EUV and soft x-ray region, it can be applied for biological cell, microscopy, lithography, etc.
Acknowledgments
This work is supported by National Natural Science Foundation of China (NSFC) (No. 61205212).
Reference and links
1. D. J. Stigliani Jr, R. Mittra, and R. G. Semonin, “Resolving power of a zone plate,” J. Opt. Soc. Am. 57(5), 610–613 (1967). [CrossRef]
2. E. Champagne, “Optical method for producing Fresnel zone plates,” Appl. Opt. 7(2), 381–382 (1968). [CrossRef] [PubMed]
3. M. D. Tipton and J. E. Dowdey, “Coded aperture imaging with on-axis Fresnel zone plates,” Opt. Eng. 12(5), 166–168 (1973). [CrossRef]
4. H. Kyuragi and T. Urisu, “Higher-order suppressed phase zone plates,” Appl. Opt. 24(8), 1139–1141 (1985). [CrossRef] [PubMed]
5. J. A. Sun and A. Cai, “Archaic focusing properties of Fresnel zone plates,” J. Opt. Soc. Am. 8(1), 33–35 (1991). [CrossRef]
6. G. Saavedra, W. D. Furlan, and J. A. Monsoriu, “Fractal zone plates,” Opt. Lett. 28(12), 971–973 (2003). [CrossRef] [PubMed]
7. J. A. Davis, S. P. Sigarlaki, J. M. Craven, and M. L. Calvo, “Fourier series analysis of fractal lenses: theory and experiments with a liquid-crystal display,” Appl. Opt. 45(6), 1187–1192 (2006). [CrossRef] [PubMed]
8. V. Ferrando, A. Calatayud, F. Giménez, W. D. Furlan, and J. A. Monsoriu, “Cantor dust zone plates,” Opt. Express 21(3), 2701–2706 (2013). [CrossRef] [PubMed]
9. L. Kipp, M. Skibowski, R. L. Johnson, R. Berndt, R. Adelung, S. Harm, and R. Seemann, “Sharper images by focusing soft X-rays with photon sieves,” Nature 414(6860), 184–188 (2001). [CrossRef] [PubMed]
10. Q. Cao and J. Jahns, “Nonparaxial model for the focusing of high-numerical-aperture photon sieves,” J. Opt. Soc. Am. A 20(6), 1005–1012 (2003). [CrossRef] [PubMed]
11. G. Andersen, “Large optical photon sieve,” Opt. Lett. 30(22), 2976–2978 (2005). [CrossRef] [PubMed]
12. F. Giménez, J. A. Monsoriu, W. D. Furlan, and A. Pons, “Fractal photon sieve,” Opt. Express 14(25), 11958–11963 (2006). [CrossRef] [PubMed]
13. H.-H. Chung, N. M. Bradman, M. R. Davidson, and P. H. Holloway, “Dual wavelength photon sieves,” Opt. Eng. (Bellingham) 47(11), 118001 (2008). [CrossRef]
14. G. Cheng, C. Hu, P. Xu, and T. Xing, “Zernike apodized photon sieves for high-resolution phase-contrast x-ray microscopy,” Opt. Lett. 35(21), 3610–3612 (2010). [CrossRef] [PubMed]
15. J. Jia and C. Xie, “Phase zone photon sieve,” Chin. Phys. B 18(1), 183–188 (2009). [CrossRef]
16. C. Zhou, X. Dong, L. Shi, C. Wang, and C. Du, “Experimental study of a multiwavelength photon sieve designed by random-area-divided approach,” Appl. Opt. 48(8), 1619–1623 (2009). [CrossRef] [PubMed]
17. C. Xie, X. Zhu, L. Shi, and M. Liu, “Spiral photon sieves apodized by digital prolate spheroidal window for the generation of hard-x-ray vortex,” Opt. Lett. 35(11), 1765–1767 (2010). [CrossRef] [PubMed]
18. A. Sabatyan and S. Mirzaie, “Efficiency-enhanced photon sieve using Gaussian/overlapping distribution of pinholes,” Appl. Opt. 50(11), 1517–1522 (2011). [CrossRef] [PubMed]
19. X. Zhao, F. Xu, J. Hu, and C. Wang, “Broadband photon sieves imaging with wavefront coding,” Opt. Express 23(13), 16812–16822 (2015). [CrossRef] [PubMed]
20. J. A. Monsoriu, A. Calatayud, L. Remón, W. D. Furlan, G. Saavedra, and P. Andrés, “Bifocal Fibonacci diffraction lenses,” IEEE Photonics J. 5(3), 3400106 (2013). [CrossRef]
21. A. Calatayud, V. Ferrando, L. Remón, W. D. Furlan, and J. A. Monsoriu, “Twin axial vortices generated by Fibonacci lenses,” Opt. Express 21(8), 10234–10239 (2013). [CrossRef] [PubMed]
22. V. Ferrando, A. Calatayud, P. Andrés, R. Torroba, W. D. Furlan, and J. A. Monsoriu, “Imaging properties of Kinoform Fibonacci lenses,” IEEE Photonics J. 6(1), 6500106 (2014). [CrossRef]
23. J. Ke, J. Zhang, and J. Zhu, “Focusing properties of modified Fibonacci photon sieve,” Chin. Opt. Lett. 13(8), 080501 (2015). [CrossRef]
24. J. Ke and J. Zhang, “Generalized Fibonacci photon sieves,” Appl. Opt. 54(24), 7278–7283 (2015). [CrossRef] [PubMed]
25. H. Dammann and E. Klotz, “Coherent optical generation and inspection of two-dimensional periodic structures,” Opt. Acta (Lond.) 24(4), 505–515 (1977). [CrossRef]
26. C. Zhou and L. Liu, “Numerical study of Dammann array illuminators,” Appl. Opt. 34(26), 5961–5969 (1995). [CrossRef] [PubMed]
27. C. Zhou, J. Jia, and L. Liu, “Circular Dammann grating,” Opt. Lett. 28(22), 2174–2176 (2003). [CrossRef] [PubMed]
28. J. Yu, C. Zhou, W. Jia, W. Cao, S. Wang, J. Ma, and H. Cao, “Three-dimensional Dammann array,” Appl. Opt. 51(10), 1619–1630 (2012). [CrossRef] [PubMed]
29. J. R. Ridenhour, “Ladder Approximations of Irrational Numbers,” Math. Mag. 59(2), 95–105 (1986). [CrossRef]
30. T. J. Osler, M. Wright, and M. Orchard, “Theon’s ladder for any root,” Int. J. Math. Educ. Sci. Technol. 36(4), 389–398 (2005). [CrossRef]
31. C. Horgan and K. Herzinger, “A further investigation of using Theon’s ladder to find roots of quadratic equations,” Int. J. Math. Educ. Sci. Technol. 45(1), 150–158 (2014). [CrossRef]
32. L. G. Hua, An Introduction to Number Theory (Science Press, 1979).