參考：https://zhuanlan.zhihu.com/p/87350161

芳香性定義參考： 休克爾規則 ：4n+2 https://baike.baidu.com/item/%E4%BC%91%E5%85%8B%E5%B0%94%E8%A7%84%E5%88%99/7496254?fr=aladdin

計算環的數量，所有類型的環都會統計

m=Chem.MolFromSmiles("*Nc1nccnc1") ri = m.GetRingInfo() print(ri.NumRings())

芳香環的數量

def Calc_AROM(mh):m = Chem.RemoveHs(mh)ring_info = m.GetRingInfo()atoms_in_rings = ring_info.AtomRings()num_aromatic_ring = 0for ring in atoms_in_rings:aromatic_atom_in_ring = 0for atom_id in ring:atom = m.GetAtomWithIdx(atom_id)if atom.GetIsAromatic():aromatic_atom_in_ring += 1if aromatic_atom_in_ring == len(ring):num_aromatic_ring += 1return num_aromatic_ringCalc_AROM(m)

總結

以上是生活随笔為你收集整理的rdkit 计算环、芳香环数的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。