文献中电流密度约10e-8,但我的结果约10e-13，差别太大，请问该如何解决

Question

go atlas

mesh 
x.mesh loc=0.0    spac=0.2 
x.mesh loc=80.0  spac=0.2  
x.mesh loc=100.0  spac=0.1  
x.mesh loc=110.0    spac=0.2 
x.mesh loc=121.0  spac=2.0  

# Y 方向从上到下：Au(100nm) → Ni(30nm) → 外延9um → 衬底650um → Ti(30nm) → Au(80nm)
y.mesh loc=0.00    spac=0.05
y.mesh loc=0.10    spac=0.05   
y.mesh loc=0.13    spac=0.02   
y.mesh loc=0.23    spac=0.1
y.mesh loc=9.13    spac=0.2    
y.mesh loc=100.0  spac=20.0
y.mesh loc=659.24  spac=50.0     
 
#-----------------------------------------------------------
# 区域顺序 从上到下：1 → 2 → 3 → 4 → 5 → 6
#-----------------------------------------------------------
region num=1 material=air x.min=100 x.max=121 y.min=0.000 y.max=0.130
region num=2 user.material=Au x.min=0 x.max=100 y.min=0.000 y.max=0.100    
region num=3 material=Ni x.min=0 x.max=100 y.min=0.100 y.max=0.130
region num=4 user.material=Ga2O3 x.min=0 x.max=121 y.min=0.13 y.max=9.13  
region num=5 user.material=Ga2O3 x.min=0 x.max=121 y.min=9.13 y.max=659.13 
region num=6 material=Ti x.min=0 x.max=121 y.min=659.13 y.max=659.16
region num=7 user.material=Au x.min=0 x.max=121 y.min=659.16 y.max=659.24

# 阳极：绑定Ni/Au电极（y=0~0.13μm，x=0~100μm），电极编号1
electrode name=anode x.min=0 x.max=100 y.min=0 y.max=0.13
# 阴极：绑定Ti/Au背电极（y=659.13~659.24μm，x=0~121μm），电极编号2
electrode name=cathode x.min=0 x.max=121 y.min=659.13 y.max=659.24
# PiN结构：P+阳极层 → N-外延层 → N+衬底

# N-外延层（y=0.23~9.13μm，x=0~121μm）
doping uniform conc=1e16 n.type region=4  
# N+衬底（y=9.13~659.13μm，x=0~121μm）
doping uniform conc=1e19 n.type region=5   

material material=Ga2O3  user.group=semiconductor \
         eg300=4.8 permittivity=10 mun=150 mup=10 nc300=3.72e18 \
         nv300=1.16e19 taun=1e-8 taup=1e-8 affinity=4.0 ni=1e-5
material material=Au user.default=gold

# 物理模型
models fermi conmob fldmob srh auger thermal lat.temp 
impact  selb
# 接触类型
contact name=anode   workfunction=5.15 
contact name=cathode    

# 温度接触
thermcontact num=3 name=anode temp=300
thermcontact num=6 name=cathode temp=300

# 第九步：求解方法（保证收敛）
method newton  maxtrap=40  climit=1e-5 
solve init
save outfile=Ga2O3_SBD_ordered.str
log outfile=IV.log
# 阶段1：0V → -1V，-0.05V/步，极致小步长建立肖特基势垒
solve vanode=0.0 vstep=-0.5 vfinal=-5 name=anode   
solve vanode=-5 vstep=-1 vfinal=-20 name=anode
log off
tonyplot IV.log

quit
![img](https://i-blog.csdnimg.cn/ask/60c8b6a2a39046ddbb2a218ce24c3c41.png "#left")

Answer 1

下午好🌅🌅🌅
本答案参考通义千问

你提供的ATLAS仿真脚本中，电流密度与文献值（约10e-8 A/cm²）相差很大，你的结果是约10e-13 A/cm²，这表明电流密度过小，可能由以下原因导致：

---

一、问题分析

1. 材料参数设置不正确

• Ga2O3 的载流子迁移率（mun=150, mup=10）：这个数值是否合理？在Ga2O3中，电子迁移率通常更高，比如1000 cm²/V·s左右。若迁移率太低，会导致电流密度下降。
• 掺杂浓度（N型外延层为1e16 cm⁻³，衬底为1e19 cm⁻³）：是否符合实际？如果掺杂浓度过低，也会导致电流不足。

2. 接触电极设置不合理

• 阳极接触（Au/Ni）的功函数（workfunction=5.15 eV）：是否匹配Ga2O3的导带位置？若肖特基势垒过高，会显著抑制电流。
• 阴极接触（Ti/Au）：是否为欧姆接触？如果是肖特基接触，也可能影响电流。

3. 网格划分不合理

• Y方向网格过于稀疏（如从0.13到9.13只用了0.2间距），可能导致电场和载流子分布计算不准确，尤其是在界面附近。
• X方向网格在100~121μm处过密（spac=2.0），可能会导致边界条件处理不当。

4. 求解器设置问题

• solve init 和 solve vanode：是否设置了足够的迭代次数或步长？如果收敛不够，结果可能不准。
• log outfile=IV.log：是否记录了正确的电流数据？

---

二、解决方案（有序列表）

1. 调整材料参数

• 修改 Ga2O3 的迁移率：

  material material=Ga2O3 user.group=semiconductor \
           eg300=4.8 permittivity=10 mun=1000 mup=100 nc300=3.72e18 \
           nv300=1.16e19 taun=1e-8 taup=1e-8 affinity=4.0 ni=1e-5
  

  重点：将电子迁移率 mun 提高到 1000 cm²/V·s，提高电流密度。

• 调整掺杂浓度：

  doping uniform conc=1e18 n.type region=4  # 外延层掺杂浓度提高
  doping uniform conc=1e20 n.type region=5  # 衬底掺杂浓度提高
  

2. 优化电极接触

• 调整阳极功函数（参考 Ga2O3 的导带位置）：

  contact name=anode workfunction=4.0   # 假设 Ga2O3 导带为 4.0 eV
  

  重点：确保 Au/Ni 接触的功函数低于 Ga2O3 的导带，形成低势垒。

• 确认阴极是否为欧姆接触：

  contact name=cathode type=ohmic
  

3. 优化网格划分

• 细化 Y 方向网格（尤其是界面区域）：

  y.mesh loc=0.00 spac=0.05
  y.mesh loc=0.10 spac=0.05
  y.mesh loc=0.13 spac=0.02   # 细化 Ni/Ga2O3 界面
  y.mesh loc=0.23 spac=0.05   # 细化 Ga2O3 外延层
  y.mesh loc=9.13 spac=0.1    # 外延层至衬底
  y.mesh loc=659.13 spac=0.1  # 衬底至 Ti/Au
  y.mesh loc=659.24 spac=0.05 # 背电极
  

• X 方向适当增加分辨率（特别是在 100~121 μm 区域）：

  x.mesh loc=100.0 spac=0.1
  x.mesh loc=121.0 spac=0.1
  

4. 优化求解器设置

• 增加 Newton 迭代次数：

  method newton maxtrap=100 climit=1e-6
  

• 使用更小的电压步长（特别是负压阶段）：

  solve vanode=0.0 vstep=-0.1 vfinal=-5 name=anode
  solve vanode=-5 vstep=-0.1 vfinal=-20 name=anode
  

5. 验证电流输出

• 确保 log 文件记录了电流数据：

  log outfile=IV.log
  

---

三、修改后的代码示例

mesh 
x.mesh loc=0.0    spac=0.2 
x.mesh loc=80.0  spac=0.2  
x.mesh loc=100.0  spac=0.1  
x.mesh loc=100.0  spac=0.1   # 细化阳极区域
x.mesh loc=121.0  spac=0.1   # 细化背电极区域

y.mesh loc=0.00    spac=0.05
y.mesh loc=0.10    spac=0.05   
y.mesh loc=0.13    spac=0.02   
y.mesh loc=0.23    spac=0.05
y.mesh loc=9.13    spac=0.1    
y.mesh loc=659.13  spac=0.1     
y.mesh loc=659.24  spac=0.05  

# ... 其他区域定义保持不变 ...

material material=Ga2O3 user.group=semiconductor \
         eg300=4.8 permittivity=10 mun=1000 mup=100 nc300=3.72e18 \
         nv300=1.16e19 taun=1e-8 taup=1e-8 affinity=4.0 ni=1e-5
material material=Au user.default=gold

# ... 其他物理模型和接触设置保持不变 ...

contact name=anode workfunction=4.0
contact name=cathode type=ohmic

method newton maxtrap=100 climit=1e-6 

solve init
save outfile=Ga2O3_SBD_ordered.str
log outfile=IV.log

solve vanode=0.0 vstep=-0.1 vfinal=-5 name=anode   
solve vanode=-5 vstep=-0.1 vfinal=-20 name=anode

log off
tonyplot IV.log

quit

---

四、总结建议

• 重点调整：材料迁移率、掺杂浓度、接触功函数、网格精度和求解器设置。
• 如果以上调整后仍无法达到目标电流密度，建议进一步检查器件结构设计是否合理，例如电极尺寸、界面质量、热效应等。

如需进一步调试，请提供完整的仿真日志文件（IV.log），以便进行更深入分析。