Influence of cerium ions on thermal bleaching of photo-darkened ytterbium-doped fibers

Abstract

We studied the influence of cerium (Ce) ions on photo-darkening (PD) behavior in ytterbium/aluminum (Yb/Al) co-doped silica fibers at room and elevated temperatures. Low levels of PD was observed for Ce co-doped Yb/Al fiber. And the Yb/Al co-doped fiber was completely bleached at ~600°C. The addition of Ce ions as co-dopant can significantly lower the initial recovery and complete bleaching temperatures. The complete recovery temperature is ~ 450°C and ~400°C for Yb/Al/Ce low fiber and Yb/Al/Ce high fiber respectively. More importantly, Ce co-doping in Yb/Al fiber also decreases the heat-induced loss.

Keywords

photo-darkening (PD) / ytterbium/aluminum (Yb/Al) fiber / cerium (Ce) ions / thermal bleaching

1 Introduction

In prospective optical networks, different modulation formats may be selectively employed depending on network size and system settings [1]. Thus, optical format conversion is an essential technology, to which has been paid much attention in recent years. However, most previous works have focused on on-off keying (OOK) formats, such as conversions from non-return-to-zero OOK (NRZ-OOK) to return-to-zero OOK (RZ-OOK) formats [2-4]. Recent studies have presented that differential phase-shift keying (DPSK) format particularly exhibits better performance than that of OOK for long-haul transmission [5]. But only a few studies have focused on such NRZ to RZ format conversion in DPSK based on cascaded second-order nonlinearities in a periodically poled lithium niobate (PPLN) waveguide [6].

Compared with conventional OOK formats, phase modulated formats are suitable for long-haul transmission. More specifically, return-to-zero differential phase-shift keying (RZ-DPSK) formats have narrower pulse and thus have a potential for multiplexing in time domain, i.e., phase-shift keying (PSK) optical time divided multiplexed (OTDM) system. In this paper, we propose a format conversion scheme from non-return-to-zero differential phase-shift keying (NRZ-DPSK) to RZ-DPSK by using phase modulators and detuning filters. Generally, in conversional NRZ-RZ conversion for OOK formats, format conversion usually requires detuning filter [7], but the application of detuning filter in DPSK system may cause demodulation of the signal [8]. But there is no report that detuning filters used in demodulated alternate mark inversion (AMI) and duobinary (DB) of the DPSK signal. By the analysis of the instinct relationships among AMI, DB and DPSK signals, numerical and simulation results of the format conversion with the help of the VPI Transmission Maker 8.5 have been showed in this paper, and the paper is organized as below: In Sect. 2, the principle of the format conversion from NRZ-DPSK to RZ-DPSK is demonstrated by exploiting the instinct relations among AMI, DB and DPSK signals. In Sect. 3, the numerical model is provided and the simulation result is shown with the help of VPI. The scheme described in this study does not require high-order nonlinearity such as four waves mixing (FWM) or cascaded second-harmonic generation (cSHG), so the high efficiency conversion may be potentially available in next generation optical network. And the scheme may be potential in practical use by mature electrical clock recovery avoiding optical-electrical-optical (O-E-O) operation.

2 Principle of conversion from NRZ-DPSK to RZ-DPSK

The instinct relations of AMI, DB and DPSK are shown in Fig. 1. DPSK signal can be demodulated to AMI and DB signals by using demodulator such as 1 bit delay interferometer [9]. And such AMI and DB signals can also build up original DPSK signal by another 1 bit delay interferometer, as long as phase information of the AMI and DB signals remain as same as the original DPSK signal [10]. If the phase information of the AMI and DB signals has been changed by wavelength conversion scheme in semiconductor optical amplifier (SOA), the binary phase-shift keying (BPSK) signal will be built instead of DPSK signal [11]. The principle of the conversion from NRZ-DPSK to RZ-DPSK is indicated in Fig. 2. Normally, NRZ-DPSK signal, which is modulated by Mach-Zehnder modulator (MZM), has intensity dips due to limited bandwidth and the voltage of drive signal. These intensity dips may be cut out in two different ways: one is by RZ pulse carving in the transmitter, the other is by high-order nonlinearity such as FWM and cSHG in the middle of the transmission spans. But our scheme is to do the RZ pulse carving “operation” in the demodulated products of the NRZ-DPSK signal. The demodulated products of NRZ-DPSK signal are NRZ-DB signal and RZ-AMI signal. If the “operation” may preserve the phase information of the original signals, the converted RZ-DB and RZ-AMI signals may build up a converted RZ-DPSK signal.

Fig.1 Relations among DPSK, AMI, DB and BPSK formats

Full size|PPT slide

Fig.2 Principle of format conversion from NRZ-DPSK to RZ-DPSK

Full size|PPT slide

The proposed setup of the format conversion from NRZ-DPSK to RZ-DPSK is shown in Fig. 3. The 20 Gbit/s NRZ-DPSK signals at 193.1 THz are modulated at the MZM based transmitter. 10% of the NRZ-DPSK signals are used as electrical clock recovery (CR) by a 10∶90 tap. Remaining 90% are demodulated by a 1 bit delay interferometer. The demodulated products are transmitted to two arms of Mach-Zehnder interferometer (MZI), and there is a phase modulator and a detuning filter in each arm. Both the same phase modulators (PMs) are driven by an electrical clock which is recovered from the CR. The filters are first-order Gauss shape with a bandwidth of 30 GHz, the center frequency of the filters is 193.180 THz, which is 80 GHz detuned from the frequency of the signals. After the MZI, a 1 bit delay interferometer is used to build up a new RZ-DPSK signal.

Fig.3 Proposed setup of format conversion from NRZ-DPSK to RZ-DPSK

Full size|PPT slide

3 Theoretical model and simulation result

The NRZ-DPSK signal field $E_{\mathrm{s}}$ can be expressed as

(1)$$E_{\mathrm{s}}=A_{\mathrm{s}}\exp (\mathrm{i}({\omega }_{\mathrm{0}}t-{\varphi }_{\mathrm{s}})),$$

where $A_{\mathrm{s}}$ is the complex amplitude of the signal, ${\omega }_{\mathrm{0}}$ is the center frequency, and ${\varphi }_{\mathrm{s}}$ is the phase information of the signal. After 1 bit delay interferometer (DI), the demodulated products can be expressed as [12]

(2)$${\mathit{E}}_{\mathrm{D}\mathrm{B}}=A_{\mathrm{s}}\exp (\mathrm{i}{\omega }_{\mathrm{0}}t)\exp \left(\mathrm{i}{\displaystyle \frac{{\varphi }_{\mathrm{s}\mathrm{1}}+{\varphi }_{\mathrm{s}\mathrm{2}}}{\mathrm{2}}}\right)\times \cos {\displaystyle \frac{{\varphi }_{\mathrm{s}\mathrm{1}}-{\varphi }_{\mathrm{s}\mathrm{2}}}{\mathrm{2}}},$$

(3)$$E_{\mathrm{A}\mathrm{M}\mathrm{I}}=A_{\mathrm{s}}\exp (\mathrm{i}{\omega }_{\mathrm{0}}t)\exp \left(\mathrm{i}{\displaystyle \frac{{\varphi }_{\mathrm{s}\mathrm{1}}+{\varphi }_{\mathrm{s}\mathrm{2}}}{\mathrm{2}}}\right)\times \sin {\displaystyle \frac{{\varphi }_{\mathrm{s}\mathrm{1}}-{\varphi }_{\mathrm{s}\mathrm{2}}}{\mathrm{2}}}.$$

$E_{\mathrm{D}\mathrm{B}},E_{\mathrm{A}\mathrm{M}\mathrm{I}}$ are the electrical fields of DB and AMI signals, respectively, ${\varphi }_{\mathrm{s}\mathrm{1}}(t),{\varphi }_{\mathrm{s}\mathrm{2}}(t)$ are the phases of the adjacent bit sequence in the two arms of the DI. The AMI and DB signals may build up original DPSK signals, because the phase information term $\pm \exp \left(\mathrm{i}{\displaystyle \frac{{\varphi }_{\mathrm{s}\mathrm{1}}(t)+{\varphi }_{\mathrm{s}\mathrm{2}}(t)}{\mathrm{2}}}\right)$ is maintained during the interference. On the other hand, by wavelength conversion scheme, only PSK signal will be built up, because only the amplitude term can be maintained, and the phase information can not.

After phase modulators, the DB and AMI signals may have a phase shift, which is related to the recovered clock of the NRZ-DPSK signal. The phase modulated DB and AMI signals $E_{\mathrm{D}\mathrm{B}}^{\prime },E_{\mathrm{A}\mathrm{M}\mathrm{I}}^{\prime }$ can be expressed as

(4)$$E_{\mathrm{D}\mathrm{B}}^{\prime }=A_{\mathrm{s}}\exp \left(\mathrm{i}{\displaystyle \frac{{\varphi }_{\mathrm{s}\mathrm{1}}+{\varphi }_{\mathrm{s}\mathrm{2}}}{\mathrm{2}}}\right)\times \cos {\displaystyle \frac{{\varphi }_{\mathrm{s}\mathrm{1}}-{\varphi }_{\mathrm{s}\mathrm{2}}}{\mathrm{2}}}\times \exp (\mathrm{i}\theta (t)),$$

(5)$$E_{\mathrm{A}\mathrm{M}\mathrm{I}}^{\prime }=A_{\mathrm{s}}\exp \left(\mathrm{i}{\displaystyle \frac{{\varphi }_{\mathrm{s}\mathrm{1}}+{\varphi }_{\mathrm{s}\mathrm{2}}}{\mathrm{2}}}\right)\times \sin {\displaystyle \frac{{\varphi }_{\mathrm{s}\mathrm{1}}-{\varphi }_{\mathrm{s}\mathrm{2}}}{\mathrm{2}}}\times \exp (\mathrm{i}\theta (t)),$$

where $\theta (t)$ is the phase shift related to the recovery clock of the NRZ-DPSK signal.

The frequency response $H(\omega )$ of the Gauss shaped filter can be expressed as

(6)$$H(\omega )=\exp \left[-{\displaystyle \frac{\mathrm{1}}{\mathrm{2}}}{\left({\displaystyle \frac{\omega -{\omega }_f}{B_{\mathrm{0}}}}\right)}^{\mathrm{2}}\right],$$

where ${\omega }_f$ is the center frequency and $B_{\mathrm{0}}$ is 3-dB bandwidth of the filter, and the related time domain response can be expressed as

(7)$$h(t)={\displaystyle \frac{B_{\mathrm{0}}}{\sqrt{\mathrm{2}\mathrm{\pi }}}}\exp \left[-{\displaystyle \frac{\mathrm{1}}{\mathrm{2}}}{(B_{\mathrm{0}}t)}^{\mathrm{2}}\right]\exp (\mathrm{i}{\omega }_{f}t).$$

The output signal fields $E_{\mathrm{D}\mathrm{B}\mathrm{o}\mathrm{u}\mathrm{t}}$ and $E_{\mathrm{A}\mathrm{M}\mathrm{I}\mathrm{o}\mathrm{u}\mathrm{t}}$ after the filter are the convolutions of the input signal field and the time domain response of the filter, which can be expressed as

(8)$$E_{\mathrm{D}\mathrm{B}\mathrm{o}\mathrm{u}\mathrm{t}}=E_{\mathrm{D}\mathrm{B}}^{\prime }\times h(t)={\displaystyle \int E_{\mathrm{D}\mathrm{B}}^{\prime }(\tau )h(t-\tau )}\mathrm{d}\tau ,$$

(9)$$E_{\mathrm{A}\mathrm{M}\mathrm{I}\mathrm{o}\mathrm{u}\mathrm{t}}=E_{\mathrm{A}\mathrm{M}\mathrm{I}}^{\prime }\times h(t)={\displaystyle \int E_{\mathrm{A}\mathrm{M}\mathrm{I}}^{\prime }(\tau )h(t-\tau )}\mathrm{d}\tau .$$

By first-order Taylor expansion, Eqs. (8) and (9) can be further expressed as [13]

(10)$$E_{\mathrm{D}\mathrm{B}\mathrm{o}\mathrm{u}\mathrm{t}}=E_{\mathrm{D}\mathrm{B}}^{\prime }\exp \left[-{\displaystyle \frac{\mathrm{1}}{\mathrm{2}}}{\left({\displaystyle \frac{{\omega }_f-{\omega }_{\mathrm{0}}+{\displaystyle \frac{\mathrm{d}\theta (t)}{\mathrm{d}t}}}{B_{\mathrm{0}}}}\right)}^{\mathrm{2}}\right],$$

(11)$$E_{\mathrm{A}\mathrm{M}\mathrm{I}\mathrm{o}\mathrm{u}\mathrm{t}}=E_{\mathrm{A}\mathrm{M}\mathrm{I}}^{\prime }\exp \left[-{\displaystyle \frac{\mathrm{1}}{\mathrm{2}}}{\left({\displaystyle \frac{{\omega }_f-{\omega }_{\mathrm{0}}+{\displaystyle \frac{\mathrm{d}\theta (t)}{\mathrm{d}t}}}{B_{\mathrm{0}}}}\right)}^{\mathrm{2}}\right].$$

From Eqs. (10) and (11), it can be seen that, first, the phase information is maintained, so $E_{\mathrm{D}\mathrm{B}\mathrm{o}\mathrm{u}\mathrm{t}}$ and $E_{\mathrm{A}\mathrm{M}\mathrm{I}\mathrm{o}\mathrm{u}\mathrm{t}}$ may build up the original DPSK signal. Second, new DB and AMI signals have a new amplitude term which comes from the phase modulators and filters. The phase modulators give the DB and AMI signals an RZ shaped clock related phase shift, and detuning filters acts as a phase shift to amplitude conversion. Because the phase shift is related to the clock, so the converted amplitude is also clock related, and this results in an NRZ-DPSK to RZ-DPSK format conversion.

It should be noticed that, using detuning filter in DPSK signal may cause demodulation while in DB and AMI signals may not. The reasons can be found in Eqs. (1)-(3). Equation (1) is the field of the DPSK signal, once there is phase change in the DPSK signal, the detuning filter will convert this phase shift to amplitude output. But in DB and AMI signals, once there is a phase change, unlike the constant amplitude of DPSK, the amplitude of the DB and AMI signals will also change due to the cosine and sine terms in Eqs. (2) and (3). Once there is a phase shift in the AMI and DB signals, the detuning filter still convert this phase shift to amplitude output, but in the mean time, the cosine or sine term is zero, the total output amplitude is zero. So using detuning filter will not cause the demodulation of the DB and AMI signals.

Numerical analysis is used to explain the theory of the scheme, so we neglect high-order term in the Taylor expansion. Next, more exactly simulation results have been shown by using commercial optical design software VPI Transmission Maker 8.5. In the simulation, no such approximation mentioned above is used.

In the following simulation, the input signal is considered as 20 Gbit/s 2³¹-1 pseudo-random binary sequence (PRBS) NRZ-DPSK signal, the recovered clock has sine pulse type, the phase modulators are just get a 180° phase change when driven by the clock signal, sample rate is set to 640 GHz to achieve high accuracy. The power, center frequency, and the time widows of the signal are set to 1 dBm, 193.1 THz and 6.4 s, respectively.

The simulation results shown in Fig. 4(a1) are the spectra of the original NRZ-DPSK signal, and Fig. 4(a2) is the spectra of the converted RZ-DPSK signal. The main band of the converted RZ-DPSK signals is smooth and broader than NRZ-DPSK signal. And the conversion efficiency is about -10 dB, which is more effective than 20 dB efficiency in FWM. Figures 4(b1) and 4(b2) are the temporal waveforms of the NRZ-DPSK signal and the demodulated products of the NRZ-DPSK signal, respectively, Figs. 4(b3) and 4(b4) are the waveforms of the electrical clock and the converted RZ-DPSK signal, respectively, Fig. 4(b5) is the waveform of the demodulated products of the RZ-DPSK signal. It can be seen from Figs. 4(b1) to 4(b5) that the format conversion from NRZ-DPSK to RZ-DPSK is obtained and the demodulated products of NRZ-DPSK and RZ-DPSK are identical. The bit pattern in Fig. 4(b2) is “10111101100111010101” and so is the bit pattern in Fig. 4(b5), Fig. 4(b2) is the NRZ form and Fig. 4(b5) is the RZ form. So the phase information is preserved. It should be noticed that the converted RZ-DPSK signal in Fig. 4(b5) is 1 bit duration later than the NRZ-DPSK signal in Fig. 4(b2). It is because every time a “delay and add” or “delay and subtract” operation happen, the first and the last of the bit sequence are lost, the first bit and the last bit do not have a “delay bit” to add or subtract. But it will not affect transmission. Figures 4(c1)-4(c4) are the eye diagrams of NRZ-DPSK signal, demodulated NRZ-DPSK signal, RZ-DPSK, and demodulated RZ-DPSK signal, respectively. According to our simulation, a clearly eye open can be expected in the format conversion from NRZ-DPSK to RZ-DPSK signal.

Fig.4 Experimental results of 20 Gbit/s NRZ-DPSK to RZ-DPSK. (a1) and (a2) Optical spectra of NRZ-DPSK and converted RZ-DPSK signal; (b1)-(b5) temporal waveforms of NRZ-DPSK, demodulated NRZ-DPSK, clock signal, RZ-DPSK, demodulated RZ-DPSK, respectively; (c1)-(c4) eye diagrams of NRZ-DPSK, demodulated NRZ-DPSK, RZ-DPSK and demodulated RZ-DPSK, respectively

Full size|PPT slide

In the conversion, high-order nonlinearity such as FWM is avoided. Although the detuning filter may induce some power loss, which is also an issue in FWM scheme. Our scheme does not suffer from the limited condition which is very common in high-order nonlinearity scheme.

In our scheme, the electrical clock is used because the electrical clock recovery is more mature and preferable in real practice, and there is no optical-electrical-optical (O-E-O) conversion in our scheme, the signal’s transmission paths are all optical paths. So the scheme is considered as all optical format conversion. On the other hand, if an optical clock is preferred, nonlinearity device such as semiconductor optical amplifier (SOA) can be used to replace the phase modulator, but the amplified spontaneous emission (ASE) and self phase modulation (SPM) in the SOA may degenerate the conversion performance.

4 Conclusion

In conclusion, a new kind of all-optical format conversion from NRZ-DPSK to RZ-DPSK signal is proposed and simulated by using phase modulators and detuning filters. The optical spectra, temporal waveforms and eye diagrams are studied in detail to confirm the successful realization of NRZ-DPSK to RZ-DPSK format conversion. The use of mature electrical clock recovery avoiding O-E-O operation may be potential in practical use.

References

[1] Richardson D J, Nilsson J, Clarkson W A. High power fiber lasers: current status and future perspectives. Journal of the Optical Society of America B, Optical Physics, 2010, 27(11): B63

[2] Paschotta R, Nilsson J, Barber P R, Caplen J E, Tropper A C, Hanna D C. Lifetime quenching in Yb-doped fibres. Optics Communications, 1997, 136(5–6): 375–378

[3] Yoo S, Basu C, Boyland A J, Sones C, Nilsson J, Sahu J K, Payne D. Photodarkening in Yb-doped aluminosilicate fibers induced by 488 nm irradiation. Optics Letters, 2007, 32(12): 1626–1628

[4] Engholm M, Norin L, Åberg D. Strong UV absorption and visible luminescence in ytterbium-doped aluminosilicate glass under UV excitation. Optics Letters, 2007, 32(22): 3352–3354

[5] Peretti R, Jurdyc A M, Jacquier B, Gonnet C, Pastouret A, Burov E, Cavani O. How do traces of thulium can explain photodarkening in Yb doped fibers? Optics Express, 2010, 18: 20455–20460

[6] Manek-Hönninger I, Boullet J, Cardinal T, Guillen F, Ermeneux S, Podgorski M, Bello Doua R, Salin F. Photodarkening and photobleaching of an ytterbium-doped silica double-clad LMA fiber. Optics Express, 2007, 15(4): 1606–1611

[7] Piccoli R, Gebavi H, Lablonde L, Cadier B, Robin T, Monteville A, Goffic O L, Landais D, Méchin D, Milanese D, Brand T, Taccheo S. Evidence of photodarkening mitigation in Yb-doped fiber lasers by low power 405 nm radiation. IEEE Photonics Technology Letters, 2014, 26(1): 50–53

[8] Guzman Chávez A D, Kir’yanov A V, Barmenkov Y O, Il'ichev N N. Reversible photo-darkening and resonant photo-bleaching of ytterbium-doped silica fiber at in-core 977-nm and 543-nm irradiation. Laser Physics Letters, 2007, 4(10): 734

[9] Piccoli R, Robin T, Brand T, Klotzbach U, Taccheo S. Effective photodarkening suppression in Yb-doped fiber lasers by visible light injection. Optics Express, 2014, 22(7): 7638–7643

[10] Gebavi H, Taccheo S, Tregoat D, Monteville A, Robin T. Photobleaching of photodarkening in ytterbium doped aluminosilicate fibers with 633 nm irradiation. Optical Materials Express, 2012, 2(9): 1286

[11] Zhao N, Xing Y B, Li J M, Liao L, Wang Y B, Peng J G, Yang L Y, Dai N L, Li H Q, Li J Y. 793 nm pump induced photo-bleaching of photo-darkened Yb³⁺-doped fibers. Optics Express, 2015, 23(19): 25272–25278

[12] Jasapara J, Andrejco M, DiGiovanni D, Windeler R.Effect of heat and H₂gas on the photo-darkening of Yb⁺³fibers. In: Proceedings of Conference on Lasers and Electro-Optics/Quantum Electronics and Laser Science Conference and Photonic Applications Systems Technologies. Long Beach, CA, USA: IEEE, 2006, CTuQ5

[13] Leich M, Jetschke S, Unger S, Kirchhof J. Temperature influence on the photodarkening kinetics in Yb-doped silica fibers. Journal of the Optical Society of America. B, Optical Physics, 2011, 28(1): 65

[14] Leich M, Röpke U, Jetschke S, Unger S, Reichel V, Kirchhof J. Non-isothermal bleaching of photodarkened Yb-doped fibers. Optics Express, 2009, 17(15): 12588–12593

[15] Engholm M, Jelger P, Laurell F, Norin L. Improved photodarkening resistivity in ytterbium-doped fiber lasers by cerium codoping. Optics Letters, 2009, 34(8): 1285–1287

[16] Sakaguchi Y, Fujimoto Y, Masuda M, Miyanaga N, Nakano H. Suppression of photo-darkening effect in Yb-doped silica glass fiber by co-doping of group 2 element. Journal of Non-Crystalline Solids, 2016, 440: 85–89

[17] Zhao N, Wang Y B, Li J M, Liu C P, Peng J G, Li H q, Dai N L, Yang L Y, Li J Y. Investigation of cerium influence on photo-darkening and photo-bleaching in Yb-doped fibers. Applied Physics A, 2016, 122(2): 75 :

[18] Koponen J, Laurila M, Söderlund M, Ponsoda J J M I, Iho A. Benchmarking and measuring photodarkening in Yb doped fibers. In: Proceedings of SPIE 7195, Fiber Lasers VI: Technology, Systems, and Applications. San Jose, California, USA: SPIE, 2009, 71950R

[19] Chen G, Xie L, Wang Y B, Zhao N, Li H Q, Jiang Z W, Peng J G, Yang L Y, Dai N L, Li J Y. Photodarkening-induced absorption and fluorescence changes in Yb fibers. Chinese Physics Letters, 2013, 30(10): 104208

[20] Cicconi M R, Neuville D R, Blanc W, Lupi J F, Vermillac M, de Ligny D. Cerium/aluminum correlation in aluminosilicate glasses and optical silica fiber preforms. Journal of Non-Crystalline Solids, 2017, 475: 85–95

[21] Lupi J F, Vermillac M, Blanc W, Mady F, Benabdesselam M, Dussardier B, Neuville D R. Steady photodarkening of thulium alumino-silicate fibers pumped at 1.07 μm: quantitative effect of lanthanum, cerium, and thulium. Optics Express, 2017, 41(12): 2771

[22] Jetschke S, Unger S, Röpke U, Kirchhof J. Photodarkening in Yb doped fibers: experimental evidence of equilibrium states depending on the pump power. Optics Express, 2007, 15(22): 14838–14843

[23] Peng X, Dong L. Temperature dependence of ytterbium-doped fiber amplifiers. Journal of the Optical Society of America B, Optical Physics, 2008, 25(1): 126

[24] Söderlund M J, Montiel i Ponsoda J J, Koplow J P, Honkanen S. Heat-induced darkening and spectral broadening in photodarkened ytterbium-doped fiber under thermal cycling. Optics Express, 2009, 17(12): 9940–9946

[25] Sylvia J, Sonja U, Anka S, Martin L, Matthias J. Role of Ce in Yb/Al laser fibers: prevention of photodarkening and thermal effects. Optics Express, 2016, 24(12): 13009

Acknowledgements

This work was financially supported by the National Natural Science Foundation of China (Grant No. 61735007), and The Youth Fund of Wuhan Donghu University (No. 2017dhzk009).

RIGHTS & PERMISSIONS

2018 Higher Education Press and Springer-Verlag GmbH Germany, part of Springer Nature