第四章-固体溶液活度理论

1第四章固体溶液活度理论（Powell,1987;江培谟，1989；Mukhopadhyayetal.,1993；PowellandHolland,1993;Spear,1995；Will,1998）最紧邻规则、长程有序与短程有序、固溶体的理想活度、活度系数、正规溶液、Margules参数、零次近似模型、简单混合物模型、似化学模型、亚正规溶液、活度系数表达式、交互固溶体、Darken二阶表达式4.1概述固体溶液是地球化学和岩石学研究中非常重要的一个方面，固溶体模型也多种多样。从固溶体发生混合的晶格结点来说，既有元素在结点内的混合，也有元素在结点之间的混合；从固溶体中组分的混合性质来说，既有理想混合，也有非理想混合；从固溶体组分的活度模型来说，既有理想活度模型，也有非理想活度模型；从溶液模型来说，既有正规溶液模型，也有亚正规溶液模型；此外，还有描述包括相变情况的活度模型。可以说，固溶体活度理论既简单，也复杂，并且是个尚未得到完满解决的问题。固体溶液（例如矿物）的成分是对它进行热力学分析的基础。由其成分，可得到矿物的某些热力学性质。但是，要使成分与矿物的热力学性质联系起来，还必须选取适当的模型。例如，化学势与成分的关系，理想溶液模型的化学势-矿物成分关系，就与非理想溶液的矿物化学势-成分关系就不同。因此，只有使用合适的模型，才能得到满意的结果。4.1.1最紧邻规则为处理液态溶液，物理化学家创立了似晶格模型的溶液理论。这个模型是把液体看作分子都排列在一定的格子里，就象晶体那样。但是，液体毕竟不是晶体。因此，此模型就叫做“似晶格模型”。正因为此模型象晶体那样处理液体，2所以目前被广泛用来研究固体溶液。这个模型不考虑溶液中的库仑力，故只适用于非极性分子的混合物，不适用于高度极性分子的混合物。在处理过程中，只考虑一个分子与它最紧邻的Z（配位数）个分子（离子）的相互作用。我们知道，矿物固溶体中离子的相互替代会造成固溶体能量的变化，所以选取合适的固溶体模型描述这种能量变化是必要的。描述固体溶液中这种能量变化的最简单的模型就是所谓的“最紧邻模型”(Nearestnearmodel)。这个模型最早由Ising(1925)所提出，因而也叫做Ising模型。Ising(1925)在描述矿物晶体内铁磁性与反铁磁性时，只考虑了最紧邻分布的原子之间的相互作用(Will,1998)。以完全有序的NaCl晶体为例。图中距离中心离子(centralatom)A为0r的是6个离子B，因而离子A是六次配位的(six-foldcoordinated)。A占据了八面体位置(octahedralsite)。距离A为20r的是次紧邻(second/nextnearestneighbours)12个离子A，距离A为30r的是8个第三级最紧邻(thirdnearestneighbours)的离子B。在固体溶液模型中，仅考虑离子A和离它最紧邻的6个离子B之间的相互作用，距离A更远（即次紧邻、第三级最紧邻）的离子不再考虑。4.1.2长程有序与短程有序所谓“长程有序”(long-rangeorder)，指在整个矿物晶体内，离子的分布呈现有序的分布。虽然不等效的结点上离子的分布不同，但是相同的结点上各种离子的分布却是有规律的。矿物的这种有规律可循的晶体化学特征，是固体溶液热力学性质的分析基础。“短程有序”(short-rangeorder)则指在特殊的“成对替代”(coupledsubstitution)情况下，矿物晶体内某些结点上离子分布的有序特点。短程有序在成对替代中很重要，也是固体溶液显示非理想混合性质的重要原因。短程有序往往是由于相互替代的离子大小、电荷不同引起的。一般来说，晶体结点内离子的替代是无序的，例如，Fe2+、Mg离子电荷相同、大小相近，它们之间的替代是无限的、无序的。但是，在某些矿物内，情况就不是这样了。例如斜长石固体溶液，K和Na可以无限制地相互替代，但是K、Na和Ca离子之间的替代就不是无限的。Na+（或K+）替代Ca2+则造成一个单位的正电荷不足，因此需要Si4+同时再取代一个Al3+，以平衡电荷。这样，实际发生的替代是NaSi=CaAl或KSi=CaAl，这种成对替代叫做“Tschermak’ssubstitution”（契尔马克替代）。契尔马克替代引起斜长石内的短程有序现象。短程有序使得中心离子外围的最紧邻离子的选择性分布，即有些离子“偏好”这个位置（结点），而另一些离子则偏好其他位置。因而，随着PT变化，短程有序使得矿物固体溶液的有序程度发生变化，同时其混合的理想程度亦发生变化。例如，可使得矿物内易于形成A-A离子对和B-B离子对，而较难以形成A-B离子对。4.1.3分子型与离子型固体溶液及交互固溶体Bradley(1962)将固体溶液的混合模型分为“分子型”和“离子型”混合模型。分子模型是最简单的混合模型。假设矿物固溶体的混合发生于不同的完整分子之3间，因而叫做“分子型混合”(molecularmixingmodel)。这种模型不认为矿物固溶体的混合发生于离子之间。例如，变质成因的石榴石固溶体，可被看作是由镁铝榴石(pyrope,Mg3Al2Si3O12)、铁铝榴石(almandine,Fe3Al2Si3O12)、锰铝榴石(spessartine,Mn3Al2Si3O12)、钙铝榴石(grossular,Mg3Al2Si3O12)分子之间的混合而形成的。但是，对于固体溶液，一般不把分子型模型作为满意的模型来看待。实际上，固溶体矿物中组分的混合基本有以下三种方式：（1）仅在晶体中的一个结点(site)上，发生两种组分的混合。例如橄榄石(Fe,Mg)2SiO4中Fe、Mg的混合；（2）仅在晶体中的一个结点上，发生两种以上的原子的混合。例如，石榴石A位置上发生Fe、Mg、Ca、Mn四种原子的混合；（3）在晶体中多种结点上，各种原子各自独立的混合。例如，石榴石固溶体M23A3)Fe(Al,Fe)(Ca,中六面体(hexahedral)A位置上Ca、Fe原子的混合，以及八面体(octahedral)M位置上Al、Fe3+原子的混合，是相互独立的。这样的固溶体叫做“交互溶液”(reciprocalsolution)或“多结点溶液”(multi-sitesolution)。这三种情况的溶液模型，描述的都是矿物晶体中不同结点上原子的独立混合，都被称为“结点混合模型”（mixing-on-sites(MOS)model），也都是所谓的“离子型固溶体”混合模型。离子型固溶体模型客观地描述了矿物固体溶液的混合情况，也是最常用的或者是隐含使用的混合模型。4.1.4关于矿物固体溶液的活度许多造岩矿物(rock-formingmineral)属于固体溶液，其混合性质明显偏离理想混合，而显示非理想混合性质，所以其活度并不能用理想混合来刻画。从宏观上来说，应该采用活度系数γ来校正其非理想性质。对于溶液中的某组分i来说，其活度往往采用表达式iiiγXa。其中，Xi称为理想混合活度(idealmixingactivity)，又称为组分的热力学摩尔分数(thermodynamicmolefraction,Powell,1978)、当量摩尔分数(equivalentmolefraction,GangulyandSaxena,1987)、构型活度(configurationalactivity,Ghiorso,1984;Price,1985)。γi是活度系数。这种关系被习惯称为溶液的“组分—活度关系”(activity-compositionrelationship)或“a—X”关系(a—Xrelationship)。一般说来，活度系数是温度、压力和溶液组成的函数。固体溶液中原子的构型(configuration)明显影响溶液中组分的理想活度(idealactivity)和真实活度(realactivity)。实际上，溶液的理想混合熵，就准确地刻画了矿物固体溶液中的长程有序现象，但短程有序则不容易刻画。这是因为，短程有序难以测量，并且没有可资利用的标准方法来描述出现短程有序的真实矿物的混合性质。由于离子局部的选择性分布，离子混合的排列数无法计算。同时，由于没有可信的物理模型来描述短程有序现象，所以出现短程有序的矿物，其混合熵及活度表达式都无法准确写出。矿物固体溶液的活度系数，就可用来描述短程有序现象。或者说，活度系数把溶液偏离理想无序混合性质的所有非理想因素（包括短程有序），都一古脑地包括进去了。需要强调的是，表达式iiiγXa是个目前已经广泛接受的习惯作法，但却是多少人为的写法。44.2活度的标准状态溶液中某组分i的化学势定义为i0iiRTlnaμμ，或者说组分i的活度定义为)RTμμexp(a0iii。从活度的定义式知道，活度的计算需要标准状态，或者说离开标准状态活度就失去了意义。在地球化学和岩石学研究中，针对不同研究对象，采用有多种标准状态。下面分别讨论不同对象的活度的标准状态。4.2.1固溶体活度的标准状态固溶体中某组分i的活度的标准状态，通常定义为所研究的温度和压力下(atthetemperatureandpressureofinterest)的纯物质i。在这样的标准状态下，其化学势为i0iiRTlnaT)(P,μX)T,(P,μ。其标准态化学势与物质成分组成无关，而成为严格的该纯物质的摩尔Gibbs自由能，即0i0iGμ。标准状态下固溶体组分的活度为1。例如，自然产出的石榴石固溶体，其中的铁铝榴石组分的化学势与标准状态下纯铁铝榴石的化学势不同，其差值就由其活度来描述。4.2.2流体活度的标准状态超临界流体(supercriticalfluid)，例如H2O、CO2、H2O—CO2等流体，其标准状态一般定义为：P=1bar、纯物质、所研究的任意温度T。这样，某溶液中组分i的化学势为i0iiRTlnaT)bar,(1μX)T,(P,μ之所以如此定义，是因为在1bar和较高的的温度下，流体接近理想气体，所以该标准状态相当于理想气体。这样的标准状态下，活度将两项因素联系了起来：一项是实际流体成分对纯流体成分的校正，另一项是对气体非理想性质的校正。流体中某组分的活度为0iiif/fa，化学势为)/fRTln(fT)bar,(1μX)T,(P,μ0ii0ii=i0iRTlnfT)bar,(1μ。其中，0if为1bar、T、纯流体组分的标准状态下的逸度，通常接近1.0。这样定义流体组分标准状态后，体积对自由能的贡献部分，就只需从标准状态到压力P积分，即iiPPrRTlnfRTlnaVdP。例如，任意温度和压力下H2O—CO2流体中H2O组分的化学势，可表达为H2O0H2OH2ORTlnfT)bar,(1X)T,(P,4.2.3水溶液中溶剂活度的标准状态5在低温水溶液地球化学系统中，溶剂一般是H2O，其标准状态定义为：所研究的任意温度T、压力P、纯H2O。较通用的标准状态是T、1bar、纯H2O。因为大多数水溶液是稀溶液，溶剂接近于纯物质，所以其活度一般接近1.0。4.2.4水溶液中溶质活度的标准状态水溶液中溶质的标准状态，一般定义为：所研究的任意温度T、压力P、假想的纯溶质。之所以称这种标准状态为假想的，是因为，当电解质溶液接近纯溶质时，溶液不再表现为理想溶液，溶质的活度并不接近1.0。不过，当溶液成为无限稀溶液时，溶质的浓度接近活度：ii0mialimm，此即Henry定律。4.3固溶体的理想活度从前面的论述中已经知道，多种组分的理想混合过程产生的理想混合熵表达为iimixingideallnXXRS。如果在矿物晶体内的某一个结点上，发生多种元素的混合（这种情况称为“结点内混合”(intrasitemixing)，或“同一结点内混合”(same-sitemixing)），其理想混合熵同样表达为iimixingideallnXXRS。如果1摩尔单位的这种矿物中共有α摩尔个这样的结点，则其理想混合熵表达为iimix