Jun 16, 2023
Polarization
Scientific Reports volume 13, Article number: 10338 (2023) Cite this article 613 Accesses Metrics details Dynamic phase-only beam shaping with a liquid crystal spatial light modulator is a powerful
Scientific Reports volume 13, Article number: 10338 (2023) Cite this article
613 Accesses
Metrics details
Dynamic phase-only beam shaping with a liquid crystal spatial light modulator is a powerful technique for tailoring the intensity profile or wave front of a beam. While shaping and controlling the light field is a highly researched topic, dynamic nonlinear beam shaping has hardly been explored so far. One potential reason is that generating the second harmonic is a degenerate process as it mixes two fields at the same frequency. To overcome this problem, we propose the use of type II phase matching as a control mechanism to distinguish between the two fields. Our experiments demonstrate that distributions of arbitrary intensity can be shaped in the frequency-converted field at the same quality as for linear beam shaping and with conversion efficiencies similar to without beam shaping. We envision this method as a milestone toward beam shaping beyond the physical limits of liquid crystal displays by facilitating dynamic phase-only beam shaping in the ultraviolet spectral range.
The first operating laser in the early 1960s1 was the dawn for many research fields in modern optics although some of their fundamental effects were already demonstrated or proposed in theory decades earlier. Holography and nonlinear optics emerged independently from each other but both fields benefited from the high coherence and high power of new light sources.
Holography is based on the interference of light waves and incorporates phase and amplitude information to go beyond photography. Dynamic phase-only beam shaping with a liquid crystal spatial light modulator (LC-SLM) is a method emerging from holography for arbitrarily controlling the intensity distribution of the beam with many applications in research2,3,4 and industry5,6,7. As this method only modulates the wave front, there are no significant losses. As a drawback, liquid crystal displays are technically limited to the visible, near-infrared and mid-infrared spectral ranges. This is not an insurmountable problem, as frequency conversion processes such as second harmonic generation or sum frequency generation are coherent processes that preserve the phase of the impinging fundamental wave. Combining nonlinear optics and holography allows for shaping the light field at the fundamental while achieving the targeted outcome at the frequency-converted field. Even though both research fields can be combined, the concept of nonlinear holography is only currently emerging.
Yariv showed decades ago that four wave mixing can be interpreted as holographic recording and reconstruction and proposed using it for the realization of real time holography8. In this process, the interaction between the fields can be interpreted as one field diffracted by the shaped pattern of another field. Meanwhile many investigations followed on the nonlinear conversion of structured light for the conservation of singularities9,10 and orbital or spin angular momentum and vortex beams11,12,13,14,15,16,17,18,19,20,21,22,23,24. Here, the physical principles of the nonlinear conversion of structured light have been well-explored and several works use beam characteristics such as polarization19,25,26, differing wavelengths15,17,18 or non-collinear geometries20 as a control mechanism for the nonlinear conversion of optical vortices. The recent review paper published by Buono and Forbes gives an overview on nonlinear optics with structured light27.
Currently, there are two major approaches to nonlinear holography: directly structuring the nonlinear crystal or imaging the plane of an LC-SLM into the crystal.
3D structuring of the nonlinear crystal leads to a modulation of the nonlinear susceptibility which shapes the wave front of the emerging light field. Such elements are called nonlinear photonic crystals as the modulation of the nonlinear susceptibility affects the beam generation and propagation28,29,30,31,32,33,34. There are demonstrations of a binary hologram in a nonlinear crystal35, a structured element combined with structured light36,37 or plasmonic metasurfaces38,39. Such 3D structured nonlinear crystals act as volume holograms or phased arrays and in theory give more degrees of freedom as a thin hologram. As their implementation is technically challenging, the freedom of design is so far strongly limited and furthermore only static solutions are possible. Such practical limitations motivate the consideration of thin holograms which are easier to realize.
In a dynamic approach, the plane of a phase-only LC-SLM is directly imaged into the nonlinear crystal40,41,42,43. It thus acts as a thin hologram and brings the SLM and the nonlinear crystal together in a common plane, where the SLM dynamically shapes the wave front of the fundamental as if the hologram was directly structured into the crystal. The phase mask on the SLM imprints a locally varying phase on the incoming light field which leads to a modification of the resulting wave front where the slope corresponds to wave vectors forming the angular spectrum. During the nonlinear interaction, all light fields which fulfill the phase matching condition interact and create a new, frequency-converted outgoing wave. If the incoming fields are not distinguishable, i.e. the fields are degenerate, all wave vectors will add up, resulting in a new outgoing wave that has an angular spectrum, containing all possible sums of incoming wave vectors. This complicates beam shaping of the outgoing, frequency-converted wave as the initially applied angular spectrum at the SLM is not conserved. Consequently, there needs to be a control mechanism which makes the incoming fields distinguishable.
Current approaches for second harmonic generation40,41,42,43 are based on a degenerate type I phase matching process and thus involve a non-collinear geometry to control the beam-shaped outcome. This imposes severe limitations in the quality and the efficiency of the achievable results. We utilize type II phase matching and demonstrate high-quality and highly-efficient nonlinear beam shaping in the second harmonic with a LC-SLM in a collinear geometry:
Basic concept of nonlinear beam shaping with linear beam shaping as reference: The applied wave front generates an angular spectrum of wave vectors which are centered around the optical axis. Their distribution results in a target image in the far field. The plane of the SLM is imaged into the nonlinear crystal. As for type II phase matching only crossed polarization states mix efficiently, the shaped light field can only mix with the unshaped one. The addition of the wave vectors results in the tailored light field with only half the initial deflection angle \(\alpha\).
Since the nonlinear conversion is only possible with phase matching, this constraint can be exploited for beam shaping. In type II phase matching, only crossed polarization states can mix and through this assignment, the two fields at the fundamental are clearly distinguishable even when working in a collinear geometry. This concept is sketched in Fig. 1. By keeping one field unshaped while the other polarization carries the full phase information, there is an unambiguous assignment of the frequency-doubled wave front for every wave vector. The experimental implementation is simple as the polarization is set diagonally before it impinges on the SLM and as the SLM is polarization-sensitive, only half of the light field is shaped. After imaging this plane into the nonlinear crystal, only crossed polarization components add up for the second harmonic field. The same concept can be transferred to sum frequency generation where the differing frequencies ensure a non-degenerate mixing process.
Based on those physical principles, we present a setup for nonlinear beam shaping with high quality while the conversion efficiency is almost defined by the crystal’s conversion efficiency without beam shaping. Besides showing experimental results, we model the chosen KTP crystal to estimate the relative conversion efficiency. Finally, we discuss optimizing the conversion efficiency and give an outlook on the applicability of the method for beam shaping in the ultraviolet spectral range.
For second harmonic generation, two fields at the fundamental frequency, \(\varvec{E_1}^{(\omega )}\) and \(\varvec{E_2}^{(\omega )}\), add up to the frequency-doubled field. From the three wave mixing process follows a quadratic relation between the fields \(\varvec{E}^{(2\omega )} \propto \varvec{E_1}^{(\omega )} \cdot \varvec{E_2}^{(\omega )}\). The phase matching condition
requires the fundamental and the second harmonic to maintain their phase relationship during propagation through the nonlinear material to avoid destructive interference and weak conversion. Here \(\varvec{k_{1/2}}^{(\omega )}\) are the wave vectors at the fundamental frequency and \(\varvec{k}^{(2\omega )}\) is the resulting wave vector at the second harmonic. We first present the experimental results before we model the nonlinear crystal to evaluate the relative conversion efficiency.
Setup for nonlinear beam shaping: the polarization-sensitive LC-SLM shapes only half of the light field. The plane of the SLM is imaged into the nonlinear KTP crystal and the unshaped and shaped field components combine during nonlinear conversion and form the target structure.
Experimentally recorded results for nonlinear beam shaping for two different crystal lengths and the initial result at the fundamental as reference. The measured conversion efficiency \(\eta\) is denoted above the individual images.
With polarization as a control mechanism, nonlinear beam shaping can be performed similar to linear beam shaping. The phase mask for wave front shaping can be calculated with the same algorithms and tools as for linear beam shaping and the setup basically consists of a half wave plate, a LC-SLM, a 4-f imaging telescope and the nonlinear crystal as Fig. 2 shows. However, as wave front shaping generates an angular spectrum of wave vectors which are centered around the optical axis, the nonlinear crystal needs to equally support vectors which slightly deviate from the optimum phase matching angle. This required angular tolerance or acceptance is one characteristic of nonlinear crystals and it is defined for an angular deviation where the conversion efficiency drops to half of its maximum value. We chose KTP as the nonlinear crystal for our application as it exhibits a high angular acceptance.
Figure 3 shows experimentally recorded profiles of the second harmonic at \(532\,\text {nm}\) and images of the fundamental at \(1064\,\text {nm}\) (Fig. 3a–d) as reference. The corresponding phase masks are calculated with the Gerchberg Saxton algorithm44. When tailoring the target image in the far field, the propagation can be mathematically described by a Fourier transform. Thus, the intensity profile of the target image effectively reflects the generated angular spectrum. Its size directly shows the required angles of the wave vectors which need to be supported during frequency conversion. Consequently, the larger the targeted image size is, the higher the angular acceptance of the nonlinear KTP crystal needs to be, to equally convert the second harmonic signal for all wave vectors.
As a shorter nonlinear crystal exhibits a higher angular acceptance, second harmonic generation is supported for a broad spectrum of wave vectors and arbitrary target structures can be shaped over a large area. We demonstrate this in Fig. 3e–h for four different target distributions with a \(2\,\text {mm}\) long KTP crystal. All structures are homogeneously converted and the quality is comparable to the results at the fundamental. The read-out at the second harmonic appears slightly degraded with respect to the fundamental. The impression of degradation probably arises from a change in speckle size with respect to the structure size. We calculate the conversion efficiency by dividing the measured power of the shaped target structure at the second harmonic by the power at the fundamental. The conversion efficiency is around \(6{-}8\,\%\) for all structures. Conversely, a \(9\,\text {mm}\) long KTP crystal promises higher conversion efficiency at lower angular acceptance. The results of the two smaller target structures for the \(9\,\text {mm}\) long crystal (Fig. 3i,j) are of the same quality as for the \(2\,\text {mm}\) crystal. The conversion efficiency is \(>30\,\%\) and this is only a little less than the initial conversion efficiency of the nonlinear crystal without beam shaping which is around \(40\,\%\). For the \(2\,\text {mm}\) crystal the values are even closer with \(8.5\,\%\) without and values around \(6{-}8\,\%\) with beam shaping. Those results demonstrate the applicability of nonlinear beam shaping in a regime of high conversion efficiency while maintaining high quality. The homogeneous conversion in the range of the initial conversion efficiency of the nonlinear crystal is due to a plateau in the conversion efficiency for small angular deflections. We further investigate this favorable effect for beam shaping in the next section. Figure 3k,l also shows the limitations of nonlinear beam shaping when working beyond this plateau. The globe and snowflake are almost cut at the borders as the required angles are not supported by phase matching. As parts of the light field are not converted, the conversion efficiency decreases. These results are shown to demonstrate the limitations outside the plateau of high conversion efficiency. It is nonetheless possible to shape a smaller target structure which is magnified with a telescope afterwards.
Normalized conversion efficiency \(\eta\) for a lateral deflection of the ordinary/extraordinary beam and projection along the two axes \(\delta \theta\) and \(\delta \phi\) for a \(2\,\text {mm}\) KTP crystal. Experimental measurements along these axes are in good agreement with the simulated curves. The refractive index ellipsoid for the biaxial KTP crystal (dimensions not to scale) shows the phase matching angles \(\theta _m = 90\,^{\circ }\) and \(\phi _m = 24.8\,^{\circ }\) with the axes for beam deflection along \(\delta \theta\) and \(\delta \phi\). The deflection of the ordinary beam along \(\delta \phi\) is indicated with two incoming fields where the extraordinary beam is polarized along the xy plane of the crystal and the deflected ordinary beam is polarized along the z-axis.
The wave vector is a function of the wavelength and the corresponding refractive index, which depends on the wavelength. Consequently, the phase matching condition in Eq. (1) is generally not fulfilled if fields with different wavelengths are propagating through a dispersive medium. Birefringence is one effect which is used for phase matching: For a biaxial crystal the refractive index is determined by three direction-dependent indices \(n_x\), \(n_y\), and \(n_z\) along the principal axes45. The effective refractive index at any position in space is given by their projection. This projection results in a refractive index ellipsoid as it is shown in Fig. 4 for a KTP crystal. The direction of the polarization determines the refractive index which is perceived by the corresponding light field.
Thus, phase matching can occur under a specific angle of incidence depending on the phase matching type. Type II phase matching is fulfilled, if two crossed polarization states mix where their average value of the refractive indices corresponds to the refractive index of the second harmonic.
KTP is a biaxial crystal which allows for type I and type II phase matching at different angles of incidence \(\theta\) and \(\phi\). As KTP exhibits a higher conversion efficiency for type II phase matching, this option is the typical choice. The phase matching angles \(\theta _m\) and \(\phi _m\) are marked in the refractive index ellipsoid in Fig. 4. The two red curves along the three mutually perpendicular planes of the coordinate system indicate the effective refractive index at the fundamental resulting for a field either polarized within the plane of incidence or perpendicular to that plane. From a perspective of the highest symmetry, the beam oscillating perpendicular to the plane of incidence is called the ordinary beam (dashed lines in Fig. 4) as the refractive index is the same for any angle while a beam with the polarization axis within the plane of incidence is called the extraordinary beam (solid lines in Fig. 4) as the effective refractive index is given by the corresponding ellipse defined by the two refractive indices of the axes spanning that plane. Likewise, the two green curves indicate the refractive index for the frequency-doubled field. As the blue curve marks the mean value of the refractive indices at the fundamental for the two polarization states, its crossing point with the green curve marks the angles for type II phase matching according to \((n_{o/e}^{(\omega )} + n_{e/o}^{(\omega )})\cdot 0.5 = n_{e/o}^{(2\omega )}\).
Wave front shaping introduces a spectrum of wave vectors centered around the optical axis. Even in a collinear geometry, the generated angular spectrum deviates from the ideal phase matching condition. This deviation causes a refractive index mismatch and affects the conversion efficiency. Due to the high angular acceptance of the chosen KTP crystal, a deflection along \(\delta \theta\) and \(\delta \phi\) only causes marginal contributions to the phase mismatch. Our application for beam shaping is different from the general definition of the angular acceptance which accounts for a deflection of the total field of the fundamental: The two mixing fields have a slightly different angle as the shaped field is deflected while the other field propagates without any modifications along the phase matching angle. The two incident beams in Fig. 4 refer to this scenario, where the ordinary beam is deflected along \(\delta \phi\), while the extraordinary beam follows the phase matching angle. In dependence of the deflection angles \(\delta \theta\) and \(\delta \phi\), this results in the phase mismatch \(\Delta \varvec{k}\):
We calculate the refractive index mismatch for a deflection along \(\delta \theta\) and \(\delta \phi\) by calculating the corresponding refractive indices46 for the wave vectors. In good approximation, the resulting wave vector at the second harmonic is generated at half of the initial deflection angle, as this is the mean value between the deflected and undeflected wave vector at the fundamental. The phase mismatch is thus given by the difference of the wave vectors projected in the direction of the frequency doubled field. The relative conversion efficiency can be deduced from the projected mismatch \(\Delta k\) by solving the following equation:
This analytical relation between the second harmonic intensity \(I^{(2\omega )}\) and the intensity at the fundamental \(I^{(\omega )}\) with respect to the crystal length L and the phase mismatch \(\Delta k\) follows when solving the nonlinear wave equation for the assumption of low depletion and the slowly varying amplitude approximation47. On this basis, we calculate the relative 2D conversion efficiency for a deviation within the plane of \(\delta \theta\) and \(\delta \phi\) in Fig. 4. The projection along \(\delta \theta\) and \(\delta \phi\) is shown in the line plots. An arrow indicates the corresponding curve of conversion for the exemplary deflection in front of the refractive index ellipsoid. According to Snell’s law and assuming the small angle approximation, the internal angle corresponds to the external angle connected via the effective refractive index which is perceived by the corresponding wave vector. The graphs in Fig. 4 show the external angle, as this is the relevant value for beam shaping.
Besides calculating the relative conversion efficiency, we perform experiments for an angular deviation introduced by the SLM for beam shaping. To evaluate different angles, we apply a blazed grating either in horizontal or vertical direction on the SLM and image this plane into the nonlinear crystal which has a length of \(2\,\text {mm}\). Similar to beam shaping, only one polarization component is shaped. The corresponding setup is sketched at the bottom of Fig. 4. The coordinate system of the refractive index ellipsoid is rotated with respect to the beam propagation axis as the beam propagates along the matching angle of the KTP crystal.
In the far field we measure the resulting power of the frequency-doubled field and divide it by the measured power of the field at the fundamental. This results in two curves for the relative conversion efficiency, either in the direction of \(\delta \theta\) or \(\delta \phi\), with respect to no deflection at \(0\,^\circ\) at the optimum phase matching angle. As either the ordinary or extraordinary polarization state can be deflected within the \(\delta \theta\) and \(\delta \phi\) plane, two different scenarios are possible and both are shown in Fig. 4. Deflecting the ordinary beam causes a highly symmetric profile of the relative conversion efficiency and is thus the favorable option for our application.
Simulation and experiment are in good agreement and show that the conversion efficiency exhibits a broad plateau. Within that range, the light is homogeneously converted, independently of the exact deflection angle. Moreover, the conversion efficiency is not significantly reduced with respect to the actual conversion efficiency without beam shaping. Thus, arbitrary intensity profiles can be shaped within that angular range without reduction in homogeneity or efficiency. The acceptance angle approximately reaches \(\pm 2\,^{\circ }\) for a \(2\,\text {mm}\) crystal and \(\pm 1\,^{\circ }\) for the \(9\,\text {mm}\) long crystal. It makes sense to compare this value with the maximum deflection angle of the SLM which is determined by \(\sin ^{-1}(\lambda /(2u))\cdot 1/M\), where u is the pixel pitch and M is the magnification between the SLM and the crystal. The two smaller images in Fig. 3 approximately acquire \(19\,\%\) of the SLM’s linear field of view, while the two bigger images acquire \(32\,\%\). This corresponds to deflection angles of approximately \(1.1\,^\circ\) and \(1.9\,^\circ\) (bearing in mind that the plane of the SLM is imaged with a factor of \(M=0.25\) into the nonlinear crystal and thus the maximum deflection angle is \(\pm 6\,^\circ\)).
At the phase matching angle \(\theta _{m}=90\,^\circ\), the crystal is non-critically phase-matched along \(\delta \theta\). This means when expanding the mismatch \(\Delta k\) in a Taylor series, the first derivative \(\partial \Delta k / \partial \theta = 0\) due to the crystal’s symmetry and thus only higher order terms with a weaker effect on the phase mismatch contribute48. Besides the control mechanism for the wave mixing process it is crucial to choose a nonlinear crystal with a high angular acceptance. We consider temperature-controlled noncritical phase matching or quasi-phase matching as another option but do not discuss this further within the scope of this paper.
Nonlinear beam shaping is a powerful tool as it enables dynamic phase-only beam shaping in new spectral ranges. To optimize this process, two parameters should be considered: While the conversion efficiency should be maximized, the angular acceptance needs to be high enough to ensure homogeneous conversion of the angular spectrum. As both parameters are mutually dependent, we discuss the proper choice of the experimental parameters to optimize the outcome. Moreover, we give an outlook on dynamic phase-only beam shaping in the ultraviolet spectral range.
When approaching a high conversion regime, it is more critical to maintain proper phase matching as the phase mismatch gains an increasing impact on the efficiency of the converted outcome. At high conversion, the angular range of tolerance narrows and this decreases the plateau of homogeneous conversion for the generated angular spectrum. Equation (3) shows the impact of the experimental parameters on the conversion efficiency while keeping the laser power constant. Both increasing the intensity and the crystal length promises higher conversion efficiency. However, both parameters also affect the angular acceptance of the nonlinear crystal.
The effect of the crystal length L is directly stated in Eq. (3). While a longer crystal increases the relative conversion with \(L^2\), it also increases the contribution of the phase mismatch as a multiplier of \(\Delta k\) in the argument of the \({{\,\textrm{sinc}\,}}\) function. We use the model of the nonlinear crystal to calculate the relative conversion of the target structures for the \(9\,\text {mm}\) crystal as they are shown in Fig. 3.
Simulated beam conversion for the \(9\,\text {mm}\) crystal. The simulated efficiency is calculated from the integrated relative conversion multiplied with the efficiency of the \(9\,\text {mm}\) crystal without beam shaping (compare Fig. 6).
The beam-shaped structures directly reflect the generated angular spectrum and thus the relative conversion can be simulated when multiplying the target image with the calculated angular conversion. Figure 5 shows the simulated results with the corresponding conversion efficiency. This value is calculated as the total measured efficiency without beam shaping multiplied by the integrated relative conversion. Working beyond the limits of phase matching results in weak conversion and the simulated as well as the measured values for the conversion efficiency drop significantly. The experimental parameters should be chosen to work within the pleateau of homogeneous conversion as otherwise not only the quality gets worse but also the efficiency drops significantly. The comparison of the simulated and experimentally measured conversion efficiency explains the strong decrease in efficiency for the two larger structures with their measured conversion efficiency marked as gray crosses in Fig. 6. This decrease can be traced back to insufficient phase matching as the structures are chosen larger than the angular acceptance. As simulation and experiment are in good agreement, proper experimental parameters can be derived on that basis even in a high conversion regime.
Measured output power vs. input power for second harmonic generation with a \(2\,\text {mm}\) and \(9\,\text {mm}\) KTP crystal without beam shaping. The crosses indicate the conversion efficiency of the structures in Fig. 3 for beam shaping, with the two intentionally cropped structures in gray.
Similar to the crystal length L, an increase in intensity also improves the conversion efficiency. It can either be increased by increasing the power of the chosen light source or by reducing the illuminated area. As the laser power is typically technically limited, we will focus on the second case. Here the plane of the SLM is imaged with a certain demagnification into the nonlinear crystal to decrease the illuminated area. It is also possible to illuminate the SLM with a smaller beam but this reduces the number of illuminated pixels and unnecessarily increases the intensity on the SLM with respect to potential damage thresholds. It is thus beneficial to design the telescope with the required demagnification. This demagnification increases the intensity in the nonlinear crystal and scales the intensity with \(1/M^2\) to gain higher conversion efficiencies. Likewise, the telescope increases the angular spectrum as the deflection angles are proportional to M. Similar to the crystal length L, this factor also affects the phase mismatch \(\Delta k\) in the argument of the \({{\,\textrm{sinc}\,}}\) function as it acts as a multiplier of the initially induced angular spectrum. Consequently, the intensity needs to be considered when modeling the phase mismatch to determine optimum parameters.
In our experimental setup we work with a telescope that decreases the image size by a factor of 0.25 and thus increases the intensity by a factor of 16. This reduction affects the deflection angles and increases them by a factor of 4 as the effective pixel size of the SLM changes. We consider that adjusting the intensity with respect to the required angular spectrum for a given crystal length is reasonable. The resulting plateau should be exactly within the required angular range—neither broader nor narrower—to ensure homogeneous conversion with maximized efficiency.
Dynamic beam shaping in the ultraviolet (UV) spectral range is highly limited as many devices absorb the UV light, including liquid crystal displays. This paper demonstrates nonlinear beam shaping from the infrared to visible spectral range. Nonetheless, other conversion processes are possible, reaching from second harmonic generation at different fundamental frequencies up to other nonlinear processes like sum frequency generation. Currently, there is research on new materials with tunable birefringence in the UV spectral range for spatial light modulators with results showing light modulation at \(303\,\text {nm}\)49. We see the potential of our method for approaching even deeper UV spectral ranges.
This section presents a few thoughts on proper nonlinear crystals and conversion processes for dynamic phase-only beam shaping in the UV range.
Ultrashort pulsed laser systems often feature nonlinear crystals for frequency conversion, for example to the second or third harmonic. High-energetic short wavelengths in the UV range of ultrashort pulsed laser systems are thus often generated with nonlinear frequency conversion from the visible or infrared range. Consequently, the initial conditions often enable a direct integration of nonlinear beam shaping into the laser setup and ideally this can be performed with conversion efficiencies close to the initial values.
To achieve high-quality results, the nonlinear crystal should have a broad angular acceptance and the conversion process needs to be non-degenerate. Sum frequency generation with an LBO crystal from \(1064\,\text {nm}\) and \(532\,\text {nm}\) to \(355\,\text {nm}\) is one option. The major benefit of this crystal is the high angular acceptance. Similarly, sum frequency generation with CLBO50 from \(1064\,\text {nm}\) and \(266\,\text {nm}\) to \(213\,\text {nm}\) should be possible. In both options the process is inherently non-degenerate as two different wavelengths mix. While the SLM can either shape the light field in the infrared or visible spectral range, only the infrared light can be shaped in the latter case as the SLM can not shape light at \(266\,\text {nm}\). With one of the options suggested here or other nonlinear processes, we see the potential of using nonlinear beam shaping to approach the UV spectral range. This enables dynamic phase-only beam shaping beyond the physical limits of the liquid crystal display.
Experimental setup for nonlinear beam shaping: A half wave plate introduces diagonal polarization on the beam impinging on the SLM. Thus, only half of the field is shaped while the crossed polarization component remains unshaped. A telescope images the plane of the SLM into the nonlinear KTP crystal. Here, the two polarization components mix and shape the targeted field. After a filter separates the second harmonic from the fundamental, a lens images the target image, where the second harmonic signal is recorded. As the shaped result is sensitive to the deflected polarization component, a half wave plate after the SLM allows for changing the deflected and non-deflected polarization component.
Nonlinear phase-only beam shaping with a liquid crystal display is possible as the applied phase information is conserved during the frequency conversion. However, as second harmonic generation is a degenerate process, the mixing waves need to be controlled by making the interacting fields distinguishable. We use the requirements of type II phase matching to unambiguously define the outcome of the second harmonic. Under this condition, only fields of crossed polarization can mix efficiently. A half wave plate, abbreviated with \(\lambda /2\) in Fig. 7 is placed in front of the SLM to generate diagonal polarization with respect to the SLM’s orientation. As the liquid crystals are polarization-sensitive, only half of the field is shaped, while the other part remains unshaped. Consequently, a shaped component can only mix with the unshaped plane wave front. In Fig. 7 the two polarization components are indicated as the ordinary and the extraordinary beam. The polarization configuration in the figure refers to where the ordinary beam is shaped. If the extraordinary beam is shaped, an additional half wave plate after the SLM helps to rotate the polarization states correspondingly. This might be beneficial when working off-axis to separate non-diffracted light as the area of quasi-homogeneous conversion is shifted away from the center of the optical axis.
A telescope images the plane of the SLM into the nonlinear crystal. This ensures a homogeneous intensity distribution within the nonlinear crystal which is required to convert all field components equally within a range of proper phase matching. In good approximation, the resulting wave vectors at the second harmonic reduce to half of the initially set angle of the shaped field component, as this is the sum of the two involved wave vectors.
The SLM is the model LSH0701010 from Hamamatsu (Hamamatsu, Japan) and has a pixel resolution of \(800\,\text {px}\times 600\,\text {px}\) with a pixel pitch of 20 μm. The beam diameter on the LC-SLM is \(14\,\text {mm}\). To increase the intensity to obtain a higher conversion efficiency, we design our telescope to demagnify the beam by \(M=0.25\). The corresponding focal lengths are \(f_1 = 200\,\text {mm}\) and \(f_2 = 50\,\text {mm}\). Before the target distribution is imaged to the far field with a \(f_3 = 200\,\text {mm}\) lens, the remaining IR light is blocked with the filter FL532-10 from Thorlabs (NJ, USA). The laser system Fuego from Time-Bandwidth Products (CA, USA) emits \(10\,\text {ps}\) pulses at a repetition frequency of \(200\,\text {kHz}\) at a wavelength of \(1064\,\text {nm}\). The power measurements were done with the power meter PM10 from Coherent (CA, USA) and PM160 from Thorlabs. We record RGB images with the camera model UI-3000SE-C-HQ from IDS (Obersulm, Germany).
Nonlinear computer-generated holography is gaining increasing attention and there is a lot of research on structuring nonlinear crystals but this process is technically challenging and static. These deficiencies can be overcome by thin holograms generated with dynamic liquid crystal displays. However, research on nonlinear beam shaping is still in its infancy. We present a high-quality and highly efficient method for nonlinear dynamic phase-only beam shaping with a simple concept: The mixing process needs to be non-degenerate by employing a control mechanism and the angular acceptance of the nonlinear crystal needs to support the angular spectrum generated by wave front shaping. Based on this, we demonstrate beam shaping at the second harmonic with conversion efficiencies close to the ones without beam shaping. This is mainly due to the fact that we work at a plateau of constant conversion efficiency for the generated angular spectrum. While working within that range, the conversion efficiency is mainly given by the conversion efficiency without beam shaping and thus this method is highly efficient. Furthermore, the quality of the frequency-converted result is of the same quality as for beam shaping at the fundamental. We see high potential in this approach as it not only enables highly-efficient and high-quality beam shaping with thin holograms but due to its simplicity, it can be easily combined with other elaborated processes. Based on nonlinear beam shaping, the spectral range of dynamic liquid crystal displays can be extended beyond the physical limits as we discuss here for the UV spectral range.
Data underlying the results presented in this paper may be obtained from the corresponding author upon reasonable request.
Mainman, T. H. Stimulated optical radiation in ruby. Nature 187, 493–494. https://doi.org/10.1038/187493a0 (1960).
Article ADS Google Scholar
Maurer, C., Jesacher, A., Bernet, S. & Ritsch-Marte, M. What spatial light modulators can do for optical microscopy. Laser Photon. Rev. 5, 81–101. https://doi.org/10.1002/lpor.200900047 (2010).
Article ADS Google Scholar
Rubinsztein-Dunlop, H. et al. Roadmap on structured light. J. Opt. 19, 013001 (2016).
Article ADS CAS Google Scholar
Accanto, N. et al. Multiplexed temporally focused light shaping for high-resolution multi-cell targeting. Optica 5, 1478–1491. https://doi.org/10.1364/OPTICA.5.001478 (2018).
Article ADS CAS Google Scholar
Sugioka, K. & Cheng, Y. Ultrafast lasers–reliable tools for advanced materials processing. Light Sci. Appl. 3, e149. https://doi.org/10.1038/lsa.2014.30 (2014).
Article ADS CAS Google Scholar
Sugioka, K. Progress in ultrafast laser processing and future prospects. Nanophotonics 6, 393–413. https://doi.org/10.1515/nanoph-2016-0004 (2017).
Article CAS Google Scholar
Orazi, L., Romoli, L., Schmidt, M. & Li, L. Ultrafast laser manufacturing: From physics to industrial applications. CIRP Ann. 70, 543–566. https://doi.org/10.1016/j.cirp.2021.05.007 (2021).
Article Google Scholar
Yariv, A. Four wave nonlinear optical mixing as real time holography. Opt. Commun. 25, 23–25. https://doi.org/10.1016/0030-4018(78)90079-2 (1978).
Article ADS Google Scholar
Abraham, N. B. & Firth, W. J. Overview of transverse effects in nonlinear-optical systems. J. Opt. Soc. Am. B 7, 951. https://doi.org/10.1364/josab.7.000951 (1990).
Article ADS CAS Google Scholar
Basistiy, I., Bazhenov, V., Soskin, M. & Vasnetsov, M. Optics of light beams with screw dislocations. Opt. Commun. 103, 422–428. https://doi.org/10.1016/0030-4018(93)90168-5 (1993).
Article ADS Google Scholar
Dholakia, K., Simpson, N. B., Padgett, M. J. & Allen, L. Second-harmonic generation and the orbital angular momentum of light. Phys. Rev. A 54, R3742–R3745. https://doi.org/10.1103/physreva.54.r3742 (1996).
Article ADS CAS PubMed Google Scholar
Courtial, J., Dholakia, K., Allen, L. & Padgett, M. J. Second-harmonic generation and the conservation of orbital angular momentum with high-order Laguerre-Gaussian modes. Phys. Rev. A 56, 4193–4196. https://doi.org/10.1103/physreva.56.4193 (1997).
Article ADS CAS Google Scholar
Mair, A., Vaziri, A., Weihs, G. & Zeilinger, A. Entanglement of the orbital angular momentum states of photons. Nature 412, 313–316. https://doi.org/10.1038/35085529 (2001).
Article ADS CAS PubMed Google Scholar
Shao, G.-H., Wu, Z.-J., Chen, J.-H., Xu, F. & Lu, Y.-Q. Nonlinear frequency conversion of fields with orbital angular momentum using quasi-phase-matching. Phys. Rev. A 88, 063827. https://doi.org/10.1103/physreva.88.063827 (2013).
Article ADS Google Scholar
Zhou, Z.-Y. et al. Generation of light with controllable spatial patterns via the sum frequency in quasi-phase matching crystals. Sci. Rep. 4, 1–5. https://doi.org/10.1038/srep05650 (2014).
Article CAS Google Scholar
Libster-Hershko, A., Trajtenberg-Mills, S. & Arie, A. Dynamic control of light beams in second harmonic generation. Opt. Lett. 40, 1944. https://doi.org/10.1364/ol.40.001944 (2015).
Article ADS PubMed Google Scholar
Li, Y., Zhou, Z.-Y., Ding, D.-S. & Shi, B.-S. Sum frequency generation with two orbital angular momentum carrying laser beams. J. Opt. Soc. Am. B 32, 407. https://doi.org/10.1364/josab.32.000407 (2015).
Article ADS CAS Google Scholar
Steinlechner, F., Hermosa, N., Pruneri, V. & Torres, J. P. Frequency conversion of structured light. Sci. Rep. 6, 1–8. https://doi.org/10.1038/srep21390 (2016).
Article CAS Google Scholar
Buono, W. T. et al. Polarization-controlled orbital angular momentum switching in nonlinear wave mixing. Opt. Lett. 43, 1439. https://doi.org/10.1364/ol.43.001439 (2018).
Article ADS CAS PubMed Google Scholar
Qiu, X., Li, F., Liu, H., Chen, X. & Chen, L. Optical vortex copier and regenerator in the Fourier domain. Photon. Res. 6, 641. https://doi.org/10.1364/prj.6.000641 (2018).
Article CAS Google Scholar
Shen, Y. et al. Optical vortices 30 years on: OAM manipulation from topological charge to multiple singularities. Light Sci. Appl. 8, 90. https://doi.org/10.1038/s41377-019-0194-2 (2019).
Article ADS CAS PubMed PubMed Central Google Scholar
Li, J. et al. Directly generating vortex beams in the second harmonic by a spirally structured fundamental wave. Chin. Opt. Lett. 19, 060005 (2021).
Article ADS Google Scholar
da Silva, B. P. et al. Spin to orbital angular momentum transfer in frequency up-conversion. Nanophotonics 11, 771–778. https://doi.org/10.1515/nanoph-2021-0493 (2021).
Article CAS Google Scholar
Hancock, S. W., Zahedpour, S. & Milchberg, H. M. Second-harmonic generation of spatiotemporal optical vortices and conservation of orbital angular momentum. Optica 8, 594. https://doi.org/10.1364/optica.422743 (2021).
Article ADS CAS Google Scholar
Wu, H.-J. et al. Radial modal transitions of Laguerre-Gauss modes during parametric up-conversion: Towards the full-field selection rule of spatial modes. Phys. Rev. A 101, 063805. https://doi.org/10.1103/physreva.101.063805 (2020).
Article ADS CAS Google Scholar
Fang, X. et al. High-dimensional orbital angular momentum multiplexing nonlinear holography. Adv. Photon. 3, 015001. https://doi.org/10.1117/1.ap.3.1.015001 (2021).
Article ADS Google Scholar
Buono, W. T. & Forbes, A. Nonlinear optics with structured light. Opto-Electron. Adv. 5, 210174–1 (2022).
Article CAS Google Scholar
Berger, V. Nonlinear photonic crystals. Phys. Rev. Lett. 81, 4136–4139. https://doi.org/10.1103/physrevlett.81.4136 (1998).
Article ADS CAS Google Scholar
Hong, X.-H., Yang, B., Zhang, C., Qin, Y.-Q. & Zhu, Y.-Y. Nonlinear volume holography for wave-front engineering. Phys. Rev. Lett. 113, 163902. https://doi.org/10.1103/physrevlett.113.163902 (2014).
Article ADS PubMed Google Scholar
Shapira, A., Naor, L. & Arie, A. Nonlinear optical holograms for spatial and spectral shaping of light waves. Sci. Bull. 60, 1403–1415. https://doi.org/10.1007/s11434-015-0855-3 (2015).
Article Google Scholar
Xu, T. et al. Three-dimensional nonlinear photonic crystal in ferroelectric barium calcium titanate. Nat. Photon. 12, 591–595. https://doi.org/10.1038/s41566-018-0225-1 (2018).
Article ADS CAS Google Scholar
Wei, D. et al. Experimental demonstration of a three-dimensional lithium niobate nonlinear photonic crystal. Nat. Photon. 12, 596–600. https://doi.org/10.1038/s41566-018-0240-2 (2018).
Article ADS CAS Google Scholar
Liu, S. et al. Nonlinear wavefront shaping with optically induced three-dimensional nonlinear photonic crystals. Nat. Commun. 10, 3208. https://doi.org/10.1038/s41467-019-11114-y (2019).
Article ADS CAS PubMed PubMed Central Google Scholar
Liu, S., Mazur, L. M., Krolikowski, W. & Sheng, Y. Nonlinear volume holography in 3d nonlinear photonic crystals. Laser Photon. Rev. 14, 2000224. https://doi.org/10.1002/lpor.202000224 (2020).
Article ADS CAS Google Scholar
Zhu, B. et al. Second-harmonic computer-generated holographic imaging through monolithic lithium niobate crystal by femtosecond laser micromachining. Opt. Lett. 45, 4132. https://doi.org/10.1364/ol.394162 (2020).
Article ADS CAS PubMed Google Scholar
Trajtenberg-Mills, S. & Arie, A. Shaping light beams in nonlinear processes using structured light and patterned crystals. Opt. Mater. Express 7, 2928. https://doi.org/10.1364/ome.7.002928 (2017).
Article ADS Google Scholar
Liu, H. & Chen, X. The manipulation of second-order nonlinear harmonic wave by structured material and structured light. J. Nonlinear Opt. Phys. Mater. 27, 1850047. https://doi.org/10.1142/s0218863518500479 (2018).
Article ADS CAS Google Scholar
Keren-Zur, S., Avayu, O., Michaeli, L. & Ellenbogen, T. Nonlinear beam shaping with plasmonic metasurfaces. ACS Photon. 3, 117–123. https://doi.org/10.1021/acsphotonics.5b00528 (2015).
Article CAS Google Scholar
Almeida, E., Bitton, O. & Prior, Y. Nonlinear metamaterials for holography. Nat. Commun. 7, 12533. https://doi.org/10.1038/ncomms12533 (2016).
Article ADS CAS PubMed PubMed Central Google Scholar
Liu, H. et al. Dynamic computer-generated nonlinear-optical holograms. Phys. Rev. A 96, 023801. https://doi.org/10.1103/physreva.96.023801 (2017).
Article ADS CAS Google Scholar
Liu, H., Zhao, X., Li, H., Zheng, Y. & Chen, X. Dynamic computer-generated nonlinear optical holograms in a non-collinear second-harmonic generation process. Opt. Lett. 43, 3236. https://doi.org/10.1364/ol.43.003236 (2018).
Article ADS CAS PubMed Google Scholar
Wu, Y., Liu, H. & Chen, X. Three-dimensional nonlinear optical holograms. Phys. Rev. A 102, 063505. https://doi.org/10.1103/physreva.102.063505 (2020).
Article ADS CAS Google Scholar
Yao, W. et al. Angle-multiplexing nonlinear holography for controllable generations of second-harmonic structured light beams. Front. Phys. 9, 751860. https://doi.org/10.3389/fphy.2021.751860 (2021).
Article Google Scholar
Gerchberg, R. W. & Saxton, W. O. A practical algorithm for the determination of phase from image and diffraction plane pictures. Optik 35, 237–246 (1972).
Google Scholar
Zhang, W., Yu, H., Wu, H. & Halasyamani, P. S. Phase-matching in nonlinear optical compounds: A materials perspective. Chem. Mater. 29, 2655–2668. https://doi.org/10.1021/acs.chemmater.7b00243 (2017).
Article CAS Google Scholar
Yao, J. Q. & Fahlen, T. S. Calculations of optimum phase match parameters for the biaxial crystal KTiOPO_4. J. Appl. Phys. 55, 65–68. https://doi.org/10.1063/1.332850 (1984).
Article ADS CAS Google Scholar
Boyd, R. W. Nonlinear Optics (Academic Press, 2020).
Google Scholar
Fan, T. Y. et al. Second harmonic generation and accurate index of refraction measurements in flux-grown KTiOPO_4. Appl. Opt. 26, 2390. https://doi.org/10.1364/ao.26.002390 (1987).
Article ADS CAS PubMed Google Scholar
Xu, Y. et al. Deep ultraviolet hydrogel based on 2D cobalt-doped titanate. Light Sci. Appl. 12, 1. https://doi.org/10.1038/s41377-022-00991-6 (2023).
Article ADS CAS PubMed PubMed Central Google Scholar
Mori, Y., Kuroda, I., Nakajima, S., Sasaki, T. & Nakai, S. New nonlinear optical crystal: Cesium lithium borate. Appl. Phys. Lett. 67, 1818–1820. https://doi.org/10.1063/1.115413 (1995).
Article ADS CAS Google Scholar
Download references
The authors gratefully acknowledge funding from the Erlangen Graduate School in Advanced Optical Technologies (SAOT) by the German Research Foundation (DFG) in the framework of the German excellence initiative and the authors acknowledge financial support by Deutsche Forschungsgemeinschaft and Friedrich-Alexander-Universität Erlangen-Nürnberg within the funding programme “Open Access Publication Funding”.
Open Access funding enabled and organized by Projekt DEAL.
Institute of Photonic Technologies, Friedrich-Alexander-Universität Erlangen-Nürnberg, Konrad-Zuse-Straße 3/5, 91052, Erlangen, Germany
Lisa Ackermann, Clemens Roider, Kristian Cvecek, Nicolas Barré, Christian Aigner & Michael Schmidt
School of Advanced Optical Technologies (SAOT), Friedrich-Alexander-Universität Erlangen-Nürnberg, Paul-Gordan-Straße 6, 91052, Erlangen, Germany
Lisa Ackermann, Kristian Cvecek, Nicolas Barré & Michael Schmidt
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
You can also search for this author in PubMed Google Scholar
L.A. and C.R. conceived the method. L.A. recorded the experimental data and performed simulations. L.A., C.R. and K.C. analyzed the data. C.A. performed provisional experiments. L.A., C.R., K.C., N.B. and M.S. contributed to the manuscript and MS acquired funding.
Correspondence to Lisa Ackermann.
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Reprints and Permissions
Ackermann, L., Roider, C., Cvecek, K. et al. Polarization-controlled nonlinear computer-generated holography. Sci Rep 13, 10338 (2023). https://doi.org/10.1038/s41598-023-37443-z
Download citation
Received: 05 April 2023
Accepted: 21 June 2023
Published: 26 June 2023
DOI: https://doi.org/10.1038/s41598-023-37443-z
Anyone you share the following link with will be able to read this content:
Sorry, a shareable link is not currently available for this article.
Provided by the Springer Nature SharedIt content-sharing initiative
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.