1_2_丙二醇水溶液玻璃化转变与结构松弛
物理化学学报(WuliHuaxueXuebao)
206[Article]
ActaPhys.-Chim.Sin.,2007,23(2):206-211February
www.whxb.pku.edu.cn
1,2-丙二醇水溶液玻璃化转变与结构松弛
高
才1,*
2
王铁军1周国燕2华泽钊2
230009;
200093)
(1合肥工业大学机械与汽车工程学院农产品与生物化工教育部重点实验室,合肥
上海理工大学低温医学与食品冷冻研究所,上海
摘要:为了考察水含量对1,2-丙二醇水溶液玻璃化转变和结构松弛参数的影响,用差示扫描量热法(DSC),测
量了5种高浓度1,2-丙二醇水溶液(60%、70%、80%、90%、100%,w)玻璃化转变区域的表观比热容.用5种降温速率(1、・・2、5、10、20Kmin-1)和10Kmin-1的升温速率获得玻璃化转变的相关参数.玻璃化转变温度分析结果表明,虽然水含量增加能从总体上降低体系的玻璃化转变温度,但与纯羟基类多元醇相比,水对1,2-丙二醇的增塑作用并不显著.结构松弛活化能计算结果表明,体系水含量的增加能明显降低结构松弛活化能.脆度分析结果表明,随着体系水含量增加,动力学脆度逐渐降低,但热力学脆度是先升高后降低,在80%左右达到最大值.结构松弛协同重排域计算结果表明,当浓度由60%增加至100%时,玻璃化转变特征长度由2.79nm增加至3.57nm.
关键词:1,2-丙二醇;
O642.3;
玻璃化转变;结构松弛;脆度;特征长度;差示扫描量热法
中图分类号:O645.16
GlassTransitionandStructureRelaxationof1,2-Propanediol
AqueousSolutions
GAOCai1,*
WANGTie-Jun1
ZHOUGuo-Yan2
HUAZe-Zhao2
(1KeyLaboratoryofFarmProductandBiochemicalEngineering,MinistryofEducationofChina,CollegeofMechanicalandAutomotiveEngineering,HefeiUniversityofTechnology,Hefei230009,P.R.China;2InstituteofCryomedicineandFood
Refrigeration,ShanghaiUniversityforScienceandTechnology,Shanghai200093,P.R.China)
Abstract:Inordertoexaminetheeffectsofwatercontentsonglasstransitionandstructurerelaxationparametersof1,2-propanediol(PD)aqueousglasses,apparentspecificheatcapacityduringglasstransitionprocessesoffivePDaqueoussolutionswithhighconcentrations(60%,70%,80%,90%,and100%,w)weredeterminedusingthedifferential
・min-1)andoneheatingrate(10K・min-1)wereusedscanningcalorimetry(DSC).Fivecoolingrates(1,2,5,10,20K
toobtaintheglasstransitionandstructurerelaxationparameters.Waterseemedtobenoobviousplasticizationeffect
“purehydroxylgroup”on1,2-propanediol(comparedtosomelinearpolyalcohol)sincetheglasstransitiontemperatures
oftheaqueoussolutionswereonlyslightlydecreasedwiththeincreasingofwatercontent.Thestructurerelaxationactivationenergieswerefoundtodecreaseasthewatercontentincreased.AcomparisonofthermodynamicfragilityCp)anddynamicfragility(Angell′sfragility)waspresentandnosimilartrendwasfoundwiththechangingofthe(Δ
component,sincetheformerincreasedfirstanddecreasedlast.Thecalculatingresultsofdynamicfragilityindicated
“stronger”thanthepuresolute.Thesizesofthat1,2-propanediolaqueoussolutionsbehavedasthoughtheywere
cooperativerearrangementregion(CRR)discussedinthetermofglasstransitioncharacteristiclengthweredetermined
sthermodynamictemperaturefluctuantmethod.ItwasfoundthatthesizesofCRRdecreasedasthewaterusingDonth’
contentincreasingin1,2-propanediolaqueoussolutions.KeyWords:
1,2-Propanediol;
Glasstransition;
Structurerelaxation;
Fragility;
Characteristiclength;
Differentialscanningcalorimetry
Received:June14,2006;Revised:August28,2006.*
Correspondingauthor.Email:gao_cai@hotmail.com;Tel:+86-13721029540.国家自然科学基金(50376040、50206013)及合肥工业大学青年学术骨干启动资金资助项目
"EditorialofficeofActaPhysico-ChimicaSinica
No.2
高
才等:1,2-丙二醇水溶液玻璃化转变与结构松弛
207
过冷液体玻璃化转变与结构松弛是物理化学领域长期存在和活跃的课题[1].由于生物活体低温玻璃化保存和冷冻干燥技术要求人们了解低温水溶液玻璃的特性[2],故针对低温水溶液玻璃的研究不仅具有科学意义,还有重要的实用价值.例如,玻璃化转变温度的测量对于低温玻璃化保存中制定合理的降温和复温程序有重要指导意义.由于玻璃态是热力学非平衡态,会经历自发地向平衡态的结构松弛过程.而结构松弛伴随的微观结构的改变会影响玻璃化低温保存的质量[3],因此结构松弛的研究对于预测生物体和生物制剂玻璃化低温保存质量有重要意义.
皿中央.
1.3温度程序
先将样品从室温骤冷至-85℃,然后以1K・
min-1降至-155℃,待热流稳定后按如下温度程序
采集数据:(1)等温1min;(2)以10K・min-1升温到
-55℃;(3)等温1min.改变从-85℃降至-155℃时
的降温速率,重复以上实验,降温速率分别取2、
5、10、20K・min-1.取和装样品相同质量的空皿(±0.01mg),骤冷至-155℃,热流稳定后,按与分析样品相同的温度程序获得基线热流.1.4
数据分析
玻璃化转变温度用Pyrissoftware5.0(Perkin-
1,2-丙二醇是生物低温玻璃化保存和药品冷冻干燥中常用的玻璃形成剂[4,5],已有文献[6-8]对中低浓度1,2-丙二醇的水溶液进行玻璃形成趋势和玻璃稳定性
考察,但迄今为止未见有关结构松弛的相关报道.
本文选择在降温过程种不至形成冰晶的五种高浓度1,2-丙二醇水溶液,利用DSC技术考察水含量对体系玻璃化转变和结构松弛参数的影响.此前,我们已对乙二醇[9-11]和丙三醇[12,13]的水溶液体系进行过类似分析.由于1,2-丙二醇分子中除了两个羟基,还有一个自由甲基,二者均被认为是能够形成玻璃的官能团,因此有必要考察二者共同作用下体系玻璃化转变和结构松弛行为是否具有某些新特征.
Elmer)从比热容曲线上读取,比热容计算采用两线
法.实验结果取5个平行样平均值.因为平行样间的误差要远远大于重复样间的误差.故未对同一样品进行重复实验.为了便于数据比较,在进行玻璃化转变温度的读取和有关结构松弛参数计算时,均采用・10Kmin-1降温后的升温的表观比热容数据.
2
2.1
结果与讨论
玻璃化转变温度(Tg)
水-溶质二元系的玻璃化转变温度随着水分含
量增加而降低的性质,称为水的“增塑作用”.表1给出几种不同浓度PD水溶液的玻璃化转变温度.
可以看出,体系的玻璃化转变温度随体系中水分含量的增加呈总体下降之势.但同分子中只有羟基的丙三醇和乙二醇,以及分子中只有甲基的二甲基亚砜相比,水对PD的增塑效果并不显著.例如,当丙三醇水溶液的浓度由100%降至60%时,其玻璃化转变温度降低约30℃[12,13].对于乙二醇和二甲基亚砜水溶液,高浓度时容易形成溶质本身的结晶,故不能直接测定体系的玻璃化转变温度,但通过对已有数据[4]的外推,发现当浓度从100%降低至60%
表1
1
1.1
实验部分
试
剂
1,2-丙二醇为分析纯试剂(上海生化试剂有限
公司),未经进一步纯化而直接使用.用二次蒸馏水(自制)按质量分数配制成60%、70%、80%、90%、100%五种浓度(分别用PD60、PD70、PD80、PD90、PD100标明).文中所有浓度均指质量浓度.环戊烷为色谱标准试剂(上海生化试剂有限公司),用于分析仪低温标定.
1.2差示扫描量热仪
Pyris-DiamondDSC(美国Perkin-Elmer公司),配备CryoFillTM液氮灌注系统.温度标定采用环戊烷的晶-晶转变点(-135.06℃,降温)和纯水(二次蒸馏)的熔点(0℃,升温)进行两点标定,标定范围为-160-0℃;热焓标定采用纯水(二次蒸馏)的融化潜热(333.88mJ・・mg-1)进行单点标定.标定速率5Kmin-1.氦气(纯度>99.999%)气氛,流量20mL・min-1.样品
皿为标准液体皿,参比侧放与样品皿相同的空皿.实验中取试样量5-8mg(精确到±0.01mg),置于样品
PD水溶液的玻璃化转变温度(K)Glasstransitiontemperatures(K)of
PDaqueoussolutions
Tg(halfCpa)162.29162.53163.23164.81166.76
Tg(Tfb)161.00160.45160.78162.47166.01
Tg(onsetc)161.27161.99161.83163.51166.00
Table1
SolutionsPD60PD70PD80PD90PD100
a
temperatureathalfoftheΔCp;blimitedfictivetemperature;
c
onsettemperature.
208
ActaPhys.-Chim.Sin.,2007
Vol.23
时玻璃化转变温度的降低在10-12℃之间.而对于
PD水溶液体系,此值只有5℃左右.从表1中还可看出,PD水溶液体系玻璃化温度随水含量的增加不
是单调降低的,在浓度低于80%时玻璃化转变温度几乎不随体系中水含量的变化而变化.这一点有别于纯羟基类玻璃形成剂(乙二醇、丙三醇)和纯甲基玻璃形成剂(二甲基亚砜).如果取极限假想温度为玻璃化转变温度,会发现PD60玻璃化转变温度比
PD70的甚至还要稍高一些.已有较多关于醇类水溶
液玻璃形成机理的探讨[14-17],但尚未达成共识.通常认为是溶质分子中的官能团与水之间的相互作用,使得溶液变得粘稠,从而使水溶液中受扩散制约的冰晶生长过程变缓,玻璃化形成趋势增强.溶质分子中的羟基可以和水分子形成氢键网络结构,从而弱化水的结晶过程.尽管碳氢化合物(烃)是憎水的,但单个自由甲基却可以改变水分子的排列结构[18,19].
其中qc表示降温速率(K・min-1),Tf为极限假想温度.
图1给出不同降温速率下升温玻璃化转变区域的表观比热容变化.显然,较低降温速率对应玻璃化转变台阶上较大吸热峰,表明降温速率越低,则相同条件下玻璃的过剩焓释放得越多,重新升温后的焓恢复过程就越明显.
利用图1的比热容数据可直接读取不同降温速率下的极限假想温度,并按式(2)计算结构松弛活化能,如图2所示.可见,随溶液浓度的增加,拟合直线的斜率呈单调递减之势.当浓度由60%增加到
图1
PD90经历不同降温速率后10K・min-1
升温玻璃化转变过程比热容
Fig.1SpecificheatcapacityofPD90measuredsamplescooledatdifferentrates
duringheatingprocessesat10K・min-1afterthe
Wilson等[20]认为孤立甲基的参与改变了水分子之间的原有结构,增加了水分子之间的相互作用力.尽管
此前有人[7]认为PD水溶液体系在玻璃化转变方面具有“简单热行为”,但表1的结果表明,PD分子和水分子之间具有更为复杂的相互作用,这种羟基和甲基的联合作用使得水对体系的增塑效果变得不明显,甚至表现出玻璃化转变温度随水含量的反常变化.
2.2结构松弛活化能
玻璃的结构松弛是一个非线性和非指数的过
程[21].其中的非线性性质是指结构松弛时间不仅和温度有关,还与结构状态有关,可用TNM(Tool-100%,ΔH/R由约1.5×104K单调增加到约2.8×104
K,增加近一倍,这和丙三醇水溶液有类似之处[13].
但同聚合物和金属玻璃等相比,1,2-丙二醇水溶液玻璃的ΔH/R是很小的.例如,PVC(聚氯乙烯)的
Narayanaswamy-Moynihan)模型[22]描述:
x!"+(1-x)!"
!(T,Tf)=!0exp
f
"
(1)
其中!0是指前因子,表示平衡态液体在无限高温度下的松弛时间;!H表示结构松弛活化能(表观活化能);R是通用气体常数;x是非线性参数,将松弛时间!分为依赖于温度(T)和依赖于假想温度(表征结构信息的温度,Tf)两部分.结构松弛活化能不仅是进行理论分析的重要参数,还是计算玻璃脆性指标的基础.由于玻璃化转变温度和降温速率有依赖关系,而松弛时间和结构松弛活化能有(1)式描述的定量关系,因此通过考察玻璃化转变温度随降温速率的变化就可求取松弛活化能.Moynihan等[22]给出计算松弛活化能的如下公式:
ΔH/R高达2.25×105K[21].在文献[21]列出的数十种玻璃
中,ΔH/R低于3.0×104K的只有丙三醇(2.6×104K),乙
图21,2-丙二醇水溶液玻璃结构松弛活化能计算CalculationoftheΔH/RofPDaqueous
solutions
dln|q|=-f(2)
Fig.2
No.2
高
才等:1,2-丙二醇水溶液玻璃化转变与结构松弛
209
二醇(1.2×104K)和LiCl(1.2×104K)三种,绝大多数玻璃的ΔH/R都在7.0×104K以上.结构松弛活化能抽象描述了结构松弛过程中松弛单元微观结构重排势垒,因而也可看成是结构重排难易程度的度量.本文的计算结果表明,尽管PD水溶液中存在复杂的氢键网络结构,但同聚合物链段参与的结构松弛相比,氢键网络的结构重排势垒微不足道.Adam-
璃化转变过程中松弛单元自由度的释放是与松弛单元的数量和种类有很大关系的.例如离子液体因有阴离子和阳离子两个次级基团参与微观结构重排,故表现出比分子液体更大的热力学脆度[29].热分析方法虽不能定量确定玻璃体结构松弛单元的种类和数量,但热力学脆度的测量结果可用于定性判断结构松弛单元的变化.由于PD水溶液中,不仅存在PD-
Gibbs理论[23]认为,结构松弛活化能与协同重排域的尺寸有关.较大的协同重排域的尺寸通常对应较大
的结构松弛活化能,即结构松弛过程较难实现.据此,可判定PD水溶液中水分含量的增加可以减小结构松弛协同重排域的尺寸,从而使结构松弛过程更容易进行.
H2O,H2O-H2O之间的氢键作用,还有PD-PD之间的氢键[30]相互作用,这种复杂的氢键网络结构使得结
构松弛单元呈现多样性.随着浓度的增加,溶质分子之间的键合几率增加,完全有可能改变松弛单元的种类和数量,从而使热力学脆度表现出如图3的变化行为.PD水溶液动力学脆度随水含量的增加而逐渐降低,和此前我们对丙三醇水溶液体系[13]的分析结果很相似.动力学脆度实际上表征了玻璃结构抵抗外界温度变化的能力,其值越小,则当玻璃体升温通过玻璃化转变区后微观结构的改变也越小.图3表明,PD水溶液玻璃的含水量越大,当玻璃升温通过玻璃化转变区后,会有更多原有结构位形得以保留.
2.3脆度
液体脆性的概念由Angell[24]提出.是玻璃物理中最重要的参数之一.根据物质的动力学性质(如粘度和松弛时间)在玻璃化转变区域随温度变化的情况下,Angell将过冷液体分为强液体和弱液体.强液体的结构对温度变化不敏感,在玻璃化转变温度以上表现出Arrhenius温度依赖关系.而弱液体则相反,表现为随着温度的升高,粘度(或松弛时间)呈快速下降趋势.用于衡量液体脆性的指标称为脆度,其定义为[25]
2.4协同重排域
玻璃结构松弛协同重排域的概念由Adam和
m=
"lg!(T)gf(3)
T=Tg
Gibbs两人提出[23].协同重排域大小(V#)可用玻璃化
3
转变特征长度$#表示(V#=$#),是表征玻璃化转变异
质性的重要指标.Sillescu[31]的综述表明,对于绝大多数玻璃而言,$#的范围是1-5nm.Donth[32,33]提出基于热力学温度波动理论的计算方法:
3%)V#=$#=kBT!(1/C!(1/Cv)=(1/Cv)glass-(1/Cv)liquid
2
≈(!Cp2p)(1+(1/4)(!Cpp)+…)
其中!(T)为结构松弛时间.上述定义式在量热方法中的具体形式[26]为
m=-
!"
g
(4)
2
(5)
式中!H是结构松弛活化能(kJ・mol-1),R是通用气体常数,Tg是玻璃化转变温度(K).m表征了玻璃化转变的动力学信息,因此也称为动力学脆度,与之对应的热力学信息则用玻璃化转变过程的比热容变化表征,称为热力学脆度[27].热力学脆度表征了升温玻璃化转变过程中自由度的释放程度,越大则表明玻璃体过渡到过冷液体后松弛单元获得的活动性就越大.统计表明,一般情况下,热力学脆度越大,则其动力学脆度也往往越大[28].这很容易理解,因为玻璃体微观结构的改变过程实际上也是松弛单元重新获得活动性的过程.图3给出PD水溶液热力学脆度和动力学脆度的变化关系.显然,二者并未表现出同增同减的变化趋势.热力学脆度是随着浓度的增加先升高后降低,浓度为80%左右达到最大值.通常,玻
(6)
图31,2-丙二醇水溶液热力学脆度和动力学脆度Thermodynamicfragility(!Cp)anddynamicfragility(m)ofPDaqueoussolutions
Fig.3
210
ActaPhys.-Chim.Sin.,2007
Vol.23
图4PD80的玻璃化转变特征长度计算
Fig.5
图5
几种玻璃结构松弛活化能和!"对照[21,34]
Fig.4TheschemeforthedeterminationoftheComparisonofstructurerelaxationactivation
energiesand!"forsomeglasses
quantitiesneededfromDSCexperimentsforPD80
其中Cv和Cp分别表示定容比热容和定压比热容(J・・・K-1g-1);!是密度(gcm-3);"T2表示协同重排域的均方温度波动(K2),"T=!T/2.5,!T是由比热容增量决定的温度跨度(K),如图4所示;p表示比热容平均值.计算中需要的各参量可以完全由DSC比热容数据获得.图4给出玻璃化转变特征长度计算所必需的各参量的示意(样品为PD80).五种溶液的#$计算结果如表
考察了相容聚合物的动力学脆度m和%$的相关性,发现m值较大的玻璃,其%$值也往往较大.对照图
3和表2的数据,可知这一规律也适合本文的情形.
3结论
利用差示扫描量热法,测量了五种浓度的PD水溶液玻璃化转变区域的表观比热容,进行玻璃化转变行为和结构松弛参数分析.结果表明,水对1,
2所示.
总体而言,体系水含量的增加能降低协同重排域的尺寸,这和上文结构松弛活化能的计算结果有很好的吻合.文献[34]列出具有代表性的28种玻璃的
%$值,其中最大的是山梨醇的3.6,最小的是PnHMA(聚乙基丙烯酸甲酯)的1.0,半数玻璃的%$值都在2.0以下.故同其他玻璃相比,PD水溶液玻璃有相对较高的%$值.为考察不同玻璃的%$和结构松弛活化能的相关性,我们选择了%$值相近的几种玻璃[21,34],与
本文研究的几种玻璃进行对照,如图5所示.可以看出,在%$相当的情况下,PD水溶液的结构松弛活化能是最低的.丙三醇和PD水溶液的%$最接近.这其中PVC比较特殊,尽管其结构松弛活化能比PD水溶液大很多,但其%$只和PD70、PD90相当.Cowie等
表2
[35]
2-丙二醇的增塑作用并不显著,尤其是当浓度低于
80%时,体系的玻璃化转变温度几乎不随体系水含量的变化而变化.水含量的增加能明显降低1,2-丙二醇水溶液的结构松弛活化能,表现为60%的1,2-
丙二醇水溶液的结构松弛活化能只有纯溶质的
50%左右.在文中选择的浓度范围内,1,2-丙二醇水
溶液溶液的动力学脆度是随着浓度的增加而增加的,但热力学脆度却表现出先升高后降级的变化趋势,在浓度为80%左右时达到最大值.协同重排域的分析结果表明,浓度越高,则体系的玻璃化转变的特征长度就越大.同其它玻璃相比,PD水溶液玻璃有相对较低的结构松弛活化能和相对较大的协同重排域.
References
1
Angell,C.A.;Ngai,K.L.;McKenna,G.B.;McMillan,P.F.;Martin,S.W.J.Appl.Phys.,2000,88(6):31132
Franks,F.Biophys.Chem.,2003,105:25
freeze-dryingandstorage:therelativeimportanceofthethermody-namicsstabilizationandglassystaterelaxationdynamics.In:Freeze-dryingofpharmaceuticalandbiologyproducts.LouisRayetal.ed.NewYork,Baseal:MareclDekker,Inc.19994
Hua,Z.Z.;Ren,H,S.Croybiologyandcryomedicinetechnology.3Pikal,M.J.Mechanismsofproteinstabilizationduring
1,2-丙二醇水溶液的玻璃化转变特征长度Glasstransitioncharacteristiclengthsof
PDaqueoussolutions
・・・・Tg/K!Cp/JK-1g-1p/JK-1g-1
Table2
SolutionsPD60PD70PD80PD90PD100
"T/K2.2561.6241.8002.0401.384
%$/nm2.7893.1103.2953.1783.565
161.00160.45160.78162.47166.01
1.1311.2021.2371.1591.047
1.9792.3712.0221.8442.170
No.2
Beijing:SciencePress,1994术.北京:科学出版社,1994]
高
才等:1,2-丙二醇水溶液玻璃化转变与结构松弛
Chem.,1981,85:123619
211
[华泽钊,任禾盛.低温生物医学技
deVisser,C.;Heuvelsland,W.J.M.J.Chem.Soc.FaradayTrans.,1978,74:1159
5Hua,Z.Z.Newtechnologyoffreezedrying.Beijing:SciencePress,20062006]
[华泽钊.冷冻干燥新技术.北京:科学出版社,
20Wilson,P.W.;Arthur,J.W.;Haymet,A.D.J.Biophys.J.,1999,
77:2850
21Hodge,I.M.J.Non-Cryst.Solids,1994,169:211
22Moynihan,C.T.;Easteal,A.J.;DeBolt,M.A.;Tucker,J.J.Am.
CeramicSoc.,1976,59:12
23Adam,G.;Gibbs,J.H.J.Chem.Phys.,1965,43:13924Angell,C.A.J.Phys.Chem.Solids,1988,49:86325Angell,C.A.J.Non-Cryst.Solids,1991,131/133:1326Angell,C.A.Science,1995,267:1924
27Bohmer,R.;Angell,C.A.Disordereffectsonrelaxational
processes.Richert,R.ed.,Berlin:Springer,19942829
Huang,D.;McKenna,G.B.J.Chem.Phys.,2001,114(13):5621MouraRamos,J.J.;Afonso,C.A.M.;Branco,L.C.J.Therm.Anal.Cal.,2003,71:659
30Crupi,V.;Longo,F.;Majolino,D.;Venuti,V.J.Molecular
Structure,2006,790(1-3):141
31Sillescu,H.J.Non-Cryst.Solids,1999,243:81.3233
Donth,E.J.Non-Cryst.Solids,1982,53:325.
Donth,E.Relaxationandthermodynamicsinpolymersglasstransition.Berlin:Akademie-Verlag,1992
34Hempel,E.;Hempel,G.;Hensel,A.;Schick,C.;Donth,E.J.Phys.
Chem.B,2000,104:246035
Cowie,J.M.G.;Harris,S.J.;GomezRibelles,L.;Meseguer,J.M.;Romero,F.;Torregrosa,C.Macromolecules,1999,32:4430
6789
Boutron,P.;Kaufmann,A.Cryobiology,1979,16(6):557Mehl,P.M.Thermochim.Acta.,1992,203:475
Jabrane,S.;Letoffe,J.M.;Claudy,P.Thermochim.Acta,1998,311:121
Gao,C.;Zhou,G.Y.;Xu,Y.;Hua,T.C.Thermochim.Acta,2005,435:38
10Gao,C.;Zhou,G.Y.;Xu,Y.;Hua,Z.Z.ActaPhys-Chim.Sin.,
2004,20(2):123
[高
才,周国燕,胥
义,华泽钊.物理化学
学报,2004,20(2):123]
11Gao,C.;Wang,W.H.;Hu,T.J.;Xu,Y.;Zhou,G.Y.;Hua,Z.Z.
ActaPhys.-Chim.Sin.,2004,20(7):701记,胥
[高才,王文华,胡桐
义,周国燕,华泽钊.物理化学学报,2004,20(7):701]
12Gao,C.;Zhou,G.Y.;Xu,Y.;Hua,Z.Z.Chin.J.Chem.Phys.,2005,18(3):457
[高
才,周国燕,胥
义,华泽钊.化学物理
学报,2005,18(3):457]
13Gao,C.;Zhou,G.Y.ActaPhys.-Chim.Sin.,2005,21(8):909
[高
才,周国燕.物理化学学报,2005,21(8):909]
14Macfarlane,D.R.Cryobiology,1987,24:18115Macfarlane,D.R.Cryobiology,1990,27:345
16Fahy,G.M.;Levy,D.I.;Ali,S.E.Cryobiology,1987,24:19617Vigier,G.;Vassoille,R.Cryobiology,1987,24:34518
Spencer,J.N.;Berger,S.K.;Powell,C.R.;Henning,B.D.J.Phys.