**1920×1080**) is set up to validate the proposed method. The experimental results indicate that the proposed method can achieve high-quality holographic femtosecond laser parallel processing with a significantly improved energy utilization rate.

## 1.

## Introduction

Femtosecond laser processing is a powerful technique which enables microfabrication in transparent materials. The advantages of femtosecond laser processing include high-spatial resolution, reduced thermal destruction, and widespread application for various materials.^{1} However, it is a drawback that such a device fabrication requires a huge number of processing points. Therefore, the throughput of femtosecond laser processing must be improved.

In recent years, with the development of spatial light modulator (SLM) technology, achieving parallel processing design by loading a computer-generated hologram (CGH) on the SLM can not only save the cost of a diffractive optical element in traditional parallel processing, but also flexibly generate the arbitrary patterns. Therefore, femtosecond laser holographic parallel processing based on SLM has attracted broad attention in research.^{2}3.4.5.^{–}^{6} However, dynamic SLM holograms based on pixel structure still face many problems for achieve the target diffraction pattern, as the zero-order diffracted light and multiorder image reproduction caused by the characteristics of the SLM pixel structure have a large impact on the target image quality. A great number of studies aiming at addressing these issues was carried out by researchers in the holographic display fields all around the world.^{7}8.9.10.^{–}^{11} For example, Palima and Daria^{7} designed a CGH which can produce a desired phase and a corrective phase, and the beam with the corrective phase can destructively interfere with the zero-order beam. Christmas et al.^{8} separated the reconstructed pattern from the zero-order light focus position and eliminated the zero-order light interference by filtering. Silvennoinen et al.^{9}^{,}^{10} used a 500-mm Fresnel lens and beam blocker to remove the unwanted zero-order diffraction from the image plane. Zhang et al.^{11} proposed a new technique to eliminate the zero-order light by adding both a linear phase and a divergent spherical phase to the precalculated phase distribution on the phase-only SLM.

In this paper, a novel method for holographic femtosecond laser parallel processing is proposed, which can not only suppress the interference of zero-order light effectively, but also improve the energy utilization rate. The paper is organized as follows. In Sec. 2, the reasons for the existence of zero-order light of SLM with a pixelated structure are demonstrated. In Sec. 3, an optimization method of hologram projection is illustrated. Finally, experimental results using the proposed method and discussions are given in Secs. 4 and 5, respectively.

## 2.

## Optical Characteristic of the Pixelated Phase-Only Spatial Light Modulator

Holographic femtosecond laser parallel processing in this paper is realized based on the pixelated structure of the SLM with a phase-only hologram. In this section, optical field characteristics of the phase-only SLM are investigated. Based on the study, the cause of zero-order light which disturbs the target pattern is explained.

Figure 1 shows the geometric schematic of an SLM pixelated structure (Holoeye, pluto NIR-2). The single pixel, whose pixel size is $a$, is square and the pixel space is $d$. ${L}_{x}$ and ${L}_{y}$ represent the width and length of the SLM, respectively. The origin of the $XOY$ coordinate system is in the center of the phase-only SLM.

The transmittance function of the SLM pixelated structure can be described as

## (1)

$$t(x,y)=\mathrm{rect}(\frac{x}{{L}_{x}},\frac{y}{{L}_{y}})[{t}_{\mathrm{ac}}(x,y)+{t}_{\mathrm{ds}}(x,y)],$$The transmittance function of coding phase hologram areas (active areas) in SLM is expressed as

## (2)

$${t}_{\mathrm{ac}}(x,y)=\mathrm{rect}(\frac{x}{a},\frac{y}{a})\otimes \left\{\frac{1}{{d}^{2}}\mathrm{comb}\right(\frac{x}{d},\frac{y}{d})\mathrm{exp}[i{\varnothing}_{\mathrm{ac}}(x,y)]\}.$$${\varnothing}_{\mathrm{ac}}(x,y)$ is the hologram coding phase information computed by the Gerchberg–Saxton algorithm,^{12} which is loaded on the active areas; $\otimes $ is the convolution operation; and $\mathrm{rect}(x/{L}_{x},y/{L}_{y})$ is the aperture function of the SLM.

The transmittance function of noncoding areas (dead space areas) in SLM is expressed as

## (3)

$${t}_{\mathrm{ds}}(x,y)=\left\{\right[\mathrm{rect}(\frac{x}{d},\frac{y}{d})-\mathrm{rect}(\frac{x}{a},\frac{y}{a})]\phantom{\rule{0ex}{0ex}}\otimes \frac{1}{{d}^{2}}\mathrm{comb}(\frac{x}{d},\frac{y}{d}\left)\right\}{A}_{\mathrm{ds}}(x,y)\mathrm{exp}[i{\varnothing}_{\mathrm{ds}}(x,y)]\},\phantom{\rule{0ex}{0ex}}$$The complex amplitude distribution $T({f}_{x},{f}_{y})$ in the reconstruction plane is the Fourier transform of $t(x,y)$, i.e., $T({f}_{x},{f}_{y})=\mathfrak{I}[t(x,y)]$, where $\mathfrak{I}$ denotes the Fourier transform.

Further, $T({f}_{x},{f}_{y})$ can be expressed as

## (4)

$$T({f}_{x},{f}_{y})={L}_{x}{L}_{y}\mathrm{sinc}({f}_{x}{L}_{x},{f}_{y}{L}_{y})\phantom{\rule{0ex}{0ex}}\otimes [{T}_{\mathrm{ac}}({f}_{x},{f}_{y})+{T}_{\mathrm{ds}}({f}_{x},{f}_{y})].$$The contribution to the complex amplitude from the active areas is

## (5)

$${T}_{\mathrm{ac}}({f}_{x},{f}_{y})={a}^{2}\mathrm{sinc}({f}_{x}a,{f}_{y}a)[\mathrm{comb}({f}_{x}d,{f}_{y}d)\phantom{\rule{0ex}{0ex}}\otimes E({f}_{x},{f}_{y})],$$The above equations indicate that the shape of pixels forms the sinc envelope and the comb function describes the replication of the reconstruction pattern $E({f}_{x},{f}_{y})$ in the output plane.

The contribution to the complex amplitude from the dead space areas is

## (7)

$${T}_{\mathrm{ds}}({f}_{x},{f}_{y})=\{[{d}^{2}\mathrm{sinc}({f}_{x}d,{f}_{y}d)\phantom{\rule{0ex}{0ex}}-{a}^{2}\mathrm{sinc}({f}_{x}a,{f}_{y}a)]\mathrm{comb}({f}_{x}d,{f}_{y}d)\}\otimes F({f}_{x},{f}_{y}),$$## (8)

$$F({f}_{x},{f}_{y})=\mathfrak{I}\{{A}_{\mathrm{ds}}(x,y)\mathrm{exp}[i{\varphi}_{\mathrm{ds}}(x,y)]\}.$$As $\text{\hspace{0.17em}}\text{\hspace{0.17em}}{L}_{x}$ and ${L}_{y}$ are large enough, the scanning function $\text{\hspace{0.17em}}\text{\hspace{0.17em}}{L}_{x}{L}_{y}\mathrm{sinc}({f}_{x}{L}_{x},{f}_{y}{L}_{y})$ here approximately equals the $\delta $ function.

According to Eq. (8), the light distribution in the center of the reconstruction plane (${f}_{x}=0,{f}_{y}=0$) can be described as

## (9)

$$T(\mathrm{0,0})={\left(\frac{a}{d}\right)}^{2}E(\mathrm{0,0})+[1-{\left(\frac{a}{d}\right)}^{2}]F(\mathrm{0,0})=\mu E(\mathrm{0,0})+(1-\mu )F(\mathrm{0,0}),$$## 3.

## Optimization Method of Hologram

## 3.1.

### Loading Digital Blazed Grating

Multiorder image reproduction exists during the processes of applying an SLM pixelated structure to achieve reproduction of the target pattern. It is located between the adjacent maximum positions of interferences with a central position of $\pm k/2d$, $k=\mathrm{1,3},5\dots $. The image reproduction is modulated by a single-slit diffraction pattern. As for the two-dimensional (2-D) case, first-order image reproduction can be obtained around the zero-order with higher energy compared with the four other diffractions. However, the efficiency of optical energy utilization is low. As in holographic parallel processing, it is only desired to get a target pattern with high energy. To realize this target, combining the former work by Tan et al.^{13} and Yu et al.,^{14} in this paper, we present a method to blaze the reproduction image to the maximum point of zero-order interface by using a digital blazed grating. In theory, the energy of the target pattern can be improved by 279.3% using the proposed method (Holoeye, pluto NIR-2).

The 2-D digital blazed grating can be characterized as

where mod is the remainder operation; $-{L}_{x}/2d\le i\le \phantom{\rule{0ex}{0ex}}{L}_{x}/2d-1$, $-{L}_{y}/2d\le j\le {L}_{y}/2d-1$; $T$ is the period of the digital blazed grating, and a digital blazed grating with different blazing angles can be obtained by varying the values of $T$; digital blazed gratings with different directions can be obtained by varying the values of $b$ and $c$. As the phase modulation range of SLM is $0\sim 2\pi $, the phase distribution loaded on SLM is## 3.2.

### Loading Divergent Spherical Phase

After the target pattern is blazed to the maximum point of the zero-order interface, it is necessary to eliminate the disturbance of the zero-order light beam spot caused by the SLM pixel structure in lens L Fourier plane before leading to the microscope processing system. Otherwise, the quality of the processing pattern will suffer a severe degradation. In order to separate the focal plane of the target pattern from the multiorder diffraction beam plane, phase information of the spherical wave ${\varphi}_{s}=-k/2r({x}^{2}+{y}^{2})$ is added to the precalculated phase distribution of the hologram plane, which is equivalent to adding a negative lens to the target pattern imaging procedure. Here, $k=2\pi /\lambda $ is the wave number in free space; $\text{\hspace{0.17em}}\text{\hspace{0.17em}}\lambda $ is the wavelength; and $r$ is the distance from the SLM to the center of the divergent spherical wave. By doing so, a reconstruction pattern can be shifted from the focal plane of the Fourier lens along the optical axis, while the zero-order illumination will remain unchanged. As shown in Fig. 2, $f$ is the focal length of the Fourier lens.

The total phase loaded on the SLM can be expressed as

The distance between the focal plane of the Fourier lens and the reconstruction plane can be expressed as

When $r>f$, a real image of the reconstruction pattern can be obtained. The zero-order light can be blocked by adding a high-pass filter in the focal plane of the Fourier lens. In consequence, a high-quality target pattern without interference can be exported and led into the microscope processing system to achieve high-quality holographic femtosecond laser parallel processing.

## 4.

## Experimental Setup

The experimental setup is shown in Fig. 3. The laser source is a mode-locked Ti:sapphire ultrafast oscillator (Coherent, Chamleon Vision-S) with a central wavelength at 800 nm, pulse duration of 75 fs, and repetition rate at 80 MHz.^{15} The power of the femtosecond laser (about 3 W) is modulated with a half-wave plate and a Glan Laser beam splitter. After passing through a shutter and beam expander, the laser beam illuminates a liquid crystal SLM to ensure the modulation effect. The SLM (Holoeye, pluto NIR-2) has $1920\times 1080$ pixels, each with a pitch of $8\text{\hspace{0.17em}}\text{\hspace{0.17em}}\mu \mathrm{m}$. In this work, the central $1080\times 1080$ pixels are used. A high-pass filter and aperture are placed at the focal plane of Lens 1 to block the zero-order light. The remainder of the modulated beam is collected by Lens 2 and focused through a microscope object (100, 0.9 NA) into the sample plane. In this work, the photoresist is SZ-2080 (provided by IESL-FORTH, Greece).

## 5.

## Result and Discussion

To validate the feasibility of the proposed method, four experiments are conducted for comparison. A Chinese character shown in Fig. 4 is selected as the target pattern to be processed. The dimension of the total pattern is $1080\times 1080$ pixels, and the dimension of the target pattern is $50\times 50$ pixels. The phase holograms in the experiments are computed by the GS algorithm. In order to eliminate random error, the hologram is loaded using the time-domain averaging method,^{16} which loads a hologram every 150 ms for a total number of 10. The target pattern is fabricated with a single exposure of 30 s and outgoing power of the SLM is about 290 MW.

In the first experiment, we first load a group of holograms on SLM without a digital blazed grating and divergent spherical wave, and the processing result is shown in Fig. 5(a). It is clear that a strong zero-order light spot exists in the center of the pattern. Processing of the pattern exhibits a gradient, in which the part close to the zero-order light spot is entirely exposed while the part away from the zero-order light spot is partly exposed. The phenomenon is in accordance with the physical characteristics.

In order to improve the energy utilization rate, in the second experiment we load a group of holograms on SLM with a digital blazed grating, and the processing result is illustrated in Fig. 5(b). Although the entire Chinese character pattern has been exposed, there is a zero-order light spot in the center of the character. The reason is that the energy of the zero-order light is too high.

In the third experiment, we also load a group of holograms on SLM with a digital blazed grating, but this time a high-pass filter located in the focal plane of Fourier lens is utilized to block the zero-order light. As shown in Fig. 5(c), the zero-order light is successfully suppressed. However, this method reduces the energy efficiency of the central beam and leads to the failure of the entire explosion in the center of the character pattern.

In the last experiment, a group of holograms on SLM with a digital blazed grating and divergent spherical wave are utilized. The results indicate that this method can not only improve the energy utilization rate, but also, by using the high-pass filter and aperture, separate the focus focal plane of the target pattern and the plane of the multiorder diffraction beam resulting from the SLM pixelated structure. As shown in Fig. 6, it can be seen that the zero-order light is suppressed and the character pattern can be entirely exposed.

## 6.

## Conclusion

In order to improve the energy utilization rate in holographic femtosecond laser parallel processing, this paper first analyzes the impact of SLM pixelated structure on phase-only hologram encoding, and a phase-only hologram containing a digital blazed grating and divergent spherical wave are designed and generated to improve the energy utilization rate and suppress the zero-order light. The experimental results indicate that the proposed method can achieve high-quality holographic femtosecond laser parallel processing with a significantly improved energy utilization rate. The proposed method can also be used to eliminate the zero-order interruption in optical systems using an SLM pixelated structure.

## Acknowledgments

This work is supported by National Science Foundation of China (Nos. 51275502 and 61377006), Anhui College of Natural Science Foundation of China (KJ2011A014), and Provincial Natural Science Foundation of China (1308085ME76).

## References

X. Donget al., “Femtosecond laser two-photon micro-/nano-fabrication and its applications,” Chin. Sci. Bull. 53(1), 2–13 (2008) (in Chinese).CSBUEF1001-6538Google Scholar

Y. Hayasakiet al., “Variable holographic femtosecond laser processing by use of a spatial light modulator,” Appl. Phys. Lett. 87(3), 031101 (2005).APPLAB0003-6951http://dx.doi.org/10.1063/1.1992668Google Scholar

H. LinB. JiaM. Gu, “Dynamic generation of Debye diffraction-limited multifocal arrays for direct laser printing nanofabrication,” Opt. Lett. 36(3), 406–408 (2011).OPLEDP0146-9592http://dx.doi.org/10.1364/OL.36.000406Google Scholar

Y. Hayasakiet al., “Experimental investigation of the closest parallel pulses in holographic femtosecond laser processing,” Appl. Phys. A 107(2), 357–362 (2012).APAMFC0947-8396http://dx.doi.org/10.1007/s00339-012-6801-1Google Scholar

S. HasegawaY. Hayasaki, “Adaptive optimization of a hologram in holographic femtosecond laser processing system,” Opt. Lett. 34(1), 22–24 (2009).OPLEDP0146-9592http://dx.doi.org/10.1364/OL.34.000022Google Scholar

S. HasegawaY. Hayasaki, “Holographic vector wave femtosecond laser processing,” Int. J. Optomechatron. 8(2), 73–88 (2014).1559-9612http://dx.doi.org/10.1080/15599612.2014.901456Google Scholar

D. PalimaV. R. Daria, “Holographic projection of arbitrary light patterns with a suppressed zero-order beam,” Appl. Opt. 46(20), 4197–4201 (2007).APOPAI0003-6935http://dx.doi.org/10.1364/AO.46.004197Google Scholar

J. ChristmasN. CollingsA. Georgiou, “Blocking zero-order in phase shift hologram generation,” UK Patent GB2438458 (2007).Google Scholar

M. Silvennoinenet al., “Parallel femtosecond laser ablation with individually controlled intensity,” Opt. Express 22(3), 2603–2608 (2014).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.22.002603Google Scholar

E. H. WallerG. von Freymann, “Multi foci with diffraction limited resolution,” Opt. Express 21(18), 21708–21713 (2013).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.21.021708Google Scholar

H. Zhanget al., “Elimination of a zero-order beam induced by a pixelated spatial light modulator for holographic projection,” Appl. Opt. 48(30), 5834–5841 (2009).APOPAI0003-6935http://dx.doi.org/10.1364/AO.48.005834Google Scholar

R. W. GerchbergW. O. Saxton, “A practical algorithm for the determination of phase from image and diffraction plane pictures,” Optik 35, 237–246 (1972).OTIKAJ0030-4026Google Scholar

K. L. Tanet al., “Dynamic holography for optical interconnections. II. Routing holograms with predictable location and intensity of each diffraction order,” J. Opt. Soc. Am. A 18(1), 205–215 (2001).JOAOD60740-3232http://dx.doi.org/10.1364/JOSAA.18.000205Google Scholar

Y. YuT. WangH. Zheng, “Optimization of optoelectronic reconstruction of phase hologram by use of digital blazed grating,” Acta Phys. Sin. 58(5), 3154–3160 (2009) (in Chinese).WLHPAR1000-3290Google Scholar

C. Zhanget al., “A rapid two-photon fabrication of tube array using an annular Fresnel lens,” Opt. Express 22(4), 3983–3990 (2014).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.22.003983Google Scholar

N. J. Jennesset al., “A versatile diffractive maskless lithography for single-shot and serial microfabrication,” Opt. Express 18(11), 11754–11762 (2010).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.18.011754Google Scholar

## Biography

**Chaowei Wang** is a postgraduate at Anhui University. He received his BS degree in automation from Anhui University in 2013. His current research interests include holographic femtosecond laser processing.

**Yahui Su** is an associate professor at Anhui University. He received his BS degree in mechatronics and his MS degree in mechanical design and theory from Hefei University of Technology in 1997 and 2002, respectively, and his PhD degree in precision instrumentation and machinery from the University of Science and Technology of China in 2007. He is the author of more than 20 journal papers. His current research interests include laser microprocessing, holography, and optoelectronic systems.

**Jinli Wang** is an engineer at Anhui Jiyuan Electric Power System Tech. Co. Ltd. He received his BS degree in measurement and control technology and instrumentation from Anhui University in 2010 and a master’s degree in measurement technique and automation equipment from Anhui University in 2013. His current research interests include holographic femtosecond laser processing.

**Chenchu Zhang** is a PhD candidate at University of Science and Technology of China. He received his BS degrees in Tianjin University in 2010. His current research interests include two-photon polymerization fabrication and laser induced nanostructures.

**Ziyu Zhang** is a postgraduate at Anhui University. He received his BS degrees in measurement and control technology and instrumentation from Anhui University in 2014. His current research interests include holographic femtosecond laser processing.

**Jiawen Li** is a lecturer at the University of Science and Technology of China. He received his BS degree from Hefei University of Technology in 2003, and PhD degrees from University of Science and Technology of China in 2011, respectively. His research interests include fundamentals of structural color, femtosecond laser fabrication, and ultrafast laser-matter interaction.