基于ATR-UV 實時濃度檢測的撲熱息痛冷卻結晶過程研究

關國強，唐凱，李紅，盧帥濤，江燕斌

（1 華南理工大學化學與化工學院，廣東 廣州 510640；2 傳熱強化與過程節能教育部重點實驗室，廣東 廣州 510640； 3 廣東肇慶星湖生物科技股份有限公司技術中心，廣東 肇慶 526040）

引 言

結晶是化學工程中重要的單元操作，廣泛應用于生物、制藥和精細化工等領域[1]。隨著高附加值功能性精細化學品的需求增加，傳統基于“黑箱模型”的試差結晶過程設計已無法滿足現代產品工程的需要。因此，基于過程分析技術（process analytical technology，PAT）的結晶過程控制與優化技術已成為結晶過程研究和應用的熱點[2-3]。

結晶是在組分過飽和度驅動下，克服分子有規化“阻力”的液固分離過程，故過飽和度對結晶成核、晶粒生長和團聚有重要的影響。由于過飽和度是組分濃度與相應溫度下平衡濃度的函數，為獲得形貌優異的晶體產品，需對結晶過程的組分濃度和溫度進行實時監控，結合粒度在線檢測，進而有效調節結晶過程操作參數。相比傳統的“先分離、后檢測”的離線濃度監測技術，如高效液相色譜（high performance liquid chromatography，HPLC）等[4]，結合衰減全反射技術（attenuated total reflection，ATR）信號采集的濃度檢測直接光譜法能實時測定結晶體系局部濃度分布，已成為PAT 濃度實時檢測的重要手段。ATR 光譜分析按其光源主要分為紅外（infrared spectrum, IR）和紫外（ultraviolet spectrum, UV）分析。相對于采用ATR-IR 的大量報道[5-16]，ATR-UV 因偵測距離與波長呈反比，故能分析更薄層表面信息，已成為新近結晶過程研究的有效方法[17]。聚焦光束反射測量技術（focused beam reflectance measurement，FBRM）是結晶過程廣泛使用的在線顆粒粒度檢測方法[18-19]。FBRM 通過測定弦長分布關聯結晶過程的晶粒尺度分布（crystal size distribution，CSD）變化、晶體成核、生長及粒間團聚等現象。采用FBRM已成為當前結晶產品控制的重要手段，其最新應用研究包括作為直接成核控制的控制傳感器[20]。

為構建“質量源于設計”（quality by design，QbD）的高效結晶分離系統，本研究在先前UV-vis法[21]研究基礎上采用更先進的ATR-UV 方法，通過近年廣泛應用于光譜分析的偏最小二乘法（partial least squares，PLS）[22]實現典型活性藥物組分——撲熱息痛（paracetamol，PCM）乙醇溶液結晶過程的實時濃度分析。進而結合FBRM 實時分析，研究不同降溫方式對結晶過程的影響。

1 實驗材料和方法

1.1 材料

PCM（CAS#103-90-2，純度≥99%）購自上海思域化工科技有限公司；分析純無水乙醇購自國藥集團化學試劑有限公司。

1.2 結晶實驗系統及主要儀器設備

PCM-乙醇體系冷卻結晶實驗系統置如圖1 所示，它由500 ml 帶攪拌夾套結晶釜、程控恒溫冷卻系統、ATR-UV 檢測系統、FBRM 分析系統及計算機數據采集系統組成。主要分析儀器見表1，分別采用ATR-UV檢測系統和FBRM實時分析結晶過程濃度和晶體粒度分布（crystal size distribution，CSD），用偏光顯微鏡觀測晶體產品的形貌，其他儀器包括樣品處理所需的真空干燥箱和分析天平等。

表1 主要儀器設備一覽表Table 1 List of main experimental instruments

1.3 ATR-UV 濃度關聯及校驗

配制不同濃度的PCM 乙醇溶液，采用ATR-UV光譜儀測定其不同溫度下光譜響應（以純乙醇為空白），見表2。

表2 ATR-UV 濃度關聯實驗條件Table 2 Experimental parameters for ATR-UV correlation

ATR-UV 光譜數據采用PLS 法進行分析，式（1）用于關聯不同溫度（T）下PCM 的濃度

式中，Ai為選定的特征波長吸光度值；模型參數a、b 和ci通過PLS 擬合獲得。

采用Umetrics 公司的SIMCA-P 11.5 分析軟件導入對已知濃度的PCM 溶液在不同溫度下測定的特征波長吸光度值，采用壓縮-中心化方法進行數據標準化后，根據霍特林T2（Hotellin T2，置信水平為0.05）分布檢驗數據點并建立回歸建模用的數據集。以最小差均方根平均偏差（root-mean-square error，RMSE）為目標函數，如式（2）所示，通過PLS 回歸式中的模型參數，使RMSE 最小。

式中，yi和?iy 分別為訓練集和校正集的濃度；n為數據點數目。

對建立的ATR-UV 在線PCM 濃度檢測模型，采用重量法進行校驗。重量法濃度測定過程如下[23-24]：

（1）取樣：采用預熱的帶砂芯過濾頭1 ml 注射器實時取樣，其中平均孔徑為1 μm；

（2）真空干燥：將溶液迅速移至稱量瓶并置于真空干燥箱中（40℃，40 kPa）干燥24 h 至恒重；

（3）稱量：在電子分析天平上測定樣品殘留固體質量。

PCM 在乙醇中的濃度（y），用式（3）計算

式中，WR為樣品殘留固體質量；WS為取樣 質量。

1.4 飽和溶解度關聯及過飽和度測定

采用文獻[25]報道的PCM 溶解度數據（y?），結合實時濃度檢測結果，即可計算結晶過程的關鍵過飽和度（S），如式（4）所示

1.5 PCM 冷卻結晶實驗

應用上述ATR-UV 濃度關聯可實時分析PCM結晶中的過飽和度變化，通過同步調節過程溫度，實現受控結晶。為此在如圖1 所示結晶系統中進行PCM 結晶實驗。

實驗具體步驟如下：

（1）母液配制：稱量85.6 g PCM 置于結晶釜中，加入300 ml 乙醇溶劑，維持攪拌200 r·min-1，60℃條件下恒溫1 h 以確保固體顆粒完全溶解；

（2）加入晶種：待母液冷卻至T0（=48℃）時加入4.7 g 13～21 μm 的晶種并穩定10 min；

（3）冷卻結晶：通過程序控制恒溫槽調控結晶溫度變化，使結晶釜內PCM 溶液冷卻結晶，同時采用ATR-UV 關聯溶液的過飽和度和通過FBRM分析結晶產物CSD，冷卻結晶持續時間（t1）為3 h。

分別采用3 種冷卻結晶溫控策略：自然冷卻（natural cooling policy，NCP）、線性冷卻（linear cooling policy，LCP）和受控冷卻（controlled cooling policy，CCP）。進行NCP 結晶時，維持結晶釜夾套中循環水溫度為固定溫度；LCP 和CCP 結晶時，控制結晶釜夾套中循環水溫度使釜內PCM 溶液溫度按式（5）冷卻至T1（=5℃）。

圖1 PCM 結晶ATR-UV 和FBRM 實時監控實驗系統Fig.1 Schematic of PCM crystallization using ATR-UV and FBRM real-time measurement

式中，量綱1 時間（τ）和溫度（Θ）分別為

2 實驗結果與討論

2.1 ATR-UV 實時濃度檢測模型的建立和校驗

圖2、圖3 分別為不同溫度和濃度的PCM 乙醇溶液ATR-UV 譜圖。由圖可知，PCM 乙醇溶液的ATR-UV 譜圖在波長（λ）255 μm 附近存在特征吸收峰。圖2 表明，當溶液溫度由60℃降低至40℃時，特征吸收峰位置由255.14 μm 變化至256.42 μm，可見溫度對ATR-UV 特征吸收峰位置的影響較小，故ATR-UV 實時濃度測定中可不對溫度進行校正。圖3 表明，當PCM 濃度從0.1686 g·g-1乙醇提高至0.3844 g·g-1乙醇時，特征吸收峰位置由251.09 μm紅移至256.42 μm。相對于溫度，PCM 濃度增加使特征吸收峰位置發生略微的紅移。

圖2 PCM 乙醇溶液在不同溫度下的ATR-UV 譜圖Fig.2 ATR-UV spectra of PCM in ethanol for varied temperatures

圖3 不同濃度PCM 乙醇溶液的ATR-UV 譜圖Fig.3 ATR-UV spectra of PCM in ethanol for varied concentrations

表3 PCM 在乙醇中的ATR-UV 濃度關聯參數Table 3 Parameters for correlating concentration of PCM in ethanol using ATR-UV spectra

為有效關聯不同溫度和濃度下PCM 乙醇溶液的ATR-UV 響應，用式（1）對包含操作范圍中特征峰位置范圍的信息進行關聯。通過如1.3 節所述方法，采用PLS 法關聯PCM 濃度，各參數關聯結果見表3，其RMSE 為0.0019。

圖4 ATR-UV實時濃度檢測結果與重量法濃度測定值的比較Fig.4 Comparison of PCM concentration measured by ATR-UV method and gravimetric method

在PCM 實際結晶操作中應用ATR-UV 進行實時濃度檢測，并在結晶開始后不同時間（例如8、21、33、40、52、75 和82 min）抽取少量的PCM乙醇溶液采用重量法測定其濃度，ATR-UV 測定結 果與重量法測定值的對比如圖4 所示。圖4 表明，采用ATR-UV 實時關聯的PCM 濃度與重量法測定結果一致，其RMSE 為0.0075，說明1.3 節所述建立的ATR-UV 濃度關聯方法可實現結晶過程PCM濃度的實時檢測。

2.2 溫度控制策略對PCM 冷卻結晶的影響

2.2.1 冷卻結晶溫控策略對晶粒產品形貌的影響 本文通過調節結晶釜夾套內循環水溫度變化實現對結晶溫度的控制，研究了NCP、LCP 和CCP 3 種冷卻結晶溫度控制策略對結晶產品質量的影響。PCM在乙醇中的結晶實驗如1.5 節所述，其中冷卻時間為180 min，終點為5℃。為便于比較溫度和濃度的變化，對不同時間的溫度及濃度進行了量綱1 化。3種不同的降溫方式結晶溫度和濃度隨時間的變化見圖5，其中PCM 的濃度采用ATR-UV 實時測定。

對于NCP 控溫，圖5 表明，由于結晶過程固含量較低，結晶熱效應不顯著，故NCP 的結晶溫度變化曲線近似于自然冷卻過程的溫度變化曲線，即：結晶初段系統與環境溫差最大而使溫度快速下降，隨后因與環境溫差的減少而降溫速率降低；結晶濃度曲線形狀與溫度曲線相似，即結晶初期濃度變化較快，隨后濃度變化逐步減少。這是由于結晶初期小粒徑晶粒數目較多有利于晶粒加速生長，從而加速了組分在液相中的消耗。隨著晶粒尺寸的增加，晶粒生長速度降低，同時NCP 中新晶粒的產生并不顯著，故組分濃度變化隨結晶進行而減少。

對于LCP 控溫，圖5 表明，體系的溫度變化具有較好的線性，表明本文采用的夾套溫控方案能較好地實現結晶過程溫度控制。在LCP 中，結晶前期的溫度變化比NCP 小，同時由于LCP 的溫度呈線性下降，LCP 結晶中段的溫度變化比NCP 更顯著。由圖5 可知，相對于NCP 過程，LCP 在結晶初期的PCM 濃度變化較小，表明LCP 結晶初段晶體成核及生長速度較NCP 慢；隨后量綱1 時間τ＞0.3（即結晶約60 min 后）時，LCP 濃度變化較逐漸比NCP 顯著，表明結晶中后期LCP 的結晶速度較 NCP 快。

圖5 不同降溫方式PCM 結晶量綱1 溫度和濃度曲線Fig.5 Dimensionless temperature (a) and concentration (b) curves of PCM crystallization

相對于LCP 的線性溫度下降，采用CCP 結晶過程的溫度變化呈現為先緩慢下降、后加速減少的變化趨勢，這種變化與NCP 過程截然相反。同時CCP 過程的濃度變化更為均勻，如圖5 所示，在3種冷卻策略中，CCP 過程的體系濃度演變曲線具有較好的線性，表明CCP 過程的結晶速率比較穩定。采用不同降溫策略所得的晶體產品形貌見圖6。偏光顯微鏡觀測發現NCP 和LCP 過程結晶產物的晶粒尺寸都較小，如圖6（a）所示，NCP 過程小尺寸晶粒出現了團聚現象；而LCP 中晶粒分散性較NCP有所改善，但發現晶粒尺寸變化較大，例如在圖6（b）中大顆粒晶體的尺寸比小顆粒大了5 倍以上。在3 種降溫策略中，采用CCP 過程的結晶產品具有 最大的晶粒尺寸和較一致的晶粒尺寸分布，如圖6（c）所示，由此說明CCP 是較好的結晶溫度控制方式。為進一步闡明CCP 能獲得更大、更均一晶體產品的原因，本文分別從過飽和度變化及晶粒CSD演變兩個方面深入說明降溫方式對結晶過程的影響。

圖6 不同降溫方式下PCM 冷卻結晶產品的顯微照片Fig.6 Micropolariscopy of produced PCM crystals using different cooling policy

2.2.2 降溫方式對結晶過程過飽和度變化影響 對特定體系的結晶過程，過飽和度是過程的推動力，其大小直接決定了結晶速度的快慢，故對結晶過程過飽和度變化的分析有助于從根本上理解結晶現象。借助ATR-UV 實時濃度關聯，可獲得不同降溫方式的過飽和度隨時間變化，如圖7 所示。

對于NCP 結晶過程，過飽和度在結晶開始約30 min 內迅速增加，而此時體系的濃度變化不大（如圖5 所示），表明此時可能處于PCM 結晶始發的誘導期；由于在高過飽和條件下可能引發晶體高速成核和生長，即類似“爆發結晶”過程，故過飽和度一旦突破其閾值將快速下降，因而在結晶后30～60 min 發現過飽和度急速下降，同時由于溶質在結晶過程的消耗導致了濃度顯著下降，這與NCP 濃度變化的測定結果一致。結晶約60 min 后的過飽和度由峰值的1.30 降低至1.10，這時結晶過程因較低的推動力而緩慢進行。

圖7 不同降溫方式下PCM 冷卻結晶過飽和度曲線Fig.7 Supersaturation curves of PCM crystallization using different cooling policy

相對于NCP，LCP前期同樣出現了過飽和度“快速增加然后迅速下降”的現象，這表明LCP 可能存在同樣的類似“爆發結晶”的過程。然而相對于NCP，LCP 的峰值過飽和度較低，這可能是由于結晶前期LCP 溫度變化比NCP 略小所致（如2.2.1 節所述），同時也正是LCP 在結晶后期具有比NCP 更大的溫度下降，導致了飽和濃度減少而使結晶約30 min 后的LCP 過飽和度比NCP 過程更高。

顯然NCP 和LCP 結晶前期出現的高過飽和度條件所產生的類似“爆發結晶”現象會產生大量的小尺寸顆粒，由此導致了圖6（a）、（b）所示較小的產品晶粒尺寸。此外，隨著結晶的進行，NCP 較LCP 因具有較小的溫度擾動和過飽和度變化而有利于晶粒團聚生長。

相比NCP 和LCP 結晶，CCP 過程過飽和度變化更為溫和，圖7 表明，CCP 結晶的過飽和度在1.08～1.13 范圍變化。特別是在結晶初期CCP 的溫度變化較小，由此使過飽和度維持在相對較低（S=1.13）的水平，這能有效地抑制NCP 和LCP出現的類似“爆發結晶”現象，從而有利于晶體的生長，獲得粒徑較大、分布較窄和晶體形貌較好的產品，如圖6（c）所示。

2.2.3 降溫方式對結晶產品CSD 的影響 進一步采用FBRM 實時分析不同降溫方式的晶粒CSD 變化，以其非加權結果識別晶粒的絕對量變化，結果如圖8 所示。由圖8 可知：（1）NCP 過程大部分晶顆是在0～30 min 內生成，此后晶粒數目（即實線CSD 曲線下的面積）隨時間增加而略微減少，而CSD 分布則變寬。這說明NCP 過程在結晶初期可能存在類似“爆發式結晶”現象，同時NCP 中存在嚴重的晶粒團聚現象，這使FBRM 將多個晶粒團聚體識別為單一大晶粒，從而測得晶粒數目隨時間而減少；（2）LCP 過程的CSD 曲線以點劃線表示，LCP在30 min 時晶粒數目因過飽和度較低而比NCP 少，隨后因LCP 具有較高的過飽和度不斷地成核結晶，由此使晶粒數目隨時間增加而增大，這與圖7 的結果一致。LCP 過程的CSD 曲線也表明晶體產品的粒度較NCP 的更小，這與圖6 觀察結果一致；（3）CCP 過程的CSD（如圖8 虛線所示）曲線下部面積隨時間推移而穩定增加，在0.5 h 時CCP 結晶過程顆粒數目明顯低于NCP 和LCP 結晶過程，隨著結晶過程的進行，CCP 結晶過程中晶體成核與生長穩定進行，使CSD 向平均粒徑增大方向移動且峰高增大，最終產品具有較大晶粒尺寸（與NCP 多晶粒團聚結果相當，而明顯大于LCP 晶粒尺寸）和更多的晶粒數量，同時CSD 分布也較NCP 和LCP 更窄。

圖8 不同降溫方式下PCM 冷卻結晶晶粒尺寸分布隨時間的變化Fig.8 Crystal size distribution in PCM crystallization using different cooling policy

不同降溫方式晶粒體積平均尺寸（d43）隨時間的變化如圖9 所示。由圖8、圖9 可知，NCP 結晶過程平均粒徑先增長后逐漸趨于穩定在較高水平，屬于小數目大粒徑分布；LCP 過程晶粒尺寸先增加，后因過飽和度較高成核增加而產生較多小尺寸晶粒，從而導致平均粒徑減小；CCP 結晶過程晶體尺寸生長較NCP 和LCP 更為穩定，其最終產品平均粒徑也最大，約為130 μm。

圖9 不同降溫方式下平均粒徑隨時間的變化Fig.9 Averaged size in PCM crystallization using different cooling policy

3 結 論

（1）基于ATR-UV 測量，采用偏最小二乘法進行濃度關聯，成功實現了乙醇-撲熱息痛體系結晶過程的實時濃度檢測，采用重量法對實時濃度檢測結果進行了校驗，其均方根平均偏差為0.0075。

（2）采用ATR-UV 實時濃度分析方法，考察了不同降溫方法（自然冷卻、線性冷卻和受控冷卻）對乙醇中撲熱息痛結晶過程的影響。結合結晶過飽和度的實時分析和FBRM 實時晶粒尺寸測量結果，發現受控冷卻結晶過程中晶體成核與生長能在更適宜獲得更大、更均一粒子的良性平衡中進行。

（3）采用實時過飽和度分析可建立基于過程分析的優化結晶控制系統，以實現質量源于設計的結晶產品調控。

符 號 說 明

A，B，C，D ——飽和溶解度關聯參數

Ai——吸光度

a, b, ci——ATR-UV 濃度關聯參數

d43——體積平均粒徑，μm

L ——晶粒尺寸，μm

n ——數據點數目

RMSE ——均方根平均偏差

S ——過飽和度

T ——溫度，℃

t ——時間，min

W ——樣品質量，kg

y ——質量濃度，g/g 溶劑

Θ ——量綱1 溫度

λ ——波長，μm

τ ——量綱1 時間

Ψ ——量綱1 濃度

上角標

* ——平衡濃度

下角標

R ——殘留量

S ——取樣量

0 ——結晶開始

1 ——結晶結束

* ——平衡濃度

[1]Wang Jingkang (王靜康).Crystallization (結晶)//Shi Jun (時鈞), Wang Jiading (汪家鼎), Yu Guocong (余國琮), Chen Minheng (陳敏恒).Handbook of Chemical Engineering (化學工程手冊) [M].Beijing: Chemical Industry Press, 1996.

[2]Yu L X, Lionberger R A, Raw A S, D'Costa R, Wu H, Hussain A S.Applications of process analytical technology to crystallization processes [J].Advanced Drug Delivery Reviews, 2004, 56 (3): 349-369.

[3]Simon L L, Pataki H, Marosi G, Meemken F, Hungerbühler K, Baiker A, Tummala S, Glennon B, Kuentz M, Steele G, Kramer H J M, Rydzak J W, Chen Z, Morris J, Kjell F, Singh R, Gani R, Gernaey K V, Louhi-Kultanen M, O’Reilly J, Sandler N, Antikainen O, Yliruusi J, Frohberg P, Ulrich J, Braatz R D, Leyssens T, von Stosch M, Oliveira R, Tan R B H, Wu H, Khan M, O’Grady D, Pandey A, Westra R, Delle-Case E, Pape D, Angelosante D, Maret Y, Steiger O, Lenner M, Abbou-Oucherif K, Nagy Z K, Litster J D, Kamaraju V K, Chiu M-S.Assessment of recent process analytical technology (PAT) trends: a multiauthor review [J].Organic Process Research & Development, 2015, 19 (1): 3-62.

[4]Palablylk i M, Onur F.The simultaneous determination of phenylephrine hydrochloride, paracetamol, chlorpheniramine maleate and dextromethorphan hydrobromide in pharmaceutical preparations [J].Chromatographia, 2007, 66 (1): 93-96.

[5]Cornel J, Mazzotti M.Estimating crystal growth rates using in situ ATR-FTIR and Raman spectroscopy in a calibration-free manner [J].Industrial & Engineering Chemistry Research, 2009, 48 (23):10740-10745.

[6]Duffy D, Barrett M, Glennon B.Novel, calibration-free strategies for supersaturation control in antisolvent crystallization processes [J].Crystal Growth & Design, 2013, 13 (8): 3321-3332.

[7]Frawley P J, Mitchell N A, O'Ciardha C T, Hutton K W.The effects of supersaturation, temperature, agitation and seed surface area on the secondary nucleation of paracetamol in ethanol solutions [J].Chemical Engineering Science, 2012, 75: 183-197.

[8]Hermanto M W, Chow P S, Tan R B H.Operating strategy to produce consistent CSD in combined antisolvent-cooling crystallization using FBRM [J].Industrial & Engineering Chemistry Research, 2012, 51 (42): 13773-13783.

[9]Mitchell N A, O'Ciardha C T, Frawley P J.Estimation of the growth kinetics for the cooling crystallisation of paracetamol and ethanol solutions [J].Journal of Crystal Growth, 2011, 328 (1): 39-49.

[10]Nagy Z K, Chew J W, Fujiwara M, Braatz R D.Comparative performance of concentration and temperature controlled batch crystallizations [J].Journal of Process Control, 2008, 18 (3/4): 399-407.

[11]Nagy Z K, Fujiwara M, Woo X Y, Braatz R D.Determination of the kinetic parameters for the crystallization of paracetamol from water using metastable zone width experiments [J].Industrial & Engineering Chemistry Research, 2008, 47 (4): 1245- 1252.

[12]O'Ciardha C T, Mitchell N A, Hutton K W, Frawley P J.Determination of the crystal growth rate of paracetamol as a function of solvent composition [J].Industrial & Engineering Chemistry Research, 2012, 51 (12): 4731-4740.

[13]Trifkovic M, Sheikhzadeh M, Rohani S.Kinetics estimation and single and multi-objective optimization of a seeded, anti-solvent, isothermal batch crystallizer [J].Industrial & Engineering Chemistry Research, 2008, 47 (5): 1586-1595.

[14]Trifkovic M, Sheikhzadeh M, Rohani S.Multivariable real-time optimal control of a cooling and antisolvent semibatch crystallization process [J].AIChE Journal, 2009, 55 (10): 2591-2602.

[15]Trifkovic M, Sheikhzadeh M, Rohani S.Determination of metastable zone width for combined anti-solvent/cooling crystallization [J].Journal of Crystal Growth, 2009, 311 (14): 3640-3650.

[16]Yu Z Q, Chow P S, Tan R B H.Application of attenuated total reflectance-Fourier transform infrared (ATR-FTIR) technique in the monitoring and control of anti-solvent crystallization [J].Industrial & Engineering Chemistry Research, 2006, 45 (1): 438-444.

[17]Saleemi A N, Rielly C D, Nagy Z K.Comparative investigation of supersaturation and automated direct nucleation control of crystal size distributions using ATR-UV/vis spectroscopy and FBRM [J].Crystal Growth & Design, 2012, 12 (4): 1792-1807.

[18]Doki N, Seki H, Takano K, Asatani H, Yokota M, Kubota N.Process control of seeded batch cooling Crystallization of the metastable α-Form glycine using an in-situ ATR-FTIR spectrometer and an in-situ FBRM particle counter [J].Crystal Growth & Design, 2004, 4 (5): 949-953.

[19]Sheikhzadeh M, Trifkovic M, Rohani S.Real-time optimal control of an anti-solvent isothermal semi-batch crystallization process [J].Chemical Engineering Science, 2008, 63 (3): 829-839.

[20]Abu Bakar M R, Nagy Z K, Saleemi A N, Rielly C D.The impact of direct nucleation control on crystal size distribution in pharmaceutical crystallization processes [J].Crystal Growth & Design, 2009, 9 (3): 1378-1384.

[21]Zhang Y, Jiang Y B, Zhang D K, Li K X, Qian Y.On-line concentration measurement for anti-solvent crystallization of β-artemether using UV-vis fiber spectroscopy [J].Journal of Crystal Growth, 2011, 314 (1): 185-189.

[22]Zhang R, Ma J, Li J, Jiang Y B, Zheng M Y.Effect of pH, temperature and solvent mole ratio on solubility of disodium 5-guanylate in water + ethanol system [J].Fluid Phase Equilibria, 2011, 303 (1): 35-39.

[23]Zhang Y, Jiang Y B, Li K X, Qian Y.Solubility of β-artemether in methanol + water and ethanol + water from (288.85 to 331.95) K [J].Journal of Chemical Engineering Data, 2009, 54 (4): 1340-1342.

[24]Samadi-Maybodi A, Hassani Nejad-Darzi S K.Simultaneous determination of paracetamol, phenylephrine hydrochloride and chlorpheniramine maleate in pharmaceutical preparations using multivariate calibration 1 [J].Spectrochimica Acta Part A: Molecular and Biomolecular Spectroscopy, 2010, 75 (4): 1270-1274.

[25]Worlitschek J, Mazzotti M.Model-based optimization of particle size distribution in batch-cooling crystallization of paracetamol [J].Crystal Growth & Design, 2004, 4 (5): 891-903.