Flexible and broadband graphene polarizer based ... - OSA Publishing

Vol. 7, No. 4 | 1 Apr 2017 | OPTICAL MATERIALS EXPRESS 1398

Flexible and broadband graphene polarizer based on surface silicon-core microfiber XIAOYING HE1,2,* AND JIACONG LIU1 1

Department of Optical Science and Engineering, Fudan University, Shanghai 200433, China Department of information, Beijing University of Technology, Beijing 100124, China *[email protected]

2

Abstract: A graphene-coated polarizer based on the surface silicon-core microfiber has been proposed and numerically analyzed. The numerical discussion shows that our proposed polarizer serves alternatively as TM-pass or TE-pass polarizers depending on the change of the radius of the silicon core. Their polarization extinction ratios can be adjusted by the structure parameters of the ellipticity, the major and minor diameter of the microfiber, and the thickness of the graphene film. When the length of the graphene film is 1.5 mm, extinction ratios of ~30 and ~32 dB can be respectively achieved at the wavelength of 1.55 μm for TMpass and TE-pass polarizers. The operation wavelength range of the TM- and TE-pass polarizer could cover ~400 and ~1000 nm, respectively. Our proposed polarizer provides a flexible way to manipulate the polarization of the light in the fiber and highly integrated optical fiber system. © 2017 Optical Society of America OCIS codes: (230.5440) Polarization-sensitive devices; (260.5430) Polarization; (060.2310) Fiber optics; (160.2290) Fiber materials.

References and links 1. 2. 3. 4. 5. 6. 7. 8. 9. 10. 11. 12. 13. 14. 15. 16. 17.

H. Sunnerud, M. Karlsson, C. Xie, and P. A. Andrekson, “Polarization-mode dispersion in high-speed fiber-optic transmission systems,” J. Lightwave Technol. 20(12), 2204–2219 (2002). D. C. Cullen, R. G. Brown, and C. R. Lowe, “Detection of immuno-complex formation via surface plasmon resonance on gold-coated diffraction gratings,” Biosensors 3(4), 211–225 (1987). M. Piliarik, J. Homola, Z. Manikova, and J. Ctyroky, “Surface plasmon resonance sensor based on a single-mode polarization-maintaining optical fiber,” Sens. Actuators B Chem. 90(1-3), 236–242 (2003). H. C. Su and L. A. Wang, “A highly efficient polarized superfluorescent fiber source for fiber-optic gyroscope applications,” IEEE Photonics Technol. Lett. 15(10), 1357–1359 (2003). T. K. Ng, M. Z. M. Khan, A. Al-Jabr, and B. S. Ooi, “Analysis of CMOS compatible Cu-based TM-pass optical polarizer,” IEEE Photonics Technol. Lett. 24(9), 724–726 (2012). Y. Ma, G. Farrel, Y. Semenova, B. Li, J. Yuan, B. Li, J. Yuan, X. Sang, B. Yan, C. Yu, and Q. Wu, “Optical microfiber-loaded surface plasmonic TE-pass polarizer,” Opt. Laser Technol. 78, 110–114 (2015). H. Esmaeilzadeh, E. Arzi, M. Mozafari, and A. Hassani, “A broadband optical fiber based inline polarizer for telecom wavelength range,” Sens. Actuators A Phys. 185, 59–65 (2012). Q. Bao, H. Zhang, B. Wang, Z. Ni, C. H. Y. X. Lim, Y. Wang, D. Y. Tang, and K. P. Loh, “Broadband graphene polarizer,” Nat. Photonics 5(7), 411–415 (2011). J. T. Kim and C. G. Choi, “Graphene-based polymer waveguide polarizer,” Opt. Express 20(4), 3556–3562 (2012). P. Mehta, N. Healy, N. F. Baril, P. J. A. Sazio, J. V. Badding, and A. C. Peacock, “Nonlinear transmission properties of hydrogenated amorphous silicon core optical fibers,” Opt. Express 18(16), 16826–16831 (2010). S. Morris, T. Hawkins, P. Foy, J. Hudson, L. Zhu, R. Stolen, R. Rice, and J. Ballato, “On loss in silicon core optical fiber,” Opt. Mater. Express 2(11), 1511–1519 (2012). D. Graham-Rowe, “Fibres get functional,” Nat. Photonics 5(2), 66–67 (2011). B. Jalali, V. Raghunathan, D. Dimitrophoulos, and O. Boyraz, “Raman-based silicon photonics,” IEEE J. Sel. Top. Quantum Electron. 12(3), 412–421 (2006). J. Ballato, T. Hawkins, P. Foy, B. Yazgan-Kokuoz, C. Mcmillen, L. Burka, S. Morris, R. Stolen, and R. Rice, “Advancements in semiconductor core optical fiber,” Opt. Fiber Technol. 16(6), 399–408 (2010). L. Yang, T. Hu, R. Hao, C. Qiu, C. Xu, H. Yu, Y. Xu, X. Jiang, Y. Li, and J. Yang, “Low-chirp high-extinctionratio modulator based on graphene-silicon waveguide,” Opt. Lett. 38(14), 2512–2515 (2013). G. W. Hanson, “Quasi-transverse electromagnetic modes supported by a graphene parallel-plate waveguide,” J. Appl. Phys. 104(8), 084314 (2008). C. Guan, L. Yuan, F. Tian, and Q. Dai, “Characteristics of near-surface-core optical fibers,” J. Lightwave Technol. 29(19), 3004–3008 (2011).

#284058 Journal © 2017

https://doi.org/10.1364/OME.7.001398 Received 3 Jan 2017; revised 19 Mar 2017; accepted 19 Mar 2017; published 27 Mar 2017

Vol. 7, No. 4 | 1 Apr 2017 | OPTICAL MATERIALS EXPRESS 1399

1. Introduction Optical polarizers are important elements for many optical systems such as fiber-optic networks, optical fiber sensors, biological sensors, optical signal process systems, and so on [1–4], because good line-polarization state lights could provide high signal-to-noise ratio and avoid signal fading. The manipulation of polarization, however, is still a challenge in the integrated optical fiber system. So far, many integratable optical polarizers [5–9] have been proposed with various kind of technology. The common technology is based on the polarization mode selectivity of the surface plasmon resonances (SPRs) [5,6], since the SPRs can only be excited by a transverse electromagnetic field perpendicular to a very thin metal surface with phase matching. However, the working spectral ranges of these optical polarizers are usually limited by their structures [6]. In order to enlarge the working spectral ranges, an in-fiber linear polarizer based on the tilted grating has been performed for TM polarization with 40dB polarization extinction ratio over 50 nm wavelength range [7]. However, its working spectral ranges cannot meet the demand in the optical fiber communication system yet. In 2011, a broadband fiber polarizer was achieved by combing the 3 mm graphene sheet with the side-polished fiber for large working spectral ranges from 800 to 1650 nm [8]. This side-polished approach, however, greatly reduced the strength and durability of the fiber, which would impose limits on the possible applications and service life. Subsequently, J. T. Kim et al., fabricated a graphene polarizer by using the planar-lightwave-circuit (PLC) and graphene [9], where the insertion losses of the TM and TE polarization light for air cladding are 4.14 and 2.18 dB/μm, and TM and TE losses for polymer cladding are 5.4 and 9.4 dB/μm, respectively. Nevertheless, such polarizer cannot integrate with the fiber system. Recently, silica-cladding silicon core optical fibers [10,11] have attracted much interest due to their infrared transparency [11], strong evanescent field, high optical confining factor, and strong optical nonlinearities, which offer great potential for many applications in sensing, biomedicine, high-speed, all-optical signal processing, supercontinuum spectrum generation [12–14]and so on. Especially, the tightly strong evanescent field of surface-core fiber [15] could enable it to combine with graphene sheet for a novel polarizer with easy-integrated, durable, and broadband performances. In this paper, a flexible and broadband graphene polarizer based on the surface silicon-core microfiber has been proposed, in which the surface silicon-core locates tangent to the cladding of the microfiber. The optical broadband polarization properties of such polarizer are investigated in detail. The impacts of the structure parameters such as ellipticity, major and minor diameter of the microfiber, the size of the silicon core, and the thickness of the graphene film on its performances are discussed. 2. Theoretical analysis and design The schematic cross section of our proposed surface silicon-core microfiber (SSCM) polarizer with surface-covered graphene film is shown in Fig. 1. A monolayer or few-layer graphene film covers upon the side of the microfiber, where the cladding of the microfiber and the graphene film locate tangent to the surface silicon-core, as shown in Fig. 1. A tight and strong evanescent field can be achieved by its small surface silicon core, which could increase the interaction between the evanescent light and the graphene film. The asymmetric dielectric structure model consists of air, graphene, silicon core, and silica cladding layers. In Fig. 1, a and b is respectively the major and minor diameter of the microfiber, and ellipticity is defined as e = b/a. The complex effective refractive index neff of such waveguide can be obtained by solving Maxwell’s equations [8]: ∇ × H = − jωε r ε 0 E + σ E

(1)

∇ × E = jωμ0 H

(2)

Vol. 7, No. 4 | 1 Apr 2017 | OPTICAL MATERIALS EXPRESS 1400

By matching the boundary conditions with the Maxwell’s equation, the dispersion relation of the graphene-microfiber waveguide could be achieved [8]: a tan ( Ca ) + a tan ( Cb ) + mπ = γ 2 d

(3)

where m is the order of the mode, and d is the diameter of the silicon core. In Eq. (3) for TM γ ε polarization modes, Ca =  1 2  ε1γ 2

-1

 γ 3ε 2 σγ 1  ; and for TE polarization mode  , Cb =  1 + ε 3γ 2  ωε 0ε1 

 γ + ωμ0σ  γ3 2 2 2 Ca =  1 , where γ 1 = k0 neff -n12 , γ 2 = k0 n22 -neff and γ 3 = k0 neff -ε 3 .  , Cb = γ γ 1 2   where, εj is the relative permittivity of region j. n1 is the refractive index of the silica cladding and n2 is the refractive index of the silicon core.

Fig. 1. The schematic cross section of the polarizer, where r is the radius of the silicon core. a and b is respectively the major and minor diameter of the microfiber. n1 = 1.4378, n2 = 3.4784.

The complex relative permittivity of the graphene ε3 is derived by [16]:

ε3 = 1 −

Im (σ )

ωε 0 d g

+i

Re (σ )

ωε 0 d g

(4)

where dg is the thickness of graphene film, σ is the complex conductivity of graphene and can be expressed by Kubo formula, including interband and intraband contributions [17]:

σ ( ω , μ c , Γ, T ) =

je2 (ω − jτ -1 )

π

 1   (ω − jτ -1 )2 



∞

0

∗  ∂f d ( ε )

ε



∂ε

−

 ∞ ∂f d ( −ε )  f d ( −ε ) − f d ( ε ) d ε − d ε  0 (ω − j 2Γ )2 − 4 (ε /  )2  ∂ε  

(5) where e is the charge of an electron,  = h / 2π is the reduced Planck’s constant, −1 (ε − μ ) k T f d ( ε ) = e c B + 1 is the Fermi-Dirac distribution, and kB is Boltzmann’s constant.

(

)

Thus, the complex effective refractive index neff of the waveguide can be obtained by combining Eq. (3),(4) and (5). The birefringence can be calculated by the real part of the effective refractive index of TE-polarized and TM-polarized modes, that is B = Re ( nTE ) − Re ( nTM ) . The attenuation constant α is related with image of the effective refractive index of such waveguide, which can be expressed as:

α = Im ( neff ) k0

(6)

The extinction ratio (ER) of such graphene-microfiber polarizer can be calculated by multiplying the difference in the attenuation constant between the TE-polarized and TMpolarized modes by the propagation distance:

Vol. 7, No. 4 | 1 Apr 2017 | OPTICAL MATERIALS EXPRESS 1401

E = (αTE − αTM ) Lg

(7)

where Lg is the graphene film length, and αTE and αTM are the losses of TE and TM modes. 3. Polarizer properties discussion (a) 0.10

with ellipse cladding b=2μm, e=0.5 with circle cladding b=4μm, e=1

r =0.152 μm

40 Extinction ratio (dB)

Birefrigence

0.08 0.06 0.04 0.02 0.00 0.0

0.1

0.2

0.3

0.4

with ellipse cladding b=2μm, e=0.5 with circle cladding b=4μm, e=1

(b) 50 30 20 10 0 -10 -20

TE-pass

-30 -40

0.5

r =0.17 μm 0.0

Radius of silicon core (μm)

(c) 0.10

0.04 0.02 0.00 -0.02

Circle cladding 0.0

0.1

0.2

50

0.3

0.4

Radius of silicon core (μm)

0.2

0.3

0.4

0.5

a=b=4μm a=b=2μm

40 Extinction ratio (dB)

Birefrigence

(d)

r =0.152 μm

0.06

0.1

Radius of silicon core (μm)

a=b=4μm a=b=2μm

0.08

TM-pass

r =0.132 μm

0.5

30 r =0.132 μm 20

TM-pass

10 0 -10 -20

r =0.17 μm

-30 -40

0.0

0.1

0.2

0.3

TE-pass 0.4

0.5

Radius of silicon core (μm)

Fig. 2. (a) Birefringence as a function of the radius of the silicon core with different ellipticity. (b) Extinction ratio versus the radius of the silicon core with different ellipticity. (c) Birefringence as a function of the radius of the silicon core with different major and minor diameter. (d) Extinction ratio versus the radius of the silicon core with different major and minor diameter.

In order to simplify the experiment and optimize the design process, it is necessary to carry out simulations of the proposed structures prior to an experimental investigation. Figure 2(a) and 2(c) show the birefringence varies with the change of the silicon core radius for different ellipticity e and different cladding size. With the change of the ellipticity from 0.5 to 1, the radius of the silicon core of ~0.152 μm related with the maximal birefringence Bmax remains unchanged, and the value of the Bmax slightly decreases. It can be seen that the ellipticity of the microfiber and the major and minor diameters of the microfiber has nothing to do with the silicon-core radius of such polarizer with the maximal birefringence, but it could affect the value of Bmax. The extinction ratio of two polarized modes is calculated by Eq. (7) with the graphene film length of 1.5 mm, which is varied with the silicon core under different ellipticity, the major and minor diameter in Fig. 2(b) and 2(d). It can be seen that, when the radius of the silicon core is low to ~0.147 μm, the fundamental TE mode losses is larger than that of fundamental TM mode. However, when the radius of the silicon core is larger than ~0.147 μm, the fundamental TM mode suffers greater loss than the fundamental TE mode. At the silicon core radius of 0.132 μm, a TM-pass polarizer has been obtained with the highest extinction ratio of ~30 dB, while a good TE-pass polarizer with 30 dB ER has been obtained at the silicon core radius of 0.17 μm. The extinction ratio of two orthogonal polarizations is a consequence of the attenuations difference between two orthogonal polarizations of the fundamental mode (TM- polarized and TE-polarized modes). So the TM-pass mode with the lowest loss in the SSCM is generated at the silicon core radius of 0.132 μm, while the TEpass mode with the lowest loss in the SSCM is realized at the silicon core radius of 0.17 μm.

Vol. 7, No. 4 | 1 Apr 2017 | OPTICAL MATERIALS EXPRESS 1402

The extinction ratios of ~30 and ~32 dB can be achieved at the wavelength of 1.55 μm for TM and TE-pass polarizer, respectively. Furthermore, the ellipticity, the major diameter, and the minor diameter of the microfiber cannot affect the polarization status of such polarizer. From Fig. 2(b) and 2(d), it is concluded that the polarization status of such polarizer is just determined by the size of the silicon core. Thus, the above phenomenon confirms that changing the silicon core of the SSCM provides a potential way to manipulate the polarization property.

TE

Graphene film

TE TM

50

(a) TM

|E|^2 (a.u.)

40 30

air

air

20 10

(b)

(c)

0 -1.0

r = 0.132μ m

-0.5

0.0

0.5

1.0

x coordinate axis (μm)

90

TM TE

80

|E|^2 (a.u.)

70 60 50 40 30 20 10 0 0.0

0.4

0.8

1.2

1.6

y coordinate axis (μm)

Fig. 3. The electric field distribution on the cross section of the SSCM polarizer with surfacecovered graphene film for TE and TM mode light at the radius of the silicon core of 0.132 μm ((a) and (b)) and 0.17 μm ((d) and (e)). (c) The amplitudes of the electric field as a function of x coordinate axis along y = 1μm at the silicon core radius of 0.132μm, a = 4 μm and b = 2 μm. (f) The intensity of the electric field as a function of y coordinate axis along x = 0 at the silicon core radius of 0.17 μm, a = 4 μm and b = 2 μm.

In order to explore deeply the reason why the polarization of the SSCM is mainly dependent on the radius of the silicon core, the electric field distributions for TM-pass and TE-pass modes at the chemical potential of 0 meV has been presented, as shown in Fig. 3. Under the silicon core radius of 0.132 μm, the electric field distributions for TE and TM modes at e = 0.5, λ = 1.55 μm are shown in Fig. 3(a) and 3(b), respectively. It can be seen that TE polarization direction is along the graphene film, and the TM polarization direction is perpendicular to the graphene film surface. Figure 3(c) is the intensity of the electric field of the core center along the x-axis at y = 1 μm, in which the Cartesian coordinate origin is at the center of the microfiber cladding. Compared to the TM mode, a large of electric field energy in TE modes is out of the silicon core to come into the graphene film by absorption. Thus, the SSCM with the silicon core radius of 0.132 μm presents as a TM-pass polarizer. Figure 3(d) and 3(e) are respectively the electric field distributions for TE and TM modes at r = 0.17 μm,

Vol. 7, No. 4 | 1 Apr 2017 | OPTICAL MATERIALS EXPRESS 1403

e = 0.5, λ = 1.55 μm. Compared the electric field distribution of the SSCM at r = 0.132 μm (Fig. 3 (a) and 3(b)) with that of r = 0.17 μm (Fig. 3(d) and 3(e)), the larger silicon core could constraint more energy in the core, which could result in the change of the attenuations difference of two orthogonal polarizations modes. Figure 3(f) is the intensity of the electric field of the core center in the SSCM at r = 0.17 μm along the y-axis at x = 0 μm, where the Cartesian coordinate origin is also at the center of the microfiber cladding. There is more energy of the TM mode coupled into the graphene film than that of the TE mode, thus achieving a TE-pass polarizer at r = 0.17 μm. 0.05

0.20

(a) 0.18

0.03

0.16

0.08 0.06

0.02

0.2

0.4

0.6

0.8

0.03

0.00

1.0

0.01

-0.02

Ellipticity

0.02

TE TM r = 0.132μm, e=1 1

0.12

TE 0.06 TM r =0.17μm, a =2μm

0.11

0.05

0.08

0.03

0.07 0.06

Ellipticity

5

6

0.01

0.06

0.055

0.05

0.050 0.04 0.045 0.03

0.040

TE TMr =0.17μm, e=1

0.02 0.6

4

0.8

1.0

0.035

1

2

3

4

5

6

Loss (dB/μm)

0.04

0.09

0.4

3

(d) 0.060

Loss (dB/μm)

0.10

0.2

2

Major diameter of the microfiber (μm)

Birefrigence

(c) 0.13

0.04

0.01

-0.01

TE TM r =0.132μm, a =2μm

0.02

Birefrigence

Birefrigence

0.10

0.02

Loss (dB/μm)

0.03

Loss (dB/μm)

0.12

0.04

Birefrigence

0.05

0.04

0.14

0.05

0.06

(b) 0.04

0.02

Major diameter of the microfiber (μm)

Fig. 4. The effect of the ellipticity, major and minor diameter on the birefringence and propagation loss for TE and TM modes. (a) and (b) for r = 0.132 μm, (c) and (d) for r = 0.17 μm.

To further discuss the impact of the ellipticity, major and minor diameter on TM-pass and TE-pass polarizer, the birefringence and propagation losses of the SSCM as a function of the ellipticity, major and minor diameter has been analyzed in Fig. 4, under keeping r = 0.132 μm and r = 0.17 μm. With the increase of the ellipticity, it is found in Fig. 4(a) and 4(c) that the birefringence at r = 0.132 μm and r = 0.17 μm all rapidly decrease, while the difference of the TE and TM modes increases at r = 0.132 μm while it decreases at r = 0.17 μm. Thus, the ellipticity should be set at 0.2 to obtain a high-ER TE-pass polarizer while it should be set at 1 to obtain a high-ER TM-pass polarizer. When the major and minor diameters increase, the birefringence at r = 0.132 μm and r = 0.17 μm increases from Fig. 4(b) and 4(d). The differences of the TE and TM modes for two polarizers at r = 0.132 μm and r = 0.17 μm both increase with the increase of the major and minor diameters. It can be concluded that the large-sized cladding in the SSCM contributes to improve the ER of the polarizer.

Vol. 7, No. 4 | 1 Apr 2017 | OPTICAL MATERIALS EXPRESS 1404

0.04

(a) 0.0383 0.0382

0.05

0.0588

0.04

0.0378 0.01

0.0377

0.0375 0.0

0.5

1.0

1.5

2.0

0.0584

0.03

0.0582 0.02

0.0580 0.0578

2.5

0.00 3.0

0.0574 0.0

0.01

TE TM r = 0.17 μm

0.0576

TE TM r =0.132 μm

0.0376

0.0586

0.5

Graphen thickness (nm)

1.0

1.5

2.0

Loss (dB/μm)

0.02

0.0379

Loss (dB/μm)

0.0380

Birefrigence

0.03

0.0381 Birefrigence

(b) 0.0592 0.0590

2.5

0.00 3.0

Graphen thickness (nm)

Fig. 5. The effect of the graphene thickness on the birefringence and propagation loss for TE and TM modes. (a) r = 0.132 μm and (b) r = 0.17 μm.

The impact of the graphene thickness on the birefringence and propagation loss is shown in Fig. 5. When the graphene layer is varied from 1 to 7, the propagation losses of the TM and TE modes increases, and on the contrary birefringence decreases. From Fig. 5, it can be observed that, in order to obtain the same high ER on polarization, the length of monolayer graphene film should be longer along the propagation direction than that of the multi-layer graphene. To be concluded, increasing the graphene layer is helpful for decreasing the length of the graphene film. The large optical absorption bandwidth is an important property for the polarizer. For analyzing the bandwidth of such SSCM polarizer, the propagation losses of the TM and TE modes have been calculated with the variation of the wavelength, as shown in Fig. 6. When the silicon core radius is 0.132 μm, the losses between the TE and TM modes dramatically increase and the difference of two polarized modes obviously decreases by varying the wavelength from 1400 to 1800 nm. On the contrary, in Fig. 6(b), by increasing the wavelength from 800 to 1750 nm, the losses of the TE and TM modes at r = 0.17 μm remarkably increase and the difference of two polarized modes also increases. For the TMpass polarizer at r = 0.132 μm, the polarization bandwidth is about 400 nm, while the TE polarization mode wave propagation bandwidth is close to ~1000 nm when the silicon core radius is set at 0.17 μm. 0.06

(b)

0.02 TE TM 0.00

0.05 0.04

0.04

Loss (dB/μm)

Loss (dB/μm)

(a)

@r=0.17μm

0.03 0.02 0.01

@r=0.132μm 1400

TE TM

1500

1600

1700

Wavelength (nm)

1800

0.00 600

800

1000

1200

1400

1600

1800

Wavelength (nm)

Fig. 6. Propagation losses of TE and TM modes against the wavelength (a) r = 0.132 μm and (b) r = 0.17 μm.

4. Summary

We proposed a graphene-coated SSCM polarizer and numerically analyzed the birefringence and losses of two polarization modes in detail. The birefringence of the SSCM has no impact on the light polarization. By adjusting the radius of the silicon core from 0.132 to 0.17 μm, we can achieve TM-pass and TE-pass polarizers. The extinction ratios of the TM-pass polarizer can be increased by enlarging the ellipticity, while the extinction ratios of the TE-pass

Vol. 7, No. 4 | 1 Apr 2017 | OPTICAL MATERIALS EXPRESS 1405

polarizer can be increased by reducing the ellipticity. The extinction ratios of the TM-pass and TE-pass polarizer can be improved by increasing the size of the cladding in the SSCM. Multi-layer the graphene is helpful for decreasing the used length of graphene film covered on the SSCM. The extinction ratios of ~30 and ~32 dB can be achieved at the wavelength of 1550 nm for TM- and TE-pass polarizers respectively, when the length of the graphene film is 1.5 mm. The operation wavelength range of the TM-pass and TE-pass polarizer could cover ~400 and ~1000 nm, respectively. Such graphene-coated SSCM polarizer can be manipulated by the physical changes on the surface silicon core without destroying the fiber structure, and therefore this polarized component has advantages of simple fabrication, easy-integration, durability, broadband and low cost, and also provide a flexible way to manipulate the polarization of the light in fiber and integrated optical fiber system. Funding

National Natural Science Foundation of China (NSFC) (No.61675046, 61205132); a Doctoral Program of Higher Education of China (20120071120023).