Fabrication and Characterization of 2D Nonlinear Structures Based on DAST Nanocrystals and SU-8 Photoresist for Terahertz Application

1. Introduction

Within the last ten years, the field of terahertz (THz) science and technology has developed greatly. Many new and promising advances in the technology for the generation, manipulation, and detection of THz radiation have been made. The most important task is the development of the new materials and technology for THz generation and detection [1]. Optical rectification in electro-optical crystals has been demonstrated as the best technique for generating subpicosecond, THz bandwidth radiation [2]. For THz application, the three most important requirements for nonlinear optical crystals are i) high nonlinear coefficient, ii) low absorption of THz irradiation by the crystal, and iii) phase matching between pump and THz waves propagations.

Several groups in the world have examined different nonlinear materials for THz emission by optical rectification. These materials include organic salts such as dimethyl amino 4-N-methylstilbazolium tosylate (DAST) [3], dielectric inorganic electro-optical materials like LiTaO₃, and semiconductors such as GaAs and GaP. Among them, DAST is the most efficient THz rectification material [5], due to its large electro-optics (EO) coefficient and its low dielectric constant ($ϵ=5.2$), giving rise to a high modulator figure of merit [5]. Indeed, according to data in [6], DAST exhibits a huge second-order nonlinear optical susceptibility ${\chi }^{\left(2\right)}=2020\pm 220$ pm/V at $\lambda =1318$ nm and its corresponding electro-optical figure of merit being $n_1^3{r}_{11}=(530\pm 60)$ pm/V at $\lambda =1313$ nm. However, phase matching between optical wave and THz radiation propagation, a major factor determining the efficiency and the bandwidth of THz emission from nonlinear crystal, is still a challenge. Actually, only few electro-optical crystals possessing a long coherence length is suitable for THz generation. Most nonlinear materials possess a coherence length in a range of few micrometers to dozens of micrometers, thus limiting the THz generation efficiency. It is necessary to find out a way to solve the problem of phase mismatch in DAST crystals.

Quasi-phase matching (QPM) is a technique which compensates the phase mismatch of nonlinear processes by creating a grating of spatially modulated nonlinear properties with a period corresponding to twice of coherence length. Typically, the QPM structures have periods in the range of few micrometers to dozens micrometers, which is quite easy to be realized today. In a quasi-phase matched structure, the interaction length between THz and pump pulses can be dramatically extended to centimeters range. For example, QPM structures were used for the generation of THz wave packets with tunable central frequency [7] or for the generation of multicycle narrow-bandwidth THz radiation based on phase synchronous optical rectification [8]. It is clear that nonlinear gratings or QPM structures may provide an efficient control of THz generation.

Conventionally, the QPM effect is obtained in a nonlinear material in which the nonlinear property is modulated between negative and positive signs (+/- QPM). However, this is not easy to realize in practice, in particular for three-dimensional (3D) QPM structures. An alternative way to achieve the same effect is to combine the nonlinear material with a linear material, corresponding to a periodical modulation of ${\chi }^{\left(2\right)}$ between zero and non-null value (+/0 QPM). The theoretical model of this particular QPM structure was presented in our recent work [9] for a general nonlinear optical effect. Very few works have been demonstrated to realize these +/0 QPM structures. For example, one of our old works [10] demonstrated a layer-by-layer direct laser writing technique to obtain a polymeric 3D QPM structure. This technique has been applied for copolymer nonlinear material requiring an additional poling technique, since the copolymer is centrosymmetric. It could be applied for patterning DAST nanocrystals as well. However, this technique requires spin-coating of different materials layer by layer and thus very time consuming, and imprecise in term of periodicity. Recently, two other groups experimentally demonstrated a direct fabrication by employing a femto-second laser based bleaching technique [11,12] in LiNbO₃ material. We also applied this technique to patterning DAST nanoscrystals structures, as shown in reference [15]. This direct technique allows to solve the problem of imprecision, but again, it does not produce a large QPM structure in a short time.

There are different methods to realize QPM structures based on DAST nonlinear material. However, DAST requires a customized approach, because DAST material is usually fabricated in the form of micro- and nanoparticles. By using the assembling of multiple layer technique, DAST/SiO₂ multilayers has been proposed for THz generation [13]. We can see that QPM elements based on DAST and SiO₂ micrometer-size layers can be effective for THz generation, but their fabrication is still a great challenge. In our recent work, we have demonstrated a way to obtain desired DAST structures by mixing DAST nanocrystals with a passive polymer matrix. Indeed, we have successfully obtained a nonlinear nanocomposite based on DAST nanocrystals embedded in a PMMA polymer matrix and showed that the nonlinear conversion efficiency doesn’t depend on the DAST crystals size [14]. We have also fabricated 2D grating structures by photobleaching DAST nanocrystals using a technique called one-photon absorption based direct laser writing [15]. The last one is very efficient but rather slow and therefore time consuming. Very recently, we proposed an alternative method by combining a mask lithography and thermal annealing to produce $+/0$ QPM DAST structure [16]. Highly efficient THz generation by this DAST nanocomposite is demonstrated upon excitation by femtosecond pulses and generation by the method of optical rectification. However, since the mask technique was used, it is not suitable to adjust the structures configurations and periodicity, and also it is not possible to produce small periods and 3D structures.

In this work, we demonstrated an alternative way to fabricate DAST/polymer large structures in a short time, which is based on filling a holographically fabricated SU-8 template by a DAST-PMMA composite. Figure 1 shows the fabrication process of desired DAST/polymer structures. First, this method consists in the preparation of the SU-8 template by an interference technique [17], which has many advantages such as rapidity, simplicity, large and uniform structures. Second, the DAST-PMMA composite is spin-coated into SU-8 template and thermally annealed to remove the solvent and to induce DAST crystallization. The final result is a linear/nonlinear structure having a controllable period and desirable configurations, which is equivalent to a $+/0$ QPM structure, and which could be useful for THz generation. We demonstrated the nonlinear properties of these fabricated structures by measuring the second-harmonic generation (SHG) signals.

2. Theory and numerical calculations of two-beam interference technique

2.1. Theory of two-beam interference

Fabrication of one-, two- and three-dimensional (1D, 2D, and 3D) micro and nano-structures is of great interest nowadays. Multi-beam interference technique allows getting different photoresist-based structures depending on the number of laser beams and their arrangement. The common advantage for all holographic techniques is a parallel fabrication process, which results in uniform structures over large area and in a short fabrication time. However since multiple laser beams should be used, the fabrication setup becomes quite complicated [18,19]. Two-beam interference technique processes many advantages because of its simplicity [17]. In particular, this technique provides high contrast between the minimal and maximal intensities of interference pattern because of identical polarization of two laser beams in interference area. Furthermore, multi-exposure of two-beam interference pattern allows to realize various structures with different configurations in multi-dimensions. In this work, we employed this two-beam interference technique to realize DAST nonlinear photonic structures.

Figure 2a illustrates the schematic setup of the two-beam interference technique. A lenses system was used to enlarge the initial laser beam and to get large and uniform working area. The extended beam is divided by a prism into two beams of the same profiles and intensities, which are then deviated by two reflecting mirrors towards the photoresist where they interfere, resulting in an intensity modulation. The electric fields of two laser beams can be written as:

$$\begin{array}{ccc}& & {\mathbf{E}}_1(\mathbf{r},t)=\mathrm{Re}\left[{\mathbf{E}}_{01}{e}^{i({\mathbf{k}}_1.\mathbf{r}-\omega t)}\right],\hfill \end{array}$$

(1)

$$\begin{array}{ccc}& & {\mathbf{E}}_2(\mathbf{r},t)=\mathrm{Re}\left[{\mathbf{E}}_{02}{e}^{i({\mathbf{k}}_2.\mathbf{r}-\omega t)}\right],\hfill \end{array}$$

(2)

where k₁ and k₂ are the wave vectors of two laser beams, r is the position vector in the overlapping space, and E₀₁ and E₀₂ are the amplitudes of electric fields. Sum of the electric field of individual plane waves is modulated (interference) at overlapping area:

$${\mathbf{E}}_T(\mathbf{r},t)={\mathbf{E}}_1(\mathbf{r},t)+{\mathbf{E}}_2(\mathbf{r},t).$$

(3)

The interference intensity distribution of the resultant wave is given by

$$I_T(\mathbf{r},t)={\langle {\mathbf{E}}_T^{*}(\mathbf{r},t).{\mathbf{E}}_T(\mathbf{r},t)\rangle }_t,$$

(4)

where ${\langle ...\rangle }_t$ represents the time average of the resultant electric field. The interference of two laser beams is a 1D pattern, whose period ($\mathsf{\Lambda }$) is calculated by

$$\mathsf{\Lambda }={\displaystyle \frac{\lambda }{2sin\theta }},$$

(5)

where $\theta$ is a half-angle between two laser beams and $\lambda$ is the wavelength of the laser beams. In our case, we assume that the polarizations of two beams are the same (vertical direction) for all values of $\theta$-angle allowing to obtain 1D structure with maximal and same contrast for any periodicity. According to the equation 5, the structure period can be adjusted from as small as $\lambda /2$, $i.e.$ few hundreds of nanometers, to as large as we want. It means that this is fully adaptable to realize the QPM structures with desired period, which is twice of the coherence length of the nonlinear crystal and usually in the range of micrometers.

As it is shown in the Figure 2a and Figure 3, the photoresist sample is fixed on a rotating holder. Both rotation angles, $\alpha$ and $\beta$, influence resultant 1D structure. Therefore, we realize a rotation of the 1D interference pattern around two axes, z by an angle $\alpha$, and y by an angle $\beta$, as seen in Figure 3b,c. The corresponding rotation matrices are:

$$R_z\left(\alpha \right)=\left[\begin{array}{ccc}cos\alpha & -sin\alpha & 0\\ sin\alpha & cos\alpha & 0\\ 0& 0& 1\end{array}\right],$$

(6)

$$R_y\left(\beta \right)=\left[\begin{array}{ccc}cos\beta & 0& sin\beta \\ 0& 1& 0\\ -sin\beta & 0& cos\beta \end{array}\right].$$

(7)

In this case, the resultant intensity distribution becomes:

$$I_{\alpha ,\beta }=2E_0^2{cos}^2[kzsin\theta sin\beta +kcos\theta cos\beta (xcos\alpha +ysin\alpha )],$$

(8)

where |E₀₁|=|E₀₂|=$E_0$.

If a multiple-exposure of the two-beam interference pattern is applied, the total intensity from all exposures is calculated by:

$$I_{\mathrm{multiple}}=\sum _i{I}_{{\alpha }_i,{\beta }_i},$$

(9)

where $i=1,2,3,..$ represent the number of exposures realized at corresponding ${\alpha }_i$ and ${\beta }_i$ angles.

2.2. Numerical calculations

By using personal MATLAB codes, we simulated the resultant structures for different exposures with different parameters. If $\alpha$ and $\beta$ are equal to zero and only one expose is applied, it is possible to get a simple 1D structure as shown in Figure 3d. Additional exposes with $\beta =0$ and $\alpha \ne 0$ provides 2D structure as shown in Figure 3e. The configuration of 2D structures depends on the number of exposures and the rotation $\alpha$–angle. For example, with only two exposures, a 2D square structure is obtained with $\alpha =0^{\circ }$ and $\alpha ={90}^{\circ }$ and a 2D hexagonal structure is obtained with $\alpha =0^{\circ }$ and $\alpha ={60}^{\circ }$. 2D quasi-periodic structures can be also created by this way with more than three exposures. One more rotation by a $\beta$–angle leads to fabrication of 3D structures, as shown in Figure 3f. Here again, the two-beam interference technique with multiple-exposure allows to obtain any 3D structure on demand [20], depending on the number of exposures and the rotation $\alpha$ and $\beta$ angles.

One of the advantages of the interference technique is that it allows to control the filling factor of fabricated structures. Indeed, manipulations with exposure time or power of the laser beam give a precise doses control. It leads to possibility to fabricate variety of structures from so-called “cylinders” (Figure 4a) to “air holes” (Figure 4b). The longer exposure time the more photoresist become cross-linked (in case of negative SU-8 photoresist), resulting in “air holes” structure. In contrast, short exposure time or less power results in separate cylinders. Controlling the filling factor is thus to get better efficiency for different QPM orders [9].

In this work, we will demonstrate the idea by simply fabricating the 2D periodic structures and characterization of their SHG effect.

3. Experimental section

3.1. Preparation of polymeric template

SU-8 2002 (Micro. Chem. Corp.) was spin-coated onto glass substrates with a thickness of $2\mu$m. The sample was soft baked at ${65}^{\circ }$C for 1 min and then at ${95}^{\circ }$C for 2 min to remove the residual solvent. For the preparation of multi-dimensional (1D, 2D and 3D) periodic structures, we have used a multiple-exposure two-beam interference method. For that, a continuous-wave (CW) UV laser (Cobalt) with $\lambda =355$ nm was used, as shown in Figure 2a. After light exposure, the sample was post baked on a hot plate at ${95}^{\circ }$C for 5 min. The sample was then developed by a SU-8 developer for 4 min, and rinsed by isopropanol and DI water. Finally, the sample was hard baked for 40 min at ${180}^{\circ }$C. This annealing of SU-8 matrix at high temperature is needed to remove residual solvents for proper environment for further DAST crystallization.

3.2. Preparation of DAST/PMMA solution

PMMA (Plexiglas® V045) was dissolved in chloroform (CL0218 Scharlau) and DAST (CAS:80969-52-4 Genolite biotek) was dissolved in methanol (322415 Aldrich). The DAST mass is equal to 2% of the PMMA mass, and DAST + PMMA is 15% of mass in chloroform. The amount of methanol was calculated from the solubility limit of DAST in methanol at room temperature [21]. Dissolved DAST was added into PMMA solution and was stirred together to achieve a DAST/PMMA solution with high homogeneity.

3.3. Mechanism of crystallization

The formation of nanocrystals in PMMA is based on DAST aggregation, upon which the size of particles is limited by the diffusion of material through the viscous polymer solution. Both solvents used for PMMA and DAST have close boiling temperatures and relatively close vapor pressures. At the beginning of the process, PMMA and DAST are dissolved in the mixture of two solvents. The formation of the PMMA matrix and DAST nanocrystals in the mixture starts with centrifugation and terminates after complete evaporation of the solvent during annealing. First, the DAST/PMMA composite was dropped on the SU-8 template and spin-coated at a low speed rate to remove extra solution from the top. It forms amorphous DAST particles distributed in PMMA. Second, during annealing at elevated temperature of ${180}^{\circ }$C for 30 min amorphous particles transform into noncentrosymmetric crystalline structure. The mechanism of the formation of DAST crystals in the PMMA matrix and the effect of the synthesis parameters on the crystallite size were considered in detail in our previous works [14,16,22].

3.4. Sample characterization

The linear and nonlinear properties of fabricated samples were characterized by different optical instruments, such as optical microscope, electron scanning microscope, as well as confocal laser scanning microscope. Figure 2b shows a laser system used to characterize both linear and nonlinear properties of DAST periodic structures. Indeed, DAST crystals can absorb and emit the light as a fluorescence emitter. For that the sample was excited by a CW 532-nm Oxxius frequency-doubled Nd-YAG laser (maximum power: 300 mW, coherence length: 300 m; pointing stability: 0.005 mrad/C). DAST crystal also possesses a strong nonlinear response. In our case, the second-harmonic generation (SHG) signal of DAST nanocrystal was demonstrated by exciting the sample by a pulsed 1064 nm laser (maximum averaged power: 80 mW; pulse duration: 1 ns, repetition rate: 24.5 kHz). These two lasers, 532 nm and 1064 nm, are perfectly aligned in front of the objective lens (OL), and were used alternatively, depending on the desired measurement. The laser beam is focused into the sample by a high numerical aperture OL. The fluorescent or SH signals were detected by the OL and sent back to the avalanche photodiode (APD), after filtering out the excitation light beam. By scanning the focusing beam (movement of the piezoelectric translator (PTZ)), we can obtain a color map of fluorescent or SH signals, which confirm the creation of periodic linear/nonlinear patterns. The emission spectra, fluorescence or SHG, were also analyzed by using a spectrometer.

4. Results and discussion

To demonstrate the idea of our proposal, we fabricated different 2D SU-8 templates by using the multiple-exposure two-beam interference technique with a period of $3\mu$m. The fabrication of 2D structure depends on the number of exposures and on sample orientation with respect to each exposure [17]. For example, two mutually perpendicular positions of the sample allowed us to fabricate 2D square structures, as shown in Figure 5b. By using a modest exposure dose (2,1 mW/cm² per interference beam and 2 seconds for each exposure), this fabricated structure is made of SU-8 pillars, as modelling by the simulation, which is shown in Figure 5a.

Figure 5c–d shows the fluorescence mapping image of the final sample, after filling the air-areas of the SU-8 template by the DAST nanocrystals. It is well seen that SU-8 pillars with a weak fluorescence signal are surrounded by DAST/PMMA composite, where DAST crystals show a fluorescence quantum yield of $14-20\%$ [23] and provide a strong fluorescence signal. By scanning the sample along transverse and longitudinal directions, we obtain complete information, $i.e.$ the structure thickness and periodicity. The sample thickness agrees well with the value measured by a profilometer, and the structure period is also very consistent with the designed structure. The fluorescence spectrum of DAST/PMMA area (Figure 6) shows typical fingerprint of DAST crystals with noncentrosymmetric “red-green” form [24].

Similar to the fluorescence experiment, when using a 1064 nm pulsed laser, we obtained the SHG mapping image of corresponding sample. Figure 7 shows the SHG wavelength and mapping image. A clear peak at 532 nm was obtained as a proof of the nonlinear property (SHG) of the DAST nanocrystals. Figure 7b shows that the SH signal is very high from DAST nanocrystals, as compared to the weak signal from SU-8 pillars. Thus, a periodic linear/nonlinear microstructure has been demonstrated by this scanning method. The advantage of this structure as compared to those fabricated by other methods [15] is that there exists a clear border between linear (SU-8 photoresist) and nonlinear (DAST nanocrystal) responses, and the structure is also very large, since it was realized by the interference method, for which the structure area depends on the diameter of the interference laser beam. Unfortunately, we cannot realize QPM measurements with these samples since the film thickness is only $2\mu$m. However, this work demonstrates that the combination of interference template and thermal annealing method is an excellent way to obtain rapidly large QPM structures. We expect to investigate this method for 3D QPM, and apply it for efficient THz generation.

5. Conclusions

In this work, we demonstrated a novel method to realize large nonlinear QPM structures for efficient nonlinear conversion. The proposed method is based first on SU-8 template preparation by a conventional interference method. Then a PMMA-DAST solution was deposited by spin coating and DAST nanocrystals were formed in the air voids of SU-8 template via the thermal annealing method. The nonlinear structure was characterized by SHG and fluorescence 3D microscopy. Mapping results evidences the formation of red crystalline form of DAST nanocrystals grown in holes of SU-8 matrix, forming well delimited periodic nonlinear domains which open interesting perspective for QPM-based nonlinear applications of these structures, such as THz emission by optical rectification.