化学工业与工程  2022, Vol. 39 Issue (4): 50-56
晶种制备的在线表征及钛液水解动力学研究
王焜1 , 路瑞芳2 , 党乐平1 , 郝琳1 , 卫宏远1     
1. 天津大学化工学院,天津 300072;
2. 钒钛资源综合利用国家重点实验室,四川 攀枝花 617000
摘要:针对硫酸氧钛水解晶种活性表征困难的问题,首次利用聚焦光束反射测量仪(FBRM)、原位拉曼光谱和动态光散射仪考察了晶种形成过程中溶液组成和粒子的变化,得到了有效的晶种表征手段, 并得出爆发成核之前为晶种制备最适宜保温时间; 考察了晶种加量和保温时间对水解反应动力学的影响,晶种量的增加会对水解速率和产率起到促进作用; 在爆发成核前取样,随着保温时间延长,后期水解速率和产率提高。水解动力学与成核-生长Avrami模型有较好的拟合度。
关键词水解    晶种制备    动力学    偏钛酸    
Online characterization of seed preparation and investigation on hydrolysis kinetics of TiOSO4 solution
WANG Kun1 , LU Ruifang2 , DANG Leping1 , HAO Lin1 , WEI Hongyuan1     
1. School of Chemical Engineering, Tianjin University, Tianjin 300072, China;
2. State Key Laboratory of Vanadium and Titanium Resources Comprehensive Utilization, Sichuan Panzhihua 617000, China
Abstract: In view of difficulties in seed characterization, Focused Beam Reflectance Measurement (FBRM), ReactRaman and Dynamic Light Scattering were applied for the first time to detect the changes in composition and particle properties of the solution during the preparation of seed, and a proper aging time of the seed crystal with rapid nucleation as the sign of critical time is obtained. Effects of volume and aging time of seed on the kinetics of hydrolysis were investigated. It turns out that hydrolysis rate and conversion are improved when the volume of seed increases or aging time of seed is prolonged before rapid nucleation. The kinetics curve fits well with Avrami nucleation-growth equation.
Keywords: hydrolysis    seed preparation    dynamics    metatitanic acid    

钛白粉在染料、化学催化剂、造纸、食品、塑料、化纤等领域有广泛的应用。工业中制备钛白粉的方法主要有硫酸法[1]和氯化法[2],而在中国,90%以上钛白颜料来自硫酸法生产。硫酸法制备钛白粉的步骤包括钛铁矿酸解、浓缩钛液水解、偏钛酸漂洗、煅烧、后处理等,其中,水解过程会影响偏钛酸的产率和粒度分布,继而影响后续煅烧产品钛白粉的产率和质量,是硫酸法制备钛白粉中重要的一步。

除了钛液质量和水解条件,晶种的数量和质量会对水解产生不容忽视的影响。根据晶种的制备方法,水解可分为自生晶种稀释水解法(布鲁门菲尔德法)和外加晶种水解法(麦克伦堡法)。目前的研究集中在自生晶种法[3-5],而外加晶种的研究较少。朱荣梅等[6]研究了外加晶种的加量对水解产率和一次团聚粒子大小及最终二氧化钛亮度的影响,得到了最佳晶种加量。由于晶种中微观粒子小,且不稳定,现缺乏系统的表征手段,对外加晶种的生成规律及晶种特性的报道几乎没有。

近年来对于钛液水解机理及动力学的研究较多[7-10],但对于水解机理还没有统一的解释。Wang等[11]探究了硫酸钛液水解机理,得到水解动力学模型,揭示了初始钛液、硫酸和水的浓度对水解的重要影响; 扈玫珑等[12]以硫酸钛为前驱物水解得到球形TiO2,探究前驱物浓度和水解温度对水解速率的影响,得到指数型水解动力学模型。但是目前关于晶种对水解过程动力学影响的研究较少。

本论文采用2种在线手段,聚焦光束反射测量仪(FBRM)和原位拉曼光谱,和一种离线手段,动态光散射仪,研究了晶种制备全过程在不同尺度的成核和生长规律; 考察了不同晶种加量和晶种保温时间对水解动力学的影响,系统解释了晶种对水解动力学的影响机理,对工业生产有一定指导意义。

1 实验部分 1.1 实验试剂与仪器

工业浓缩钛液(总钛含量为196 g·L-1, F值为1.90, 铁钛比为0.31,钛白粉厂提供); 氢氧化钠(分析纯,质量分数≥96%,购于天津江天化工试剂有限公司); 纯水通过优普纯水制造系统(UPY-Ⅱ-103,四川优普超纯科技有限公司)提纯。

晶种制备在平行反应器Easymax102(Mettler Toledo)中进行; 0.5 L夹套式四口玻璃结晶器用来进行水解反应。温度由数控恒温油浴箱(SC-15B,宁波新芝生物科技股份有限公司)来控制,机械搅拌保证一定速度下的均匀传质和传热。

1.2 实验方法

将工业钛液加入到氢氧化钠溶液中,95 ℃保温一定时间得到水解晶种(本论文中将加料结束后晶种在95 ℃下反应的时间定义为晶种保温时间)。

将一定量晶种加入到预热至95 ℃的水解钛液中,保温10 min,然后升温至沸腾,一沸至钛液变灰,停止加热,熟化30 min后再次沸腾,直到水解结束,得到锐钛型偏钛酸固体。

1.3 分析方法 1.3.1 晶种成核生长

晶种制备过程用原位拉曼(ReactRaman 785, Mettler Toledo)进行在线监测(波长785 nm,曝光时间为2 s,扫描次数为10次,扫描间隔为30 s),得到晶种浓度实时变化的定性曲线; FBRM(Particle Track G400, Mettler Toledo)在线表征溶液中微米级晶种粒子数目变化情况; 在晶种制备不同时间进行取样,利用动态光散射仪(DynaPro NanoStar, Wyatt)在25 ℃表征纳米级微观粒子分布情况。

1.3.2 水解动力学

在水解二沸的不同时间取样,过滤,用容量法测定不同晶种加量下、不同时间的水解率,绘制转化率曲线。

2 结果与讨论 2.1 晶种成核与生长

晶种为水解沉淀提供规则的结晶中心,晶种的质量、数量会影响偏钛酸的粒度分布和产率,从而对钛白粉的颜料性能、光学性能和产率产生影响。由于晶种的成核生长过程变化快、粒子粒度较小,难以检测,现有关于晶种的报道较少。本研究通过控制晶种制备过程,并利用在线测量手段,对晶种的化学成分、浓度变化及粒度变化,进行了系统研究。

拉曼光谱由于分析的快捷性和准确性,近年来被越来越多地应用于Ti(Ⅳ)在硫酸-水溶液体系中的表征分析[13, 14]。原位拉曼可以在无需取样的情况下,利用物质的拉曼散射现象,实时测量晶种溶液的拉曼图谱,包括溶液中溶剂、溶质及析出固体的图谱,其中出峰位置可定性分析溶液组成,峰强度可表征对应物质或基团的实时浓度。拉曼图谱的平滑、基线校准、峰高计算等操作在仪器自带软件iC RamanTM统一进行。

图 1,由于钛液中存370~380 g·L-1硫酸,所以晶种溶液的拉曼图谱中存在980 cm-1处的SO42-特征峰和890和1 050 cm-1处HSO4-的特征峰,590 cm-1处的峰同时为SO42-和HSO4-的特征峰[13]。仪器的激光发射接收探头的部件蓝宝石会在417和749 cm-1处出峰,分析时应引起注意。晶种溶液在152、515和636 cm-1处出峰,表现出与锐钛型偏钛酸一致的特征峰[15],可见晶种制备发生的主要化学反应为式(1):

$ \mathrm{TiOSO}_{4}+\mathrm{H}_{2} \mathrm{O} \rightarrow \mathrm{H}_{2} \mathrm{TiO}_{3} \downarrow+\mathrm{H}_{2} \mathrm{SO}_{4} $ (1)
图 1 晶种溶液与偏钛酸固体、硫酸溶液的拉曼谱 Fig.1 Normalized Raman spectra of seed and H2TiO3, H2SO4 solution

152 cm-1处特征峰表示锐钛型偏钛酸中Ti—O键的弯曲振动[16],其峰强度随时间的变化可以较好地反映出晶种浓度的变化趋势。如图 2所示,以加料结束时刻为0时刻,晶种制备可分为3个阶段:2 200 s之前为第1阶段,偏钛酸浓度接近直线上升,水解匀速进行; 2 200~2 250 s为第2阶段,偏钛酸浓度快速降低; 2 250 s之后为第3阶段,偏钛酸浓度又开始上升。

图 2 晶种制备过程152 cm-1特征峰相对强度随时间的变化 Fig.2 Variation of relative intensity of peak at 152 cm-1 with time during preparation of seed

图 3表示晶种制备过程中大于1 μm的粒子个数的变化情况。1 800 s之前,1 μm以上的总粒数较少,并保持平缓上升,结合图 3可知,1 800 s之前主要生成1 μm以下的小微粒; 1 800 s之后,1 μm以上粒子数发生爆发式增长,表明反应在不断生成小粒子的同时,小粒子结合为1 μm以上的粒子。2 200 s时粒子数目到达峰值,并开始减少,结合图 2,这时偏钛酸浓度也同步降低,由此可推测溶液中生成的晶种发生了溶解。

图 3 晶种制备过程大于1 μm粒子数的变化 Fig.3 Variation of numbers of particle(> 1 μm) with time during preparation of seed

晶种生成过程中,溶液中不断发生着沉淀-溶解竞争与平衡,由式(1)可知,随着晶种生成,酸度也在不断增大,偏钛酸的溶解度随硫酸含量发生变化[17],2 200 s时溶解度大于晶种浓度,晶种溶解; 晶种的溶解推动反应正向进行,2 250 s之后,粒子数又开始增加, 与图 3中晶种浓度升高相对应,晶种开始逐渐析出。2 400 s之后,20 μm以下的粒子减少,大于20 μm的粒子数增加,溶液中发生团聚现象。

由于仪器的检测限制,FBRM测得的溶液中均为1 μm以上粒子数目的变化。为了在纳米尺度对晶种进行研究,本研究采用动态光散射仪对其晶种制备过程进行表征。动态光散射仪在测量纳米颗粒粒径和粒径分布时,相比扫描电子显微镜(SEM),不需对样品进行离心等分离处理,具有更高的可靠度和准确性[18]。动态光散射仪利用粒子的布朗运动对散射光信号的影响(如图 4所示),得到粒子扩散常数。假定粒子是以相同速率扩散的球体,那么它的大小可定义为流体力学半径(本论文中简称半径)。利用Stokes-Einstein方程[式(2)]可以计算出粒子的流体力学半径Rh

$ R_{\mathrm{h}}=\frac{k_{\mathrm{B}} T}{6 \pi \eta D_{\mathrm{T}}} $ (2)
图 4 动态光散射仪原理示意图 Fig.4 Schematic diagram of dynamic light scattering

式(2)中:kB是玻尔兹曼常数,T是绝对温度(K),η是溶液的黏滞系数,DT为平移扩散系数。本研究利用动态光散射仪的优势,测量了晶种溶液中纳米级别粒子的半径及其分布情况。

根据Sekhar Sathyamoorthy等的研究[19],在沉淀过程中微粒尺径主要分布在3个区间,分别为晶粒(crystallites)、初级聚集体(primary agglomerates)和团聚体(aggregates)。4~8 nm的晶粒在晶桥(主要为SO42-)的作用下结合为60~100 nm的初级聚集体,这种结合比较牢固,不易被物理剪切等作用力破坏; 初级聚集体进而在物理吸附作用下聚集为1~2 μm的稳定性较差的团聚体。

在晶种生成过程不同时间取样,得到的粒子分布符合Sekhar Sathyamoorthy等的研究。如图 5所示,晶种晶粒的半径在2.3~4.2 nm之间波动,随着保温时间的延长,其体积分数从10.7%下降至1.8%;同时,初级聚集体半径由43.6 nm增大到98.0 nm,体积分数也从48.8%增大到95.1%;团聚体的半径在前30 min内在896.1~1 024.0 nm之间波动,体积分数由40.4%下降到3.4%,第40 min时,半径增大为7 385 nm,同时,其体积比例为3.1%。可以推断,随着水解的进行,溶液过饱和度逐渐增大,在成核的同时,晶粒在过饱和度的推动下形成初级聚集体。同时,由于表面电荷等的作用,团聚体之间的吸引力减弱,团聚体破裂减少。晶种中起作用的重要粒子为43.6~98.0 nm之间的初级团聚体,其平均粒径随时间增大。

图 5 晶种制备过程中纳米级粒子分布变化 Fig.5 Change of distribution of particle (in nanometer scale) size with time during preparation of seed

综合图 2图 3图 5数据可以推断,在晶种制备的前30 min为潜伏期,进行缓慢而匀速的成核,同时晶粒之间形成的初级团聚体在水解中起到重要的诱导作用。30 min之后发生爆发成核、溶解和团聚,晶种数量和质量变得不可控,所以建议在爆发成核之前完成晶种的取样。在工业生产中,由于实际工艺参数的变化,具体爆发成核时间会有变化,本研究为工业生产提供了表征手段和参数的参考。

2.2 晶种加量对水解动力学的影响

晶种保温10 min后,分别以体积为水解钛液体积的1%、2%、3%和4%的晶种诱导水解与不加晶种直接水解的水解率曲线进行对照。

图 6,转化率曲线都呈对数型增长,不加晶种转化率最低,为66.4%,随着晶种加量的增加,最终的转化率有所提高。晶种加量为1%时,最终转化率只有72.6%;当晶种加量为4%时,最终转化率达到98.4%。水解前期产率随时间增长较快,随着反应的进行,转化率的增速放缓,直至水解结束。在晶种加量较小时,前期水解速率较慢,水解不充分; 随着晶种的增加,前期水解速率有所提高,最终产率也随之增大。

图 6 晶种加量及水解时间对产率的影响 Fig.6 Effect of volume of seed and hydrolysis time on hydrolysis conversion

当不加晶种或晶种量过少时,溶液中没有足够的结晶诱导中心来加快沉淀和反应,导致前期水解速率降低,转化率减小。反之,当晶种量增大时,晶种对二次成核起到诱导作用,降低偏钛酸生长活化能,因此加快前期水解速率,提高转化率。

水解速率由成核和生长同时控制,可由成核-生长Avrami方程进行拟合。Avrami方程表示如式(3):

$ -\ln \left(1-\frac{\alpha}{100}\right)=k t^{n} $ (3)

式(3)中:α为转化率,t为时间,k为反应速率常数,n为Avrami指数,与形核机理和生长维数有关。将式(3)两边同时取对数得到式(4):

$ \lg \left[-\ln \left(1-\frac{\alpha}{100}\right)\right]=\lg k+n \lg t $ (4)

图 7,lg[-ln(1-α/100)]与lgt呈直线关系,由直线斜率可求的Avrami指数n,截距可求解得到反应速率常数k。求解得到动力学参数如表 1所示。

图 7 不同晶种加量下lg[-ln(1-α/100)]和lgt的关系 Fig.7 lg[-ln(1-α/100)] vs. lgt for different seed volume
表 1 不同晶种加量下Avrami方程的拟合参数 Table 1 Fitted parameter of Avrami model for different seed volume
晶种加量/% 转化率α/% 斜率n 截距lgk k/min-1 R2
0 66.4 0.917 -2.099 0.007 96 0.984
1 72.6 0.845 -1.866 0.013 61 0.993
2 86.8 0.859 -1.734 0.018 45 0.996
3 93.6 0.836 -1.495 0.031 99 0.966
4 98.4 0.891 -1.365 0.043 15 0.980

结果显示,拟合相关系数均大于0.96,说明Avrami模型可以较好地反映晶种加量对水解沉淀过程的影响。当不加晶种、晶种体积分数为1%、2%、3%和4%时,Avrami指数分别为0.917、0.845、0.859、0.836和0.891,反应速率常数对应分别为0.007 96, 0.013 61, 0.018 45, 0.031 99和0.043 15。Avrami指数在0.83~0.92区间内波动,表明晶种加量对水解成核生长模式有较小的影响。随着晶种量的增大,反应速率常数依次增大,说明晶种加量对水解动力学产生了一定的影响。

2.3 晶种保温时间对水解动力学的影响

为了研究不同性质的晶种对水解动力学的影响规律,本研究在晶种保温10、20、30和40 min后分别取样水解钛液体积的3%,加入水解反应中诱导沉淀。

图 8,在不同的晶种保温时间下,转化率曲线都呈对数型分布。晶种保温时间在30 min以内时,水解前期水解速率相近,随着晶种保温时间的延长,水解后期水解速率有所增大,最终转化率也相应提高,当晶种保温10 min时,最终转化率为93.6%,当保温时间延长到30 min时,转化率增大到98.7%。当晶种保温40 min时,前期水解速率较高,但后期水解速率降低,最终转化率降低为87.4%。

图 8 晶种保温时间及水解时间对产率的影响 Fig.8 Effect of aging time of seed and hydrolysis time on hydrolysis conversion

保温30 min以内时,晶种中主要发生缓慢成核,延长晶种保温时间相当于增加晶种量,增加水解结晶的诱导中心,因此提高了水解速率和最终转化率,又因为起诱导作用的初级聚集体粒径有所增加,所以与增加晶种加量引起的动力学变化不完全相同。在晶种保温30 min之后发生了爆发成核、溶解与团聚,导致晶种粒子大小不均匀,质量下降,降低了最终转化率。

对转化率曲线与Avrami模型进行拟合,得到lg[-ln(1-α/100)]与lgt的线性关系,如图 9所示。相关动力学参数系数如表 2所示。

图 9 不同晶种保温时间下lg[-ln(1-α/100)]和lgt的关系 Fig.9 lg[-ln(1-α/100)] vs lgt for different aging time of seed
表 2 不同晶种保温时间下Avrami方程的拟合参数 Table 2 Fitted parameter of Avrami model for different aging time of seed
晶种保温时间/min 转化率/% 斜率n 截距lgk k/min-1 R2
10 93.6 0.836 -1.495 0.031 99 0.966
20 97.8 0.986 -1.758 0.017 46 0.994
30 98.7 1.077 -1.909 0.012 33 0.998
40 87.4 0.589 -1.023 0.094 84 0.919

结果显示,拟合相关系数均大于0.91,说明Avrami模型对晶种保温时间对水解沉淀过程的影响有较好的拟合性。晶种保温时间为10、20和30 min,Avrami指数分别为0.836、0.986、1.077,都接近于1,表现出水解有着类似的成核生长模式; 反应速率常数分别为0.031 99、0.017 46和0.012 33,依次减小,说明晶种保温时间对水解动力学产生了一定的影响。当晶种保温时间为40 min时,Avrami指数为0.589,说明其成核生长模式发生了变化; 反应速率常数增大为0.094 84。

3 结论

研究了晶种制备过程中晶种浓度、微米级别粒子数目、纳米微观粒度分布的变化情况,结果显示,晶种制备前期缓慢成核,到达一定过饱和度时发生爆发成核、溶解和团聚,影响晶种质量。可用在线拉曼、FBRM 2种在线手段结合确定晶种爆发成核时间。增加晶种量,水解前期速率增大,转化率提高; 在爆发成核之前提高保温时间,水解后期速率增大,水解率提高。爆发成核之后的晶种会降低水解后期速率和水解率,改变水解的成核生长模式,影响水解质量。

致谢: 感谢攀钢集团研究院有限公司对本研究的资金和技术支持。

参考文献
[1]
HIXSON A W, PLECHNER W W. Hydrated titanium oxide. thermal precipitation from titanium sulfate solutions[J]. Industrial & Engineering Chemistry, 1933, 25(3): 262-274.
[2]
CHENG H, MA J, ZHAO Z, et al. Hydrothermal preparation of uniform nanosize rutile and anatase particles[J]. Chemistry of Materials, 1995, 7(4): 663-671. DOI:10.1021/cm00052a010
[3]
HAO L, WEI H. On-line investigation of anatase precipitation from titanyl sulphate solution[J]. Chemical Engineering Research and Design, 2010, 88(9): 1264-1271. DOI:10.1016/j.cherd.2010.01.008
[4]
TIAN C. Internal influences of hydrolysis conditions on rutile TiO2 pigment production via short sulfate process[J]. Materials Research Bulletin, 2018, 103: 83-88. DOI:10.1016/j.materresbull.2018.03.025
[5]
TIAN C. Effects of structural factors of hydrated TiO2 on rutile TiO2 pigment preparation via short sulfate process[J]. Scientific Reports, 2020. DOI:10.1038/s41598-020-64976-4
[6]
朱容梅, 陈葵, 朱家文, 等. 外加晶种对钛白粉水解过程及其亮度的影响[J]. 无机盐工业, 2019, 51(2): 30-34.
ZHU Rongmei, CHEN Kui, ZHU Jiawen, et al. Influence of crystal seeds amount on hydrolysis process of rutile titanium dioxide and its brightness[J]. Inorganic Chemicals Industry, 2019, 51(2): 30-34. (in Chinese)
[7]
ZENG F, LUO D, ZHANG Z, et al. Study on the behavior of sulfur in hydrolysis process of titanyl sulfate solution[J]. Journal of Alloys and Compounds, 2016, 670: 249-257. DOI:10.1016/j.jallcom.2016.02.031
[8]
TIAN M, LIU Y, ZHAO W, et al. Study on the grain size control of metatitanic acid in a mixture acid system based on Arrhenius and Boltzmann fitting[J]. RSC Advances, 2020, 10(2): 1055-1065. DOI:10.1039/C9RA08503C
[9]
TIAN M, LIU Y, WANG W, et al. Mechanism of synthesis of anatase TiO2pigment from low concentration of titanyl sulfuric-chloric acid solution under hydrothermal hydrolysis[J]. Journal of the Chinese Chemical Society, 2020, 67(2): 277-287. DOI:10.1002/jccs.201900071
[10]
PICHAVANT A, PROVOST E, FVRST W, et al. Determination of the temperature dependence of Titanium(Ⅳ) hydrolysis and complexation constants in aqueous sulfuric or chlorhydric solutions[J]. The Journal of Chemical Thermodynamics, 2019, 131: 184-191. DOI:10.1016/j.jct.2018.11.001
[11]
WANG W, LIU Y, XUE T, et al. Mechanism and kinetics of titanium hydrolysis in concentrated titanyl sulfate solution based on infrared and Raman spectra[J]. Chemical Engineering Science, 2015, 134: 196-204. DOI:10.1016/j.ces.2015.05.006
[12]
扈玫珑, 徐盛明, 白晨光, 等. 水解制备球形TiO2及其水解过程动力学[J]. 物理化学学报, 2009, 25(8): 1511-1516.
HU Meilong, XU Shengming, BAI Chenguang, et al. Preparation of spherical TiO2 and kinetics of its hydrolysis process[J]. Acta Physico-Chimica Sinica, 2009, 25(8): 1511-1516. DOI:10.3866/PKU.WHXB20090730 (in Chinese)
[13]
BAILLON F, PROVOST E, FVRST W. Study of titanium(Ⅳ) speciation in sulphuric acid solutions by FT-Raman spectrometry[J]. Journal of Molecular Liquids, 2008, 143(1): 8-12. DOI:10.1016/j.molliq.2008.04.010
[14]
SZILÁGYI I, KÖNIGSBERGER E, MAY P M. Characterization of chemical speciation of titanyl sulfate solutions for production of titanium dioxide precipitates[J]. Inorganic Chemistry, 2009, 48(5): 2200-2204. DOI:10.1021/ic801722r
[15]
吴健春, 任亚平, 王斌. 利用拉曼光谱定量分析金红石型二氧化钛[J]. 分析试验室, 2012, 31(12): 100-103.
WU Jianchun, REN Yaping, WANG Bin. Quantitative analysis of the content of rutile TiO2 by Raman spectrum[J]. Chinese Journal of Analysis Laboratory, 2012, 31(12): 100-103. (in Chinese)
[16]
MA W, LU Z, ZHANG M. Investigation of structural transformations in nanophase titanium dioxide by Raman spectroscopy[J]. Applied Physics A, 1998, 66(6): 621-627. DOI:10.1007/s003390050723
[17]
蒋银花, 殷恒波, 邱凤仙, 等. 酸在制备金红石型纳米TiO2中的作用[J]. 化工新型材料, 2006, 34(9): 46-48, 51.
JIANG Yinhua, YIN Hengbo, QIU Fengxian, et al. Effect of acid on the preparation of rutile TiO2 nanoparticles[J]. New Chemical Materials, 2006, 34(9): 46-48, 51. (in Chinese)
[18]
陈哲敏, 胡朋兵, 孟庆强. 动态光散射及电子显微镜纳米颗粒测量方法的比较研究[J]. 光散射学报, 2015, 27(1): 54-58.
CHEN Zhemin, HU Pengbing, MENG Qingqiang. Comparison study on DLS and SEM methods of measuring nanoparticle size[J]. The Journal of Light Scattering, 2015, 27(1): 54-58. (in Chinese)
[19]
SATHYAMOORTHY S, MOGGRIDGE G D, HOUNSLOW M J. Particle formation during anatase precipitation of seeded titanyl sulfate solution[J]. Crystal Growth & Design, 2001, 1(2): 123-129.