基于经验势分子动力学模拟的N719染料/二氧化钛界面在室温离子液体溶剂中的振动光谱
2020年1月25日
采用室温离子液体（RTILs）对染料敏化太阳能电池进行溶剂化处理；通过优化的经验势函数，应用分子动力学模拟计算其振动特性。将所得振动光谱与实验结果及从头算分子动力学结果进行比较；多种经验势函数的光谱表明，离子液体的部分电荷参数化如何影响振动光谱的预测。
本方案旨在探讨经验势方法在典型染料敏化太阳能电池体系的结构和振动特性合理预测方面的适用性。这一点至关重要，因为经验势方法具有显著降低计算成本的优势，并可与多种采样方法相结合。要进行分子动力学模拟，请打开 DL_POLY 模拟软件文件，并使用 Lopez 等人已充分验证的力场。
通过经验势函数对室温离子液体构型进行弛豫。然后在弛豫过程中采用包含二氧化钛迁移性的Matsui-Akaogi力场来模拟锐钛矿结构。在NVT系综下，设定场文件中的优化参数，采用共轭梯度最小化方法进行几何优化，相对终止梯度为0.0001，模拟时间为15皮秒，时间步长为1飞秒，温度为300开尔文。
在终端中运行 DL_POLY，请输入 DL_POLY.X。在控制文件中，为 NVT 系综指定 Nose-Hoover 方法，并选择每 1 飞秒输出一次位置和速度轨迹。对于锐钛矿表面，需确保其沿实验室坐标系的 x 和 y 轴具有周期性，从而在室温离子液体中形成一对平行表面，其中 x 轴为 21 埃，y 轴为 21 埃。
确保包含显式溶剂的染料敏化太阳能电池体系由827个原子组成。对于液泡内的体系，应包含347个原子。在Lennard-Jones参数的场文件中，应用Lorentz-Berthelot混合规则，并在场文件的非键相互作用选项卡下方部分，填入Lennard-Jones半径的算术平均值和Lennard-Jones势阱深度的几何平均值，用于经验力场。
为处理长程静电相互作用，采用埃瓦尔德（Ewald）方法，设置非键截断长度为 10 埃，并在控制文件中将埃瓦尔德方法的实空间衰减参数设为非键截断长度的约 3.14 倍。然后在控制文件中选择适当的埃瓦尔德波矢数量，以确保埃瓦尔德求和计算的相对容差达到 1×10⁻⁵。通过一系列势能计算，将重命名为 config 的 REVCON 文件作为输入，持续运行直至输出结果中的体系压力收敛至百分之几的范围内，以确定合适的非键截断长度。
从历史文件中，在终端使用 Python 脚本 diatom velocity separate.py 提取每一步的 x、y 和 z 方向速度。为计算速度自相关函数，使用 vacf151005。
py 并在终端中输入 classical_dye_autocore.sh。将针对所有硅藻计算速度自相关函数。
为了从分子动力学中计算光谱，使用染料原子速度自相关函数的质加权傅里叶变换，运行 Python MWPS.py 并在终端中输入 run_all_4.sh。
将计算质量加权功率谱。为了便于相互比较，在针对室温离子液体的帝国势基分子动力学模拟中评估部分电荷集的变化时，需准备一份文献电荷表格，并将该表格输入DL-POLY的场文件格式中。为计算Mulliken室温离子液体电荷，需对从从头算分子动力学轨迹中选取的四个时间点的电荷进行平均，重新归一化后，制备文献电荷表格，并将该表格输入DL-POLY的场文件格式中。
为了进行EHT分析，采用分子模拟环境软件包（molecular operating environment software package）中实现的方法，对从头算分子动力学轨迹中的四个点取平均值。重新归一化文献中的电荷并制表，并将该场数据转换为DL-POLY的场文件格式。此处展示了分子动力学模拟15皮秒后，四种不同部分电荷集的代表性结合模式。
如所示，根据文献推导的电荷可见与表面质子存在显著的氢键相互作用。然而，三种由从头算分子动力学推导的电荷集并未表现出与表面质子如此强的库仑相互作用。事实上，与通过不同方式从从头算分子动力学采样得到的部分电荷相比，文献推导情况下的部分电荷绝对值较小，导致其相对于室温离子液体中较大绝对值的部分电荷，电荷屏蔽效应较弱。
有趣的是，基于Mulliken方法得到的电荷分布显示，染料分子呈现出一定程度的持续弯曲构型，从而与锐钛矿表面的桥接氧原子形成显著的氢键。相较于从头算分子动力学获得的更高质量电荷拟合结果，使得N719的结合模式更为真实，并与基于PEV的Born-Oppenheimer分子动力学模拟结果一致，后者包含了Grimme-D3色散修正。质量加权的从头算分子动力学速度自相关函数光谱结果显示，在300至400 cm⁻¹波数范围内出现了一组光谱峰。
在经典光谱中，文献导出电荷的主要峰位于 600 和 800 厘米处，Mulliken 电荷位于 525 和 800 厘米处，EHT 电荷位于 675、810 和 900 厘米处，Hirshfeld 电荷集则位于 650、800 和 900 厘米处。务必检查二氧化钛几何结构在所有方向上的初始构型。为进行优化，可使用力场来评估分子动力学预测的振动光谱中的经验势能。
从头算分子动力学是一种有效且前景广阔的策略。分子动力学的使用在模拟辅助染料敏化太阳能电池原型方面也可能具有重要意义
查看完整文字稿并访问数千部科学视频
概述
本文介绍了一种利用基于经验势的分子动力学（MD）模拟来预测染料敏化太阳能电池（DSC）界面结构和振动特性的详细实验方案，重点研究N719染料在锐钛矿相TiO2表面的吸附及其被室温离子液体（RTIL）溶剂化的体系。该研究评估了RTIL不同部分电荷参数化方案对振动光谱预测结果的影响，并将经验分子动力学模拟结果与从头算分子动力学（AIMD）模拟及实验数据进行了比较。
关键研究要素
科学领域
- 物理化学
- 计算化学
- 材料科学
- 分子模拟
背景
- 准确预测吸附染料的光活性金属氧化物表面的振动光谱具有挑战性。
- 染料敏化太阳能电池（DSCs）在太阳能转换方面具有前景，其性能受界面特性影响。
- 室温离子液体（RTILs）被用作DSCs中的电解质，影响染料与表面之间的相互作用。
- 经验势分子动力学（MD）具有计算效率高的优势，但需要谨慎的参数化和验证。
研究目的
- 评估经验势分子动力学（MD）在预测DSC界面振动光谱和结构特性方面的准确性。
- 研究不同RTIL分电荷集对振动光谱预测结果的影响。
- 将经验MD结果与AIMD及实验光谱进行比较，以验证其有效性。
使用的方法
- 使用经验证的力场（RTIL 采用 Lopez 等人的力场，锐钛矿 TiO2 采用 Matsui-Akaogi 力场），在 DL_POLY 中进行分子动力学模拟。
- 在 300 K 下进行几何优化和 NVT 系综模拟。
- 对 Lennard-Jones 参数采用 Lorentz-Berthelot 混合规则，对长程静电相互作用采用 Ewald 方法。
- 提取原子速度并计算速度自相关函数。
- 对质量加权的速度自相关函数进行傅里叶变换，以获得振动光谱。
- 利用不同来源的部分电荷集（来自 AIMD 和文献的文献值、Mulliken、EHT、Hirshfeld 电荷）进行光谱比较。
主要结果
- 采用不同RTIL电荷集的经验分子动力学模拟在DSC界面处产生不同的振动光谱和结合模式。
- 文献推导的电荷表现出与表面质子形成强烈的氢键作用，而AIMD推导的电荷则产生更真实的结合模式，并减少电荷屏蔽效应。
- 经验分子动力学模拟得到的振动光谱主峰频率因电荷集不同而异，其中AIMD光谱在300–400 cm-1区域显示出一组峰簇。
- 对经验力场进行仔细拟合并验证后，可合理预测DSC的振动特性。
结论
- 当经验势分子动力学经过适当参数化和验证后，可有效预测染料敏化太阳能电池界面的振动光谱。
- 室温离子液体的部分电荷参数化显著影响预测的界面性质和光谱。
- 将经验分子动力学与从头算分子动力学及实验验证相结合，可提高光活性界面建模的可靠性。
常见问题
行业重要性概述
准确预测染料/金属氧化物界面的振动光谱对于降低能源和光电子研发早期阶段材料的风险至关重要。本方案展示了在经过严格参数化和验证后，基于经验势的分子动力学方法如何以远低于从头计算方法的计算成本，提供对界面性质的预测性洞察。此类能力使得新型光活性材料及界面的快速假设检验和项目筛选成为可能。
生物制药研发中的战略应用
早期发现与靶点验证
- 能够研究与功能材料性能相关的界面分子相互作用及结合模式。
- 通过阐明电荷参数化对预测光谱和界面结构的影响，支持机制层面的风险评估。
- 有助于提升材料筛选的预测可靠性，并为器件原型开发的早期阶段提供决策依据。
筛选与检测方法开发
- 为染料-表面及电解质组合的下游筛选制备经过验证的计算模型。
- 标准化模拟输出结果，以实现与实验光谱的可重复性和交叉比较。
- 支持参数集的可扩展性评估，助力高通量计算筛选工作流程。
转化与临床前研究
- 将计算预测与实验和从头算（ab initio）基准对齐，支持从模拟到实验室验证的转化连续性。
- 基于定量光谱输出，为候选材料的风险调整型推进提供框架。
- 在进行高成本实验放大之前，降低界面设计中机制不确定性的风险。
流程与工作流整合
该模拟方案适用于从发现到临床前研究的连续过程，架起了计算假设检验与先进材料实验验证之间的桥梁。
- 发现生物学：在分子水平上支持对界面电荷效应和结合基序的假设验证。
- 筛选：提供可重复的定量振动光谱，用于参数集的比较评估。
- 分析：提供质量加权的光谱读出结果和速度自相关分析，以支持可靠的数据驱动决策。
- 转化研究：通过将实证预测与从头计算和实验数据进行基准比对，确保研究连续性。
- 企业级复用：建立可复用的模拟工作流程，可适配于其他染料、表面和电解质体系。
运营与企业影响
- 科学价值：提高界面设计中的预测置信度，降低机制上的模糊性。
- 操作价值：实现标准化、可扩展且可重复的计算工作流程。
- 战略价值：提升资本使用效率，支持材料研发管线中更优的立项/终止决策。
- 组合影响：促进基于风险调整的候选材料优先排序与推进。
实施注意事项
- 需要具备分子模拟、力场参数化和光谱分析方面的专业知识。
- 需要能够访问高性能计算基础设施以及经过验证的模拟软件（例如 DL_POLY）。
- 要求跨团队对参数集和数据格式进行严格的标准化。
- 针对新体系的适应需要谨慎分配电荷，并通过实验或从头计算数据进行验证。
- 局限性包括对力场质量的敏感性以及持续基准测试的需求。
