结晶过程理解：从实验室到工业放大

结晶是原料药生产的最后一步，也是决定产品质量的关键一步：晶型、粒径分布、晶习（形貌）、纯度、流动性——这些都在结晶过程中确定，且下游几乎无法修正。

结晶的两个基本过程

	成核（nucleation）	生长（growth）
过程	形成新的晶核	晶核长大
驱动力	过饱和度	过饱和度
对过饱和度的敏感性	极强（指数关系）	较弱（近似线性）
结果	决定晶体数量	决定晶体大小

控制粒径的核心逻辑：成核与生长的相对速率决定了最终粒径。过饱和度高 → 大量成核 → 得到很多小晶体；过饱和度低 → 少量成核 + 充分生长 → 得到少量大晶体。因此控制粒径的关键是控制过饱和度的产生速率。

介稳区：结晶控制的核心概念

# 溶解度曲线与超溶解度曲线之间的区域 = 介稳区（MSZW）
#
#   浓度
#    ↑
#    │        不稳区（自发成核）
#    │  ────────────────────  ← 超溶解度曲线
#    │        介稳区
#    │        （不自发成核，但晶种能生长）
#    │  ────────────────────  ← 溶解度曲线
#    │        欠饱和区（晶体溶解）
#    └──────────────────────→ 温度
#
# 工艺控制策略：
#   1) 冷却到介稳区（不自发成核）
#   2) 加入晶种（提供生长表面）
#   3) 控制冷却速率，让体系保持在介稳区内
#      → 全部过饱和度用于生长，而非二次成核
#   4) 得到粒径均匀的大晶体
#
# 如果冷却太快冲进不稳区：
#   → 大量自发成核 → 细小晶体 → 过滤困难、流动性差

关键工艺参数

参数	影响
过饱和度控制方式	冷却 / 反溶剂 / 蒸发 / 反应结晶
晶种	加入时机、加入量、晶种粒径与晶型
冷却/加料曲线	线性 vs 立方曲线；控制过饱和度恒定
搅拌	影响传质、二次成核、晶体破碎
溶剂选择	影响溶解度曲线、晶习、晶型
杂质	可能抑制生长或改变晶习
老化/熟化时间	Ostwald 熟化，改善粒径分布

晶种策略是工业结晶最重要的控制手段：加入正确晶型、合适粒径的晶种，在合适的过饱和度下加入，能同时控制晶型、粒径和批间一致性。不加晶种的自发结晶几乎无法保证重现性。

立方冷却曲线

# 为什么不用线性冷却？
#
# 结晶初期晶体表面积小，生长慢，
# 如果按线性冷却，过饱和度会快速累积 → 二次成核
#
# 立方冷却曲线：初期慢、后期快
#   T(t) = T_start - (T_start - T_end) × (t/t_total)³
#
# 目的：让过饱和度在整个过程中保持大致恒定
#      → 全部用于生长，避免二次成核

import numpy as np

def cubic_cooling_profile(T_start, T_end, total_time, n_points=100):
    t = np.linspace(0, total_time, n_points)
    T = T_start - (T_start - T_end) * (t / total_time) ** 3
    return t, T

# 反溶剂结晶同理：加料速率应先慢后快

# 更先进：基于在线浓度测量的反馈控制
#   实时测浓度 → 计算当前过饱和度 → 调整冷却/加料速率
#   保持过饱和度在设定值

过程分析技术（PAT）

技术	测什么
FBRM（聚焦光束反射测量）	颗粒计数与弦长分布，实时监测成核与生长
ATR-FTIR / 拉曼	溶液浓度 → 计算过饱和度；也可监测晶型
PVM（颗粒视觉测量）	原位显微图像，看晶习与聚集
在线 PXRD	晶型监测
浊度	简单的成核检测

PAT 让结晶从「黑箱操作」变成「可观测可控制的过程」：知道什么时候成核、过饱和度多少、粒径如何演变，才能真正理解并优化工艺。

放大时的挑战

• 混合不均匀：这是放大最主要的问题。大反应釜中，反溶剂加入点附近的局部过饱和度远高于平均值，导致局部大量成核。实验室小烧瓶中混合几乎瞬时，这个问题不存在。
• 传热速率下降：散热面积/体积比下降，冷却速率受限，实际温度曲线可能偏离设定。
• 搅拌剪切：大规模搅拌的剪切分布不均，可能造成晶体破碎或磨损，产生细粉。
• 停留时间分布：连续结晶中的返混问题。
• 放大准则的选择：等搅拌功率/体积？等叶端速度？等混合时间？不同准则会给出不同的转速，需要根据主导因素选择——对结晶通常关注混合时间与局部过饱和度。

常见问题与对策

问题	可能原因	对策
晶体太细，过滤慢	过饱和度过高、二次成核	加晶种、放慢冷却、优化搅拌
晶习不好（针状）	溶剂、杂质影响	换溶剂、加晶习调节剂
批间粒径不一致	成核不受控	晶种策略
得到错误晶型	动力学控制、缺少晶种	加正确晶型晶种、延长熟化
包裹母液（纯度低）	生长太快	降低过饱和度、延长生长时间
结块 / 聚集	过饱和度高、搅拌不足	优化搅拌与过饱和度

建模与优化

# 群体平衡模型（PBM）是结晶建模的标准框架
#
# 描述晶体粒径分布随时间的演变：
#   ∂n/∂t + ∂(G·n)/∂L = B - D
#
#   n = 粒径分布函数
#   G = 生长速率（过饱和度的函数）
#   B = 成核速率（过饱和度的强函数）
#   D = 消失项（溶解、破碎）
#
# 结合物料衡算与能量衡算，可以模拟整个结晶过程
#
# 用途：
#   - 预测不同工艺参数下的粒径分布
#   - 优化冷却/加料曲线
#   - 指导放大
#
# 参数（成核与生长动力学常数）需从实验数据拟合
# → PAT 数据是拟合这些参数的关键来源

延伸资源

• 多晶型：392《晶型预测 CSP》、393《多晶型 Polymorphism》；盐型与共晶：394《盐型选择 Salt Selection》、395《共晶 Cocrystal》；
• 工艺优化：397《工艺路线优化》；溶解度：396《溶解度与溶出预测》。