Search

Article

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

Spatiotemporal evolution and spectral character of second harmonic generation in optical microresonator

Xu Xin Jin Xue-Ying Hu Xiao-Hong Huang Xin-Ning

Citation:

Spatiotemporal evolution and spectral character of second harmonic generation in optical microresonator

Xu Xin, Jin Xue-Ying, Hu Xiao-Hong, Huang Xin-Ning
PDF
HTML
Get Citation
  • With the consideration of the second and the third order nonlinear effect, the Lugiato-Lefeve equation which describes the field evolution of the fundamental frequency wave and the second harmonic wave is introduced. Based on the Lugiato-Lefeve equation, the generation of the second harmonic wave in the SiN microresonator is analyzed, and the effect of the each parameter on the dual field is studied. Simulation results indicate that the stable field of the fundamental frequency wave is of flat top pulse, and the field of the second harmonic wave is of sinusoidal distribution. When the detuning parameter increases, the power of the dual wave inside the microresonator oscillates, and the stable power weakens, the stable light field is periodically varied. Moreover, the chaos emerges as detuning parameter becomes large. The stable field can be generated in the microresonator with the weak pump power. However, because of the high pump power, the dispersion and nonlinear effect are enhanced, resulting in the periodic light field. Furthermore, the oscillation of the dual power curve is aggravated, as the pump power increases. In addition, the turning patterns can be observed by choosing the special dimension of microresonator. Theoretical analysis results are significant for studying the generation of the second harmonic wave in the microresonator.
      Corresponding author: Jin Xue-Ying, xyjin007@hfut.edu.cn
    • Funds: Project supported by the Young Scientists Fund of the National Natural Science Foundation of China (Grant No. 51705121) and the Open Research Fund of the State Key Laboratory of Transient Optics and Photonics, Chinese Academy of Sciences (Grant No. SKLST201612)
    [1]

    邢书剑, 张福民, 曹士英, 王高文, 曲兴华 2013 物理学报 62 170603Google Scholar

    Xing S J, Zhang F M, Cao S Y, Wang G W, Qu X H 2013 Acta Phys. Sin. 62 170603Google Scholar

    [2]

    Niering M, Holzwarth R, Reichert J, Pokasov P, Udem T, Weitz M, Hansch T W, Lemonde P, Santarelli G, Abgrall M, Laurent P, Salomon C, Clairon A 2000 Phys. Rev. Lett. 84 5496Google Scholar

    [3]

    孟飞, 曹士英, 蔡岳, 王贵重, 曹建平, 李天初, 方占军 2011 物理学报 60 100601Google Scholar

    Meng F, Cao S Y, Cai Y, Wang G Z, Cao J P, Li T C, Fang Z J 2011 Acta Phys. Sin. 60 100601Google Scholar

    [4]

    Swann W C, McFerran J J, Coddington I, Newbury N R, Hartl I, Fermann M E, Westbrook P S, Nicholson J W, Feder K S, Langrock C, Fejer M M 2006 Opt. Lett. 31 3074

    [5]

    Washburn B R, Fox R W, Newbury N R, Nicholson J W, Feder K, Westbrook P S, Jorgensen C G 2004 Opt. Express 12 4999Google Scholar

    [6]

    韩海年, 魏志义, 赵刚 2012 物理 41 249

    Han H N, Wei Z Y, Zhao G 2012 Physics 41 249

    [7]

    Herr T, Brasch V, Jost J D, Wang C Y, Kondratiev N M, Gorodetsky M L, Kippenberg T J 2013 Nat. Photon. 8 145

    [8]

    Song Q H 2019 Sci. China-Phys. Mech. Astron. 62 074231Google Scholar

    [9]

    Zhang X Y, Cao Q T, Wang Z, Liu Y X, Qiu C W, Yang L, Gong Q H, Xiao Y F 2019 Nat. Photon. 13 21Google Scholar

    [10]

    Godey C, Balakireva I V, Aurélien C, Chembo Y K 2014 Phys. Rev. A 89 722

    [11]

    Huang S W, Zhou H, Yang J, McMillan J F, Matsko A, Yu M, Kwong D L, Maleki L, Wong C W 2014 Phys. Rev. Lett. 114 053901

    [12]

    Wang W Q, Lu Z Z, Zhang W F, Chu Sai T, Little B E, Wang L R, Xie X P, Liu M L, Yang Q H, Wang L, Zhao J G, Wang G X, Sun Q B, Liu Y S, Wang Y S, Zhao W 2018 Opt. Lett. 43 2002Google Scholar

    [13]

    Xue X X, Xuan Y, Liu Y, Wang P H, Chen S, Wang J, Leaird D E, Qi M H, Weiner A M 2015 Nat. Photon. 9 594Google Scholar

    [14]

    Coillet A, Balakireva I, Henriet R, Saleh K, Larger L, Dudley J M, Menyuk C R, Chembo Y K 2013 IEEE Photon. J. 5 6100409Google Scholar

    [15]

    Matsko A B, Liang W, Savchenkov A A, Maleki L 2013 Opt. Lett. 38 525Google Scholar

    [16]

    Coen S, Randle H G, Sylvestre T, Erkintalo M 2013 Opt. Lett. 38 37Google Scholar

    [17]

    Xue X X, Leo F, Xuan Y, Villegas J A J, Wang P H, Leaird D E, Erkintalo M, Qi M H, Weiner Andrew M 2017 Light: Sci. Appl. 6 e16253Google Scholar

    [18]

    Lin J T, Xu Y X, Ni J L, Wang M, Fang Z W, Qiao L L, Fang W, Cheng Y 2016 Phys. Rev. Appl. 6 014002Google Scholar

    [19]

    Chembo Y K, Menyuk C R, 2014 Phys. Rev. A 87 053852

    [20]

    Lin J T, Xu Y X, Fang Z W, Wang M, Song J X, Wang N W, Qiao L L, Fang W, Cheng Y 2015 Sci. China-Phys. Mech. Astron. 58 114209Google Scholar

    [21]

    Coen S, Haelterman M 2001 Opt. Lett. 26 39Google Scholar

    [22]

    And T, Boccaletti A, Grebogi C, Lai Y C, Mancini H, Maza D 2000 Phys. Rep. 329 2000

    [23]

    Akhmediev N, Pelinovsky E 2010 Eur. Phys. J. Spec. Top. 185 1Google Scholar

  • 图 1  SiN微腔结构示意图

    Figure 1.  Structure of the SiN microresonator.

    图 2  SiN光学微腔中二次谐波的产生 (δ1 = 0, Ein = 1 W1/2) (a) 倍频光场在腔内的演化过程; (b) 基频光功率随光在腔内循环次数的变化曲线; (c) 倍频光功率随光在腔内循环次数的变化曲线; (d)基频场的稳定分布; (e)倍频光场的稳定分布; (f)基频光的光谱; (g) 基频光的光谱

    Figure 2.  Second harmonic waves generates inside the SiN microresonator (δ1 = 0, Ein = 1 W1/2): (a) Evolution of the second harmonic waves; (b) curves of the fundamental frequency waves power and (c) the second frequency waves power vary with the round trip number; (d) stationary distribution of the fundamental frequency waves light field and (e) the second frequency waves light field; (f) spectra of the fundamental frequency waves and (g) the second frequency waves.

    图 3  发生频率失谐后, 微腔内光场变化情况 (δ1 = 0.1, Ein = 1 W1/2) (a) 基频光功率和(b)倍频光功率随光在腔内循环次数的变化曲线; 光场稳定后, (c)基频光场和(d)倍频光场随时间的演化; (e) 基频光场和(f)倍频广场的稳定分布

    Figure 3.  Light field evolution in the microresonator with the frequency detuning (δ1 = 0.1 and Ein = 1 W1/2): (a) Curves of the fundamental frequency waves power and (b) the second frequency waves power vary with the round trip number; evolution of (c) the fundamental frequency waves (d) the second frequency waves after the light fields are stable; stationary distribution of (e) the fundamental frequency waves light field and (f) the second frequency waves light field.

    图 4  失谐参量δ1取不同值时, 微腔内基频光和倍频光功率变化曲线 (a) 0.02 ≤ δ1 ≤ 0.08时, 基频光功率变化曲线; (b) 0.2 ≤ δ1 ≤ 0.8时, 基频光功率变化曲线; (c) 0.02 ≤ δ1 ≤ 0.08时, 倍频光功率变化曲线; (d) 0.2 ≤ δ1 ≤ 0.8时, 倍频光功率变化曲线

    Figure 4.  Influence of the frequency detuning δ1 on the power change curves: (a) Curves of the power variation for the fundamental frequency waves, 0.02 ≤ δ1 ≤ 0.08; (b) curves of the power variation for the fundamental frequency waves, 0.2 ≤ δ1 ≤ 0.8; (c) curves of the power variation for the second harmonic waves, 0.02 ≤ δ1 ≤ 0.08; (d) curves of the power variation for the second harmonic waves, 0.2 ≤ δ1 ≤ 0.8.

    图 5  失谐参量δ1 = 1时, 腔内光场稳定后基频光和倍频光的光场演化 (a) 基频光场的演化; (b) 倍频光场的演化; (c) 基频光光场变化周期内, 光在腔内每循环100次, 绘制其波形; (d) 倍频光光场变化周期内, 光在腔内每循环100次, 绘制其波形

    Figure 5.  Stable evolution of the dual light fields when δ1 = 1: (a) Evolution of the fundamental frequency waves; (b) evolution of the second harmonic waves; (c) intensity profiles of the fundamental frequency waves at six different moments within a period, the waveforms are plotted every hundred times; (d) intensity profiles of the second harmonic waves at six different moments within a period, the waveforms are plotted every hundred times.

    图 6  失谐参量δ1取值过大时, 腔内光场的混沌状态 (a) 基频光光场的混沌状态; (b) 倍频光光场的混沌状态; (c) 某一时刻基频光场的分布; (d) 某一时刻倍频光场的分布; (e) 与图(c)中光场对应的基频光光谱; (f) 与图(d)中光场对应的倍频光光谱

    Figure 6.  Chaos inside the microresonator, when the value of detuning parameter is too large: (a) Chaos of the fundamental frequency waves; (b) chaos of the second harmonic waves; (c) intensity profile of the fundamental frequency waves in a moment; (d) intensity profile of the second harmonic waves in a moment; (e) spectrum of the fundamental frequency waves; (f) spectrum of the second harmonic waves.

    图 7  抽运功率对腔内光功率的影响 (a)腔内基频光功率的变化情况; (b)腔内倍频光功率的变化情况

    Figure 7.  Influence of the pump power on the power change curves: (a) Power variation for the fundamental frequency waves; (b) power variation for the second harmonic waves.

    图 8  Ein = 100 W1/2时, 微腔内光场的演化 (a) 基频光功率的变化曲线; (b) 倍频光功率的变化曲线; (c) 光场稳定后, 腔内基频光光场的周期性演化; (d) 光场稳定后, 腔内倍频光光场的周期性演化

    Figure 8.  Evolution of the light field in the microresonator at Ein = 100 W1/2: (a) Curve of power variation for the fundamental frequency waves; (b) curve of power variation for the second harmonic waves; (c) periodic evolution of the fundamental frequency waves; (d) periodic evolution of the second harmonic waves.

    图 9  基频光和倍频光的功率变化曲线 Ein = 800 W1/2时(a)基频光和(b)倍频光的功率变化; Ein = 1000 W1/2时(c)基频光和(d) 倍频光的功率变化; Ein = 1200 W1/2时(e)基频光和(f)倍频光的功率变化

    Figure 9.  Curves of the power variation for the fundamental frequency waves and the second harmonic waves: Power variation for (a) the fundamental frequency waves and (b) the second harmonic waves at Ein = 800 W1/2; power variation for (c) the fundamental frequency waves and (d) the second harmonic waves at Ein = 1000 W1/2; power variation for (e) the fundamental frequency waves and (f) the second harmonic waves at Ein = 1200 W1/2.

    图 10  微腔FSR = 300 GHz时, 微腔内出现“图灵环”(Ein = 100 W1/2, δ1 = 0.1) (a)基频光功率变化曲线; (b)倍频光功率变化曲线; (c) 腔内光场稳定后, 基频光光场随时间的演化; (d) 腔内光场稳定后, 倍频光光场随时间的演化; (e) 腔内光场稳定后, 基频光的光谱; (f) 腔内光场稳定后, 倍频光的光谱

    Figure 10.  Turning patterns in the microresonator, when FSR = 300 GHz (Ein = 100 W1/2, δ1 = 0.1): (a) Curves of the power variation for the fundamental frequency waves; (b) curves of the power variation for the second harmonic waves; (c) evolution of the fundamental frequency waves; (d) evolution of the second harmonic waves; (e)spectra of the fundamental harmonic waves; (f) spectra of the second harmonic waves.

    图 11  微腔FSR = 298 GHz时, 稳定后的光场分布及光谱 (a)基频光的光场分布; (b) 倍频光的光场分布; (c) 基频光光谱; (d) 倍频光光谱

    Figure 11.  Stable intensity profile and spectra of the fundamental frequency waves for FSR = 298 GHz: (a) Intensity profile of the fundamental frequency waves; (c) spectrum of the fundamental frequency waves; (b) intensity profile of second harmonic waves; (d) spectrum of second harmonic waves.

    图 12  微腔FSR = 295 GHz时, 稳定后的基频光和倍频光光场分布及光谱 (a) 基频光光场随时间的演化; (b) 倍频光光场随时间的演化; (c) 基频光场的瞬时分布; (d) 倍频光场的瞬时分布; (e) 基频光场的光谱; (f) 倍频光场的光谱

    Figure 12.  Stable intensity profile and spectra of the fundamental frequency waves for FSR = 295 GH: (a) Evolution of the fundamental frequency waves; (b) evolution of the second frequency waves; (c) intensity profile of the fundamental frequency waves; (d) spectrum of the fundamental frequency waves; (e) intensity profile of the second harmonic waves; (f) spectrum of the second harmonic waves.

  • [1]

    邢书剑, 张福民, 曹士英, 王高文, 曲兴华 2013 物理学报 62 170603Google Scholar

    Xing S J, Zhang F M, Cao S Y, Wang G W, Qu X H 2013 Acta Phys. Sin. 62 170603Google Scholar

    [2]

    Niering M, Holzwarth R, Reichert J, Pokasov P, Udem T, Weitz M, Hansch T W, Lemonde P, Santarelli G, Abgrall M, Laurent P, Salomon C, Clairon A 2000 Phys. Rev. Lett. 84 5496Google Scholar

    [3]

    孟飞, 曹士英, 蔡岳, 王贵重, 曹建平, 李天初, 方占军 2011 物理学报 60 100601Google Scholar

    Meng F, Cao S Y, Cai Y, Wang G Z, Cao J P, Li T C, Fang Z J 2011 Acta Phys. Sin. 60 100601Google Scholar

    [4]

    Swann W C, McFerran J J, Coddington I, Newbury N R, Hartl I, Fermann M E, Westbrook P S, Nicholson J W, Feder K S, Langrock C, Fejer M M 2006 Opt. Lett. 31 3074

    [5]

    Washburn B R, Fox R W, Newbury N R, Nicholson J W, Feder K, Westbrook P S, Jorgensen C G 2004 Opt. Express 12 4999Google Scholar

    [6]

    韩海年, 魏志义, 赵刚 2012 物理 41 249

    Han H N, Wei Z Y, Zhao G 2012 Physics 41 249

    [7]

    Herr T, Brasch V, Jost J D, Wang C Y, Kondratiev N M, Gorodetsky M L, Kippenberg T J 2013 Nat. Photon. 8 145

    [8]

    Song Q H 2019 Sci. China-Phys. Mech. Astron. 62 074231Google Scholar

    [9]

    Zhang X Y, Cao Q T, Wang Z, Liu Y X, Qiu C W, Yang L, Gong Q H, Xiao Y F 2019 Nat. Photon. 13 21Google Scholar

    [10]

    Godey C, Balakireva I V, Aurélien C, Chembo Y K 2014 Phys. Rev. A 89 722

    [11]

    Huang S W, Zhou H, Yang J, McMillan J F, Matsko A, Yu M, Kwong D L, Maleki L, Wong C W 2014 Phys. Rev. Lett. 114 053901

    [12]

    Wang W Q, Lu Z Z, Zhang W F, Chu Sai T, Little B E, Wang L R, Xie X P, Liu M L, Yang Q H, Wang L, Zhao J G, Wang G X, Sun Q B, Liu Y S, Wang Y S, Zhao W 2018 Opt. Lett. 43 2002Google Scholar

    [13]

    Xue X X, Xuan Y, Liu Y, Wang P H, Chen S, Wang J, Leaird D E, Qi M H, Weiner A M 2015 Nat. Photon. 9 594Google Scholar

    [14]

    Coillet A, Balakireva I, Henriet R, Saleh K, Larger L, Dudley J M, Menyuk C R, Chembo Y K 2013 IEEE Photon. J. 5 6100409Google Scholar

    [15]

    Matsko A B, Liang W, Savchenkov A A, Maleki L 2013 Opt. Lett. 38 525Google Scholar

    [16]

    Coen S, Randle H G, Sylvestre T, Erkintalo M 2013 Opt. Lett. 38 37Google Scholar

    [17]

    Xue X X, Leo F, Xuan Y, Villegas J A J, Wang P H, Leaird D E, Erkintalo M, Qi M H, Weiner Andrew M 2017 Light: Sci. Appl. 6 e16253Google Scholar

    [18]

    Lin J T, Xu Y X, Ni J L, Wang M, Fang Z W, Qiao L L, Fang W, Cheng Y 2016 Phys. Rev. Appl. 6 014002Google Scholar

    [19]

    Chembo Y K, Menyuk C R, 2014 Phys. Rev. A 87 053852

    [20]

    Lin J T, Xu Y X, Fang Z W, Wang M, Song J X, Wang N W, Qiao L L, Fang W, Cheng Y 2015 Sci. China-Phys. Mech. Astron. 58 114209Google Scholar

    [21]

    Coen S, Haelterman M 2001 Opt. Lett. 26 39Google Scholar

    [22]

    And T, Boccaletti A, Grebogi C, Lai Y C, Mancini H, Maza D 2000 Phys. Rep. 329 2000

    [23]

    Akhmediev N, Pelinovsky E 2010 Eur. Phys. J. Spec. Top. 185 1Google Scholar

  • [1] Hu Sheng-Run, Ji Xue-Qiang, Wang Jin-Jin, Yan Jie-Yun, Zhang Tian-Yue, Li Pei-Gang. Ultralow switching threshold optical bistable devices based on epsilon-near-zero Ga2O3-SiC-Ag multilayer structures. Acta Physica Sinica, 2024, 73(5): 054201. doi: 10.7498/aps.73.20231534
    [2] Zhang Xiao-Li, Wang Qing-Wei, Yao Wen-Xiu, Shi Shao-Ping, Zheng Li-Ang, Tian Long, Wang Ya-Jun, Chen Li-Rong, Li Wei, Zheng Yao-Hui. Influence of thermal lens effect on second harmonic process in semi-monolithic cavity scheme. Acta Physica Sinica, 2022, 71(18): 184203. doi: 10.7498/aps.71.20220575
    [3] Xu Fan, Zhao Yan, Wu Yu-Hang, Wang Wen-Chi, Jin Xue-Ying. Stability and non-linear dynamic analysis of Kerr optical frequencycombs in dual-coupled microcavities with high-order dispersion. Acta Physica Sinica, 2022, 71(18): 184204. doi: 10.7498/aps.71.20220691
    [4] Guo Qi-Qi, Chen Yi-Hang. Enhanced nonlinear optical effects based on strong coupling between epsilon-near-zero mode and gap surface plasmons. Acta Physica Sinica, 2021, 70(18): 187303. doi: 10.7498/aps.70.20210290
    [5] Li Hai-Peng, Zhou Jia-Sheng, Ji Wei, Yang Zi-Qiang, Ding Hui-Min, Zhang Zi-Tao, Shen Xiao-Peng, Han Kui. Effect of edge on nonlinear optical property of graphene quantum dots. Acta Physica Sinica, 2021, 70(5): 057801. doi: 10.7498/aps.70.20201643
    [6] Bai Rui-Xue, Yang Jue-Han, Wei Da-Hai, Wei Zhong-Ming. Research progress of low-dimensional semiconductor materials in field of nonlinear optics. Acta Physica Sinica, 2020, 69(18): 184211. doi: 10.7498/aps.69.20200206
    [7] Xu Xin, Jin Xue-Ying, Gao Hao-Ran, Cheng Jie, Lu Yang, Chen Dong, Yu Lian-Dong. Analysis of frequency tuning process of dual coupled optical microcavities. Acta Physica Sinica, 2020, 69(18): 184207. doi: 10.7498/aps.69.20200530
    [8] Guan Xiao-Tong, Fu Wen-Jie, Lu Dun, Yang Tong-Bin, Yan Yang, Yuan Xue-Song. Design of second harmonic terahertz gyrotron cavity based on double confocal waveguide. Acta Physica Sinica, 2020, 69(6): 068401. doi: 10.7498/aps.69.20191222
    [9] Zeng Zhou-Xiao-Song, Wang Xiao, Pan An-Lian. Second harmonic generation of two-dimensional layered materials: characterization, signal modulation and enhancement. Acta Physica Sinica, 2020, 69(18): 184210. doi: 10.7498/aps.69.20200452
    [10] Gu Hong-Ming, Huang Yong-Qing, Wang Huan-Huan, Wu Gang, Duan Xiao-Feng, Liu Kai, Ren Xiao-Min. Theoretical analysis of new optical microcavity. Acta Physica Sinica, 2018, 67(14): 144201. doi: 10.7498/aps.67.20180067
    [11] Deng Jun-Hong, Li Gui-Xin. Nonlinear photonic metasurfaces. Acta Physica Sinica, 2017, 66(14): 147803. doi: 10.7498/aps.66.147803
    [12] Du Jin-Jin, Li Wen-Fang, Wen Rui-Juan, Li Gang, Zhang Tian-Cai. Precision measurement of resonate frequency and the effective cavity length of the high finesse optical micro-cavity. Acta Physica Sinica, 2013, 62(19): 194203. doi: 10.7498/aps.62.194203
    [13] Cheng Zheng-Fu, Long Xiao-Xia, Zheng Rui-Lun. Influence of temperature on the Bose condensation of photons and excitons in optic microcavity. Acta Physica Sinica, 2010, 59(12): 8377-8384. doi: 10.7498/aps.59.8377
    [14] Yang Guang, Chen Zheng-Hao. Large optical nonlinearities in Ag-doped BaTiO3 nanocomposite films. Acta Physica Sinica, 2007, 56(2): 1182-1187. doi: 10.7498/aps.56.1182
    [15] Huang Xiao-Ming, Tao Li-Min, Guo Ya-Hui, Gao Yun, Wang Chuan-Kui. Theoretical studies of nonlinear optical properties of a novel double-conjugated-segment molecule. Acta Physica Sinica, 2007, 56(5): 2570-2576. doi: 10.7498/aps.56.2570
    [16] Chen Liang, Liang Chang-Hong, Dang Xiao-Jie. Second-harmonic generation in nonlinear left-handed metamaterials. Acta Physica Sinica, 2007, 56(11): 6398-6402. doi: 10.7498/aps.56.6398
    [17] Liang Xiao-Rui, Zhao Bo, Zhou Zhi-Hua. Ab initio study on the second-order nonlinear optical properties of some coumarin derivatives. Acta Physica Sinica, 2006, 55(2): 723-728. doi: 10.7498/aps.55.723
    [18] Zhang Ming-Xin, Wu Ke-Chen, Liu Cai-Ping, Wei Yong-Qin. Computational study on the exchange-correlation function in density functional theory and optical nonlinearity of transition-metal complexes. Acta Physica Sinica, 2005, 54(4): 1762-1770. doi: 10.7498/aps.54.1762
    [19] Liu Tao, Zhang Tian-Cai, Wang Jun-Min, Peng Kun-Chi. Optical dipole trap in a high-finesse micro-cavity. Acta Physica Sinica, 2004, 53(5): 1346-1351. doi: 10.7498/aps.53.1346
    [20] Ni Pei-Gen, Ma Bo-Qin, Cheng Bing-Ying, Zhang Dao-Zhong. Two-dimensional periodically poled LiNbO3 nonlinear photonic crystal. Acta Physica Sinica, 2003, 52(8): 1925-1928. doi: 10.7498/aps.52.1925
Metrics
  • Abstract views:  9149
  • PDF Downloads:  162
  • Cited By: 0
Publishing process
  • Received Date:  26 August 2019
  • Accepted Date:  24 September 2019
  • Available Online:  01 January 2020
  • Published Online:  20 January 2020

/

返回文章
返回