外延层厚度与自掺杂过渡区分析器｜工艺参数分析工具
这个工具解决什么问题
同一道工序不同机台的数据要不要合并算，是个口径问题。这工具把选择权交给你。
外延厚度按「生长速率 x 时间」算，但决定击穿特性的是扣除自掺杂过渡区之后的净厚度。 本工具回归生长速率的活化能、定时长、算片内均匀性、用 SIMS 剖面拟合过渡区宽度、 算净外延厚度，并判断击穿电压是雪崩限制还是穿通限制。
输入数据格式
数据文件 epi_growth_rate.csv，列顺序为：温度_摄氏、生长速率_微米每分钟、1overT_每K、lnGR。
示例首行：1000.0,0.2858,0.00078545,-1.25233
示例输出（真实运行结果节选）
==============================================================================
外延层厚度与自掺杂过渡区分析报告
==============================================================================
一、生长速率-温度回归（Arrhenius：ln GR = ln GR_ref - Ea/kB x (1/T - 1/Tref)）
温度(度C) 1/T(1/K) 实测速率(微米/分) 模型速率 相对偏差
----------------------------------------------------------------------
1000 0.00078545 0.2858 0.2874 +0.54%
1050 0.00075577 0.4928 0.4916 -0.25%
1100 0.00072825 0.8200 0.8087 -1.38%
1150 0.00070267 1.2683 1.2846 +1.29%
1180 0.00068816 1.6734 1.6701 -0.20%
拟合活化能 Ea = 1.559 eV（拟合优度 R^2 = 0.99981）
Ea 的物理意义：它反映生长是「表面反应控制」还是「质量输运控制」。
Ea > 1.0 eV 说明低温段由表面反应控制（对温度极敏感，1% 温度偏差就有可观厚度偏差）；
高温段 Ea 会明显下降并转入输运控制 —— 所以外推温度不要跨越这个转折点。
二、工艺温度 1120 度C 下的生长方案
模型生长速率 0.9770 微米/分钟
目标厚度 3.00 微米
所需生长时间 3.07 分钟（184 秒）
温度灵敏度 0.94 %/度C（1 度温度偏差带来的厚度偏差）
三、片内 5 点厚度均匀性
中心 3.021 / 上 3.048 / 下 2.987 / 左 3.010 / 右 3.032 微米
平均厚度 3.0196 微米；极差 0.0610 微米（
一句话
厚度合格、击穿电压不合格 —— 多半是自掺杂过渡区吃掉了有效外延厚度。
运行环境与快速开始
python epi_autodoping_transition_analyzer.py
输入数据格式（换成你的数据即可）
· epi_growth_rate.csv：温度 / 生长速率 / 1overT / lnGR（用于活化能回归） · epi_sims_profile.csv：深度 / 载流子浓度（用于过渡区拟合） 脚本内 FIVE_POINT 为片内 5 点厚度，顶部常量可改工艺温度、目标厚度与外延掺杂浓度。
输出说明
生长速率-温度回归与活化能、工艺温度下的生长方案与温度灵敏度、片内均匀性、 过渡区拟合与净外延厚度、击穿电压判定（雪崩/穿通）与目标电压所需厚度、结论与建议。
算法口径（关键公式）
· Arrhenius：GR(T) = GR_ref * exp( -Ea/kB * (1/T - 1/Tref) )，Ea 由 lnGR 对 (1/T - 1/Tref) 回归求出 · 自掺杂剖面：C(x) = Cb + (Cs - Cb) * erfc(x/W)，在 (Cs, W) 上做相对残差最小二乘 · 过渡区宽度：浓度降到本底 10 倍处的深度 · 击穿：耗尽宽度 W = eps*Ec/(q*Nd)；净厚度不足则按穿通公式 BV = Ec*t - q*Nd*t^2/(2*eps)
判定标准说明
自检包含活化能回归接近真值（示例 1.56 eV）、过渡区 W 接近真值、片内极差小于 5%、 示例落在穿通限制区且击穿电压低于雪崩极限。
常见问题
· 为什么强调「净厚度」？ 过渡区浓度渐变，既不完整耐压又抬高有效掺杂。 · 活化能有什么用？ 大于 1 eV 说明生长由表面反应控制，对温度极敏感（示例每度约 0.9%）。 · 怎么缩短过渡区？ 降低埋层表面浓度或优化氢预烘焙，均需 SIMS 复测确认。
可调参数与计算口径
脚本顶部集中了主要阈值与常量：REPORT=os.path.join(BASE, "示例输出报告.txt")、FIVE_POINT=[3.021, 3.048, 2.987, 3.010, 3.032]、GEN_SEED=20261006。改完直接重跑即可，不必改动计算逻辑。
会生成哪些文件
运行后除控制台输出外，还会生成：epi_growth_rate.csv、epi_sims_profile.csv、示例输出报告.txt。
脚本里的主要处理步骤
gr_at()：Arrhenius 生长速率模型（微米/分钟）。fit_arrhenius()：按 ln GR 对 (1/T - 1/Tref) 做最小二乘，返回 (Ea, ln GR_ref, R2)。fit_transition()：在 (Cs, W) 上做相对残差最小二乘，返回 (Cs, W, 相对R2)。
以上步骤在脚本中是分开的函数，可以单独调用或按需替换其中一段，不必整体重写。
工具中出现的字段与提示
- Arrhenius 生长速率模型（微米/分钟）。
- 温度_摄氏
- 生长速率_微米每分钟
- 1overT_每K
- 深度_微米
- 载流子浓度_每立方厘米
- 在 (Cs, W) 上做相对残差最小二乘，返回 (Cs, W, 相对R2)。
- 外延层厚度与自掺杂过渡区分析报告
- 拟合活化能 Ea = %.3f eV（拟合优度 R^2 = %.5f）
- Ea 的物理意义：它反映生长是「表面反应控制」还是「质量输运控制」。
- 高温段 Ea 会明显下降并转入输运控制 —— 所以外推温度不要跨越这个转折点。
- 二、工艺温度 %.0f 度C 下的生长方案
- 模型生长速率 %.4f 微米/分钟
- 目标厚度 %.2f 微米
- 所需生长时间 %.2f 分钟（%.0f 秒）
- 温度灵敏度 %.2f %%/度C（1 度温度偏差带来的厚度偏差）
- 三、片内 5 点厚度均匀性
- 平均厚度 %.4f 微米；极差 %.4f 微米（极差/均值 = %.2f%%）
其余字段以界面与示例报告为准。
下载与文件清单
压缩包内包含：主程序脚本、示例数据文件、示例输出报告、中文说明文档，以及同主题 AI 提示词一套。解压后按说明运行即可，示例数据可直接试跑。
主程序文件：epi_autodoping_transition_analyzer.py。
本站工具免费下载，无需注册。工具为 Python 脚本，Python 3.8 及以上环境可直接运行，核心功能尽量不依赖第三方库。
常见问题
输入数据要准备成什么样？
按 示例首行中的列顺序准备即可，共 4 列：温度_摄氏、生长速率_微米每分钟、1overT_每K、lnGR。列为空或格式不符时结果会失真，建议先用示例数据试跑一遍。
阈值可以改成我们自己的标准吗？
可以，脚本顶部已把主要阈值与常量集中声明，按自己产线的标准修改后重跑即可，不需要改计算逻辑。
附赠的提示词怎么用？
把提示词与自己的数据一起交给大模型，可以让它按同一套口径帮你改写分析流程或扩展报告字段。
示例数据能直接跑吗？
可以。示例数据与主程序在同一压缩包内，解压后直接运行即可看到完整输出，确认无误再替换成自己的台账。
同分类的其它工具
所属分类 半导体工艺参数分析工具包，也可回到 工具资源包下载页 按分类查找。
