化学学报 上一篇    下一篇

研究论文

基于多尺度的专业大语言模型预测有机化合物键解离能

刘嘉乐, 周佳*   

  1. 哈尔滨工业大学(深圳)理学院,城乡水资源与水环境全国重点实验室,深圳 518000
  • 投稿日期:2026-06-30
  • 通讯作者: *E-mail: jiazhou@hit.edu.cn ; Tel.: +86 15704617427
  • 基金资助:
    城乡水资源与水环境全国重点实验室(哈尔滨工业大学)自主课题(课题编号:2025TS36).

Specialized Large Language Models with Multiscale Fusion for Predicting Bond Dissociation Energies of Organic Compounds

Jiale Liu, Jia Zhou*   

  1. State Key Lab of Urban-rural Water Resource and Environment, School of Science, Harbin Institute of Technology, Shenzhen 518055, China
  • Received:2026-06-30

键解离能(Bond Dissociation Energy, BDE)是表征化学键强度的基础热力学参数,决定有机分子的反应活性与热力学行为。传统实验测量与高精度量子化学计算耗时费力、资源开销大,难以满足大规模分子筛选的需求。尽管数据驱动的机器学习方法已成为高效替代方案,但现有方法多依赖分子图或分子指纹等单一尺度表征,忽略了键类型、官能团环境与碎片结构等自然语言描述中蕴含的化学语义信息。本研究提出一种整合分子结构序列与化学自然语言文本的多尺度融合框架,用于预测有机化合物的键解离能。框架包含三个模块:预训练分子编码器(ChemBERTa或MolFormer)从SMILES序列提取结构特征,大语言模型(DeepSeek-7B、Qwen2-7B或Mistral-7B)编码化学文本语义,交叉注意力融合层实现两者的自适应交互。采用低秩适配(LoRA)微调,仅更新约0.1%的参数。在包含531244条BDE记录、65540个有机分子的数据集上系统评估了6种模型组合,最优的ChemBERTa+Qwen2-7B组合在验证集上MAE为0.98 kcal mol-1、决定系数为0.9830。在数据发布方公开划分下重训并评测,官方测试集MAE为0.98 kcal mol-1,仍高于同划分下以键为中心的图神经网络报告值(0.61-0.64 kcal mol-1),差距集中于含杂原子的极性键。规模对照(Qwen2-0.5B/1.5B/7B:1.15/1.04/0.98)与输出头对照(纯文本配置下回归头1.39、生成式1.22 kcal mol-1)表明,分子编码器的融合而非输出头形式决定精度。结果表明,多尺度融合有效利用分子拓扑与化学语义的互补优势,为在不引入三维构型与量子化学特征的前提下开展大规模BDE预测提供了一条通用技术路线。

关键词: 键解离能, 大语言模型, 多尺度融合, 分子编码器, 低秩适配, 交叉注意力

Bond dissociation energy (BDE) is a fundamental thermodynamic parameter that characterizes chemical bond strength and governs the reactivity and thermal behavior of organic molecules. Traditional experimental measurements and high-level quantum chemical calculations are time-consuming and resource-intensive, thus impractical for large-scale molecular screening. Although data-driven machine learning methods have emerged as efficient alternatives, existing approaches rely predominantly on single-modal representations such as molecular graphs or fingerprints, overlooking the chemical semantics embedded in natural language descriptions of bond types, functional group environments, and fragment structures. Here we propose a multimodal fusion framework that integrates molecular structure sequences with chemical natural language text for BDE prediction. The framework comprises three modules: a pre-trained molecular encoder (ChemBERTa or MolFormer) extracting structural features from SMILES, a large language model (DeepSeek-7B, Qwen2-7B or Mistral-7B) encoding chemical text semantics, and a cross-attention-based fusion layer that enables adaptive interaction between them. Low-Rank Adaptation (LoRA) fine-tuning updated only approximately 0.1% of the parameters. Six model combinations were evaluated on a dataset of 531244 BDE records from 65540 organic molecules; the optimal ChemBERTa+Qwen2-7B combination reached a mean absolute error (MAE) of 0.98 kcal mol-1 and a coefficient of determination (R2) of 0.9830 on the validation set. On the publicly released split, the model attains an MAE of 0.98 kcal mol-1 on the official test set, still above the 0.61-0.64 kcal mol-1 reported for bond-centred graph neural networks, with the gap concentrated on heteroatom-containing polar bonds. Scale ablation (Qwen2-0.5B/1.5B/7B: 1.15/1.04/0.98 kcal mol-1) and an output-head ablation (regression head 1.39 vs. generative 1.22 kcal mol-1) show that the fusion of molecular representations, not the output-head form, governs accuracy. The framework thus exploits the complementary strengths of molecular topology and chemical semantics, offering a general route to large-scale BDE prediction without three-dimensional geometries or quantum-chemical features.

Key words: Bond dissociation energy, Large language model, Multiscale fusion, Molecular encoder, Low-rank adaptation, Cross-attention