In recent years, with the developments of the laser technology, femtosecond laser technology is becoming more and more consummate as a novel technology. Femtosecond laser pulse is also becoming into a powerful tool for microfabrication and micro-machining of various multi-functional structures in dielectric materials through multi-photon absorption because of its high-quality and damage-free processing. Up to now, many high-quality material processing techniques have been achieved by using femtosecond laser pulses with the methods of directly writing [1-10] and holographic fabrication [11-22], such as waveguide , special diffractive optical elements (DOE) [4-10], micro-gratings [11-15], and photonic crystals [16-20]. Because multiphoton nonlinear effects play a major role in this process, the resulting change in refractive index or cavity formation can be highly localized only in the focal volume where the fluence is above a certain material dependent threshold, which makes it possible to micro-fabricate devices inside the bulk of transparent materials. These structures were usually fabricated with a focused beam and written dot by dot by translation of the sample with respect to the focal point.
Compared to directly writing technology by femtosecond laser pulses, holographic lithography has been considered as a more effective method for fabricating periodic structures because it can be controlled easily by the number of the beam, angles between every two beams, energies of the laser beams, and most importantly, only one pulse needed for holographic fabrication. Especially, holographic lithography is considered to be the most effective method for the fabrication of the photonic crystals. And Cai
There are also many groups who have attempted to fabricate the periodic structures in photosensitive transparent materials or on the surface of the silica glass and the metal film by interfered multiple femtosecond laser pulses. Especially, it is very easy to fabricate the periodic structures inside of the photosensitive transparent materials with the aid of diffractive beam-splitter (DBS) [18-20]. However it is limited by the angles between two beams and the energy of the pulses. It is therefore difficult to fabricate microstructure with smaller period, especially in the materials with big band gap, such as silica glass. So there are also some researching groups [13, 16-17] focusing on the realization of periodic structures on the surface of the silica glass by a single shot of two or three femtosecond laser pulses originating from one pulse by the beam splitters.
In this chapter, we have reviewed the progress of fabrication of the periodic structures on the surface or inside of the transparent silica glass by a single shot of several (two, three, four) femtosecond laser pulses. When a single shot of two pulses interfered with each other, there will be one dimensional grating structures being formed inside of the photosensitive transparent materials or on the surface of the materials. When a single shot of three coplanar pulses interfered with each other, the one-dimensional M-shape grating can be formed on the surface of the silica glass. However, when a single shot of noncoplanar three or four pulses interfered with each other, two-dimensional periodic microstructure can be obtained, which distributed as a hexagonal lattice or tetragonal lattice.
2. Holographic fabrication by a single shot of two femtosecond pulses
2.1. The formation of the common grating and the extraordinary grating on the surface of the silica glass
The experimental setup for the holographic fabrication of microgratings by a single shot of two interfered femtosecond laser pulses is depicted in Fig. 2 (a). A regeneratively amplified Ti:sapphire laser (Coherent. Co.) with a central wavelength of 800 nm, pulse duration of 120 fs, and pulse repetition of 1–1000 Hz was used. A single femtosecond laser pulse with a beam diameter of 6 mm could be selected and split into two pulses that were then redirected at controllable incident angles on the surface of the fused silica glass (K9 glass) which is transparent for the laser with a wavelength of 800 nm. These split pulses generating from the single pulse were focused on the glass surface by two lenses with focal lengths of 20 cm to give a spot size of ~50 μm at the focal plane. Attenuators could be used to obtain a proper energy of the pulses for controlling the results of the holographic fabrication. The two pulses could be adjusted both spatially and temporarily by the optical delay device perfectly. The angle between two pulse and the pulse energies could be controlled easily.
When we set the angle between the split beams as 40o and the pulse energy as 45 μJ, the fabricated grating is shown in the Fig. 2 (b-e) by optical microscopy and atomic force microscopy (AFM) respectively. From the Fig. 2 (d), we can obtain the period of the fabricated grating is about 1.06 μm which agree well with the calculated result according to the formula d=λ/[2sin(θ/2)] in which θ is colliding angle (θ=40o in this experiment) between two incident beams and λ is the incident laser wavelength of 800
However, when we set the angle θ as 20o and the pulse energy as 65 μJ, the fabricated grating is shown in the Fig. 3. Not only did we get the ordinary grating whose periods accorded with the theoretic equation d=λ/[2sin(θ/2)], but also obtained the extraordinary grating [15, 23] whose period is a half of the ordinary grating. The extraordinary grating formed at the middle of each bulge of the ordinary grating as shown in Fig. 3 (a-b), so the period of the extraordinary grating is a half of the ordinary grating’s observed easily from Fig. 3(c) which also shows that the depth of the extraordinary grating is a half of the ordinary grating’s nearly.
In the experiments, with the decrease of the incident energy in the same angle between two beams, the modulation depths of the extraordinary grating are decreasing gradually. With the increase of the angle between two beams with the same energy, the modulation depths of the extraordinary grating are also decreasing gradually. At last the modulation will be vanishing from the central part to the edge of the ordinary grating gradually.
The formation of the modulation grating could be attributed to the higher-order modulation arising from second-harmonic generation (SHG) when the femtosecond laser pulse was incident to the surface of silica glass. As a rule, because of the inversion symmetry of the silica glass, there should be no second order nonlinearity in the silicon glass. However, a layer of plasma could be formed on the surface of the glass when the pulse incident to the sample because of the ultra high electrical field of the femtosecond laser in a time given by the duration of the laser pulse[24,25]. Therefore, when the femtosecond laser pulse with high intensity is incident on this thin plasma layer, there will be the higher harmonic generation because of the electrons quivering nonlinearly, and the SHG can reach to 2% of the fundamental laser . The second-harmonic radiation is emitted in the direction of the reflected fundamental laser direction as depicted in Fig. 2 (d) . So the period of the modulation grating is a half of the common grating because of the same of the angle between the two beams but a half of the wavelength. Although the intensity of the SHG is much smaller than the fundamental laser, the modulated depth of the silica glass can reach to a correspondingly large depth because it is much easier for the silica glass to be ionized by a mechanism of single or two photons absorption than multiphoton absorption.
2.2. Fabrication of multiple layers of grating
The multiple gratings were written inside of the soda–lime glass one after another by focusing the beam 1 and beam 2 into special position inside of the glass sample with the recording plane. Without loss of generality, three layers of micrograting could be recorded at depths of 200, 400, and 600 μm, respectively. The grating at depth of 600μm should be encoded firstly. Then the sample was translated along the
Because the micrograting was formed inside of the samples around the focal point without damage to the surface or other parts of the sample, multiple layers of grating can be recorded successfully. And the images of the fabricated multiple layers of microgratings could be read out by beam 2 and taken over three different recording planes as depicted in Fig. 4 (b) and (c). The schematics on the right illustrate the relative positions of the sample, three layers of grating, the recording plane, and the objective. The experimental results demonstrated that the readout image of the first layer (not shown here) and the second layer (as depicted in Fig. 4 (b)) of the recorded grating consisted of well-defined straight bright lines alternating with black lines. However, for the readout of the third layer of the micrograting, some of the straight lines became curved as shown in Fig. 4 (c). This aberration may be caused by the wave-front distortion of the incident beams due to the accumulation of nonlinear effects when the focused high peak-power pulse propagated through the sample. The longer the optical path inside the sample was, the more severe the wave-front distortion and the more obvious the resulting grating aberration.
2.3. Fabrication of micrograting inside the glass doped Au nanoparticles
Noble metal nanoparticle-contained glasses exhibit large third-order nonlinear susceptibility and ultrafast nonlinear response due to the local field effect near surface Plasmon resonance and quantum size effect [27, 28]. In recent years, many studies have been carried out on the fabrication of nanoparticle-doped glasses [29-31]. Shiliang Qu
A typical silicate glass is composed of 70SiO2.20Na2O.10CaO doped with 0.1Au2O (mol%). Reagent grade SiO2, CaCO3, Na2CO3, and AuCl3.HCl.4H2O were used as starting materials. An approximately 40g batch was mixed and placed into a platinum crucible. Melting was carried out in an electric furnace at 1550 °C for 1 h. The glass sample was obtained by quenching the melt to room temperature. The sample thus obtained was transparent and colorless. The annealed sample was cut and polished, and then subjected to successive experiments for femtosecond laser.
The used laser system and the experimental setup are the same as the used setup as shown in Fig. 2 (a). The two incident beams were first focused onto the front surface of Au2O-doped glass to optimize the incident pulse energy. In the case of sufficiently high energy, the two coherent beams can induce periodic ablation, forming a grating in the glass. Herein, we reduced the incident pulses’ energy to a certain lower level at which the two coherent beams cannot directly induce periodic ablation on the surface of the glass. Then the sample was moved 50μm to the lens in the direction of angular bisector of two incident light paths and made the laser pulses be focused inside the glass. After irradiation by a single shot of two interfered pulses, the micrograting can be recorded inside the glass, which can not be observed immediately because there is no obvious changes in the focusing place, however, after heat treating the sample at 550 °C for 1 h, the formed micrograting can be observed because of the Au nanoparticle precipitation. Such grating was constituted by the laser-heating induced Au nanoparticle precipitation in the Au2O-doped glass.
In Qu’s experiments, this lower pulse energy was selected to be 30 and 38 μJ for comparison, and the colliding angle θ between the two incident beams was fixed at 45°. The fabricated micrograting is shown in Fig.5 (a-d). The period d of the obtained gratings was about ~1μm which was agreeing well with the value calculated from the colliding angle θ and the laser wavelength λ according to the formula d=λ/[2sin(θ/2)].
The absorption spectrum of the grating was measured by a spectrophotometer (JASCOV-570) as shown in Fig. 5 (e) (line a). Apparently, a weak peak occurs around 508 nm, which is induced by the surface plasmon resonance of Au nanoparticles in the glass. The Au nanoparticles with 3 nm average size in the glass were observed in the grid of the fabricated micrograting from a transmitted electronic microscopy (TEM). However, if the glass sample was irradiated only by the interfered pulses but not heat treated, neither could absorption peak in the range of 500–600 nm be observed in the absorption spectrum as shown in Fig. 5 (e) (line b), nor could Au nanoparticles be observed by the TEM. This indicates that the Au nanoparticles can be precipitated in the periodic one-dimensional arrays in the glass through the irradiation of two coherent beams with the aid of heat treatment.
2.4. Fabrication of multiple gratings by a single shot of two interfered femtosecond laser pulses
The interference of two laser beams can create microstructures with one-dimensional periodic patterns in certain materials due to the periodic modulation of the laser intensity with a period scale of the order of the laser wavelength. As stated as above, there has been many groups focusing the fabrication of microgratings in glasses [11, 12], thin films , polymers, and inorganic–organic hybrid materials  by use of two interfered femtosecond laser pulses in a single shot. However, there is a disadvantage for micrograting fabrication with this method: Only a single micrograting can be formed for one pulse. Shiliang Qu
In the experimental setup as shown in Fig. 6 (a) which is similar with the Fig. 2 (a) and Fig. 4 (a), a regeneratively amplified Ti:sapphire laser (Spectra-Physics) with a wavelength of 800 nm, pulse duration of 120 fs, and pulse repetition of 1–1000 Hz was used. A single laser pulse with a beam diameter of 8 mm was selected and split into two beams that were then redirected at approximately equal incident angles on a silicate glass surface. The two beams were focused on the glass surface by two lenses with focal lengths of 10 cm (L1) and 20 cm (L2) to yield spot sizes of 40μm and 80μm, respectively.
The colliding angle θ between the two beams was fixed at 40o. After the optical paths were adjusted to realize perfect overlap of the two beams both spatially and temporally, the surface of the glass was adjusted to be approximately normal to the perpendicular bisector of the two incident beams, so that the glass surface became the laser interfering plane. A mask used for laser beam modulation was placed in the optical path in which lens L1 is located. The mask consists of three equilaterally and triangularly arrayed apertures, whose diameters are 2.5 mm and the spaces between them are 3.5 mm as shown in Fig. 6(b).
After irradiation by one single shot of two interfered pulses, the multi-micrograting comprises three microgratings as shown in Fig. 6(c), which has a high fidelity to the configuration of the mask used. The multi-micrograting was formed through periodic ablation resulting from the interference of the reference laser beam with the three beams caused by the mask. And the period of the formed microgratings is also agree well with the theoretical expectation of the common grating d=λ/[2sin(θ/2)]. This means that a multi-micrograting comprising even more microgratings and configurations can also be one-off written by changing only the mask structure. However, multi-micrograting formation can be realized only when the two interfered beams are overlapped at an appropriate position out of their focus on the front surface of the glass.
3. Holographic fabrication by a single shot of multiple femtosecond pulses
3.1. The formation of the M-shaped micrograting by three coplanar interfered beams
As stated above, two interfered femtosecond pulses’ interference can induce one and two-dimensional periodic structures by single and double-exposure techniques, respectively [11−15, 40]. However, usually the second pulse could not overlap completely with the microstructure formed by the first pulse in the double-exposure technique due to the rather small size of the focal spot. Here, we will show that the fabrication of M-shape gratings with controllable modulation depth could be realized by adding the third beam into a two-beam interference system. The experimental results show that the depth ratio between neighbor grooves can be conveniently controlled by changing the pulse energy of the third beam. Morphology characterizations of as-fabricated periodic M-shape gratings with a period of 2
The schematic experimental setup for fabricating M-shape grating is shown in Fig.7. The used ultrafast pulses (with a pulse width of 120 fs, central wavelength of 800 nm, and repetition rate of 1–1000 Hz) are produced by a Ti: sapphire regenerative amplified laser system (Coherent Inc). A laser pulse with a diameter of
The optical microscope images of the gratings formed on the silica glass at different pulse energies of L2 (
It is obviously shown that the structure is formed by periodically arranging the M-shape units with a size of about 2.6 μm. The modulation depths of the deeper grooves and shallower grooves of the M-shape grating are ~500 nm and ~240 nm, respectively. The deeper grooves and the shallower grooves of the M-shape grating each have a period of
The simulated intensity distributions in the interfering region formed by three coplanar interfering beams are shown in Figs. 9 (a), (b) and (c), which correspond to the middle beam (L2) pulse energies of 20 μJ, 50 μJ and 100 μJ, respectively. In Figs. 9 (a), (b) and (c), the upper images show the three-dimensional patterns and the nether curves display the two-dimensional cross-section patterns. From the cross-section images of Figs. 9 (a-c), we can see clearly that the periodic intensity patterns are formed by the periodic arrangement of inversely M-shape structures, which is attributed to stronger intensity peaks and weaker ones arrayed with the same period of about 2.6 μm alternately and periodically. These results are consistent with the AFM results as shown in Fig. 8. We notice that with the energy of the middle pulses increasing, the stronger intensity peaks in the interfering region increase, while the weaker intensity peaks decrease. It indicates that the intensity ratio of the stronger peak and the weaker one increases. If laser pulses with such a periodic inversely M-shape intensity distribution irradiate the materials, the stronger intensity peaks will lead to the formation of deeper grooves, and the weaker intensity peaks will induce the formation of shallower ones. As a result, the M-shape surface structures will be resultantly formed. This fabricating technology for the formation of M-shape gratings provides a fabricating method for special gratings with special use in industrial applications, such as the calibration for 3D reconstruction in computer vision application. And the fabricated M-shape grating can also be used in microfluidic chip devices as transport channels with different flowing speeds.
3.2. The formation of the triangular lattice by three noncoplanar interfered beams
When a single shot of three coplanar pulses interfered with each other, the M-shaped gratings could be fabricated as above. However, when a single shot of three non-coplanar pulses interfered with each other as depicted in Fig. 10 (a) by adjusting two time-delay for obtaining perfect overlap of the three pulses both spatially and temporarily, the two-dimensional periodic microstructure have been obtained, which distributed as a hexagonal lattice as shown in Fig. 10 (b) and agreed well with the simulated results [see Fig. 10 (c)] [16, 17]. In experiments, the geometric angles were kept as θ=300 and ϕ=350 respectively as shown in Fig. 10 (a). When we set the pulse energy as 75 μJ, we can obtain two-dimensional periodic hexagonal lattice of microholes [see Fig. 10 (d-f)], however, when we set the pulse energy as 30 μJ, the two dimensional periodic microstructures present doughnut orbicular platform [see Fig. 10 (g-i)]. The period in the direction of line “
The different microstructures in our experiments were attributed to the formation of plasma and molten liquid at different pulse energy levels. Generally it is hard to form plasma and molten liquid on the surface of silica glass by a laser with a wavelength of 800 nm because of the bigger band gap. However, the intensity at the focal point of the femtosecond laser beam where three beams interfere together could reach to 100 TW/cm2 nearly. Such a high-energy influence within the focal volume ionized the silica glass quickly through the combined action of the avalanche and multiphoton processes [17, 35]. A layer of plasma formed on the surface of the silica glass at the time of the laser pulse duration. While the intensity decreased to be lower than a certain value, which can be called the ionized threshold as shown in Fig. 11 (a), the plasma vanished and the molten liquid of the material appeared. With the decrease of the intensity to be lower than a certain value, which could be called the molten threshold [see Fig. 11 (a)], the molten liquid of the material disappeared, so there would be a doughnut molten liquid formed in every enhanced spot of the interfered field [see Fig. 11 (b)]
When the pulse energy was set as 75 μJ, a layer of modulated plasma was formed by the interfered field with hexagonal lattice depicted in Fig. 10 (c) after the anterior part of the pulse was incident to the surface of the silica glass. Higher intensity induced a relatively larger plasma area in every enhanced spot of the interfered field. The subsequent posterior part of the pulse removed the plasma very swiftly because of the high light pressure originating from the reflection. Therefore, the microholes formed on the surface of the silica glass, as depicted in Fig. 10 (d-f). When the pulse energy decreased to 30 μJ, there was just a layer of the plasma with smaller areas than that of 75 μJ in the center of every enhanced spot in the interfered field after the anterior part of the pulse was incident to the surface. The
subsequent posterior part of the pulse also removed the central plasma very swiftly by means of light pressure, so there is a tiny hole formed in the center of the enhanced spot. However, in the region of the molten liquid, there were two distinct interaction components because of Marangoni effect, thermocapillary and chemicapillary, which resulted from the thermal potential of a temperature gradient and the chemical potential of a compositional gradient, respectively [17, 36-37]. The thermocapillary force moved the molten material outward from the center, while the chemicapillary force moved the molten material toward the center [38-39] as depicted in Fig. 10 (b). When the chemicapillary force dominated, a platform formed in the center of the spot. The combined action of the light pressure to the plasma and the chemicapillary force to the molten liquid induced a periodic orbicular platform on the surface of the silica glass as depicted in Fig. 10 (g-i).
In order to verify our explanation on the different induced microstructures, we set the pulses energy as 50 μJ, and the resulted microstructure is shown in Fig. 12 (a). Because of the Gaussian type intensity distribution of the pulse, three regions are selected from the central part (I), outer part (II) and the edge part (III) of the microstructure as depicted in Fig. 12 (a). The corresponding three dimensional images of the selected areas are shown in Fig. 12 [(b)-(d)], respectively. For the region I, the intensity is high enough to make the light pressure dominate, therefore the periodic microholes formed [see Fig. 12 (b)]. In contrast, in the region II, where the intensity is comparatively low, the chemicapillary force and the light pressure dominate, so that the microstructure present a orbicular platform [see Fig. 12 (c)]. However, in the selected area III, the intensity is very low and there is no plasma layer but just a layer of liquid formed on the surface of the silica glass. In this case, if the chemicapillary force dominates, the microcones can be observed. As shown in Fig. 12 (d), just several microcones (indicated by arrows) formed in the edge of the interfered districts.
3.3. The formation of the tetragonal lattice by four noncoplanar interfered beams
The interference of two beams creates a one-dimensional (1D) periodic pattern. By increasing the number of beams [16-20] or the double-exposure techniques , in principle, two-dimensional (2D) and three-dimensional (3D) periodic patterns can be designed. Although, as stated above, the interference of the three beams [13, 16-17] can create a 1D  or 2D [16, 17] periodic pattern on the surface of the materials of inside the transparent materials, the complicated optical setup is required for the interference of multiple laser beams, and its precise adjustment is difficult. Kondo
The optical setup used for the present experiments is shown in Fig. 13(a). Briefly, a DBS (G1023A or G1025A; MEMS Optical Inc.) divides the input laser beam into several beams, and the beams are collected on the sample by two lenses. Temporal overlap of the divided pulses is achieved without adjusting the optical path lengths. Each beam was made to be parallel or slightly focused by the adjustment of the distance between the two lenses. Slight focusing increased laser power density and helped to make the MPA efficient. The beams meant to form interference were selected by an aperture array, which is placed between the two lenses. Negative photoresist SU-8 (Microlithography Chemical Corp.) was used as an initial material for the fabrication. The absorption spectrum of SU-8 indicates that one-photon absorption is negligible at an 800-nm wavelength. Consequently, it is expected that photopolymerization, if occurring, is due to a multiphoton reaction. The layer of SU-8 was spin-coated on a coverglass plate having a thickness of about 4 μm, and prebaked before exposure to fs pulses. The interference angles θair (the angle between the main optical axis and the other beams in air) applied in the experiments were measured to be 33.6°, 21.9°, and 10.8°.
Femtosecond pulses from a Ti:sapphire regenerative amplifier (wavelength of 800 nm, pulse duration of 150 fs, repetition rate of 1 kHz) were used for experiments. By using this method, periodic structures were fabricated. Periodic structures fabricated with an interference angle of 33.6° are presented in Fig. 14. In this figure, (a) and (b) show SEM images of the same sample from different perspectives. The oblique view shown in (b) clearly demonstrates that the periodic structure consists of high-aspect-ratio rods. The height and diameter of the rods are about 4 μm and 0.6 μm, respectively; thus, an aspect ratio of about 7 was achieved. It should be noted that the height of the obtained structure is not limited by the coherent length of the pulse. A higher structure could be obtained if the stiffness of the material allowed the rods to withstand capillary forces during the development procedure.
Figure 14 (c) provides a close-up view of the rods fabricated with an interference angle at 21.9°. To obtain this image, the coverglass containing the fabricated structures was deliberately broken and a small fragment was observed. As seen, the rods are slightly bellows-shaped. Distinct ring-like features repeat periodically along each rod with a period about 0.3μm. This could be attributed as the result of interference between the incident and reflected beams at the resist–coverglass interface.
In conclusion, we have reviewed the fabrications of the one-dimensional and two-dimensional periodic microstructures on the surface or inside of materials by multiphoton absorption using a single shot of two or multiple interfered femtosecond laser pulses.
Firstly, we have introduced the fabrication of the microgratings by a single shot of two interfered femtosecond laser pulses. When two interfered pulses overlapped on the surface of the silica glass, not only did we get the ordinary grating whose periods accorded with the theoretic equation d=λ/[2sin(θ/2)], but also obtained the extraordinary grating whose period is a half of the ordinary grating. The formation of the modulation grating could be attributed to the higher-order modulation arising from second-harmonic generation (SHG) when the femtosecond laser pulse was incident to the surface of silica glass. The multiple layers of the microgratings have been fabricated in the different depths of the silica glass sample by a single shot of two interfered femtosecond laser pulses, and the results show that the fabricated gratings can be read out by one of the recorded beams very easily. The noble metal nanoparticles consisted microgratings have also been realized in silicate glasses by two interfered femtosecond laser pulses with the aid of heat treatment because the noble metal nanoparticles in silica glass can be precipitated after the irradiation of the femtosecond laser and the successive heat treatment. At the same time, the multiple gratings can also be realized on the surface of the silica glass samples by a single shot of two interfered femtosecond laser pulses with the aid of a mask, which is very significative for enhancing the processing efficiency of the fabricated microgratings.
Secondly, we have also introduced the fabrication of the 1-D or 2-D periodic microstructures by a single shot of multiple interfered femtosecond laser pulses. When a single shot of three coplanar interfered femtosecond laser pulses interfered with each other on the surface of the silica glass, M-shaped gratings can be formed and the morphologies of the M-shaped gratings can also be modulated by tuning the incident pulse energy. However, when a single shot of noncoplanar three or four pulses interfered with each other, two-dimensional periodic microstructure can be obtained, which distributed as a hexagonal lattice. Different morphologies of the induced structures such as microvoid, orbicular platform and nanotip, could be formed with the changes of the incident pulse energy. The fabrication of the 2D tetragonal lattice have also been fabricated in Negative photoresist SU-8 by four noncoplanar interfered pulses originating from the single pulse with the aid of the special DBS (diffractive beam splitter). Although the special DBS provided an easy way for realizing multiple beams interference, it also had some shortages, such as difficulty for fabricating microstructure with a smaller period, especially in the materials with big bandgap, such as silica glass.
This work was supported by the National Science Foundation of China (NSFC: 10904027; 61108018), the China Postdoctoral Science Foundation (AUGA41001348) and the Heilongjiang Province Postdoctoral Science Foundation (AUGA1100074), and development program for outstanding young teachers in Harbin Institute of Technology, HITQNJS. 2009. 033.