Tunable diffraction-free array in nonlinear photonic crystal.

Diffraction-free beams have attracted increasing research interests because of their unique performances and broad applications in various fields. Although many methods have been developed to produce such beams, it is still challenging to realize a tunable non-diffracting beam. Here, we report the generation of a tunable diffraction-free array through second-harmonic generation in a nonlinear photonic crystal, i.e., a 2D periodically-poled LiTaO3 crystal. In such a crystal, the second-harmonic wave is engineered by properly designing the domain structure based on the Huygens-Fresnel principle. The characteristics of the generated diffraction-free array including its period, propagation length, and wavelength can be tuned by simply changing the input wavelength. Our observation not only enriches the diffraction-free optics, but also has potential applications for photolithography and imaging.

Diffraction, originating from the Helmholtz equation, has long been considered as a universal characteristic of all classical waves. However, Durnin et al. reported in 1987 an exact diffraction-free mode solution of the Helmholtz equation12, which has a transverse intensity distribution independent of the propagation distance. The first experimental demonstration was a nearly non-diffracting Bessel beam1. Since then, investigations of such diffraction-free beams and their applications in metrology3, nonlinear optics4, atomic optics5, optical micro-manipulation678, medical imaging9, electron microscope10, and wireless optical communications11 have become an active research area. Besides this Bessel beam, other non-diffracting solutions1213 including the Airy beam1415 were also discovered. So far, the diffraction-free beams are mainly generated through linear optical methods, such as Fabry-Pérot interferometer16, spatial light modulator17, holographic process18, diffractive phase elements19, axicon20 and surface plasmon polariton (SPP)2122. The experiments using nonlinear optical techniques are less reported232425. In most of these methods, the performances of the generated diffraction-free beam, such as its wavelength, beam size, and propagating length, are fixed in the devices. The few tunable methods require certain complicated instruments like spatial light modulators17. In this Letter, we propose and demonstrate a novel nonlinear optical method to produce a tunable diffraction-free array of beams in a single nonlinear photonic crystal, i.e., a two-dimensional (2D) periodically-poled LiTaO3 (PPLT) crystal.

PPLT crystals have been extensively investigated because they can realize highly-efficient frequency conversions through the quasi-phase-matching (QPM) technique26. Since the concept of nonlinear photonic crystal, i.e., 2D PPLT crystal, was proposed by Berger27 in 1998, numerous interesting phenomena have been discovered such as non-collinear second-harmonic generation (SHG)28, nonlinear Čerenkov radiations293031, and nonlinear Talbot self-imaging32. Recently, domain engineering in nonlinear photonic crystals for spatial light modulation attracts an increasing research interest. Scientists have developed various domain structures to realize dual-focused second-harmonic (SH) spots33, conical SHG3435, optical orbital angular momentum states3637, beam shaping3839, and superfocusing40. By utilizing the domain-engineering method based on the Huygens-Fresnel principle, we can design the desired domain structure to realize certain tunable diffraction-free SH array in a single PPLT chip. Our results not only extend the concept of diffraction-free optics, but also open a door for broader applications of non-diffracting beams in photolithography and imaging.

Results

Theory

The idea is intrigued by the free-space Bessel beam, in which the diffraction-free field can be decomposed into plane-wave components with wave vectors on a cone1. Similarly, two plane waves can form a diffraction-free array with a cosine transverse profile, which has been demonstrated in SPPs21. Such cosine beam can be considered as the 2D counterpart of the Bessel beam21. Here, we produce a tunable diffraction-free array through a SHG process in a 2D PPLT crystal as shown in Fig. 1a. The coupled-wave equation can be written as

Figure 1

Experimental configuration.

The experimental setup is shown in (a). The input laser propagates along the y-axis of the squarely-poled LiTaO3 crystal. The generated SH pattern is recorded by a CCD camera. The diffraction-free cosine beam in (b) can be decomposed into two components (c).

where E1 and E2 are the electrical fields of the fundamental wave and SH wave, respectively; k2 is the wave vector of the SH wave; K is the coupling coefficient; and t(x, y) is the structural function of the PPLT crystal. We seek a diffraction-free cosine solution for the SH wave (Fig. 1b), which can be written as

Here, A is a constant, k2x and k2y are the x and y components of the SH wave vectors, respectively. The electrical field distribution along z direction is uniform in this case. The cosine beam described by Eq. (2) can be decomposed into two plane-wave components (Fig. 1c). Obviously, such solution has a transverse intensity profile independent of the propagation direction y, which represents a diffraction-free SH array. Based on Eqs (1) and (2), we use the nonlinear Huygens-Fresnel principle33 to engineer the domain structure t(x, y) for the realization of such beam (see Methods for the details). As well known, it is impossible to experimentally realize an ideal Bessel beam because it carries infinite energy1. Alternatively, one can generate a Gaussian-Bessel beam (i.e. the Bessel solution modulated by a Gaussian envelope) in experiment, which preserves the diffraction-free properties in the paraxial approximation41. Because the non-diffracting solution in Eq. (2) suffers from the same problem, we introduce the cosine-Gaussian beam21 in the experiment.

The cosine beam consisting of two SH components can be understood through another view point, i.e., QPM. In a typical QPM configuration, SHG can be greatly enhanced by using the reciprocal vectors in a PPLT crystal to compensate for the phase mismatch between the fundamental wave and the SH wave26. Luckily, there exist abundant non-collinear reciprocal vectors in a 2D domain structure, which can realize non-collinear QPM SHG2728 as shown in Fig. 2a. For instance, in a squarely-poled LiTaO3 crystal as shown in Fig. 1a, the reciprocal vectors (Fig. 2b) are defined by

Figure 2

Noncollinear QPM.

The diffraction-free cosine beam can be understood from the non-collinear QPM configuration (a). The reciprocal vectors G and G can produce two non-collinear SH waves, which can be considered as the decomposed components of the cosine beam. The reciprocal vectors in a squarely-poled PPLT crystal are shown in (b).

where the subscripts m and n denote the orders of the reciprocal vector along the longitudinal and transverse directions, respectively. is the period of the domain structure. The QPM condition under non-collinear configuration (Fig. 2a) requires

where k1 and k2 are the wave vectors of the fundamental and SH waves, respectively. Interestingly, Gm,n and its mirror-symmetrical vector Gm,−n can simultaneously generate two SH waves as shown in Fig. 2a, which can be considered as the decomposed components of the cosine beam in Eq. (2). They interfere with each other and result in a diffraction-free SH array as shown in Fig. 1b and c. Considering that the two SH components of the cosine beam have an in-between angle of 2θ (decided by Equation (4)), the transverse profile of the SH intensity (Fig. 1b) can be easily deduced to be

Obviously, the period of the obtained SH array is

From Fig. 1c, the propagation distance of the diffraction-free beam can be written as

where D is the diameter of the SH component. Unlike other diffraction-free schemes which are usually fixed by the sample structures, a given PPLT structure can actually generate diffraction-free arrays with varied periods and propagation distances by involving different Gm,n. This can be easily realized by changing the input wavelength or tuning the operation temperature. In our scheme, the prerequisite to realize diffraction-free is to satisfy the non-collinear QPM condition. If not fulfilled, one cannot obtain the diffraction-free array because no non-collinear SH beams are efficiently generated.

Experimental demonstration of the tunable diffraction-free array

The PPLT crystal for the generation of a diffraction-free SH array is designed to have a squarely-poled structure with a period of Λ = 5.5 μm (Fig. 1a). A Ti:Sapphire femtosecond laser serves as the input fundamental field, which can be continuously tuned from 690 nm to 1050 nm in wavelength. This fundamental beam is first reshaped to produce a near-parallel beam. Then, it travels along the y direction with its polarization parallel to the z-axis of the crystal (Fig. 1a). The coordinate system is set according to the crystal axis. Under this experimental configuration, the involved nonlinear optical coefficient d33 is the biggest one in the LiTaO3 crystal. A short-pass filter is placed after the crystal to block the fundamental field. The SH patterns near the PPLT crystal are magnified by a 100 objective lens with a N.A. = 0.7 and then recorded by a CCD camera. By moving the objective lens along the y direction, we can investigate the diffractive characteristics of the SH patterns.

The input laser is first set to be 906 nm. At this wavelength, the non-collinear SHG can be phase-matched with the reciprocal vectors G1,3 and G1,−3 (Fig. 2). At the output face of the PPLT crystal, one can observe a SH array with a period of 0.92 μm as shown in Fig. 3a, which is well consistent with the theoretical period of 0.917 μm from Eq. (6) with |n| = 3. The corresponding numerical simulation based on the Huygens-Fresnel principle is shown in Fig. 3d. The period of the simulated SH pattern is 0.93 μm. The small deviation may result from that the dispersion relation of the LiTaO3 crystal used in the calculation does not perfectly match our sample. By moving the objective lens along the y direction, we can record the SH patterns at different observation planes. Figure 4a shows the measured evolution of the SH carpet within y = 200 μm, which clearly presents the diffraction-free performance. As propagating along the y direction, the intensity of the SH array decreases because of the cosine-Gaussian mode; however, the array period does not change. The numerical simulation in Fig. 4b also confirms such diffraction-free behavior. For simplicity, we have assumed a plane-wave illumination in simulations, which cannot predict the attenuation of the SH intensity along the propagation direction in Fig. 4a. As shown in Fig. 4c and d, the standard deviation from the theoretical diffraction-free array increases from 1.37 × 10−4 to 2.29 × 10−4 as the experimentally generated cosine-Gaussian beam propagates from y = 25 μm to y = 198 μm. It should be noted that the SH array presents such diffraction-free performance only near the center of the whole picture (within the area of ~100 μm × 10  μm in our experiment) because of the Gaussian modulation of the Bessel solution.

Figure 3

Diffraction-free SH arrays generated by different fundamental wavelengths.

The measured (a–c) and simulated (d–f) cross sections of the diffraction-free SH arrays at certain observation planes. The periods of the array in the experiment are 0.92 μm, 1.38 μm, and 2.75 μm at the fundamental wavelengths of 906 nm (a), 928 nm (b) and 944 nm (c), respectively, which are well consistent with the corresponding numerical simulations.

Figure 4

Diffraction-free “carpets”.

Experimental (a) and theoretical (b) diffraction-free “carpets” along the propagation direction are obtained at a 906 nm input laser. (c) and (d) show the deviations of the measured beam profile from an ideal non-diffracting array at propagation distances of 25 μm and 198 μm, respectively.

Next, we change the input laser wavelength to 928 nm and 944 nm, respectively, to tune the diffraction-free array. The fundamental beam power is kept at 50 mW for all the wavelengths. As shown in Fig. 3b and c, the SH arrays change dramatically comparing to the pattern excited by a 906 nm fundamental beam (Fig. 3a). The period of the SH pattern at the pump wavelength of λ = 928 nm is 1.38 μm (Fig. 3b), which is one quarter of the domain period. When further increasing the input wavelength to λ = 944 nm, the array period becomes 2.75 μm (Fig. 3c). The dependence of the period of the diffraction-free array on the wavelength originates from the involvement of different reciprocal vectors, and therefore G1,2/G1,−2 (Fig. 2) at 928 nm and G1,1/G1,−1 (Fig. 2) at 944 nm correspond to 1/4 and 1/2 of the domain period, respectively, according to Eq. (6). The numerical simulations based on Huygens-Fresnel principle for these two cases are shown in Fig. 3e and f, which are well in agreement with the experimental results. Our measurement shows that the non-diffracting SH pattern can be observed at a distance of up to 4.5 mm away from the sample with an input wavelength of 944 nm. This is slightly shorter than the theoretically predicted propagation length of 4.8 mm from Eq. (7) with D = 400 μm. From the experimental images (Fig. 3a–c), one can easily see that the intensity of the SH pattern increases when the laser is tuned to a longer wavelength, which is mainly caused by the more effective nonlinear coefficient. Usually, a lower-order Gm,n has a higher effective nonlinear coefficient, which can realize a higher SHG efficiency. At a non-QPM wavelength, we can hardly observe a diffraction-free array pattern in the experiment because the prerequisite condition has been broken.

In principle, for each pair of Gm,n and Gm,−n in a 2D PPLT crystal, one can always find a suitable wavelength to satisfy the non-collinear QPM condition and then to generate the diffraction-free array. However, this might not be realizable in experiment because (1) the high-order reciprocal vector may have an effective nonlinear coefficient which is too small to efficiently generate the SH waves; (2) the above prediction is only valid under the paraxial approximation, which rules out its applicability to the reciprocal vectors with big subscript n. To design a practical PPLT crystal for generating such non-diffracting arrays, it is important to suppress the collinear SHG process because the unwanted background could completely ruin the diffraction-free beam pattern. Usually, the input wavelength should be chosen as far as possible away from the QPM wavelength for the collinear SHG process.

Discussion

In conclusion, we have presented the generation of a tunable diffraction-free array, i.e. cosine-Gaussian beam, through non-collinear QPM SHG processes in a nonlinear photonic crystal. Beyond the previously demonstrated techniques, this work has extended the generation of diffraction-free beam in two fronts. First, the introduction of SHG produces a non-diffraction array at a shorter wavelength, which can have potential applications in photolithography and optical imaging. The resolution of the beam is improved by a factor of 2 due to SHG. Second, one can easily tune several characteristics of the diffraction-free cosine-Gaussian beam in a single chip, which makes it more convenient to utilize it in integrated photonic devices. Here, we have demonstrated the wavelength- and range-tunable non-diffracting array in a PPLT crystal by using varied input wavelength. Actually, there are more tools to modulate and optimize the generated array in the PPLT crystals. For example, by utilizing the excellent thermo-optical, electro-optical and acoustic-optical performances of the LiTaO3 crystal, the diffraction-free array can be modulated by changing the operation temperature, and applying an electrical or acoustic field. The performance of the non-diffracting array can be further improved through additional domain engineering techniques. For instance, it has been experimentally shown that chirped, ring-shaped and quasi-periodic structures can greatly enhance the tunable properties of the QPM processes4243. One can utilize these structures to further tune and modify the non-diffracting array. Most importantly, such nonlinear photonic crystals provide a useful integrated platform to manipulate the propagations of diffraction-free beam arrays and other spatial light beams for their potential applications in photolithography and optical imaging.

Methods

Analytical expressions of the fundamental field E 1 and the SH field E 2 for numerical simulations

The numerical simulations are performed by using the Huygens-Fresnel principle33, in which each part of the crystal is considered as a point source which emits the SH wave. There is a π phase-shift between the SH waves generated from positive and negative domains44. The input fundamental beam propagates along the y axis of the crystal (Fig. 1). Under the slowly-varying-envelope approximation, the evolutions of the fundamental field E1 and the SH field E2 in the squarely-poled LiTaO3 crystal can be described by

where is the coupling coefficient with d33 being the nonlinear coefficient of the crystal. is the phase mismatch between the fundamental wave and the SH wave. From Eq. (8), one can write the difference equations as

Here, Δx and Δy are the space steps. Then, we apply the finite difference method45 to calculate SHG process in the 2D PPLT crystal. It should be noted that the above equations are valid in the paraxial approximation.

Additional Information

How to cite this article: Liu, D. et al. Tunable diffraction-free array in nonlinear photonic crystal. Sci. Rep. 7, 40856; doi: 10.1038/srep40856 (2017).

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.