从原子尺度提示碳酸盐熔盐比热容增强机制
聚光太阳能（CSP）电站要靠熔盐把白天的热量储存起来、供夜间发电。但熔盐自身比热容偏低，储热密度一直是瓶颈。实验发现：往熔盐里只需加入 1% 的纳米颗粒，比热容就能提升近 20% —— 可“为什么”始终没人说清。最近，分子动力学模拟给出了原子级的答案。
背景：1% 纳米颗粒，比热容提升 20%
熔盐是 CSP 热储能（TES）的主力介质——热稳定性好、成本低、工作温区宽。实验上，向 Li₂CO₃-K₂CO₃（62:38）共晶盐中掺入 1 wt% 纳米颗粒后：SiO₂+19%、Al₂O₃ +21%、MgO +22%。研究者推测，纳米颗粒表面会生长出“盐基枝状纳米结构”，这些致密的纳米级结构域具有更高的比热容。但高温纳米尺度实验表征极其困难，机理长期停留在猜测阶段。
建模：两类模拟微结构对比
用 MedeA 构建两类模拟域，其余条件完全一致：
块体：44 Å × 44 Å × 44 Å 立方盒，10,800 原子，周期边界、体积可调；
纳米结构：缩小尺寸的同种盐，固定边界以规避周期自相互作用。
随后用 MedeA LAMMPS 求解：Lennard-Jones 势叠加长程库仑作用，键角用简谐势、二面角用 CVFF 力场；块体走 NVE 平衡后转 NPT，纳米结构走 NVT，时间步长 0.5 fs，温度范围 800–900 K。
关键构造：把 Li–O、K–O、Na–O 键长统一压缩至 1.792 Å，用以模拟纳米结构中的致密堆积。
图1 块体与纳米尺度结构：(a)(b)(c) 块体 Li₂CO₃ / K₂CO₃ / Na₂CO₃；(d)(e)(f) 对应的纳米结构
先验证：块体数据比热容数据验证
820 K 下，MD 给出三种碳酸盐的密度分别为 1.64 / 1.97 / 2.13 g/cm³，与高温实验值吻合；比热容分别为 2.01 / 1.06 / 1.39 kJ/(kg·°C)，与 NIST Shomate 方程预测值（2.18 / 1.16 / 1.46）及 Na₂CO₃ 实测值 1.43 高度一致——说明力场与流程可靠，计算准确。
图2（左）块体密度 MD 值；图3（右）块体比热容 MD 值与 NIST Shomate 方程对比
关键结果：致密化让比热容最高提升 24.1%
把离子键长压缩到 1.792 Å 后，三种熔盐的比热容同步跃升——Li₂CO₃ +19.4%（2.401 J/(g·°C)）、K₂CO₃ +24.1%（1.318）、Na₂CO₃ +23.3%（1.714），与实验掺杂观测到的 +19%~+21% 高度吻合。而只缩小尺寸、不压缩键长时，提升并不显著——证明增强源于“致密化”本身，而非单纯的纳米尺寸效应。
图4 块体 vs 纳米尺度比热容：短键长（致密化）带来 19.4% / 24.1% / 23.3% 的提升
机理：库仑能升高是机制
键长缩短 → 阳离子与氧的距离减小 → 由库仑定律 F = Kq₁q₂/r，库仑能显著增大。更强的静电相互作用限制了原子振动幅度、增强了非简谐性，体系需要吸收更多能量才能激发振动，于是储存更多热能、比热容升高。
完整链条：纳米颗粒表面 –OH → 吸引阳离子 → CO₃²⁻ 跟进 → 离子间距缩短 → 库仑能升高 → 振动谱扩展 → 比热容增强。
图5 短/长键长下的库仑能对比：短键长对应更高库仑能，验证强静电作用是增强的微观起源
本案例为纳米颗粒掺杂提升熔盐储热密度提供了原子级解释。MedeA + LAMMPS 可作为熔盐纳米流体比热容预测与机理研究的高效工作流，并可推广至氟盐、氯盐、硝酸盐等其它熔盐体系。
文献来源：Q. Altwarah et al. Enhanced Specific Heat Capacity in Molten Salt Nanofluids (2025).
推荐阅读案例
