Open eva-not opened 1 year ago
Hi, sorry for the delay. This looks pretty tricky - thanks for all the info, it really helps track stuff down.
My first idea is to capture errors in the function call to QCEngine and re-raise the error with extra info about the fragment. That wouldn't fix it, but it would provide more information about the specific fragment that's causing the failure. Since geomeTRIC is already providing a decent amount of information about the failure, I don't think it's reasonable to ask it to provide more info ("and here is the mapped SMILES of the full molecule:" etc.)
A more robust solution would probably be https://github.com/openforcefield/openff-bespokefit/issues/295; I figure you shouldn't normally have linear triplets around a double bond, so this is just a funky chemistry forced by driving the torsion, so skipping that one particular structure is justified.
We recently released a new version (0.2.3) with a couple of improvements related to torsion drives (parallelization #277 and re-starting failed jobs #279) but if I understand the issue correctly, neither of these really address the root cause. It still might be worth a shot.
Hi, thank you very much for the reply. Unfortunately I can't skip this particular structure as it's one of the few I would like to optimise. I will probably have to resort to manual torsion scans.
Hi, I am trying to get bespoke torsions for the molecule
CNC(=O)[C@H](CCC\C=C/CCC[C@H](NC(C)=O)C(=O)NC)NC(C)=O
. However, the QC generation step exits with the following message:Is there a way to see which atoms these are and which fragment is causing the error? I suspect it's atoms around the C=C double bond, since I ran bespokefit on the same molecule with a single bond instead
CNC([C@@H](NC(C)=O)CCCCCCCC[C@@H](C(NC)=O)NC(C)=O)=O
and the calculations finished fine (I got a .offxml file).I also tested bespokefit on a smaller molecule with a C=C double bond
O=C(O)/C=C/C(O)=O
and I got a similar error:Bespokefit is ran via the CLI:
The workflow file is: factory.json
Any help with this would be greatly appreciated!
Software versions
conda list
?Output of
conda list