### Python 实现计算溶液 pH 值的算法 以下是一个基于 Python 的解决方案，该方案可以根据输入的离子和盐类信息动态计算溶液的 pH 值。程序的核心逻辑在于解析输入数据、构建化学平衡模型并利用数值方法求解氢离子浓度。 --- #### 数据结构设计 为了便于管理复杂的离子体系及其对应的电离常数，采用字典存储相关信息： - `ions` 字典保存每种离子的最大带电量及各级电离平衡常数。 - `salts` 列表记录每种盐的具体组成及其浓度。 --- #### 数学建模 根据质子条件（Proton Condition），建立关于氢离子浓度 `[H+]` 的代数方程组。对于复杂体系中的多种酸碱成分，需综合考虑所有可能的贡献项，并通过迭代法求数值解。 --- #### 程序实现 ```python import math from scipy.optimize import fsolve def calculate_pH(ions, salts): """ 计算溶液的pH值 :param ions: dict, 存储离子的相关属性 {name: (max_charge, pKa_list)} :param salts: list of tuples, [(salt_name, concentration)] :return: float, 溶液的pH值 """ # 初始化变量 ion_concentrations = {} # 构造初始离子浓度分布 for salt in salts: name, conc = salt components = parse_salt(name) # 将盐名分解为构成它的离子 for component in components: if component not in ion_concentrations: ion_concentrations[component] = 0.0 ion_concentrations[component] += conc * abs(components[component]) # 根据摩尔比例分配浓度 # 定义目标函数：质子条件方程 def proton_condition(h_plus): total_h_plus = h_plus oh_minus = 1e-14 / h_plus # Kw = [H+][OH-] for ion, properties in ions.items(): max_charge, pkas = properties current_ion_conc = ion_concentrations.get(ion, 0.0) if max_charge > 0: # 酸性离子 for i in range(max_charge): ka_i = 10 ** (-pkas[i]) numerator = ka_i * h_plus denominator = sum((ka_j * h_plus)**j for j in range(i + 2)) total_h_plus -= current_ion_conc * numerator / denominator elif max_charge < 0: # 碱性离子 for i in range(abs(max_charge)): kb_i = 10 ** (-14 + pkas[i]) # Kb from Ka and Kw relation numerator = kb_i * oh_minus denominator = sum((kb_j * oh_minus)**j for j in range(i + 2)) total_h_plus += current_ion_conc * numerator / denominator return total_h_plus - h_plus # 使用fsolve求解[H+] initial_guess = 1e-7 # 默认假设中性环境下的[H+]=1e-7 M solution = fsolve(proton_condition, initial_guess) # 返回pH值 ph_value = -math.log10(solution[0]) return round(ph_value, 2) def parse_salt(salt_name): """ 解析盐的名字，返回其中包含的离子及其系数 :param salt_name: str, 盐的名称形式如 "(NH4)(CO3)" :return: dict, 含有的离子及其系数 {ion_name: coefficient} """ result = {} parts = salt_name.strip('()').split(')') for part in filter(None, parts): # 过滤掉空字符串 sub_parts = part.split('(') ion = sub_parts[-1].strip() coeff_str = ''.join(sub_parts[:-1]).replace('(', '') try: coeff = int(coeff_str or '1') # 如果没有显式的系数，默认为1 except ValueError: raise ValueError(f"无法解析盐 '{salt_name}' 中的系数部分") if ion in result: result[ion] += coeff else: result[ion] = coeff return result # 测试案例 if __name__ == "__main__": # 输入定义 ions_data = { "NH4": (1, [9.25]), # NH4+: 单级电离，pKa=9.25 "CO3": (-2, [6.37, 10.25]) # CO3^2-: 双级电离，pKa=[6.37, 10.25] } salts_data = [ ("(NH4)(CO3)", 0.1), ("(NH4)4(CO3)(PO4)(H)", 0.2), # 更复杂的例子可扩展支持更多离子 ("(NH4)2(PO4)(H)", 0.1) ] # 执行计算 final_ph = calculate_pH(ions_data, salts_data) print(final_ph) ``` --- #### 关键点解释 1. **质子条件**：核心思想是确保溶液中所有的质子来源与消耗达到平衡状态。 2. **数值求解**：由于实际问题往往涉及非线性方程，因此选用 SciPy 提供的 `fsolve` 函数寻找满足条件的 `[H+]` 值。 3. **灵活性**：代码允许用户自定义任意数量的离子和盐类组合，适应性强。 ---