搜索

x

留言板

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

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

二维双金属铁磁半导体CrMoI6的电子结构与稳定性

黄玉昊 张贵涛 王如倩 陈乾 王金兰

引用本文:
Citation:

二维双金属铁磁半导体CrMoI6的电子结构与稳定性

黄玉昊, 张贵涛, 王如倩, 陈乾, 王金兰

Electronic structure and stability of two-dimensional bimetallic ferromagnetic semiconductor CrMoI6

Huang Yu-Hao, Zhang Gui-Tao, Wang Ru-Qian, Chen Qian, Wang Jin-Lan
PDF
HTML
导出引用
  • 二维磁性半导体由于兼具磁性、半导体性和特殊的二维结构而受到人们的广泛关注, 为纳米级自旋电子和光电子器件的研发应用和相关的基础理论研究提供了新的思路和平台. 基于第一性原理计算, 在对一系列二维双金属碘化物CrTMI6的交换能进行初步筛选的基础上, 选出了具有铁磁性的CrMoI6单层结构. 进一步计算表明, CrMoI6单层的电子能带结构展示出理想的半导体特性, 计算得到的带隙约为1.7 eV, 而且还具有很大的磁各向异性能(741.3 μeV/TM). 通过基于海森伯模型的蒙特卡罗模拟预测这一材料的居里温度达到92 K, 约为CrI3单层的2倍. 此外, 分子动力学和声子谱的计算还证明了其良好的热稳定性和动力学稳定性. 这类可以通过合金化方法合成的磁性过渡金属卤化物将进一步拓展二维磁性材料家族及其在自旋电子学器件领域的应用.
    Two-dimensional magnetic semiconductors have received extensive attention due to their combination of magnetism, semi-conductivity and special two-dimensional structures, which also provide a new idea and platform for developing the nanometer spintronic and optoelectronic devices and also for conducting the related basic theoretical research. However, in addition to the common problems of two-dimensional magnetic semiconductor materials, such as volume manufacturing and environmental stability, the two-dimensional magnetic semiconductor materials have the unique difficulty, i.e. low Curie temperature, which makes it difficult to maintain ferromagnetic coupling at higher temperature. For example, the Curie temperature of the existing CrI3 monolayer is lower than 45 K, while that of the Cr2Ge2Te6 double-layer is only 20 K, which is far lower than the room temperature. Therefore, how to improve the Curie temperature of two-dimensional magnetic semiconductor materials through various approaches is one of the important issues that need to be resolved in this field of research. Based on the first-principles calculations, the exchange energies of a series of two-dimensional bimetallic iodides CrTMI6 (TM denotes transition metal elements in fourth and fifth rows) constructed from the lattice of CrI3 monolayer are preliminarily calculated and screened. Structures are fully relaxed until the force and the energy are converged to 0.01 eV/Å and 10–6 eV, respectively, and the ferromagnetic CrMoI6 monolayer is selected. Further calculations show that the band structure of the CrMoI6 monolayer exhibits ideal semiconductor characteristics with a band gap of about 1.7 eV. At the same time, theoretical calculations with considering the spin-orbit coupling show that the CrMoI6 monolayer has a considerable magnetic anisotropy (741.3 μeV/TM), and its easy axis is perpendicular to the two-dimensional plane. Monte Carlo simulation based on the Heisenberg model predicts that the Curie temperature of CrMoI6 monolayer reaches 92 K, which is about twice that of the CrI3 monolayer. The molecular dynamics and phonon spectrum calculations also prove that it has both thermal and kinetic stability. In addition, under the condition of applying compressive and tensile strain, its ferromagnetic coupling shows strong stability. This kind of magnetic transition metal halide which can be synthesized by alloying will further expand the family of two-dimensional magnetic materials and their applications in the field of spintronic devices.
      通信作者: 陈乾, qc119@seu.edu.cn ; 王金兰, jlwang@seu.edu.cn
    • 基金项目: 国家自然科学基金(批准号: 21973011)和国家重点研发计划(批准号: 2017YFA0403600)资助的课题
      Corresponding author: Chen Qian, qc119@seu.edu.cn ; Wang Jin-Lan, jlwang@seu.edu.cn
    • Funds: Project supported by the National Natural Science Foundation of China (Grant No. 21973011) and the National Key R&D Program of China (Grant No. 2017YFA0403600).
    [1]

    Novoselov K S, Mishchenko A, Carvalho A, Castro Neto A H 2016 Science 353 aac9439Google Scholar

    [2]

    Wang C, Xia K, Wang H, Liang X, Yin Z, Zhang Y 2019 Adv. Mater. 31 1801072Google Scholar

    [3]

    Gao Q, Zhang W, Shi Z, Yang L, Tang Y 2019 Adv. Mater. 31 1802880Google Scholar

    [4]

    Wang B, Zhang X, Zhang Y, Yuan S, Guo Y, Dong S, Wang J 2020 Mater. Horiz. 7 1623Google Scholar

    [5]

    Chen Q, Wang R, Huang Z, Yuan S, Wang H, Ma L, Wang J 2021 Nanoscale 13 6024Google Scholar

    [6]

    Huang C, Feng J, Wu F, Ahmed D, Huang B, Xiang H, Deng K, Kan E 2018 J. Am. Chem. Soc. 140 11519Google Scholar

    [7]

    Huang B, Clark G, Navarro-Moratalla E, Klein D R, Cheng R, Seyler K L, Zhong D, Schmidgall E, McGuire M A, Cobden D H 2017 Nature 546 270Google Scholar

    [8]

    Gong C, Li L, Li Z, Ji H, Stern A, Xia Y, Cao T, Bao W, Wang C, Wang Y 2017 Nature 546 265Google Scholar

    [9]

    Bonilla M, Kolekar S, Ma Y, Diaz H C, Kalappattil V, Das R, Eggers T, Gutierrez H R, Phan M H, Batzill M 2018 Nat. Nanotechnol. 13 289Google Scholar

    [10]

    Feng Y P, Shen L, Yang M, Wang A, Zeng M, Wu Q, Chintalapati S, Chang C R 2017 Wiley Interdiscip. Rev. Comput. Mol. Sci. 7 e1313Google Scholar

    [11]

    Li X, Yang J 2016 Natl. Sci. Rev. 3 365Google Scholar

    [12]

    Xue F, Hou Y, Wang Z, Wu R 2019 Phys. Rev. B 100 224429Google Scholar

    [13]

    Guan J, Huang C, Deng K, Kan E 2019 J. Phys. Chem. C 123 10114Google Scholar

    [14]

    Kresse G, Furthmüller J 1996 Phys. Rev. B 54 11169Google Scholar

    [15]

    Perdew J P, Burke K, Ernzerhof M 1996 Phys. Rev. Lett. 77 3865Google Scholar

    [16]

    Blöchl P E 1994 Phys. Rev. B 50 17953Google Scholar

    [17]

    Monkhorst H J, Pack J D 1976 Phys. Rev. B 13 5188Google Scholar

    [18]

    Joyce G 1967 Phys. Rev. 155 478Google Scholar

    [19]

    Baroni S, De Gironcoli S, Dal Corso A, Giannozzi P 2001 Rev. Mod. Phys. 73 515Google Scholar

    [20]

    Ryee S, Han M J 2018 J. Phys. Condens. Matter 30 275802Google Scholar

    [21]

    Wang L, Maxisch T, Ceder G 2006 Phys. Rev. B 73 195107Google Scholar

    [22]

    Jain A, Hautier G, Ong S P, Moore C J, Fischer C C, Persson K A, Ceder G 2011 Phys. Rev. B 84 045115Google Scholar

    [23]

    Kumagai Y, Harada K, Akamatsu H, Matsuzaki K, Oba F 2017 Phys. Rev. Appl. 8 014015Google Scholar

    [24]

    Solovyev I, Dederichs P, Anisimov V 1994 Phys. Rev. B 50 16861Google Scholar

    [25]

    Parry-Jones A, Weightman P, Andrews P 1979 J. Phys. C:Solid State Phys. 12 1587Google Scholar

    [26]

    Martyna G J, Klein M L, Tuckerman M 1992 J. Chem. Phys. 97 2635Google Scholar

    [27]

    Goodenough J B 1955 Phys. Rev. 100 564Google Scholar

    [28]

    Kanamori J 1959 J. Phys. Chem. Solids 10 87Google Scholar

    [29]

    Anderson P W 1959 Phys. Rev. 115 2Google Scholar

    [30]

    Zhang W, Qu Q, Zhu P, Lam C 2015 J. Mater. Chem. C 3 12457Google Scholar

  • 图 1  (a)二维等价合金化合物CrTMI6中TM元素的选取示意图; (b)高通量计算中考虑的单层CrTMI6的铁磁构型和两种反铁磁构型; (c)使用GGA+U方法前后单层CrTMI6材料的交换能. 其中Eex1 = EAFM1EFM, Eex2 = EAFM2EFM; E(U)ex1E(U)ex2为使用GGA+U后的相应值

    Fig. 1.  (a) Schematic diagram of the selection of TM element in the two-dimensional equivalent alloy compound CrTMI6; (b) ferromagnetic configuration and two antiferromagnetic configurations of monolayer CrTMI6 considered in high-throughput calculations; (c) exchange energy of CrTMI6 monolayer before and after GGA+U method applied, where Eex1 = EAFM1EFM, Eex2 = EAFM2EFM; E(U)ex1 and E(U)ex2 are the corresponding values using GGA+U.

    图 2  (a)单层CrMoI6可能存在的6种不同构型; (b) C1—C6各构型相对C2构型平均到每个化学式的能量差和Mo-Mo间的平均距离; (c)单层CrMoI6的声子谱; (d)时长10 ps, 温度300 K的分子动力学模拟下, 系统温度和总能量的变化以及模拟结束时的结构

    Fig. 2.  (a) Six possible configurations of CrMoI6 monolayer that may exist; (b) energy difference of C1−C6 configuration relative to C per chemical formula and the average bond distance between Mo-Mo; (c) phonon spectrum of CrMoI6 monolayer; (d) molecular dynamics simulation with duration of 10 ps, 300 K and the variation curves of system temperature and total energy, as well as the structure at the end of simulation.

    图 3  (a)单层CrMoI6的6种反铁磁构型示意图; (b)不同反铁磁构型相对于铁磁构型平均到每个化学式的能量差; (c)单层CrMoI6的能带结构示意图, 红色实线表示自旋向上通道, 蓝色实线表示自旋向下通道; (d) Cr, Mo合金体系的轨道演化及能级差示意图; (e)单层CrMoI6的自旋向上状态的分轨道能带结构, Cr的$ {\text{e}}_{\text{g}} $${\text{t}}_{\text{2g}}$轨道分别用浅蓝色和深蓝色表示, Mo的$ {\text{e}}_{\text{g}} $${\text{t}}_{\text{2g}}$轨道分别用浅红色和深红色表示, 灰色表示Cr, Mo和I原子的s轨道以及p轨道; (f)使用HSE06杂化泛函计算得到的单层CrMoI6的能带结构图

    Fig. 3.  (a) Schematic diagram of six antiferromagnetic configurations of CrMoI6 monolayerand (b) energy difference per formula of different antiferromagnetic configurations. (c) Band structure of CrMoI6 monolayer. The solid red and blue lines represent the spin-up and spin-down channels, respectively. (d) schematic diagram of orbital evolution and energy level difference of Cr and Mo alloy systems. (e) Band structure of CrMoI6 monolayer in the spin-up state. The $ {\text{e}}_{\text{g}} $ and ${\text{t}}_{\text{2g}}$ orbitals of Cr and Mo atoms are shown in light blue, dark blue, light red and dark red, respectively. The s and p orbitals of Cr, Mo, and I atoms are shown in gray. (f) Band structure of CrMoI6 monolayer calculated with HSE06 functional.

    图 4  (a) θφ的示意角度及xy轴在结构平面上的方向; (b)全空间内的磁各向异性能分布; (c)垂直于y轴方向的xz平面内的磁各向异性能分布极坐标图; (d)交换相互作用J1J6的示意图; (e)单层CrMoI6的磁矩及磁化率随温度的变化曲线; (f)单层CrI3的磁矩及磁化率随温度的变化曲线

    Fig. 4.  (a) Schematic angles of θ and φ and the direction of the xy axis on the structure plane; (b) distribution of magnetic anisotropy energy in the whole space; (c) polar coordinate diagram of magnetic anisotropic energy distribution in the xz plane perpendicular to the y axis; (d) schematic diagram of exchange interaction J1J6; The magnetic moment and susceptibility of (e) CrMoI6 monolayer and (f) CrI3 monolayer as a function of temperature.

    图 5  (a)双轴应变下CrMoI6单层的交换能以及总能量的变化; (b)施加拉伸及压缩应变时Cr—I—Cr和Mo—I—Mo的键长与键角随应变大小的变化; (c)施加拉伸及压缩应变时Cr—I—Mo的键角以及其中Cr—I和Mo—I的键长随应变大小的变化; (d)单轴应变下的交换能随应变大小的变化

    Fig. 5.  (a) Variation of exchange energy and total energy of CrMoI6 monolayer under biaxial strain; (b) variation of bond lengths and bond angles of Cr—I—Cr and Mo—I—Mo with strain under tensile and compressive strain; (c) variation of the bond angles of Cr—I—Mo and the bond lengths of Cr—I and Mo—I as a function of strain under tensile and compressive strain; (d) variation of exchange energy with strain under uniaxial strain.

    表 1  单层CrTMI6的2×2超胞的交换能

    Table 1.  Exchange energy of 2×2 supercells of CrTMI6 monolayer.

    材料交换能/eV
    Eex1Eex2E(U)ex1E(U)ex2
    CrMoI6219.40315.84108.58196.34
    CrZnI662.1462.0633.5234.48
    CrScI619.5319.4923.1323.00
    CrRhI617.9717.9425.9925.91
    CrYI610.6510.6613.9013.88
    CrTiI672.12146.71150.49–82.80
    CrRuI6–101.1566.0614.58–8.83
    CrCoI625.8425.8880.72–235.76
    CrMnI6–150.41–150.73150.10157.28
    CrAgI6–29.974.7421.69–32.62
    CrCdI6–20.45–52.78–39.279.51
    CrCuI6–52.43–132.37–7.68–65.38
    CrNiI6–194.93–231.25–51.18–45.46
    CrPdI6–313.25–359.92–24.18–67.84
    CrZrI6–290.64–385.452.73–77.70
    CrFeI6–292.42–399.20–62.40–121.16
    CrTcI6–98.88–146.83–301.14–1104.20
    CrNbI665.53337.29–698.92–681.61
    CrVI6219.46211.92–776.31–718.52
    下载: 导出CSV

    表 2  过渡金属元素TM所使用的U

    Table 2.  Value of U used for the transition metal element TM.

    材料U/eV
    V3.1[22]
    Cr3.5[22]
    Mn3.9[22]
    Fe4.0[22]
    Co3.4[22]
    Ni6.0[22]
    Cu4.0[22]
    Zn5.0[23]
    Nb2.0[24]
    Mo3.5[22]
    Tc2.3[24]
    Ru2.4[24]
    Rh2.8[24]
    Pd3.3[24]
    Ag5.1[25]
    Cd6.0[25]
    下载: 导出CSV
  • [1]

    Novoselov K S, Mishchenko A, Carvalho A, Castro Neto A H 2016 Science 353 aac9439Google Scholar

    [2]

    Wang C, Xia K, Wang H, Liang X, Yin Z, Zhang Y 2019 Adv. Mater. 31 1801072Google Scholar

    [3]

    Gao Q, Zhang W, Shi Z, Yang L, Tang Y 2019 Adv. Mater. 31 1802880Google Scholar

    [4]

    Wang B, Zhang X, Zhang Y, Yuan S, Guo Y, Dong S, Wang J 2020 Mater. Horiz. 7 1623Google Scholar

    [5]

    Chen Q, Wang R, Huang Z, Yuan S, Wang H, Ma L, Wang J 2021 Nanoscale 13 6024Google Scholar

    [6]

    Huang C, Feng J, Wu F, Ahmed D, Huang B, Xiang H, Deng K, Kan E 2018 J. Am. Chem. Soc. 140 11519Google Scholar

    [7]

    Huang B, Clark G, Navarro-Moratalla E, Klein D R, Cheng R, Seyler K L, Zhong D, Schmidgall E, McGuire M A, Cobden D H 2017 Nature 546 270Google Scholar

    [8]

    Gong C, Li L, Li Z, Ji H, Stern A, Xia Y, Cao T, Bao W, Wang C, Wang Y 2017 Nature 546 265Google Scholar

    [9]

    Bonilla M, Kolekar S, Ma Y, Diaz H C, Kalappattil V, Das R, Eggers T, Gutierrez H R, Phan M H, Batzill M 2018 Nat. Nanotechnol. 13 289Google Scholar

    [10]

    Feng Y P, Shen L, Yang M, Wang A, Zeng M, Wu Q, Chintalapati S, Chang C R 2017 Wiley Interdiscip. Rev. Comput. Mol. Sci. 7 e1313Google Scholar

    [11]

    Li X, Yang J 2016 Natl. Sci. Rev. 3 365Google Scholar

    [12]

    Xue F, Hou Y, Wang Z, Wu R 2019 Phys. Rev. B 100 224429Google Scholar

    [13]

    Guan J, Huang C, Deng K, Kan E 2019 J. Phys. Chem. C 123 10114Google Scholar

    [14]

    Kresse G, Furthmüller J 1996 Phys. Rev. B 54 11169Google Scholar

    [15]

    Perdew J P, Burke K, Ernzerhof M 1996 Phys. Rev. Lett. 77 3865Google Scholar

    [16]

    Blöchl P E 1994 Phys. Rev. B 50 17953Google Scholar

    [17]

    Monkhorst H J, Pack J D 1976 Phys. Rev. B 13 5188Google Scholar

    [18]

    Joyce G 1967 Phys. Rev. 155 478Google Scholar

    [19]

    Baroni S, De Gironcoli S, Dal Corso A, Giannozzi P 2001 Rev. Mod. Phys. 73 515Google Scholar

    [20]

    Ryee S, Han M J 2018 J. Phys. Condens. Matter 30 275802Google Scholar

    [21]

    Wang L, Maxisch T, Ceder G 2006 Phys. Rev. B 73 195107Google Scholar

    [22]

    Jain A, Hautier G, Ong S P, Moore C J, Fischer C C, Persson K A, Ceder G 2011 Phys. Rev. B 84 045115Google Scholar

    [23]

    Kumagai Y, Harada K, Akamatsu H, Matsuzaki K, Oba F 2017 Phys. Rev. Appl. 8 014015Google Scholar

    [24]

    Solovyev I, Dederichs P, Anisimov V 1994 Phys. Rev. B 50 16861Google Scholar

    [25]

    Parry-Jones A, Weightman P, Andrews P 1979 J. Phys. C:Solid State Phys. 12 1587Google Scholar

    [26]

    Martyna G J, Klein M L, Tuckerman M 1992 J. Chem. Phys. 97 2635Google Scholar

    [27]

    Goodenough J B 1955 Phys. Rev. 100 564Google Scholar

    [28]

    Kanamori J 1959 J. Phys. Chem. Solids 10 87Google Scholar

    [29]

    Anderson P W 1959 Phys. Rev. 115 2Google Scholar

    [30]

    Zhang W, Qu Q, Zhu P, Lam C 2015 J. Mater. Chem. C 3 12457Google Scholar

  • [1] 高金玮, 陈璐, 李旭洪, 史俊勤, 曹腾飞, 范晓丽. 具有磁弹耦合的本征多铁半导体: 单分子层MoTeX (X = F, Cl, Br, I). 物理学报, 2024, 73(19): 197501. doi: 10.7498/aps.73.20240829
    [2] 刘冰心, 李宗良. CrO2单层: 一种兼具高居里温度和半金属特性的二维铁磁体. 物理学报, 2024, 73(10): 106102. doi: 10.7498/aps.73.20240246
    [3] 谢玲凤, 董金瓯, 赵雪芹, 杨巧林, 宁凡龙. In掺杂对磁性半导体Li1.05(Zn0.925, Mn0.075)As中铁磁序的调控. 物理学报, 2024, 73(8): 087501. doi: 10.7498/aps.73.20231949
    [4] 樊晓筝, 李怡莲, 吴怡, 陈俊彩, 徐国亮, 安义鹏. 二维磁性半导体笼目晶格Nb3Cl8单层的磁性及自旋电子输运性质. 物理学报, 2023, 72(24): 247503. doi: 10.7498/aps.72.20231163
    [5] 孙敬淇, 吴绪才, 阙志雄, 张卫兵. 基于材料组分信息的高居里温度铁磁材料预测. 物理学报, 2023, 72(18): 180202. doi: 10.7498/aps.72.20230382
    [6] 宋蕊, 王必利, 冯凯, 王黎, 梁丹丹. 二维VOBr2单层的结构畸变及其磁性和铁电性. 物理学报, 2022, 71(3): 037101. doi: 10.7498/aps.71.20211516
    [7] 宋蕊. 二维VOBr2单层的结构畸变及其磁性和铁电性研究. 物理学报, 2021, (): . doi: 10.7498/aps.70.20211516
    [8] 张浩杰, 张茹菲, 傅立承, 顾轶伦, 智国翔, 董金瓯, 赵雪芹, 宁凡龙. 一种具有“1111”型结构的新型稀磁半导体(La1–xSrx)(Zn1–xMnx)SbO. 物理学报, 2021, 70(10): 107501. doi: 10.7498/aps.70.20201966
    [9] 蒋小红, 秦泗晨, 幸子越, 邹星宇, 邓一帆, 王伟, 王琳. 二维磁性材料的物性研究及性能调控. 物理学报, 2021, 70(12): 127801. doi: 10.7498/aps.70.20202146
    [10] 刘雨亭, 贺文宇, 刘军伟, 邵启明. 二维材料中贝里曲率诱导的磁性响应. 物理学报, 2021, 70(12): 127303. doi: 10.7498/aps.70.20202132
    [11] 王海宇, 刘英杰, 寻璐璐, 李竞, 杨晴, 田祺云, 聂天晓, 赵巍胜. 大面积二维磁性材料的制备及居里温度调控. 物理学报, 2021, 70(12): 127301. doi: 10.7498/aps.70.20210223
    [12] 文林, 胡爱元. 双二次交换作用和各向异性对反铁磁体相变温度的影响. 物理学报, 2020, 69(10): 107501. doi: 10.7498/aps.69.20200077
    [13] 魏钟鸣, 夏建白. 低维半导体偏振光探测器研究进展. 物理学报, 2019, 68(16): 163201. doi: 10.7498/aps.68.20191002
    [14] 肖嘉星, 鲁军, 朱礼军, 赵建华. 垂直磁各向异性L10-Mn1.67Ga超薄膜分子束外延生长与磁性研究. 物理学报, 2016, 65(11): 118105. doi: 10.7498/aps.65.118105
    [15] 邹文琴, 路忠林, 王申, 刘圆, 陆路, 郦莉, 张凤鸣, 都有为. Mn和N共掺ZnO稀磁半导体薄膜的研究. 物理学报, 2009, 58(8): 5763-5767. doi: 10.7498/aps.58.5763
    [16] 张继业, 骆 军, 梁敬魁, 纪丽娜, 刘延辉, 李静波, 饶光辉. 赝二元固溶体TbGa1-xGex(0≤x≤0.4)的结构与磁性. 物理学报, 2008, 57(10): 6482-6487. doi: 10.7498/aps.57.6482
    [17] 王叶安, 秦福文, 吴东江, 吴爱民, 徐 茵, 顾 彪. 基于电子回旋共振-等离子体增强金属有机物化学气相沉积技术生长GaMnN稀磁半导体的研究. 物理学报, 2008, 57(1): 508-513. doi: 10.7498/aps.57.508
    [18] 刘喜斌, 沈保根. Mn5Ge2.7M0.3 (M=Ga,Al,Sn) 化合物的磁性和磁熵变. 物理学报, 2005, 54(12): 5884-5889. doi: 10.7498/aps.54.5884
    [19] 李锐鹏, 王 劼, 李红红, 郭玉献, 王 锋, 胡志伟. 软x射线磁性圆二色吸收谱研究铁单晶薄膜的面内磁各向异性. 物理学报, 2005, 54(8): 3851-3855. doi: 10.7498/aps.54.3851
    [20] 陈伟, 钟伟, 潘成福, 常虹, 都有为. La0.8-xCa0.2MnO3纳米颗粒的居里温度与磁热效应. 物理学报, 2001, 50(2): 319-323. doi: 10.7498/aps.50.319
计量
  • 文章访问数:  6546
  • PDF下载量:  240
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-05-19
  • 修回日期:  2021-06-16
  • 上网日期:  2021-10-13
  • 刊出日期:  2021-10-20

/

返回文章
返回