各位老师，我想请教一下怎么正确地生成盐类分子的三维坐标，我用RDkit和pybel生成的原子重叠，阴离子构型也不合理（以[Li+].F[P-](F)(F)(F)(F)F为例）：

收藏

libo371324 发表于 2022-12-29 10:40
我想问一下这个模块目前默认的产生坐标的设置是useExpTorsionAnglePrefs=True, useBasicKnowledge=True， ...

抱歉啊，刚看见，不知道你是否已经解决。
按照我之前的解决方案，应该将你的代码改为
“AllChem.EmbedMolecule(a, useRandomCoords=True)”，应该可以避免原子重叠。
至于openbabel的pybel的安装依照“冰释之川”的回答即可。

libo371324 发表于 2022-12-30 09:30
我昨天弄了一天的openbabel包，但是一直不成功，from openbabel import pybel，一直导入失败，到现在也没 ...

直接conda装就行了 conda install openbabel

冰释之川 发表于 2022-12-29 21:13
换pybel 试试

我昨天弄了一天的openbabel包，但是一直不成功，from openbabel import pybel，一直导入失败，到现在也没有下载成功，我按照网站上下载安装，也是不成功，能给一下详细的下载教程吗

libo371324 发表于 2022-12-29 10:40
我想问一下这个模块目前默认的产生坐标的设置是useExpTorsionAnglePrefs=True, useBasicKnowledge=True， ...

换pybel 试试
RDKit_3D_mols_gen.py (2.7 KB, 下载次数 Times of downloads: 61)

2022-12-29 21:13 上传 Uploaded

点击下载
Click to download

wangxubin 发表于 2021-11-9 21:51
谢谢社长，我试过直接用openbabel产生结构，对于复杂分子，存在生成极不合理的分子结构的情况，比如多个 ...

我想问一下这个模块目前默认的产生坐标的设置是useExpTorsionAnglePrefs=True, useBasicKnowledge=True，但是我在使用的过程中出现了那种极不合理的分子结构，所有的原子都重合了，即a= AllChem.EmbedMolecule(),等于-1的那种情况，我师弟给我说是生成失败了，这种情况应该怎么办，换用open babel这个软件包吗，还是说rdkit有其他的构象生成的方法。

sobereva 发表于 2021-11-8 22:31
把openbabel搞出来的结构，以及进一步用openbabel做MMFF94力场优化后的结构都截个图。如果后者就已经足够合 ...

谢谢社长，我试过直接用openbabel产生结构，对于复杂分子，存在生成极不合理的分子结构的情况，比如多个原子重叠。
目前已通过下述方法实现生成、优化分子结构。
AllChem.EmbedMolecule(mol, useRandomCoords=True)
AllChem.MMFFOptimizeMolecule(mol)

把openbabel搞出来的结构，以及进一步用openbabel做MMFF94力场优化后的结构都截个图。如果后者就已经足够合理，就用它就完了

wangxubin 发表于 2021-11-8 11:05
都有xyz文件了，直接用openbabel转换更方便一些。
obabel -ixyz test.xyz -omol -O test.mol

嗯呢 这不是想着少用一个工具算一个么

心向暖阳 发表于 2021-11-8 09:36
插楼问一下哈，xyz转成mol有自带的函数么

都有xyz文件了，直接用openbabel转换更方便一些。
obabel -ixyz test.xyz -omol -O test.mol

插楼问一下哈，xyz转成mol有自带的函数么