我们在日常生活中遇到的大多数材料在受热时都遵循一个简单的规则:它们会发生混合。如果你搅拌油和水,它们会分离,但如果你加热某些塑料或聚合物的混合物,热量通常能提供足够的能量来克服它们天然的混合阻力,使它们融合成为一种单一、均匀的物质。这种行为如此普遍,以至于科学家长期以来一直认为这是复杂液体的一种默认状态,而发生混合时的温度被称为上临界溶液温度。然而,大自然还有一些奇招。有些特定的液体表现出截然相反的行为:它们在低温时保持混合,但在加热时却突然分离成不同的层。这种违反直觉的现象被称为下临界溶液温度。更为罕见的是,有些材料在进一步加热后会决定再次混合,从而形成一个闭环,即物质在低温和高温下都是分离的,而在中间温度下则是均匀的。几十年来,科学家们一直认为这些异常行为需要特殊的成分,例如能够彼此形成并断开特定化学键的分子,或者不同组分在尺寸或形状上存在结构失配的系统。
麻省理工学院的一个研究小组现在挑战了这一长期以来的观点,他们证明了这些复杂的温度响应可以源于最基础、最对称的构建模块。通过计算机模拟,研究人员模拟了多组分聚合物液体,其中每一个粒子在结构和尺寸上都是完全相同的,且粒子之间仅以一种简单的、不随温度变化的排斥力相互作用。在聚合物科学领域,这类完美的对称且无粘性的系统此前被预测总是会以标准方式运行:即加热时发生混合。然而,研究人员发现,当这些模型的密度足够低时,加热这一简单的动作会导致它们发生分离,在某些情况下,甚至会导致它们先分离然后再重新混合,从而形成那种难以捉取的闭环模式。这一发现的关键不在于特殊的化学键或结构失配,而在于粒子在稀疏环境中如何排列的一个基本属性。
要理解其原理,可以将此想象为拥挤的房间与近乎空旷的房间之别。在稠密液体中,粒子排列得非常紧密,不断地与邻居碰撞,创造了一个混沌但稳定的环境,其中温度变化的影响是可预测的。但在稀疏液体中,粒子彼此相距较远,粒子之间的空间变得至关重要。研究人员发现,在这些低密度条件下,随着温度的变化,粒子寻找彼此的方式发生了剧烈变化。在系统处于低温时,由于粒子因排斥作用而避开邻居,因此在相互作用范围内找到相邻单体的概率趋于零。随着温度升高,粒子获得能量,这增加了它们在有效相互作用范围内能触及的邻居数量。这种有效邻居数量的增加导致不同类型粒子之间的相互作用强度发生了非线性变化。这种相互作用强度并非仅仅随着热量的增加而减弱,而是首先增强(导致混合物分离),然后在极高温度下最终再次减弱(使得混合物重新融合)。
该团队通过创建两种类型的聚合物系统数字模型测试了这一想法:一种是两种不同类型长链的混合物,另一种是嵌段共聚物的熔体(由两个不同片段连接而成的链)。他们确保模拟中的每一个粒子在尺寸和形状上都是完全相同的,且粒子之间唯一的力是简单的、不随温度变化的排斥力。他们在各种密度和温度下运行了这些模拟。在高密度下,模型表现得完全符合预期,即加热时发生混合。但当他们将密度降低到特定阈值时,行为发生了逆转。在这些稀疏模型中,系统在低温时保持混合,随着升温而分离成不同的相,并且在处理具有较长键的嵌段共聚物时,甚至在极高温度下重新混合。研究人员证实,这种行为是由每个粒子在温度变化时所能触及的邻居数量的变化驱动的,这一因素只有在液体足够稀薄、使得这些强烈的空间相关性变得占主导地位时才会显现。
这一发现表明,现实世界材料中看到的奇异温度响应并不总是需要复杂的化学性质或结构差异。相反,它们可能是液体粗粒化模型在“有效”粒子密度较低时的一种自然结果。在这些模型中,单个粒子代表了一大群真实的原子,因此模拟中的低密度并不意味着真实的材料是气体;它仅仅意味着该模型是以特定的方式捕捉大型分子团簇的行为。研究人员指出,这种机制可以解释为什么一些聚合物共混物和溶液即使在看起来化学性质简单且对称的情况下,也会表现出下临界溶液温度或闭环相图。通过识别这一基本机制,这项研究为观察对温度做出反应的“智能”材料的行为提供了一个新的视角,暗示着调节这些属性的能力可能比我们之前认为的更加广泛地存在于现有材料之中。这项工作并不声称解决了聚合物相行为的所有谜团,但它确实揭示了:只要密度适中,最简单、最对称的系统也能产生最复杂的的热响应。
技术摘要:非缔合全对称多组分聚合物系统中的 LCST 与闭环相行为
问题陈述
在聚合物科学和液态理论中,低临界溶解温度(LCST)以及闭环相图(即在低温和高温下均发生相分离,并由上临界溶解温度 UCST 限制)的出现,传统上被归因于特定的机制。这些机制包括:
- 特定缔合作用: 存在随加热而断裂的方向性键(如氢键)。
- 结构不对称性: 组件之间在单体尺寸、形状或状态方程属性上的差异。
- 温度依赖性势能: 显式设计的、随温度变化的相互作用势。
因此,“全对称”系统——即组件在结构上完全相同、体积分数相等,且通过温度无关的纯排斥势进行相互作用的系统——在理论上被预测仅会表现出 UCST 行为(随加热而混合)。本研究旨在解决的问题是:在这些缺乏上述传统驱动因素的最简单粗粒化(CG)模型中(具体而言是无缔合、全对称的多组分液体,通过任意纯排斥势进行相互作用),是否也能涌现出 LCST 和闭环相行为。
方法论
作者结合了解析理论和分子动力学(MD)模拟,采用了双重研究方法:
理论框架:
- 研究重点在于单体为点状粒子、通过纯排斥且温度无关的非键合势 Unonbonded=aijv(rij) 进行相互作用的粗粒化(CG)模型。
- 作者推导了有效 Flory-Huggins 参数 χe 的温度依赖性,χe 决定了相行为。χe 被表示为能量失配 α=aAB−axx 的幂次的交错和,并由有效配位数分布 z 的累积量进行加权。
- 理论分析了两种单体数密度(ρ)的极限状态:
- 高 ρ(平均场)状态: 此时位置相关性较弱,单体间距离的分布在很大程度上与温度无关。
- 低 ρ(稀薄气体状)状态: 此时会出现强烈的成分间位置相关性。作者将系统近似为由不相连的 CG 单体组成的系统,推导了配位数 z∞ 以及截断半径附近的相互作用势 v(r) 的温度依赖性。
模拟:
- 作者使用 LAMMPS 软件在 NVT 系综下模拟了全对称二元共混物和二嵌段共聚物(diBCP)熔体。
- 模型利用谐振键合势和纯排斥谐振非键合势(这在耗散粒子动力学 DPD 中很常见)。
- 模拟在一定的温度范围(T∈[0.5,5])和密度范围(ρ∈[0.5,5])内进行。
- 相行为通过测量归一化峰值散射强度 S(qpeak) 并利用从导出的理论定义直接计算 χe(T) 来表征。
主要贡献与结果
发现了一种新的理论机制:
作者证明了,在全对称、非缔合系统中,LCST 和闭环相行为可以仅由于 CG 液体模型在低密度下的基本性质 而产生。
- 高 ρ 状态: 在高密度下,有效配位数 z∞ 对温度的依赖性很弱。因此,χe∝1/T,导致标准的 UCST 行为(随加热而趋向无序)。
- 低 ρ 状态: 在低密度下,强烈的单体-单体位置相关性开始发展。有效配位数 z∞ 表现出强烈的温度依赖性,其标度关系为 z∞∝T1/m(其中 m 取决于势函数的形状,对于谐振势通常 m=2)。
- 由此产生的 χe(T): 在低 ρ 极限下,Flory-Huggins 参数的线性项满足 χe(1)∝z∞(T)/T∝T1/m−1。对于 m=2,χe∝T1/2。这创造了一个 χe 随温度升高而增加的机制,从而驱动相分离(LCST)。随着温度进一步升高进入高 T 极限,系统会回归到 χe∝1/T 的状态,从而可能产生闭环相图。
模拟验证:
- 高密度模型 (ρ=5): 对二元共混物和 diBCP 熔体的模拟证实了理论预测的 UCST 行为,即 χe 随温度单调下降。
- 低密度模型 (ρ=0.7):
- 具有短键长(r0=1.0)的模型表现出明显的 LCST 行为:在加热时发生相分离,这与预测的低 T 状态下 χe∝T1/2 的标度关系一致。
- 具有长键长(r0=1.5)的模型表现出 闭环相行为:系统在加热时有序化(LCST),随后在更高温度下趋于无序(UCST),展示了理论所预测的再入转变(reentrant transition)。
- 模拟确认,从 UCST 到 LCST/闭环行为的转变发生在临界密度阈值(对于所用势函数,约为 ρ≈1.0)之下。
意义与主张
本文声称已经发现了一种导致多组分液体复杂温度响应的基本机制,这种机制并不依赖于特定的化学相互作用(缔合作用)或结构不对称性。
- 普遍性: 该机制源于在低单体数密度下出现的强单体-单体位置相关性,这是 任何 CG 液体模型固有的属性。
- 简洁性: 该效应在最简单的、广泛使用的 CG 模型(全对称、纯排斥、温度无关势)中被观察到,而这些模型此前被认为仅会表现出 UCST 行为。
- 对现有材料的意义: 作者指出,这种机制可能有助于解释许多现有聚合物材料(当它们被建模或表现为低密度 CG 系统时)中 LCST 和闭环相行为的出现。他们提到,CG 模型中的低 ρ 并不一定意味着真实的系统是气态的,因为一个 CG “单体”代表了一组真实的原子(例如在稀释聚合物溶液或特定的嵌段共聚物熔体中)。
这项工作得出结论:即使是最“简单”的模型也能够表现出“复杂”的行为,这挑战了以往认为需要不对称性或特定缔合作用才能解释非 UCST 相图的理论假设。作者对于该机制能否直接映射到特定的实验系统保持谨慎,指出确定该机制适用于哪些真实世界的精确范围,需要进行能够保留温度响应特性的系统性粗粒化研究,这也是未来的研究课题。
每周获取最佳 condensed matter 论文。
受到斯坦福、剑桥和法国科学院研究人员的信赖。
请查收邮箱确认订阅。
出了点问题,再试一次?
无垃圾邮件,随时退订。