搜索

x

留言板

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

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

MoS2中S原子空位形成的非绝热动力学研究

王月 马杰

引用本文:
Citation:

MoS2中S原子空位形成的非绝热动力学研究

王月, 马杰

Non-adiabatic dynamic study of S vacancy formation in MoS2

Wang Yue, Ma Jie
PDF
HTML
导出引用
  • 缺陷是半导体领域中最核心的问题. 采用含时密度泛函方法, 模拟了S原子脱离MoS2晶格形成空位缺陷过程中的电子动力学行为, 发现该过程中存在显著的非绝热效应. 非绝热效应导致S原子需要消耗更多能量以脱离晶格形成空位缺陷. 随着S原子的初始动能增大, 其脱离晶格形成空位的能量势垒也持续增大, 并且在初始动能达到22 eV附近时发生了阶跃式的增长. 这是由朗道-齐纳电子跃迁和能级间库仑作用共同导致的. 非绝热效应还改变了脱离晶格的S原子上电荷的轨道分布, 以及晶格中缺陷附近的电荷分布. 此外, 还发现该过程中自旋轨道耦合十分重要, 必须被考虑. 本文阐明了MoS2中S原子空位的形成机制, 尤其是电子非绝热动力学的重要作用, 为进一步研究缺陷对材料物理性质的调控提供了理论基础.
    Defect is one of the central issues in semiconductors. MoS2 is sensitive to irradiation and can be damaged by electron beams, accompanied with the generation of sulfur vacancies. However, the dynamics for the defect generation process is still unclear. In this work, we employ the time-dependent density functional theory to simulate the process of a sulfur atom sputtering from the MoS2, producing a sulfur vacancy defect in the lattice. We find that there exists a strong non-adiabatic effect in the process.During the formation of the sulfur vacancy, there exist electron transitions which can be described by the Landau-Zener model. As the sulfur atom leaves away from the lattice, two energy levels from the valence bands rise up and one energy level from the conduction band falls down. When the spin-orbit coupling (SOC) is not considered, those levels do not couple with each other. However, when the SOC is taken into account, electrons can transit between those levels. The transition probability is related to the kinetic energy of the sputtered sulfur atom. As the kinetic energy of the sulfur atom increases, the non-adiabatic electron transitions are enhanced. The evolution of the energy levels is also strongly dependent on the kinetic energy of the sputtered sulfur atom, which is induced by the non-adiabatic electron transition. It is worth noting that the SOC plays a key role in sputtering sulfur atoms, although the system produces no magnetic moments in the whole process.The non-adiabatic effect enhances the energy barrier of the sulfur sputtering. As the initial kinetic energy of the sputtered sulfur atom increases, the energy barrier increases, and exhibits a jump around the initial kinetic energy of ~22 eV, which can be explained by the non-adiabatic electron occupation and the Coulomb repulsion. Beside the energy barrier, the non-adiabatic effect also modifies the charge distribution. When the kinetic energy of the sputtered sulfur atom is relatively low, more electrons occupy the pz level; when the kinetic energy is relatively high, more electrons occupy the pxy level instead. The sputtered sulfur atom always carries a bit more electrons, leaving holes around the vacancy defect.Our work reveals the dynamics of the sulfur sputtering and vacancy formation in MoS2, particularly the non-adiabatic effect in the process. It builds the theoretical foundation for defect engineering.
      通信作者: 马杰, majie@bit.edu.cn
    • 基金项目: 国家自然科学基金(批准号: 12074032)资助的课题.
      Corresponding author: Ma Jie, majie@bit.edu.cn
    • Funds: Project supported by the National Natural Science Foundation of China (Grant No. 12074032).
    [1]

    Bozheyev F, Ellmer K 2022 J. Mater. Chem. A 10 9327Google Scholar

    [2]

    Khatua D P, Singh A, Gurung S, Khan S, Tanwar M, Kumar R, Jayabalan J 2022 J. Phys. Condens. Matter. 34 155401Google Scholar

    [3]

    Ganatra R, Zhang Q 2014 ACS Nano 8 4074Google Scholar

    [4]

    Singh E, Singh P, Kim K S, Yeom G Y, Nalwa H S 2019 ACS Appl. Mater. Interfaces 11 11061Google Scholar

    [5]

    Nawz T, Safdar A, Hussain M, Sung Lee D, Siyar M 2020 Crystals 10 902Google Scholar

    [6]

    Sanikop R, Sudakar C 2019 ACS Appl. Nano Mater. 3 576Google Scholar

    [7]

    Islam M M, Dev D, Krishnaprasad A, Tetard L, Roy T 2020 Sci. Rep. 10 21870Google Scholar

    [8]

    Wu F, Tian H, Shen Y, Hou Z, Ren J, Gou G Y, Sun Y B, Yang Y, Ren T L 2022 Nature 603 259Google Scholar

    [9]

    Cao Y 2021 ACS Nano 15 11014Google Scholar

    [10]

    Zhao G Y, Deng H, Nathaniel T, Michael G, Abdellah L, Peng Q, Yan J A, Wang C D, Lan Y C 2019 Appl. Sci. 9 678Google Scholar

    [11]

    Pandey M, Rasmussen F A, Kuhar K, Olsen T, Jacobsen K W, Thygesen K S 2016 Nano Lett. 16 2234Google Scholar

    [12]

    Wang Y Q, Deng L J, Wei Q L, Wan Y, Liu Z, Lu X, Li Y, Bi L, Zhang L, Lu H P, Chen H Y, Zhou P H, Zhang L B, Cheng Y C, Zhao X X, Ye Y, Huang W, Pennycook S J, Loh K P, Peng B 2020 Nano Lett. 20 2129Google Scholar

    [13]

    Taffelli A, Dirè S, Quaranta A, Pancheri L 2021 Sensors 21 2758Google Scholar

    [14]

    Sorkin V, Zhou H, Yu Z G, Ang K W, Zhang Y W 2022 Sci. Rep. 12 18001Google Scholar

    [15]

    Wang Y, Liu S, Li Q, Quhe R, Yang C, Guo Y, Zhang X, Pan Y, Li J, Zhang H, Xu L, Shi B, Tang H, Li Y, Yang J, Zhang Z, Xiao L, Pan F, Lu J 2021 Rep. Prog. Phys. 84 056501Google Scholar

    [16]

    Cheng Y, Wan R, Li L, Liu Z, Yan S, Li L, Wang J, Gao Y 2022 J. Opt. 24 2040Google Scholar

    [17]

    Cheng Y, Song H, Wu H, Zhang P, Tang Z, Lu S 2020 Chem. Asian J. 15 3123Google Scholar

    [18]

    Xie Y, Liang F, Chi S, Wang D, Zhong K, Yu H, Zhang H, Chen Y, Wang J 2020 ACS Appl. Mater. Interfaces 12 7351Google Scholar

    [19]

    Yang Z, Zhu L, Lü C, Zhang R, Wang H, Wang J, Zhang Q 2021 Mater. Chem. Front. 5 5880Google Scholar

    [20]

    Komsa H P, Kotakoski J, Kurasch S, Lehtinen O, Kaiser U, Krasheninnikov A V 2012 Phys. Rev. Lett. 109 035503Google Scholar

    [21]

    Susi T, Meyer J C, Kotakoski J 2019 Nat. Rev. Phys. 1 397Google Scholar

    [22]

    Yoshimura A, Lamparski M, Kharche N, Meunier V 2018 Nanoscale 10 2388Google Scholar

    [23]

    Ugurlu O, Haus J, Gunawan A A, Thomas M G, Maheshwari S, Tsapatsis M, Mkhoyan K A 2011 Phys. Rev. B 83 113408Google Scholar

    [24]

    Lehnert T, Lehtinen O, Algara-Siller G, Kaiser U 2017 Appl. Phys. Lett. 110 033106Google Scholar

    [25]

    Tai K L, Huang C W, Cai R F, Huang G M, Tseng Y T, Chen J, Wu W W 2020 Small 16 1905516Google Scholar

    [26]

    Meyer J C, Eder F, Kurasch S, Skakalova V, Kotakoski J, Park H J, Roth S, Chuvilin A, Eyhusen S, Gred B, Krasheninnikov A V, Kaiser U 2012 Phys. Rev. Lett. 108 196102Google Scholar

    [27]

    Kretschmer S, Lehnert T, Kaiser U, Krasheninnikov A V 2020 Nano Lett. 20 2865Google Scholar

    [28]

    Runge E, Gross E K 1984 Phys. Rev. Lett. 52 997Google Scholar

    [29]

    Jia W, Fu J, Cao Z, Wang L, Chi X, Gao W, Wang L W 2013 J. Comput. Phys. 251 102Google Scholar

    [30]

    Jia W, Cao Z, Wang L, Fu J, Chi X, Gao W, Wang L W 2013 Comput. Phys. Commun. 184 9Google Scholar

    [31]

    Hamann D R, Schlüter M, Chiang C 1979 Phys. Rev. Lett. 43 1494Google Scholar

    [32]

    Hamann D R 2013 Phys. Rev. B 88 085117Google Scholar

    [33]

    Shen Y T, Xu T, Tan X D, He L B, Yin K B, Wan N, Sun L T 2018 Adv. Mater. 30 1705954Google Scholar

    [34]

    Zener C 1932 Proc. R. Soc. London, Ser. A 137 696Google Scholar

    [35]

    Olson A J, Wang S J, Niffenegger R J, Li C H, Greene C H, Chen Y P 2014 Phys. Rev. A 90 013616Google Scholar

    [36]

    Wang L W 2002 Phys. Rev. B 65 153410Google Scholar

  • 图 1  (a) 5 × 5的完整MoS2超胞; (b) 发生弹性碰撞后, S原子获得一定的初始速度v脱离晶格, 其速度方向垂直于晶格; (c) S原子脱离晶格后运动到距离晶格为d的位置. 图中黄色球与紫色球分别代表S原子与Mo原子, 高亮的原子为脱离晶格的S原子

    Fig. 1.  (a) The 5 × 5 supercell of MoS2; (b) after the elastic collision with electrons, the S atom gets an initial velocity v, whose direction is perpendicular to the lattice; (c) the moving distance d of the S atom after sputtering from the lattice. Yellow and purple balls represent the sulfur and molybdenum atoms, respectively. The highlighted ball represents the sputtered S atom.

    图 2  (a) 当S原子以不同初始动能脱离晶格时, 系统势能Ep随S原子运动距离d的变化关系, 其中初始时刻完整MoS2晶格的情况为势能零点; (b) S原子脱离晶格所需越过的能量势垒与初始动能的关系

    Fig. 2.  (a) For the sulfur atom with various initial kinetic energies, the energy Ep as a function of the atom distance d, where the energy of the perfect MoS2 lattice is set to zero; (b) energy barriers for the S sputtering process as a function of the initial kinetic energy.

    图 3  系统能级随S原子运动距离d的变化情况. 能级上的红色气泡表示该能级被电子占据, 气泡大小表示电子占据数. 插图为箭头所指能级的电荷分布 (a) 不考虑SOC; (b) 考虑SOC且脱离晶格的S原子的初始动能较低(Ek = 15.2 eV); (c) 考虑SOC且脱离晶格的S原子的初始动能较高(Ek = 30.4 eV)

    Fig. 3.  Energy levels of the system as a function of the atom distance d. Red bubble on an energy level indicates that the energy level is occupied by electrons, and the size of the bubble indicates the number of electrons occupied. The insets show the charge densities of the corresponding states: (a) Without SOC; (b) with SOC and a low initial kinetic energy of the S atom (Ek = 15.2 eV); (c) with SOC and a high initial kinetic energy of the S atom (Ek = 30.4 eV).

    图 4  (a) 当脱离晶格的S原子初始动能较低时(Ek = 15.2 eV), 相关能级上电子占据数的变化; (b) 当脱离晶格的S原子初始动能较高时(Ek = 30.4 eV), 相关能级上电子占据数的变化; (c) 原子运动距离d = 3 Å时, 高能量能级上的电子占据数随S原子初始动能$ {E}_{{\rm{k}}} $的变化

    Fig. 4.  (a) When the initial kinetic energy of the sputtered S atom is low (Ek = 15.2 eV), the electron occupations on the relevant states as functions of the atom distance; (b) when the initial kinetic energy of the sputtered S atom is high (Ek = 30.4 eV), the electron occupations on the relevant states as functions of the atom distance; (c) electron occupation on the high-energy state as a function of the initial kinetic energy $ {E}_{{\rm{k}}} $ of the sputtered S atom, when the atom distance d = 3 Å.

    图 5  当脱离晶格原子的运动距离d取不同值时, 非绝热模拟与绝热模拟所得的电荷密度之差(蓝色和黄色等值面分别代表电荷的减少和增加) (a)—(c) 低速条件(Ek = 15.2 eV)下, (a) d = 1.75 Å, (b) d = 2.25 Å, (c) d = 5.00 Å; (d)—(f) 高速条件(Ek = 30.4 eV)下, (d) d = 1.75 Å, (e) d = 2.25 Å, (f) d = 5.00 Å

    Fig. 5.  Charge density difference between the non-adiabatic and adiabatic results at different atom distances d (Blue and yellow isosurfaces represent the charge depletion and accumulation respectively): (a)–(c) When the initial kinetic energy of the sputtered S atom is low (Ek = 15.2 eV), (a) d = 1.75 Å, (b) d = 2.25 Å, (c) d = 5.00 Å; (d)–(f) when the initial kinetic energy of the sputtered S atom is high (Ek = 30.4 eV), (d) d = 1.75 Å, (e) d = 2.25 Å, (f) d = 5.00 Å.

    图 6  不同初始动能下, 非绝热模拟中脱离晶格的S原子上的电荷数随d的关系

    Fig. 6.  Number of charges on the sputtered S atom as a function of d with different initial kinetic energies.

  • [1]

    Bozheyev F, Ellmer K 2022 J. Mater. Chem. A 10 9327Google Scholar

    [2]

    Khatua D P, Singh A, Gurung S, Khan S, Tanwar M, Kumar R, Jayabalan J 2022 J. Phys. Condens. Matter. 34 155401Google Scholar

    [3]

    Ganatra R, Zhang Q 2014 ACS Nano 8 4074Google Scholar

    [4]

    Singh E, Singh P, Kim K S, Yeom G Y, Nalwa H S 2019 ACS Appl. Mater. Interfaces 11 11061Google Scholar

    [5]

    Nawz T, Safdar A, Hussain M, Sung Lee D, Siyar M 2020 Crystals 10 902Google Scholar

    [6]

    Sanikop R, Sudakar C 2019 ACS Appl. Nano Mater. 3 576Google Scholar

    [7]

    Islam M M, Dev D, Krishnaprasad A, Tetard L, Roy T 2020 Sci. Rep. 10 21870Google Scholar

    [8]

    Wu F, Tian H, Shen Y, Hou Z, Ren J, Gou G Y, Sun Y B, Yang Y, Ren T L 2022 Nature 603 259Google Scholar

    [9]

    Cao Y 2021 ACS Nano 15 11014Google Scholar

    [10]

    Zhao G Y, Deng H, Nathaniel T, Michael G, Abdellah L, Peng Q, Yan J A, Wang C D, Lan Y C 2019 Appl. Sci. 9 678Google Scholar

    [11]

    Pandey M, Rasmussen F A, Kuhar K, Olsen T, Jacobsen K W, Thygesen K S 2016 Nano Lett. 16 2234Google Scholar

    [12]

    Wang Y Q, Deng L J, Wei Q L, Wan Y, Liu Z, Lu X, Li Y, Bi L, Zhang L, Lu H P, Chen H Y, Zhou P H, Zhang L B, Cheng Y C, Zhao X X, Ye Y, Huang W, Pennycook S J, Loh K P, Peng B 2020 Nano Lett. 20 2129Google Scholar

    [13]

    Taffelli A, Dirè S, Quaranta A, Pancheri L 2021 Sensors 21 2758Google Scholar

    [14]

    Sorkin V, Zhou H, Yu Z G, Ang K W, Zhang Y W 2022 Sci. Rep. 12 18001Google Scholar

    [15]

    Wang Y, Liu S, Li Q, Quhe R, Yang C, Guo Y, Zhang X, Pan Y, Li J, Zhang H, Xu L, Shi B, Tang H, Li Y, Yang J, Zhang Z, Xiao L, Pan F, Lu J 2021 Rep. Prog. Phys. 84 056501Google Scholar

    [16]

    Cheng Y, Wan R, Li L, Liu Z, Yan S, Li L, Wang J, Gao Y 2022 J. Opt. 24 2040Google Scholar

    [17]

    Cheng Y, Song H, Wu H, Zhang P, Tang Z, Lu S 2020 Chem. Asian J. 15 3123Google Scholar

    [18]

    Xie Y, Liang F, Chi S, Wang D, Zhong K, Yu H, Zhang H, Chen Y, Wang J 2020 ACS Appl. Mater. Interfaces 12 7351Google Scholar

    [19]

    Yang Z, Zhu L, Lü C, Zhang R, Wang H, Wang J, Zhang Q 2021 Mater. Chem. Front. 5 5880Google Scholar

    [20]

    Komsa H P, Kotakoski J, Kurasch S, Lehtinen O, Kaiser U, Krasheninnikov A V 2012 Phys. Rev. Lett. 109 035503Google Scholar

    [21]

    Susi T, Meyer J C, Kotakoski J 2019 Nat. Rev. Phys. 1 397Google Scholar

    [22]

    Yoshimura A, Lamparski M, Kharche N, Meunier V 2018 Nanoscale 10 2388Google Scholar

    [23]

    Ugurlu O, Haus J, Gunawan A A, Thomas M G, Maheshwari S, Tsapatsis M, Mkhoyan K A 2011 Phys. Rev. B 83 113408Google Scholar

    [24]

    Lehnert T, Lehtinen O, Algara-Siller G, Kaiser U 2017 Appl. Phys. Lett. 110 033106Google Scholar

    [25]

    Tai K L, Huang C W, Cai R F, Huang G M, Tseng Y T, Chen J, Wu W W 2020 Small 16 1905516Google Scholar

    [26]

    Meyer J C, Eder F, Kurasch S, Skakalova V, Kotakoski J, Park H J, Roth S, Chuvilin A, Eyhusen S, Gred B, Krasheninnikov A V, Kaiser U 2012 Phys. Rev. Lett. 108 196102Google Scholar

    [27]

    Kretschmer S, Lehnert T, Kaiser U, Krasheninnikov A V 2020 Nano Lett. 20 2865Google Scholar

    [28]

    Runge E, Gross E K 1984 Phys. Rev. Lett. 52 997Google Scholar

    [29]

    Jia W, Fu J, Cao Z, Wang L, Chi X, Gao W, Wang L W 2013 J. Comput. Phys. 251 102Google Scholar

    [30]

    Jia W, Cao Z, Wang L, Fu J, Chi X, Gao W, Wang L W 2013 Comput. Phys. Commun. 184 9Google Scholar

    [31]

    Hamann D R, Schlüter M, Chiang C 1979 Phys. Rev. Lett. 43 1494Google Scholar

    [32]

    Hamann D R 2013 Phys. Rev. B 88 085117Google Scholar

    [33]

    Shen Y T, Xu T, Tan X D, He L B, Yin K B, Wan N, Sun L T 2018 Adv. Mater. 30 1705954Google Scholar

    [34]

    Zener C 1932 Proc. R. Soc. London, Ser. A 137 696Google Scholar

    [35]

    Olson A J, Wang S J, Niffenegger R J, Li C H, Greene C H, Chen Y P 2014 Phys. Rev. A 90 013616Google Scholar

    [36]

    Wang L W 2002 Phys. Rev. B 65 153410Google Scholar

  • [1] 杨静, 冯少蓉, 张涛, 牛旭平, 王荣, 李敏, 于润升, 曹兴忠, 王宝义. B位空位补偿型钐掺杂PZT(54/46)陶瓷中的缺陷分析及其对压电性能的影响. 物理学报, 2024, 73(7): 077701. doi: 10.7498/aps.73.20231872
    [2] 段聪, 刘俊杰, 陈永杰, 左慧玲, 董健生, 欧阳钢. MoS2/SiO2界面黏附性能的尺寸和温度效应. 物理学报, 2024, 73(5): 056801. doi: 10.7498/aps.73.20231648
    [3] 刘东静, 周福, 陈帅阳, 胡志亮. 氮化镓/石墨烯/碳化硅异质界面热输运特性的分子动力学研究. 物理学报, 2023, 72(15): 157901. doi: 10.7498/aps.72.20230537
    [4] 武鹏, 谈论, 李炜, 曹立伟, 赵俊博, 曲尧, 李昂. 大面积单层二硫化钼的制备及其光电性能. 物理学报, 2023, 72(11): 118101. doi: 10.7498/aps.72.20230273
    [5] 王婉玉, 石凯熙, 李金华, 楚学影, 方铉, 匡尚奇, 徐国华. MoO3覆盖层对MoS2基光伏型光电探测器性能的影响. 物理学报, 2023, 72(14): 147301. doi: 10.7498/aps.72.20230464
    [6] 孔宇晗, 王蓉, 徐明生. CuPc/MoS2范德瓦耳斯异质结荧光特性. 物理学报, 2022, 71(12): 128103. doi: 10.7498/aps.71.20220132
    [7] 刘东静, 王韶铭, 杨平. 石墨烯/碳化硅异质界面热学特性的分子动力学模拟. 物理学报, 2021, 70(18): 187302. doi: 10.7498/aps.70.20210613
    [8] 邓文, 汪礼胜, 刘嘉宁, 余雪玲, 陈凤翔. 光电协控多层MoS2记忆晶体管的阻变行为与机理研究. 物理学报, 2021, 70(21): 217302. doi: 10.7498/aps.70.20210750
    [9] 李晓克, 冯伟. 非绝热分子动力学的量子路径模拟. 物理学报, 2017, 66(15): 153101. doi: 10.7498/aps.66.153101
    [10] 兰生, 李焜, 高新昀. 基于分子动力学的石墨炔纳米带空位缺陷的导热特性. 物理学报, 2017, 66(13): 136801. doi: 10.7498/aps.66.136801
    [11] 张红, 尹海峰, 张开彪, 林家和. 基于含时密度泛函理论的表面等离激元研究进展. 物理学报, 2015, 64(7): 077303. doi: 10.7498/aps.64.077303
    [12] 彭亚晶, 蒋艳雪. 分子空位缺陷对环三亚甲基三硝胺含能材料几何结构、电子结构及振动特性的影响. 物理学报, 2015, 64(24): 243102. doi: 10.7498/aps.64.243102
    [13] 王志萍, 朱云, 吴亚敏, 张秀梅. 质子与羟基碰撞的含时密度泛函理论研究. 物理学报, 2014, 63(2): 023401. doi: 10.7498/aps.63.023401
    [14] 王志萍, 陈健, 吴寿煜, 吴亚敏. 碳分子线C5在激光场中的含时密度泛函理论研究. 物理学报, 2013, 62(12): 123302. doi: 10.7498/aps.62.123302
    [15] 杨平, 王晓亮, 李培, 王欢, 张立强, 谢方伟. 氮掺杂和空位对石墨烯纳米带热导率影响的分子动力学模拟. 物理学报, 2012, 61(7): 076501. doi: 10.7498/aps.61.076501
    [16] 李海宏, 刘文, 刘德胜. 理论计算中电势能零点的选取对电荷注入的影响. 物理学报, 2011, 60(9): 097201. doi: 10.7498/aps.60.097201
    [17] 袁剑辉, 程玉民, 张振华. 空位结构缺陷对C纳米管弹性性质的影响. 物理学报, 2009, 58(4): 2578-2584. doi: 10.7498/aps.58.2578
    [18] 李海宏, 李冬梅, 刘 文, 李 元, 刘晓静, 刘德胜, 解士杰. 金属/掺杂聚合物/金属结构中载流子的注入与输运. 物理学报, 2008, 57(2): 1117-1122. doi: 10.7498/aps.57.1117
    [19] 崔 磊, 顾 斌, 滕玉永, 胡永金, 赵 江, 曾祥华. 脉冲激光偏振方向对氮分子高次谐波的影响--基于含时密度泛函理论的模拟. 物理学报, 2006, 55(9): 4691-4694. doi: 10.7498/aps.55.4691
    [20] 张超, 吕海峰, 张庆瑜. 低能Pt原子与Pt(111)表面相互作用的分子动力学模拟. 物理学报, 2002, 51(10): 2329-2334. doi: 10.7498/aps.51.2329
计量
  • 文章访问数:  2569
  • PDF下载量:  353
  • 被引次数: 0
出版历程
  • 收稿日期:  2023-05-15
  • 修回日期:  2023-07-13
  • 上网日期:  2023-09-05
  • 刊出日期:  2023-11-20

/

返回文章
返回