金属配位平衡形态分析软件

软件由AI辅助编程完成，可能存在错误，可以联系本人探讨。点击最后阅读原文访问软件界面。推荐使用chrome浏览器桌面功能显示或类似能模仿台式机界面功能

金属配位平衡形态分析仪

逐级稳定常数 · 累积常数 · 分布分数 · 平均配位数 · 物料守恒数值求解

中 文 使 用 说 明 书

配套文件：金属配位平衡形态分析.html

一、软件概述

本软件是一个单文件网页工具，用于分析金属离子 M 与单齿配体 L 的逐级配位平衡M + nL ⇌ MLₙ。它把配位化学中几个容易算错、画错的量全部交给计算机完成，并把抽象的平衡关系变成可以直接读数的图表。

软件做四件事：

·把逐级稳定常数 Kₙ 自动换算为累积稳定常数 βₙ（βₙ = K₁K₂…Kₙ）。

·绘制各组分的分布分数 αᵢ 随 pL 变化的曲线（α–pL 图）。

·绘制平均配位数 n̄ 随 pL 变化的形成曲线（n̄–pL 图）。

·在给定的金属总浓度 C_M 与配体总浓度 C_L 下，数值求解真正的游离配体浓度 [L]，并给出完整的平衡组成表。

1.1 与其他工具的区别

多数教材与课件只画「示意图」或只给出某几个 pL 点的 α 值，而本软件的核心是第 4 项——它不假定 [L] 已知，而是通过物料守恒方程C_L = [L] + C_M·n̄ 反解出真实的 [L]。这一点很关键：配位平衡中加入的配体只有一部分是游离的，其余被金属结合，直接用 C_L 当 [L] 会算错。

1.2 四个核心概念

二、运行环境与启动

三、界面总览

此外，左栏底部有「公式说明」折叠区，点击可展开全部核心公式。

四、输入详解

4.1 五个内置体系预设

点击预设按钮会一次性载入该体系的名称、配位数、逐级常数与默认浓度。蓝色高亮的为当前选中的体系。

4.2 最高配位数 N

下拉选择 1 ~ 8 级，体系生成 M、ML₁ … ML_N 共 N+1 个物种。切换 N 时，原有已填的 log K 会被保留，新增的级数按 3 − 0.5i 的默认值补齐（这是一个合理的逐级递减初值，需按实际体系修改）。

4.3 浓度输入方式

C_M 与 C_L 均采用「尾数 × 10^指数」两段式输入：

·尾数框：0.01 ~ 9.99，步长 0.1

·指数下拉：10⁰ ~ 10⁻⁸

例如输入 1.0 × 10⁻³ 即 1.0×10⁻³ mol/L。两个浓度任一项修改后，所有图表与结果立即重算。

4.4 逐级常数输入

第 n 行对应反应 MLₙ₋₁ + L ⇌ MLₙ，填的是 log Kₙ（以 10 为底）。

·可以为负值（表示该级配合物不稳定，如 Cd–Cl⁻ 的 K₃、K₄）。

·留空或非法输入时，软件显示黄条提示「存在未填写或非法的 log K 值」，补全后自动恢复计算。

·表格每行左侧的色块与分布分数图中该组分的曲线颜色一致，便于对照。

五、输出详解

5.1 常数换算表

按 βₙ = K₁K₂…Kₙ 逐行累积相乘，同时给出 Kₙ、log βₙ 与 βₙ。注意 β₀ = 1（对应游离金属 M）。

5.2 平衡组成区（核心输出）

5.3 分布分数图（α–pL）

5.4 平均配位数形成曲线（n̄–pL）

纵轴为 n̄，范围 0 ~ N。曲线呈 S 形（单调递增）：pL 越小（配体越多）n̄ 越接近 N；pL 越大（配体越少）n̄ 越接近 0。图上以圆点标出当前工作点 n̄ 的数值。

5.5 辅助功能

·恢复当前体系默认值：把参数还原为所选预设的初始值（快捷键式按钮）。

·复制结果摘要：把体系、常数、浓度、pL、[L]、n̄ 与全部物种的 α 和浓度复制为纯文本，可直接粘贴进实验报告或课件。

六、模型与算法原理

6.1 基本定义

逐级稳定常数：Kₙ = [MLₙ] / ([MLₙ₋₁]·[L])

累积稳定常数：βₙ = K₁·K₂·…·Kₙ = [MLₙ] / ([M]·[L]ⁿ)，β₀ = 1

分布分数：αᵢ = βᵢ[L]ⁱ / Σⱼ βⱼ[L]ʲ，满足 Σαᵢ = 1

平均配位数：n̄ = Σ i·αᵢ

pL = −log[L]

6.2 数值稳定化处理

直接计算 βᵢ[L]ⁱ 在 i 较大、[L] 很小时容易下溢。软件采用「先取对数、再指数归一化」的写法：

tᵢ = log βᵢ + i·log[L]

m= max(tᵢ)

αᵢ = exp(tᵢ − m) / Σⱼ exp(tⱼ − m)

由于最大项对应 exp(0) = 1，分母恒不小于 1，从根本上避免了溢出与下溢。

6.3 求解游离配体浓度 [L]

这是本软件的关键步骤。配体总量一部分游离、一部分被金属结合，满足：

C_L = [L] + C_M · n̄([L])

移项构造目标函数：

g(log[L]) = C_L − 10^(log[L]) − C_M · n̄(log[L])

g 是关于 log[L] 的单调递减函数，因此用二分法求解 g = 0：

·搜索区间 log[L] ∈ [−25, 12]，先以 0.25 为步长粗扫定位变号区间；

·再二分迭代至多 300 次，收敛判据为 |g| < 10⁻¹⁵(C_L + C_M) 或区间宽度 < 10⁻¹³；

·解得 log[L] 后回代，得到全部 αᵢ、n̄ 与各物种平衡浓度。

若在整个区间内 g 不变号（例如 C_L 极小或常数异常），软件返回无解并显示黄条提示，属正常结果而非故障。

6.4 pL 显示范围自动适配

autoRange 函数在 pL ∈ [−2, 22] 上以 960 步扫描，找出「至少有 2 个组分 α ≥ 5%」的 pL 区间，再向两侧各扩展 1.5 个单位并取整；若区间不足 6 个单位则补足到 6。这样能自动把图形聚焦在真正发生形态转化的区域，避免曲线挤在一角。

6.5 一个可直接验证的规律

相邻两个组分的分布分数相等时，有简洁的解析结果：

αᵢ = αᵢ₊₁⟺βᵢ[L]ⁱ = βᵢ₊₁[L]ⁱ⁺¹⟺[L] = βᵢ/βᵢ₊₁ = 1/Kᵢ₊₁

即：pL = log Kᵢ₊₁

换言之，第 i 条与第 i+1 条曲线的交点，横坐标恰好等于 log Kᵢ₊₁。这是一个可以直接在图上读数验证的规律，也是检验常数输入是否正确的好办法。说明书示例 4 给出了五个体系的逐项验证，差值均为 0.0000%。

七、内置常数表（含累积常数）

注意 Ni²⁺–NH₃ 的 log β₆ = 8.31 小于 log β₅ = 8.34，Cd²⁺–Cl⁻ 的 log β₃、log β₄也依次下降——这是因为对应的逐级 Kₙ < 1。累积常数并非必须单调递增，详见示例 5。

八、计算示例

以下所有数值均由本软件的算法复算得出，可直接用于核对软件读数。

示例 1　Cu²⁺–NH₃ 默认条件的完整求解

条件：C_M = 1.0×10⁻³ mol/L，C_L = 2.0×10⁻³ mol/L，log K = [4.04, 3.43, 2.80, 1.48]。

物料守恒：C_L = [L] + C_M·n̄

解得：pL = 3.393，[L] = 4.048×10⁻⁴ mol/L，n̄ = 1.5952

残差 = 9.5×10⁻¹⁷ mol/L（收敛良好）

读法：加入的配体总量 C_L = 2.0×10⁻³，但游离的只有 4.048×10⁻⁴，约 80% 被金属结合。体系以 ML₁ 与 ML₂ 为主（合计 63.3%），平均每个铜结合 1.60 个 NH₃——注意 n̄ 不是整数，也不等于优势组分的级数。

示例 2　五个预设体系的求解结果对比

对比可见：Cd²⁺–CN⁻ 与 Zn²⁺–en 的常数很大（log β₄ = 17.92、log β₃ = 13.89），配体过量不多时 n̄ 就已达 2.5 ~ 2.8，几乎全部金属以高配位型体存在；而 Ni²⁺–NH₃ 常数小（log K₁ 仅 2.72），即使 C_L/C_M = 3，n̄ 也只有 1.05，游离金属仍占 35%。

示例 3　配体浓度对形态的影响（Cu²⁺–NH₃，C_M = 1×10⁻³）

规律：随配体增加，pL 下降、n̄ 单调上升，优势物种依次从 M 迁移到 ML₂、ML₃。但 n̄ 上升越来越慢——继续加大配体，n̄ 趋近于上限 N = 4 而永远不会达到（因为总有少量低配位型体存在）。这正是形成曲线呈 S 形的原因。

示例 4　验证交点规律：αᵢ = αᵢ₊₁ 处 pL = log Kᵢ₊₁

五个体系、全部相邻交点的预测值与软件实测值偏差均为 0.0000%，可作为软件正确性的自检手段，也适合作为课堂练习：让学生先在图上读出交点 pL，再与 log K 对照。

示例 5　逐级常数小于 1 的反常情形（Cd²⁺–Cl⁻）

该体系 log K₃ = −0.20、log K₄ = −0.70，即 K₃ = 0.631、K₄ = 0.1995，均小于 1。后果是累积常数不再单调：

log β：0.00 → 1.98 → 3.60 → 3.40 → 2.70

β₂ = 3.981×10³>β₃ = 2.512×10³>β₄ = 5.012×10²

由交点规律可知：ML₃ 要超过 ML₂ 需 [L] > β₂/β₃ = 1.585 mol/L（pL < −0.20），ML₄ 要超过 ML₃ 需 [L] > β₃/β₄ = 5.012 mol/L（pL < −0.70）。也就是说，ML₃、ML₄ 只在配体浓度大于 1 mol/L 的极端条件下才占优势——这类高配位型体在水溶液中实际上很难成为主要形态。本软件把 pL 下限定为 −2，正是为了把这段区域纳入显示范围。

示例 6　边界情形

注意最后一列：即使配体浓度加到 1 mol/L（是金属的 1000 倍），n̄ 也只有 3.793 而非 4。这是配位平衡的本质——任何有限浓度下都不可能达到理论上限。

九、教学应用建议

·概念引入：先用 Cu²⁺–NH₃，让学生观察 α 图上 5 条曲线的此消彼长，建立「形态分布是连续的、随 [L] 变化」而非「非此即彼」的认识。

·逐级递减的直观化：指出 Cu–NH₃ 的 log K 从 4.04 逐级降到 1.48，对应图上各曲线峰值依次降低、宽度变窄。

·交点规律练习（最推荐）：让学生用鼠标悬停读出相邻曲线交点的 pL，再与常数表中的 log K 对照，验证 pL = log Kₙ。示例 4 给出标准答案。

·n̄ 的物理意义：强调 n̄ 是加权平均、可以不是整数；演示示例 1 中优势物种是 ML₂ 而 n̄ = 1.595 的来源。

·物料守恒观念：对比 C_L 与解得的 [L]，说明「加入量 ≠ 平衡浓度」，这是学生最容易出错的地方。

·反常体系讨论：用 Cd²⁺–Cl⁻ 讨论 Kₙ < 1 的情形，说明累积常数不一定单调、高配位型体未必是主要形态。

·参数探究：固定 C_M，逐步加大 C_L（示例 3），让学生自己总结 n̄ 的 S 形变化与饱和趋势。

十、注意事项与已知限制

10.1 模型口径

·只考虑单齿配体 L 与 MLₙ 型单核配合物，不考虑螯合环外的其他配体、羟合/水解型体、多核配合物与沉淀。

·输入的是浓度常数而非热力学常数，未计入活度系数与离子强度效应。内置常数对应 25 °C、离子强度 ≈ 0 的条件。

·不考虑配体的质子化（如 NH₃ + H⁺ ⇌ NH₄⁺）与 pH 影响。实际氨体系必须在缓冲条件下讨论，游离 [NH₃] 随 pH 变化很大；本软件的 [L] 应理解为「有效游离配体浓度」。

·未考虑竞争配位与掩蔽剂；C_L 指可与 M 配位的该配体总量。

·逐级常数相互独立输入，软件不校验其热力学自洽性（理论上应满足统计效应导致的逐级递减趋势，但允许输入任意值，包括负值）。

10.2 数值与显示

·求解区间为 log[L] ∈ [−25, 12]。若真实解落在该区间之外（C_L 极小或极大），软件会显示无解提示。

·分布分数图纵轴固定为 0 ~ 1，n̄ 图纵轴为 0 ~ N；两者共用同一 pL 横轴范围。

·「自动适配 pL 范围」默认勾选。若手动缩放后想固定视角，可取消勾选；取消后切换体系时范围不会自动更新，需重新勾选。

·浓度输入的两段式（尾数 × 10^指数）中，指数下拉仅提供 10⁰ ~ 10⁻⁸，不支持正指数（浓度大于 1 mol/L 的场景）。

·数值显示规则：绝对值 ≥ 10⁵ 或 < 10⁻⁴ 时用科学记数法（保留 3 位有效数字），否则保留 4 位有效数字。

10.3 与教材数值的出入

不同教材、数据库的稳定常数因测定方法、离子强度、温度不同常有出入，例如 Cu²⁺–NH₃ 的 log β₄ 文献值在 11.8 ~ 13.0 之间均有报道。本软件内置常数取自 Martell & Smith《Critical Stability Constants》（页脚已注明），内部自洽，适用于教学演示与趋势分析。若需与特定教材或实验数据严格对齐，请直接在第 ② 区手动修改 log Kₙ。

本说明书依据软件实际算法与界面撰写；所有示例数值均由该算法复算，可用于核对软件读数。