rdkit.Chem中的化学式与分子量计算

RDKit是一种开源化学信息学软件包，用于分子建模和化学信息学。在RDKit中，可以轻松地计算化学分子的化学式和分子量。下面是使用RDKit计算化学式和分子量的示例：

from rdkit import Chem

# 1. 计算化学式
# 输入一个SMILES字符串
smiles = 'CC(=O)O'
mol = Chem.MolFromSmiles(smiles)

# 使用当前的Mol对象计算分子的化学式
formula = Chem.rdMolDescriptors.CalcMolFormula(mol)
print(f"Molecular formula: {formula}")

# 2. 计算分子量
# 使用当前的Mol对象计算分子的分子量
mw = Chem.rdMolDescriptors.CalcExactMolWt(mol)
print(f"Molecular weight: {mw}")

# 3. 计算多个分子的化学式和分子量
# 输入一个SMILES字符串列表
smiles_list = ['CC(=O)O', 'C1=CC=CC=C1']
molecules = [Chem.MolFromSmiles(smiles) for smiles in smiles_list]

# 使用循环计算每个分子的化学式和分子量
for i, mol in enumerate(molecules):
    formula = Chem.rdMolDescriptors.CalcMolFormula(mol)
    mw = Chem.rdMolDescriptors.CalcExactMolWt(mol)
    print(f"Molecular formula for molecule {i+1}: {formula}")
    print(f"Molecular weight for molecule {i+1}: {mw}")
    print('-' * 30)

在上面的示例中，我们首先使用给定的SMILES字符串创建一个Mol对象。然后，我们使用rdMolDescriptors模块中的CalcMolFormula函数计算该分子的化学式，并使用CalcExactMolWt函数计算分子的分子量。如果需要计算多个分子的化学式和分子量，我们可以将多个SMILES字符串转换为Mol对象，然后通过循环对每个分子进行计算。

RDKit非常方便易用，可以用于各种化学信息学任务，包括计算化学式和分子量。上述示例提供了在RDKit中计算化学式和分子量的简单用法。