化学水浴法合成硫化铅探测器的红外响应研究

杨冬, 申钧, 高恺聪, 冷重钱, 聂长斌, 张之胜

杨冬, 申钧, 高恺聪, 冷重钱, 聂长斌, 张之胜. 化学水浴法合成硫化铅探测器的红外响应研究[J]. 红外技术, 2023, 45(6): 559-566.
引用本文: 杨冬, 申钧, 高恺聪, 冷重钱, 聂长斌, 张之胜. 化学水浴法合成硫化铅探测器的红外响应研究[J]. 红外技术, 2023, 45(6): 559-566.
YANG Dong, SHEN Jun, GAO Kaicong, LENG Chongqian, NIE Changbin, ZHANG Zhisheng. Infrared Response of Lead Sulfide Detector Synthesized from Chemical Bath Deposition[J]. Infrared Technology , 2023, 45(6): 559-566.
Citation: YANG Dong, SHEN Jun, GAO Kaicong, LENG Chongqian, NIE Changbin, ZHANG Zhisheng. Infrared Response of Lead Sulfide Detector Synthesized from Chemical Bath Deposition[J]. Infrared Technology , 2023, 45(6): 559-566.

化学水浴法合成硫化铅探测器的红外响应研究

基金项目: 

北京市工程技术研究中心基金项目 BG0231

重庆市自然科学基金面上项目 cstc2021jcyj-msxmX0389

详细信息
    作者简介:

    杨冬(1996-),男,硕士研究生,主要研究方向为碳基光电探测器。E-mail: yangdong20@mails.ucas.ac.cn

    通讯作者:

    张之胜(1978-),男,工程师,主要研究方向为碳基电子和光电器件。E-mail: zhangzhisheng@cigit.ac.cn

  • 中图分类号: TN362

Infrared Response of Lead Sulfide Detector Synthesized from Chemical Bath Deposition

  • 摘要: 硫化铅探测器具有短波红外高灵敏度、低俄歇噪声等优点,其中化学水浴法合成的硫化铅薄膜可与CMOS半导体工艺兼容,有利于实现低成本高性能的面阵探测器。然而,目前对化学水浴法合成硫化铅探测器的研究主要集中在较大尺寸的单元探测器。本文基于化学水浴法合成硫化铅薄膜,利用离子束刻蚀工艺,制备了10~200 μm尺寸的硫化铅光电探测器,研究了器件光电性能随电阻、长宽比、线宽等参数的变化。结果表明,随着尺寸的减小,硫化铅光电探测器的响应度逐渐增加,在1550 nm短波红外光的照射下,10 μm级器件的响应度达到了51.68 A/W,约为200 μm级器件的123倍,且在可见光和2.7 μm红外波长下也具有良好的宽波段光电响应。本文研究的微米尺寸探测器件可为硫化铅探测器研究提供一定的支撑。
    Abstract: A lead sulfide detector has the advantages of high short-wave infrared sensitivity and low auger noise. The lead sulfide film synthesized by chemical bath deposition can be compatible with the CMOS semiconductor process, which is beneficial for identifying low-cost and high-performance surface array detectors. However, the current research on lead sulfide detectors synthesized by chemical bath deposition primarily focuses on the larger unit detectors. In this study, the synthesis of lead sulfide film was based on chemical bath deposition. Lead sulfide photodetectors of 10–200 μm size were prepared using an ion beam etching process and the photoelectric performance of the device was studied in terms of resistance, aspect ratio, line width, and other parameters. The results showed that the responsivity of the PbS photodetector increased gradually as the size decreased. Under the irradiation of 1550 nm short-wave infrared light, the responsivity of the 10 µm photodetector was 51.68 A/W, which was approximately 123 times the responsivity of the 200 μm photodetector. Moreover, the PbS photodetector also had a good wide-range photoelectric response at visible light and 2.7 μm infrared wavelengths. The micron-size detector of this study can provide support for the research of lead sulfide detectors.
  • 随着遥感器技术的进步,空间遥感器相机向着多谱段、大幅宽、实时传输等方向不断发展。获取多个谱段的景物信息对资源普查、矿物勘测、农作物监测等具有极高的应用前景[1]。多光谱遥感相机具有较宽谱段的成像能力,美国于2008年发射Geoeye-1高分辨率卫星,运行于高度681 km,轨道倾角98°的太阳同步轨道。其全色谱段地面像元分辨率为0.41 m,包括4个经典波段(红、绿、蓝和近红外)多光谱的地面像元分辨率为1.65 m[2]。国内于2020年成功发射资源三号03星,运行于轨道高度506 km的太阳同步轨道,地面分辨率下视2.1 m,前后视2.5 m,多光谱5.8 m,覆盖宽度单载荷大于50 km,立体成像大于45 km[3]

    本文所介绍的空间相机运行于近地圆轨道,相比于太阳同步轨道,全寿命周期内β角(轨道面与太阳光矢量之间的夹角)变化很大,在-58.4°~58.4°之间,相机安装于卫星舱体外,造成相机每个面都有可能受照[4]。相机在轨可实现全天时工作,工作期间姿态多变且伴随着大角度的俯仰与侧摆,轨道外热流非常复杂,热流扰动大。相机拥有可见光与红外双成像通道,在外热流变化复杂的空间环境下,相机需要对以下方面开展散热设计:可见光焦面需要在较低的温度水平下工作,以减小暗电流和热噪声对成像质量的影响[5];低温红外焦面所配备的压缩制冷机功耗大,且工作时间长。为保证其在轨期间正常进行制冷工作,需将其产生的热量快速排散。

    为满足该空间相机内热源的散热要求,本文对该相机所受到的外热流进行分析并根据空间环境进行了散热设计。针对热设计开展热仿真计算,验证了相机散热设计的正确性与合理性。

    该相机采用同轴三反光学系统,主体构型示意图与主要内热源分布情况见图 1,内热源由可见光组件、红外组件和制冷机组件等组成。

    图  1  相机主要结构及内热源分布示意图
    Figure  1.  Schematic diagram of the main structure of the camera and the distribution of internal heat sources

    相机内热源中红外组件与制冷机组件需要长期工作。其中制冷机组件包含制冷机热端部分与制冷机压缩机部分,压缩机部分安装在相机本体外,固定在卫星舱板上。相机工作热耗219 W,待机热耗为110 W,温控指标如表 1所示。

    表  1  内热源温控指标
    Table  1.  Temperature index of the internal heat source 
    Component Operating temperature
    Visible components -5~40
    Infrared components -25~55
    Refrigerator compressor -15~30
    refrigerator hot end ≤20
    下载: 导出CSV 
    | 显示表格

    该相机运行于近地圆轨道,β角在-58.4°~58.4°之间变化。相机在轨运行姿态为:阳照区非工作时段-Z对日,工作期间调整为+Z对地,阴影区长期+Z对地,+Z向为相机光线入射方向的反向,+X向为速度矢量方向,+Y向由右手坐标系确定。其中,在阳照区成像过程中相机入光口方向需要与太阳光矢量方向成最小30°规避角,避免太阳光入射相机入光口影响成像。

    根据相机飞行轨道参数分析轨道外热流大小和变化规律,是确定相机散热面方位和定义极端工况的重要输入条件[6]。利用热分析软件计算出相机在冬至时刻运行姿态下随β角变化的平均到达外热流,如图 2所示。

    图  2  平均到达外热流随β角变化趋势
    Figure  2.  Available average heat flows with different β

    相机+Z面、+X面和-X面平均到达外热流变化趋势相同,都随着β角绝对值增大呈下降的趋势。-Y面平均到达外热流随β角增大呈下降趋势,+Y面外热流变化与-Y面相反,±Y面平均到达外热流以β角0°呈对称趋势,且变化较大。

    相机成像过程中,由于规避角影响,使相机±X和±Y面受到太阳直射,相较于正常姿态其外热流会发生较大变化。由于相机在不同方向摆动时,-X面与+X面,-Y面与+Y面到达外热流变化趋势相同,故以单向摆动为例,在极限β角(0°与58.4°)下计算±X与±Y面平均到达外热流,如图 3图 4

    图  3  β角0°时±X、±Y平均到达外热流
    Figure  3.  Available average heat flow when β=0° Time/s
    图  4  β角58.4°时±X、±Y平均到达外热流
    Figure  4.  Available average heat flow when β=60°

    ±Y面在β角为58.4°时受阳光直射时间长,平均到达外热流大,相较于β角为0°下的±X面,热环境恶劣,不利于内热源的散热设计。而±X面虽然均会受到阳光直射,但在相机在轨运行中存在单一侧面长时间背阳的情况,相较于±Y面,到达外热流小,使内热源散热压力较小。根据相机结构布局与到达外热流情况,将相机散热面设置在±X侧。

    针对以上相机情况,结合相机所处空间环境与内热源特性,总结以下散热设计难点:

    相机所处空间环境外热流变化复杂。在相机处于工作姿态时,除入光口外相机各面均会受到阳光直射。由于相机在轨期间姿态多变,使其受到的外热流变化复杂。

    相机具有可见光与红外双成像通道,结构复杂,制冷机压缩机部分位于相机本体外,内热源位置分散。部分内热源在轨期间处于长期工作状态,发热量大,且相机所处空间环境复杂,不利于大功率热源散热。

    针对以上相机散热设计的难点,同时兼顾相机内热源的热排散与热补偿,采用被动热控和主动热控相结合的热设计方案。

    为保证相机内热源产生的热量有效排散,需要对其进行被动热控设计。由于无法在相机本体上开设散热面,热设计中额外设置了辐射冷板进行散热。辐射冷板表面喷涂KS-ZA白漆来增大散热效果,内热源通过槽道热管将热量传导至辐射冷板,并最终辐射至外部空间。冷板与卫星舱体通过聚酰亚胺隔热垫与钛合金螺钉进行隔热安装,冷板与热管,热管与热管工装面之间填充导热填料,增大热量传导效率[6]

    相机各部分内热源工作模式与工作指标各不相同,需要设置多个辐射冷板进行散热设计。可见光组件热耗65 W,红外组件总热耗24 W,二者温度指标相近,且可见光组件单轨工作60 s,散热压力较小,使用一块辐射冷板对可见光组件与红外组件同时散热。可见光组件通过两根热管与辐射冷板相连,三部分红外组件各通过一根热管与辐射冷板相连,辐射冷板设置在相机-X侧。制冷机压缩机部分长期工作,其待机功耗为38.5 W,成像工作时功耗为45.5 W。使用一块辐射冷板对压缩机散热,压缩机通过一根热管与辐射冷板相连,辐射冷板设置在相机+X侧。制冷机热端部分是相机内部最大热源,其待机功耗为71.5 W,成像工作时功耗为84.5 W。由于制冷机热端长期工作,且功耗大,而相机在成像过程中±X面均有可能受到太阳直射,单一辐射冷板不利于热量排散,并且由于单侧设置辐射冷板受到的外热流大,使辐射冷板面积过大,不利于相机的结构设计。采用双辐射冷板耦合散热的方案,可以利用成像过程中相机一侧背阳,到达外热流较小的辐射冷板提升制冷机热端的散热效率,同时减小了辐射冷板面积,有利于相机的结构设计。辐射冷板耦合设计如图 5所示,制冷机热端通过两根热管将热量分别引至±X侧辐射冷板,两侧的辐射冷板通过一根U型热管连接。

    图  5  制冷机热端散热方案示意图
    Figure  5.  Schematic diagram of the heat dissipation for refrigerator hotend

    为使短期工作热源与长期工作热源均满足温度指标与合理规划热控资源,需要精确计算辐射冷板面积并代入仿真计算进行校核,最终确定辐射冷板的面积[6]。根据外热流与辐射冷板涂层属性,计算辐射冷板吸收热量Q1为:

    $$ Q_1=Q_{\text {内 }}+\alpha \times\left(q_1+q_2\right) \times A+\varepsilon \times q_3 \times A $$ (1)

    式中:Q为内热源功耗;α为辐射冷板太阳吸收率;q1为太阳直射的到达外热流密度;q2为阳光反照的到达外热流密度;ε为辐射冷板表面发射率;q3为地球红外的到达外热流密度;A为辐射冷板面积。

    辐射冷板散出的热量Q2为:

    $$ Q_2=A \times \varepsilon \times \sigma \times T^4 $$ (2)

    式中:σ为斯特藩-玻尔兹曼常数,其值为5.67×10-8 W/(m2·K4)[7]

    当热平衡时,有:

    $$ Q_1=Q_2 $$ (3)

    通过仿真计算得出,可见光组件与红外组件辐射冷板面积为0.19 m2,制冷机压缩机辐射冷板面积为0.26 m2,制冷机热端辐射冷板两面各0.25 m2,总面积0.5 m2

    主动热控是在被动热控的基础上进一步提高并维持内热源温度水平的热设计方法,采用补偿加热的方式控制内热源温度,使其在低温工况下或待机状态下也满足温控指标[8-9]。补偿加热回路采用趋势控温法进行闭环控温,当温度低于温度指标时,回路连接,加热片工作,当温度达到温度指标时,回路断开,加热片停止工作。加热回路设置在各个辐射冷板上,保证内热源温度指标的同时也可以确保热管处于正常工作温区内。加热片使用聚酰亚胺薄膜型加热片,热敏电阻测温,由相机下位机采集控温。

    可见光组件内热源包含两部分,其中一部分通过其安装框架将热量传导至电箱基板,通过安装在电箱基板上的热管将热量传导至辐射冷板,最终辐射至外部空间。另一部分通过一个铝扩热板将热量集中传导至电箱外壳,通过热管将热量传导至辐射冷板并辐射至外部空间。红外组件上的发热元件与红外组件电箱外壳直接接触,接触界面填充导热填料增强导热,电箱外壳通过热管将热量传导至辐射冷板,最终将热量辐射至外部空间。制冷机热端产生的热耗直接由安装在其表面的两根热管传导至相机两侧辐射冷板,两侧的辐射冷板同时将热量辐射至外部空间。制冷机压缩机安装在相机本体外,通过铝合金安装架固定在卫星舱板上。压缩机产生的热耗传导至压缩机上压盖,接触面之间填充导热填料增强导热,通过安装在上压盖上的热管传导至辐射冷板,并辐射至外部空间。相机各内热源传热路径如图 6所示,图中点划线为传导导热路径,实线为热管传热,虚线为辐射传热路径。

    图  6  内热源传热路径
    Figure  6.  Heat transfer path of the internal heat source

    导热热阻分为材料自身的传导热阻与接触面之间的接触传热热阻,计算公式分别为:

    $$ {R_{\text{c}}} = \frac{\delta }{{{A_{\text{c}}} \cdot \lambda }} $$ (4)
    $$ {R_{{\text{ct}}}} = \frac{1}{{K \cdot {A_{{\text{ct}}}}}} $$ (5)

    式中:Rc为材料自身的传导热阻;δ为热量传导方向的长度;Ac为横截面面积;λ为材料导热率;Rct为接触传热热阻;K为接触传热系数;Act为接触面积[10]

    根据上述热设计方案,利用有限元软件对相机进行建模,如图 7所示。壳单元划分,单元厚度按照等效厚度计算,采用热耦合的方式对模型简化,共建立了209个热耦合,划分了16437个单元。

    图  7  相机有限元模型
    Figure  7.  Finite element model of camera

    仿真计算中根据卫星太阳位置、β角和散热面涂层表面属性定义高温与低温两个极端工况,如表 2所示。

    表  2  相机工况定义
    Table  2.  Definition of camera working conditions
    Solar position β Heat sink surface coating surface properties
    High temperature working condition Winter solstice α/ε=0.23/0.92(End period)
    Low temperature working condition Summer Solstice 60° α/ε=0.13/0.92(Initial period)
    下载: 导出CSV 
    | 显示表格

    根据工况定义,利用仿真软件计算相机在各工况下的温度结果,单轨结束时热源温度分布云图与热源温度变化曲线如图 89所示。

    图  8  高温工况温度分布云图与温度变化曲线
    Figure  8.  High temperature distribution clouds and temperature variation curves
    图  9  低温工况温度分布云图与温度变化曲线
    Figure  9.  Low temperature distribution clouds and temperature variation curves

    图 89可以看出,相机在温度稳定后高温工况下可见光组件温度为-1.1℃~12.9℃,红外组件温度为3℃~10.5℃,制冷机热端温度为-5℃~0.3℃,制冷机压缩机温度为-1℃~21.3℃。相机在低温工况下可见光组件温度为-1.9℃~9℃,红外组件温度为1.7℃~6.7℃,制冷机热端温度为-12℃~-10.8℃,制冷机压缩机温度为-11.3℃~5.4℃。均满足热控指标要求。

    本文为解决相机在外热流复杂的空间环境下相机大功率热源的散热问题,通过分析相机各个侧面到达外热流情况与相机内热源工作模式,采用在相机两侧设置辐射冷板的散热方式,并使用热管将两侧辐射冷板耦合的方案,增大热量的排散效率,减小了辐射冷板的面积。最终,通过热仿真计算验证了设计的合理性。仿真结果表面:可见光组件温度为-1.9℃~12.9℃,红外组件温度为1.7℃~10.5℃,制冷机热端温度为-12℃~0.3℃,制冷机压缩机温度为-11.3℃~21.3℃,均满足温控指标要求,解决了复杂外热流下相机大功率热源的散热问题。

  • 图  1   PbS光电探测器三维结构示意图与工作原理图及制备过程

    Figure  1.   Three-dimensional structure, working principle and preparation process of PbS photodetector

    图  2   PbS光电探测器的微观形貌图

    Figure  2.   Microstructure images of PbS photodetector

    图  3   不同尺寸的器件的光电流和响应度与波长之间的关系

    Figure  3.   Photocurrent and responsivity of photodetectors with different sizes versus wavelength

    图  4   1550 nm光源下的光电流响应曲线

    Figure  4.   Photocurrent response curves under 1550 nm light source

    图  5   探测器件电阻与尺寸的关系

    Figure  5.   Resistance versus size of detectors

    图  6   响应度和光电流与电阻的关系

    Figure  6.   Responsivity and photocurrent versus resistance

    图  7   不同尺寸探测器在1550 nm入射光下的光电性能与光功率的关系

    Figure  7.   Photoelectric properties of photodetectors with different sizes versus light power under incident light at 1550 nm

    图  8   不同尺寸探测器在2.7 μm入射光下的光电性能与光功率的关系

    Figure  8.   Photoelectric properties of photodetectors with different sizes versus light power under incident light at 2.7 μm

    图  9   探测器的光电性能与尺寸之间的关系

    Figure  9.   Photoelectric properties of photodetectors versus sizes

    表  1   CBD法制备PbS所需原料

    Table  1   The raw materials required for PbS prepared by CBD method

    Name of raw material Chemical formula Use of raw material
    Lead acetate Pb(CH3COO)2 The reactant, which provides the lead source
    Sodium hydroxide NaOH The reactants, which provides an alkaline environment
    Thiourea SC(NH2)2 The reactants, which provides the sulfur source
    Sodium citrate Na3C6H5O7 The complexing agent, which controls reaction rate
    下载: 导出CSV
  • [1] 胡伟达, 李庆, 陈效双, 等. 具有变革性特征的红外光电探测器[J], 物理学报, 2019, 68(12): 7-41. https://www.cnki.com.cn/Article/CJFDTOTAL-WLXB201912001.htm

    HU W D, LI Q, CHEN X S, et al. Recent progress on advanced infrared photodetectors[J]. Acta Physica Sinica, 2019, 68(12): 7-41. https://www.cnki.com.cn/Article/CJFDTOTAL-WLXB201912001.htm

    [2]

    SHI L L, LIANG Q B, WANG W Y, et al. Research progress in organic photomultiplication photodetectors[J]. Nanomaterials (Basel), 2018, 8(9): 713. DOI: 10.3390/nano8090713

    [3] 李国伟. PbS红外探测器的制备和性能研究[D]. 成都: 电子科技大学, 2011.

    LI G W. Study on Preparation and Performance of PbS Infrared Detector[D]. Chengdu: University of Electronic Science and Technology of China, 2011.

    [4] 丁瑞军, 杨建荣, 何力, 等. 碲镉汞红外焦平面器件技术进展[J]. 红外与激光工程, 2020, 49(1): 93-99. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ202001010.htm

    DING R J, YANG J R, HE L, et al. Development of technologies for HgCdTe IRFPA[J]. Infrared and Laser Engineering, 2020, 49(1): 93-99. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ202001010.htm

    [5] 李雪, 邵秀梅, 李淘, 等. 短波红外InGaAs焦平面探测器研究进展[J]. 红外与激光工程, 2020, 49(1): 64-71. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ202001006.htm

    LI X, SHAO X M, LI T, et al. Developments of short-wave infrared InGaAs focal plane detectors[J]. Infrared and Laser Engineering, 2020, 49(1): 64-71. https://www.cnki.com.cn/Article/CJFDTOTAL-HWYJ202001006.htm

    [6] 周连军, 韩福忠, 白丕绩, 等. 高温碲镉汞中波红外探测器的国内外进展[J]. 红外技术, 2017, 39(2): 116-124. http://hwjs.nvir.cn/article/id/hwjs201702002

    ZHOU L J, HAN F Z, BAI P J, et al. Review of HOT MW infrared detector using MCT technology[J]. Infrared Technology, 2017, 39(2): 116-124. http://hwjs.nvir.cn/article/id/hwjs201702002

    [7]

    DIEZHANDINO J, VERGARA G, PEREZ G, et al. Monolithic integration of spectrally selective uncooled lead selenide detectors for low cost applications[J]. Applied Physics Letters, 2003, 83(14): 2751-2753. DOI: 10.1063/1.1615314

    [8]

    VERGARA G, LINARES HERRERO R, GUTíERREZ ÁLVAREZ R, et al. 80×80 VPD PbSe: the first uncooled MWIR FPA monolithically integrated with a Si-CMOS ROIC[C]//Pro. of SPIE, 2013, 8704: 87041M.

    [9]

    WENG B B, QIU J J, ZHAO L H, et al. Recent development on the uncooled mid-infrared PbSe detectors with high detectivity[C]// Pro. of SPIE, 2013, 8993: 899311.

    [10]

    KONSTANTATOS G, CLIFFORD J, LEVINA L, et al. Sensitive solution-processed visible-wavelength photodetectors[J]. Nature Photonics, 2007, 1(9): 531-534. DOI: 10.1038/nphoton.2007.147

    [11]

    MI L F, CHANG Y J, ZHANG Y, et al. Hybrid perovskite exchange of PbS quantum dots for fast and high-detectivity visible-near-infrared photo-detectors[J]. Journal of Materials Chemistry C, 2020, 8(23): 7812-7819. DOI: 10.1039/D0TC01373K

    [12]

    YAN S Y, YANG Q, FENG S L, et al. Effect of air atmosphere sensitization on formation of PbSe p-n junctions for high-performance photodetectors[J]. Journal of Electronic Materials, 2020, 49(8): 4929-4935. DOI: 10.1007/s11664-020-08215-6

    [13]

    DONG R, BI C, DONG Q F, et al. An Ultraviolet-to-NIR broad spectral nanocomposite photodetector with gain[J]. Advanced Optical Materials, 2014, 2(6): 549-554. DOI: 10.1002/adom.201400023

    [14]

    KUFER D, LASANTA T, BERNECHEA M, et al. Interface engineering in hybrid quantum dot-2D phototransistors[J]. Acs Photonics, 2016, 3(7): 1324-1330.

    [15]

    SONG X X, ZHANG Y T, ZHANG H T, et al. Graphene and PbS quantum dot hybrid vertical phototransistor[J]. Nanotechnology, 2017, 28(14): 145201.

    [16]

    KONSTANTATOS G. SARGENT E H. PbS colloidal quantum dot photoconductive photodetectors: Transport, traps, and gain[J]. Applied Physics Letters, 2007, 91(17): 173505.

    [17]

    LAUER R B. WILLIAMS F. Photoelectronic properties of graded composition crystals of Ii-Vi semiconductors[J]. Journal of Applied Physics, 1971, 42(7): 2904-2910.

    [18]

    BUSCEMA M, ISLAND J O, GROENENDIJK D J, et al. Photocurrent generation with two-dimensional van der Waals semiconductors[J]. Chemical Society Reviews, 2015, 44(11): 3691-3718.

    [19]

    KONSTANTATOS G, BADIOLI M, GAUDREAU L, et al. Hybrid graphene-quantum dot phototransistors with ultrahigh gain[J]. Nat Nanotechnol, 2012, 7(6): 363-368.

    [20]

    KONSTANTATOS G. SARGENT E H. Solution-processed quantum dot photodetectors[C]//Proceedings of IEEE, 2009, 97(10): 1666-1683.

图(9)  /  表(1)
计量
  • 文章访问数:  350
  • HTML全文浏览量:  52
  • PDF下载量:  65
  • 被引次数: 0
出版历程
  • 收稿日期:  2022-05-10
  • 修回日期:  2022-06-06
  • 刊出日期:  2023-06-19

目录

/

返回文章
返回