如何捕获RDKit中Chem.MolFromSmiles的错误信息？

问题原因 #

RDKit的Chem.MolFromSmiles()函数在解析失败时不会主动抛出Python异常，只会返回None，同时将错误信息输出到标准错误流（stderr）。这就是你添加try-except后无法捕获错误的核心原因——函数根本没触发异常。

解决方法 #

下面提供两种可靠的捕获错误信息的方式：

方法1：利用RDKit内置的错误信息获取函数

直接检查MolFromSmiles的返回值，若为None，调用rdBase.GetLastError()获取详细错误：

from rdkit import Chem
from rdkit.rdBase import GetLastError

# 尝试解析SMILES
smiles = 'OCC1OC(C(C(C1O)O)O)[C]1(C)(CO)CC(=O)C=C(C1CCC(=O)C)C'
mol = Chem.MolFromSmiles(smiles)

if mol is None:
    # 获取并打印错误信息
    error_msg = GetLastError()
    print(f"解析失败：{error_msg}")
else:
    # 解析成功，处理结构（Jupyter环境中直接显示）
    display(mol)

方法2：将RDKit日志重定向到Python logging模块

适合需要统一日志管理的场景，把RDKit的错误日志转发到Python的logging系统：

import logging
from rdkit import Chem
from rdkit.rdBase import LogToPythonLogger

# 配置Python日志，只输出ERROR级别信息
logging.basicConfig(level=logging.ERROR, format='%(message)s')
# 把RDKit的日志重定向到Python logger
LogToPythonLogger()

smiles = 'OCC1OC(C(C(C1O)O)O)[C]1(C)(CO)CC(=O)C=C(C1CCC(=O)C)C'
mol = Chem.MolFromSmiles(smiles)

if mol is None:
    print("无法解析SMILES字符串，详细错误见上方日志")
else:
    display(mol)

额外提示：你的SMILES本身存在语法错误 #

报错提示的"Explicit valence for atom # 11 C, 5, is greater than permitted"，是因为SMILES里的[C]1(C)(CO)...部分，该碳原子的成键数达到了5，超出了碳的最大允许价（4）。如果是想表示碳正离子，应该写成[C+]（此时成键数为3，符合规则），修正后的SMILES如下：

OCC1OC(C(C(C1O)O)O)[C+]1(C)(CO)CC(=O)C=C(C1CCC(=O)C)C

内容的提问来源于stack exchange，提问作者user10913305