Closed entropybit closed 1 year ago
Hi @entropybit ,
The error reporting is not great at that stage so its not always clear what has happened, there should be within the optimisation folder a json file called result.json
which should have more information could you post that here it should help diagnose what's gone wrong.
Sure, I assume you meant workflow_result.json:
{
"version": "2.1.1",
"creation_date": "2023_03_07",
"modified_date": "2023_03_07",
"input_molecule": {
"atoms": [
{
"atomic_number": 6,
"atom_index": 0,
"atom_name": "C1",
"formal_charge": 0,
"aromatic": false,
"stereochemistry": null,
"bonds": [
1,
2,
3,
4
],
"aim": {
"volume": null,
"charge": null,
"c8": null
},
"dipole": null,
"quadrupole": null,
"cloud_pen": null
},
{
"atomic_number": 8,
"atom_index": 1,
"atom_name": "O1",
"formal_charge": 0,
"aromatic": false,
"stereochemistry": null,
"bonds": [
0,
5
],
"aim": {
"volume": null,
"charge": null,
"c8": null
},
"dipole": null,
"quadrupole": null,
"cloud_pen": null
},
{
"atomic_number": 1,
"atom_index": 2,
"atom_name": "H1",
"formal_charge": 0,
"aromatic": false,
"stereochemistry": null,
"bonds": [
0
],
"aim": {
"volume": null,
"charge": null,
"c8": null
},
"dipole": null,
"quadrupole": null,
"cloud_pen": null
},
{
"atomic_number": 1,
"atom_index": 3,
"atom_name": "H2",
"formal_charge": 0,
"aromatic": false,
"stereochemistry": null,
"bonds": [
0
],
"aim": {
"volume": null,
"charge": null,
"c8": null
},
"dipole": null,
"quadrupole": null,
"cloud_pen": null
},
{
"atomic_number": 1,
"atom_index": 4,
"atom_name": "H3",
"formal_charge": 0,
"aromatic": false,
"stereochemistry": null,
"bonds": [
0
],
"aim": {
"volume": null,
"charge": null,
"c8": null
},
"dipole": null,
"quadrupole": null,
"cloud_pen": null
},
{
"atomic_number": 1,
"atom_index": 5,
"atom_name": "H4",
"formal_charge": 0,
"aromatic": false,
"stereochemistry": null,
"bonds": [
1
],
"aim": {
"volume": null,
"charge": null,
"c8": null
},
"dipole": null,
"quadrupole": null,
"cloud_pen": null
}
],
"bonds": [
{
"atom1_index": 0,
"atom2_index": 1,
"bond_order": 1.0,
"aromatic": false,
"stereochemistry": null
},
{
"atom1_index": 0,
"atom2_index": 2,
"bond_order": 1.0,
"aromatic": false,
"stereochemistry": null
},
{
"atom1_index": 0,
"atom2_index": 3,
"bond_order": 1.0,
"aromatic": false,
"stereochemistry": null
},
{
"atom1_index": 0,
"atom2_index": 4,
"bond_order": 1.0,
"aromatic": false,
"stereochemistry": null
},
{
"atom1_index": 1,
"atom2_index": 5,
"bond_order": 1.0,
"aromatic": false,
"stereochemistry": null
}
],
"coordinates": [
-0.35075347726940603,
-0.005072983373261349,
-0.018028132595706226,
0.9643704239531457,
-0.36240238530810037,
-0.2601199990033363,
-0.5984566510608368,
0.061649323412126845,
1.0523734303936758,
-0.9862545465164104,
-0.8156752796030018,
-0.4626034355081824,
-0.6578568673402398,
0.9268983092448698,
-0.5010770555240878,
1.6289511182337482,
0.19460301562736884,
0.18945519223763826
],
"multiplicity": 1,
"name": "methanol",
"provenance": {
"creator": "QUBEKit",
"version": "2.1.1",
"routine": [
"QUBEKit.ligand.from_rdkit",
"QUBEKit.ligand.from_smiles",
"CO"
]
},
"extra_sites": {
"sites": {}
},
"BondForce": {
"type": "HarmonicBondForce",
"parameters": null
},
"AngleForce": {
"type": "HarmonicAngleForce",
"parameters": null
},
"TorsionForce": {
"type": "PeriodicTorsionForce",
"parameters": null,
"ordering": "default"
},
"ImproperTorsionForce": {
"type": "PeriodicImproperTorsionForce",
"parameters": null
},
"RBTorsionForce": {
"type": "RBProperTorsionForce",
"parameters": null
},
"ImproperRBTorsionForce": {
"type": "RBImproperTorsionForce",
"parameters": null
},
"NonbondedForce": {
"type": "NonbondedForce",
"parameters": null,
"coulomb14scale": 0.8333333333,
"lj14scale": 0.5,
"combination": "amber"
},
"chargemol_coords": null,
"hessian": null,
"qm_scans": null,
"wbo": null
},
"qc_spec": {
"program": "gaussian",
"basis": "6-311++G(d,p)",
"method": "wb97x-d",
"td_settings": null
},
"local_resources": {
"cores": 96,
"memory": 200
},
"current_molecule": {
"atoms": [
{
"atomic_number": 6,
"atom_index": 0,
"atom_name": "C1",
"formal_charge": 0,
"aromatic": false,
"stereochemistry": null,
"bonds": [
1,
2,
3,
4
],
"aim": {
"volume": null,
"charge": null,
"c8": null
},
"dipole": null,
"quadrupole": null,
"cloud_pen": null
},
{
"atomic_number": 8,
"atom_index": 1,
"atom_name": "O1",
"formal_charge": 0,
"aromatic": false,
"stereochemistry": null,
"bonds": [
0,
5
],
"aim": {
"volume": null,
"charge": null,
"c8": null
},
"dipole": null,
"quadrupole": null,
"cloud_pen": null
},
{
"atomic_number": 1,
"atom_index": 2,
"atom_name": "H1",
"formal_charge": 0,
"aromatic": false,
"stereochemistry": null,
"bonds": [
0
],
"aim": {
"volume": null,
"charge": null,
"c8": null
},
"dipole": null,
"quadrupole": null,
"cloud_pen": null
},
{
"atomic_number": 1,
"atom_index": 3,
"atom_name": "H2",
"formal_charge": 0,
"aromatic": false,
"stereochemistry": null,
"bonds": [
0
],
"aim": {
"volume": null,
"charge": null,
"c8": null
},
"dipole": null,
"quadrupole": null,
"cloud_pen": null
},
{
"atomic_number": 1,
"atom_index": 4,
"atom_name": "H3",
"formal_charge": 0,
"aromatic": false,
"stereochemistry": null,
"bonds": [
0
],
"aim": {
"volume": null,
"charge": null,
"c8": null
},
"dipole": null,
"quadrupole": null,
"cloud_pen": null
},
{
"atomic_number": 1,
"atom_index": 5,
"atom_name": "H4",
"formal_charge": 0,
"aromatic": false,
"stereochemistry": null,
"bonds": [
1
],
"aim": {
"volume": null,
"charge": null,
"c8": null
},
"dipole": null,
"quadrupole": null,
"cloud_pen": null
}
],
"bonds": [
{
"atom1_index": 0,
"atom2_index": 1,
"bond_order": 1.0,
"aromatic": false,
"stereochemistry": null
},
{
"atom1_index": 0,
"atom2_index": 2,
"bond_order": 1.0,
"aromatic": false,
"stereochemistry": null
},
{
"atom1_index": 0,
"atom2_index": 3,
"bond_order": 1.0,
"aromatic": false,
"stereochemistry": null
},
{
"atom1_index": 0,
"atom2_index": 4,
"bond_order": 1.0,
"aromatic": false,
"stereochemistry": null
},
{
"atom1_index": 1,
"atom2_index": 5,
"bond_order": 1.0,
"aromatic": false,
"stereochemistry": null
}
],
"coordinates": [
-0.35075347726940603,
-0.005072983373261349,
-0.018028132595706226,
0.9643704239531457,
-0.36240238530810037,
-0.2601199990033363,
-0.5984566510608368,
0.061649323412126845,
1.0523734303936758,
-0.9862545465164104,
-0.8156752796030018,
-0.4626034355081824,
-0.6578568673402398,
0.9268983092448698,
-0.5010770555240878,
1.6289511182337482,
0.19460301562736884,
0.18945519223763826
],
"multiplicity": 1,
"name": "methanol",
"provenance": {
"creator": "QUBEKit",
"version": "2.1.1",
"routine": [
"QUBEKit.ligand.from_rdkit",
"QUBEKit.ligand.from_smiles",
"CO"
]
},
"extra_sites": {
"sites": {}
},
"BondForce": {
"type": "HarmonicBondForce",
"parameters": [
{
"type": "HarmonicBondParameter",
"attributes": [],
"parameter_eval": [],
"atoms": [
0,
1
],
"length": 0.14293303934379997,
"k": 277982.70612434536
},
{
"type": "HarmonicBondParameter",
"attributes": [],
"parameter_eval": [],
"atoms": [
0,
2
],
"length": 0.10939105249969999,
"k": 315503.50521639985
},
{
"type": "HarmonicBondParameter",
"attributes": [],
"parameter_eval": [],
"atoms": [
0,
3
],
"length": 0.10939105249969999,
"k": 315503.50521639985
},
{
"type": "HarmonicBondParameter",
"attributes": [],
"parameter_eval": [],
"atoms": [
0,
4
],
"length": 0.10939105249969999,
"k": 315503.50521639985
},
{
"type": "HarmonicBondParameter",
"attributes": [],
"parameter_eval": [],
"atoms": [
1,
5
],
"length": 0.09713231822139,
"k": 464991.4883167736
}
]
},
"AngleForce": {
"type": "HarmonicAngleForce",
"parameters": [
{
"type": "HarmonicAngleParameter",
"attributes": [],
"parameter_eval": [],
"atoms": [
0,
1,
5
],
"k": 562.7562748394744,
"angle": 1.9249208206429511
},
{
"type": "HarmonicAngleParameter",
"attributes": [],
"parameter_eval": [],
"atoms": [
1,
0,
2
],
"k": 415.19503141569464,
"angle": 1.9836878140552876
},
{
"type": "HarmonicAngleParameter",
"attributes": [],
"parameter_eval": [],
"atoms": [
1,
0,
3
],
"k": 415.19503141569464,
"angle": 1.9836878140552876
},
{
"type": "HarmonicAngleParameter",
"attributes": [],
"parameter_eval": [],
"atoms": [
1,
0,
4
],
"k": 415.19503141569464,
"angle": 1.9836878140552876
},
{
"type": "HarmonicAngleParameter",
"attributes": [],
"parameter_eval": [],
"atoms": [
2,
0,
3
],
"k": 278.45479940981784,
"angle": 1.9948080932715693
},
{
"type": "HarmonicAngleParameter",
"attributes": [],
"parameter_eval": [],
"atoms": [
2,
0,
4
],
"k": 278.45479940981784,
"angle": 1.9948080932715693
},
{
"type": "HarmonicAngleParameter",
"attributes": [],
"parameter_eval": [],
"atoms": [
3,
0,
4
],
"k": 278.45479940981784,
"angle": 1.9948080932715693
}
]
},
"TorsionForce": {
"type": "PeriodicTorsionForce",
"parameters": [
{
"type": "PeriodicTorsionParameter",
"attributes": [],
"parameter_eval": [],
"atoms": [
2,
0,
1,
5
],
"k1": 0,
"k2": 0,
"k3": 0.9743050933247757,
"k4": 0,
"periodicity1": 1,
"periodicity2": 2,
"periodicity3": 3,
"periodicity4": 4,
"phase1": 0,
"phase2": 3.141592653589793,
"phase3": 0.0,
"phase4": 3.141592653589793
},
{
"type": "PeriodicTorsionParameter",
"attributes": [],
"parameter_eval": [],
"atoms": [
3,
0,
1,
5
],
"k1": 0,
"k2": 0,
"k3": 0.9743050933247757,
"k4": 0,
"periodicity1": 1,
"periodicity2": 2,
"periodicity3": 3,
"periodicity4": 4,
"phase1": 0,
"phase2": 3.141592653589793,
"phase3": 0.0,
"phase4": 3.141592653589793
},
{
"type": "PeriodicTorsionParameter",
"attributes": [],
"parameter_eval": [],
"atoms": [
4,
0,
1,
5
],
"k1": 0,
"k2": 0,
"k3": 0.9743050933247757,
"k4": 0,
"periodicity1": 1,
"periodicity2": 2,
"periodicity3": 3,
"periodicity4": 4,
"phase1": 0,
"phase2": 3.141592653589793,
"phase3": 0.0,
"phase4": 3.141592653589793
}
],
"ordering": "smirnoff"
},
"ImproperTorsionForce": {
"type": "PeriodicImproperTorsionForce",
"parameters": []
},
"RBTorsionForce": {
"type": "RBProperTorsionForce",
"parameters": []
},
"ImproperRBTorsionForce": {
"type": "RBImproperTorsionForce",
"parameters": []
},
"NonbondedForce": {
"type": "NonbondedForce",
"parameters": [
{
"type": "LennardJones612",
"attributes": [],
"parameter_eval": [],
"atoms": [
0
],
"charge": 0.11669999999999998,
"epsilon": 0.4577296,
"sigma": 0.33996695084235345
},
{
"type": "LennardJones612",
"attributes": [],
"parameter_eval": [],
"atoms": [
1
],
"charge": -0.5988,
"epsilon": 0.8803136,
"sigma": 0.30664733878390477
},
{
"type": "LennardJones612",
"attributes": [],
"parameter_eval": [],
"atoms": [
2
],
"charge": 0.02869999999999999,
"epsilon": 0.06568879999999999,
"sigma": 0.24713530441213008
},
{
"type": "LennardJones612",
"attributes": [],
"parameter_eval": [],
"atoms": [
3
],
"charge": 0.02869999999999999,
"epsilon": 0.06568879999999999,
"sigma": 0.24713530441213008
},
{
"type": "LennardJones612",
"attributes": [],
"parameter_eval": [],
"atoms": [
4
],
"charge": 0.02869999999999999,
"epsilon": 0.06568879999999999,
"sigma": 0.24713530441213008
},
{
"type": "LennardJones612",
"attributes": [],
"parameter_eval": [],
"atoms": [
5
],
"charge": 0.396,
"epsilon": 0.0002204968,
"sigma": 0.05345392308842035
}
],
"coulomb14scale": 0.8333333333,
"lj14scale": 0.5,
"combination": "amber"
},
"chargemol_coords": null,
"hessian": null,
"qm_scans": null,
"wbo": null
},
"results": {
"parametrisation": {
"stage": "OpenFF",
"stage_settings": {
"type": "OpenFF",
"force_field": "openff_unconstrained-1.3.0.offxml"
},
"status": "done",
"error": null
},
"optimisation": {
"stage": "Optimiser",
"stage_settings": {
"type": "Optimiser",
"pre_optimisation_method": "gfn2xtb",
"seed_conformers": 40,
"maxiter": 350,
"convergence_criteria": "GAU_TIGHT"
},
"status": "error",
"error": [
" File \"/work/scratch/b_mayer/miniconda3/envs/qubekit/lib/python3.10/site-packages/qubekit/workflow/workflow.py\", line 335, in _run_stage\n result_mol = stage.run(\n",
" File \"/work/scratch/b_mayer/miniconda3/envs/qubekit/lib/python3.10/site-packages/qubekit/workflow/helper_stages.py\", line 144, in run\n return self._run_qm_opt(\n",
" File \"/work/scratch/b_mayer/miniconda3/envs/qubekit/lib/python3.10/site-packages/qubekit/workflow/helper_stages.py\", line 206, in _run_qm_opt\n raise GeometryOptimisationError(\n"
]
},
"hessian": {
"stage": "Hessian",
"stage_settings": {
"type": "Hessian"
},
"status": "waiting",
"error": null
},
"charges": {
"stage": "DDECCharges",
"stage_settings": {
"type": "DDECCharges",
"solvent_settings": {
"program": "gaussian",
"solver_type": "IPCM",
"epsilon": 4.0,
"volume_contour": 0.0004
},
"program": "gaussian",
"basis": null,
"method": null,
"td_settings": null,
"ddec_version": 6
},
"status": "waiting",
"error": null
},
"virtual_sites": {
"stage": null,
"stage_settings": null,
"status": "skip",
"error": null
},
"non_bonded": {
"stage": "LennardJones612",
"stage_settings": {
"type": "LennardJones612",
"lj_on_polar_h": true,
"free_parameters": {
"H": {
"v_free": 7.6,
"b_free": 6.5,
"r_free": 1.738
},
"X": {
"v_free": 7.6,
"b_free": 6.5,
"r_free": 1.083
},
"C": {
"v_free": 34.4,
"b_free": 46.6,
"r_free": 2.008
},
"N": {
"v_free": 25.9,
"b_free": 24.2,
"r_free": 1.765
},
"O": {
"v_free": 22.1,
"b_free": 15.6,
"r_free": 1.499
}
},
"alpha": 1.0,
"beta": 0.0
},
"status": "waiting",
"error": null
},
"bonded_parameters": {
"stage": "ModSeminario",
"stage_settings": {
"type": "ModSeminario",
"vibrational_scaling": 1.0
},
"status": "waiting",
"error": null
},
"torsion_scanner": {
"stage": "TorsionScan1D",
"stage_settings": {
"type": "TorsionScan1D",
"special_torsions": [],
"default_scan_range": [
-165,
180
],
"avoided_torsions": [
{
"smirks": "[*:1]-[CH3:2]"
},
{
"smirks": "[*:1]-[NH2:2]"
}
],
"torsion_driver": {
"type": "torsiondriver",
"n_workers": 1,
"grid_spacing": 15,
"energy_decrease_thresh": null,
"energy_upper_limit": null,
"starting_conformations": 4
}
},
"status": "waiting",
"error": null
},
"torsion_optimisation": {
"stage": "ForceBalanceFitting",
"stage_settings": {
"type": "ForceBalanceFitting",
"penalty_type": "L1",
"job_type": "optimize",
"max_iterations": 10,
"convergence_step_criteria": 0.01,
"convergence_objective_criteria": 0.01,
"convergence_gradient_criteria": 0.01,
"n_criteria": 1,
"eig_lowerbound": 0.01,
"finite_difference_h": 0.01,
"penalty_additive": 0.1,
"constrain_charge": false,
"initial_trust_radius": -0.25,
"minimum_trust_radius": 0.05,
"error_tolerance": 1.0,
"adaptive_factor": 0.2,
"adaptive_damping": 1.0,
"normalize_weights": false,
"extras": {},
"priors": {
"Proper_k": 6.0
},
"targets": {
"TorsionProfile_OpenMM": {
"target_name": "TorsionProfile_OpenMM",
"description": "Relaxed energy and RMSD fitting for torsion drives only.",
"writelevel": 2,
"keywords": {
"pdb": "molecule.pdb",
"coords": "scan.xyz"
},
"openmm_platform": "Reference",
"energy_denom": 1.0,
"energy_upper": 10.0,
"attenuate": false,
"restrain_k": 1.0
}
}
},
"status": "waiting",
"error": null
}
}
}
No there should be an output from qcengine in a folder called optimisation
, the json file should contain the raw output from the QM software.
Okay I have simply appended a zip of the entire results folder here:
Thanks very much! It looks like our most recent version on conda-forge (2.1.1) is missing a fix to work with Gaussian16. Ill look into updating this but for now, it should work if you install the main branch from github. These commands should help
Remove qubekit only from conda
conda uninstall qubekit --force
clone github code
git clone https://github.com/qubekit/QUBEKit.git
move into package
cd QUBEKit
install package
pip install .
Thanks a lot, I did all of this and just submitted a job with the fresh install. Let's see if it works this time ^^. If so I'll close it afterwards.
I think you should really update your requirements.txt. openmm should be in it I guess, and after doing what you just proposed it first failed with module not found chemper ^^'.
To be continued:
qubekit.utils.exceptions.WorkFlowExecutionError: The workflow stopped unexpectedly
due to the following error at stage: torsion_optimisation
Thanks for testing this out, it looks like our new fragmentation changes caused some issues, the main branch should now work provided you use a protocol that does not use virtual sites. If you do need v-sites you will need a specific branch of the openff-toolkit, the 10-x branch which you can install via pip or from github for now.
Hmh, I did exactly that and after mistakenly uninstalling and installing qubekit again I had to do the following things
pip install click
pip install qcengine
pip install torsiondrive
pip install forcebalance
conda install ndcctools -c conda-forge
pip install xtb
XTB i installed after this
(qubekit) [b_mayer@logc0002 fabi_ligands]$ qubekit run -sm CO -n methanol
Warning: Unable to load toolkit 'OpenEye Toolkit'. The Open Force Field Toolkit does not require the OpenEye Toolkits, and can use RDKit/AmberTools instead. However, if you have a valid license for the OpenEye Toolkits, consider installing them for faster performance and additional file format support: https://docs.eyesopen.com/toolkits/python/quickstart-python/linuxosx.html OpenEye offers free Toolkit licenses for academics: https://www.eyesopen.com/academic-licensing
Traceback (most recent call last):
File "/work/scratch/b_mayer/miniconda3/envs/qubekit/bin/qubekit", line 8, in <module>
sys.exit(cli())
File "/work/scratch/b_mayer/miniconda3/envs/qubekit/lib/python3.10/site-packages/click/core.py", line 1130, in __call__
return self.main(*args, **kwargs)
File "/work/scratch/b_mayer/miniconda3/envs/qubekit/lib/python3.10/site-packages/click/core.py", line 1055, in main
rv = self.invoke(ctx)
File "/work/scratch/b_mayer/miniconda3/envs/qubekit/lib/python3.10/site-packages/click/core.py", line 1657, in invoke
return _process_result(sub_ctx.command.invoke(sub_ctx))
File "/work/scratch/b_mayer/miniconda3/envs/qubekit/lib/python3.10/site-packages/click/core.py", line 1404, in invoke
return ctx.invoke(self.callback, **ctx.params)
File "/work/scratch/b_mayer/miniconda3/envs/qubekit/lib/python3.10/site-packages/click/core.py", line 760, in invoke
return __callback(*args, **kwargs)
File "/work/scratch/b_mayer/miniconda3/envs/qubekit/lib/python3.10/site-packages/qubekit/cli/run.py", line 136, in run
workflow.new_workflow(molecule=molecule, skip_stages=skip_stages, end=end)
File "/work/scratch/b_mayer/miniconda3/envs/qubekit/lib/python3.10/site-packages/qubekit/workflow/workflow.py", line 272, in new_workflow
return self._run_workflow(
File "/work/scratch/b_mayer/miniconda3/envs/qubekit/lib/python3.10/site-packages/qubekit/workflow/workflow.py", line 294, in _run_workflow
self.validate_workflow(workflow=workflow, molecule=molecule)
File "/work/scratch/b_mayer/miniconda3/envs/qubekit/lib/python3.10/site-packages/qubekit/workflow/workflow.py", line 122, in validate_workflow
pre_spec.validate_specification()
File "/work/scratch/b_mayer/miniconda3/envs/qubekit/lib/python3.10/site-packages/qubekit/utils/datastructures.py", line 151, in validate_specification
self.validate_program()
File "/work/scratch/b_mayer/miniconda3/envs/qubekit/lib/python3.10/site-packages/qubekit/utils/datastructures.py", line 100, in validate_program
raise SpecificationError(
qubekit.utils.exceptions.SpecificationError: The program xtb is not available, available programs are {'gaussian', 'rdkit'}
Now it seems to me it is not using gaussian/g16 anymore. But with xtb it also doesn't work ^^
(qubekit) [b_mayer@logc0002 fabi_ligands]$ qubekit run -sm CO -n methanol
Warning: Unable to load toolkit 'OpenEye Toolkit'. The Open Force Field Toolkit does not require the OpenEye Toolkits, and can use RDKit/AmberTools instead. However, if you have a valid license for the OpenEye Toolkits, consider installing them for faster performance and additional file format support: https://docs.eyesopen.com/toolkits/python/quickstart-python/linuxosx.html OpenEye offers free Toolkit licenses for academics: https://www.eyesopen.com/academic-licensing
If QUBEKit ever breaks or you would like to view timings and loads of other info, view the log file.
Our documentation (README.md) also contains help on handling the various commands for QUBEKit.
Fragmenting molecule using the WBOFragmenter
Traceback (most recent call last):
File "/work/scratch/b_mayer/miniconda3/envs/qubekit/lib/python3.10/site-packages/qubekit/workflow/workflow.py", line 348, in _run_stage
result_mol = stage.run(
File "/work/scratch/b_mayer/miniconda3/envs/qubekit/lib/python3.10/site-packages/qubekit/fragmentation/wbo_fragmenter.py", line 71, in run
return self._run(molecule, *args, **kwargs)
File "/work/scratch/b_mayer/miniconda3/envs/qubekit/lib/python3.10/site-packages/qubekit/fragmentation/wbo_fragmenter.py", line 79, in _run
fragmentation_result = self.fragment(
File "/work/scratch/b_mayer/miniconda3/envs/qubekit/lib/python3.10/site-packages/openff/fragmenter/fragment.py", line 916, in fragment
result = self._fragment(molecule, target_bond_smarts)
File "/work/scratch/b_mayer/miniconda3/envs/qubekit/lib/python3.10/site-packages/openff/fragmenter/fragment.py", line 1012, in _fragment
molecule = assign_elf10_am1_bond_orders(
File "/work/scratch/b_mayer/miniconda3/envs/qubekit/lib/python3.10/site-packages/openff/fragmenter/chemi.py", line 41, in assign_elf10_am1_bond_orders
molecule = _generate_conformers(molecule, max_confs, rms_threshold)
File "/work/scratch/b_mayer/miniconda3/envs/qubekit/lib/python3.10/site-packages/openff/fragmenter/chemi.py", line 98, in _generate_conformers
canonical_molecule.generate_conformers(
File "/work/scratch/b_mayer/miniconda3/envs/qubekit/lib/python3.10/site-packages/openff/toolkit/topology/molecule.py", line 3030, in generate_conformers
return toolkit_registry.call(
File "/work/scratch/b_mayer/miniconda3/envs/qubekit/lib/python3.10/site-packages/openff/toolkit/utils/toolkit_registry.py", line 380, in call
raise ValueError(msg)
ValueError: No registered toolkits can provide the capability "generate_conformers" for args "(Molecule with name '' and SMILES '[H][O][C]([H])([H])[H]',)" and kwargs "{'n_conformers': 800, 'rms_cutoff': <Quantity(1.0, 'angstrom')>, 'clear_existing': True, 'make_carboxylic_acids_cis': True}"
Available toolkits are: [ToolkitWrapper around The RDKit version 2022.09.1, ToolkitWrapper around AmberTools version 22.0, ToolkitWrapper around Built-in Toolkit version None]
ToolkitWrapper around The RDKit version 2022.09.1 <class 'AttributeError'> : Neither Quantity object nor its magnitude (1.0) has attribute 'value_in_unit'
The above exception was the direct cause of the following exception:
Traceback (most recent call last):
File "/work/scratch/b_mayer/miniconda3/envs/qubekit/bin/qubekit", line 8, in <module>
sys.exit(cli())
File "/work/scratch/b_mayer/miniconda3/envs/qubekit/lib/python3.10/site-packages/click/core.py", line 1130, in __call__
return self.main(*args, **kwargs)
File "/work/scratch/b_mayer/miniconda3/envs/qubekit/lib/python3.10/site-packages/click/core.py", line 1055, in main
rv = self.invoke(ctx)
File "/work/scratch/b_mayer/miniconda3/envs/qubekit/lib/python3.10/site-packages/click/core.py", line 1657, in invoke
return _process_result(sub_ctx.command.invoke(sub_ctx))
File "/work/scratch/b_mayer/miniconda3/envs/qubekit/lib/python3.10/site-packages/click/core.py", line 1404, in invoke
return ctx.invoke(self.callback, **ctx.params)
File "/work/scratch/b_mayer/miniconda3/envs/qubekit/lib/python3.10/site-packages/click/core.py", line 760, in invoke
return __callback(*args, **kwargs)
File "/work/scratch/b_mayer/miniconda3/envs/qubekit/lib/python3.10/site-packages/qubekit/cli/run.py", line 136, in run
workflow.new_workflow(molecule=molecule, skip_stages=skip_stages, end=end)
File "/work/scratch/b_mayer/miniconda3/envs/qubekit/lib/python3.10/site-packages/qubekit/workflow/workflow.py", line 272, in new_workflow
return self._run_workflow(
File "/work/scratch/b_mayer/miniconda3/envs/qubekit/lib/python3.10/site-packages/qubekit/workflow/workflow.py", line 310, in _run_workflow
molecule = self._run_stage(
File "/work/scratch/b_mayer/miniconda3/envs/qubekit/lib/python3.10/site-packages/qubekit/workflow/workflow.py", line 372, in _run_stage
raise WorkFlowExecutionError(
qubekit.utils.exceptions.WorkFlowExecutionError: The workflow stopped unexpectedly due to the following error at stage: fragmentation
I have appended you an zip of the entire generated output folder again. methanol3.zip
Hi @entropybit sorry for the delay. For the first part, there is a pre-optimisation stage in the workflow which uses a cheaper method to hopefully speed up the full QM optimisation. By default, this runs xtb-GFN2 but can be set to something else like ANI2-x. For the second part it looks like when installing manually you picked up a newer version of the openff-fragmenter
for now you will need openff-toolkit=0.10.7
and openff-fragmenter=0.1.2
we should have a new conda package soon which makes this easier.
Holy crap
the first time it worked, thanks a lot. And don't worry about the delay, although I have to say I was very near to write something like "bump" here. I will see if I produce new errors with my other molecules and if not I will close this. Also I'll probably never touch that conda environment again until my research is done and possibly make a few copies of it ^^'.
Fantastic glad it finally works thanks for sticking with it and helping me get these issues sorted out, feel free to tag me in any other issues or help so I pick it up faster!
Hmh, well I got a similiar result for my molecule that doesn't have S in it as I wrote in the seconde issue. This:
Warning (not error because allow_undefined_stereo=True): Unable to make OFFMol from RDMol: RDMol has unspecified stereochemistry. Undefined chiral centers are:
- Atom C (index 7)
Stereochemistry for atom 37 flipped from R to None
Stereochemistry for atom 37 flipped from R to S
If QUBEKit ever breaks or you would like to view timings and loads of other info, view the log file.
Our documentation (README.md) also contains help on handling the various commands for QUBEKit.
Fragmenting molecule using the WBOFragmenter
Traceback (most recent call last):
File "/work/scratch/b_mayer/miniconda3/envs/qubekit/lib/python3.10/site-packages/qubekit/workflow/workflow.py", line 348, in _run_stage
result_mol = stage.run(
File "/work/scratch/b_mayer/miniconda3/envs/qubekit/lib/python3.10/site-packages/qubekit/fragmentation/wbo_fragmenter.py", line 71, in run
return self._run(molecule, *args, **kwargs)
File "/work/scratch/b_mayer/miniconda3/envs/qubekit/lib/python3.10/site-packages/qubekit/fragmentation/wbo_fragmenter.py", line 83, in _run
depict_fragmentation_result(
File "/work/scratch/b_mayer/miniconda3/envs/qubekit/lib/python3.10/site-packages/openff/fragmenter/depiction.py", line 307, in depict_fragmentation_result
depict_fragments(
File "/work/scratch/b_mayer/miniconda3/envs/qubekit/lib/python3.10/site-packages/openff/fragmenter/depiction.py", line 332, in depict_fragments
header_svg = _oe_render_parent(parent, [*fragments])
File "/work/scratch/b_mayer/miniconda3/envs/qubekit/lib/python3.10/site-packages/openff/fragmenter/depiction.py", line 92, in _oe_render_parent
oe_parent = parent.to_openeye()
File "/work/scratch/b_mayer/miniconda3/envs/qubekit/lib/python3.10/site-packages/openff/toolkit/utils/base_wrapper.py", line 50, in wrapped_function
raise ToolkitUnavailableException(msg)
openff.toolkit.utils.exceptions.ToolkitUnavailableException: This function requires the OpenEye Toolkit toolkit
The above exception was the direct cause of the following exception:
Traceback (most recent call last):
File "/work/scratch/b_mayer/miniconda3/envs/qubekit/bin/qubekit", line 8, in <module>
sys.exit(cli())
File "/work/scratch/b_mayer/miniconda3/envs/qubekit/lib/python3.10/site-packages/click/core.py", line 1130, in __call__
return self.main(*args, **kwargs)
File "/work/scratch/b_mayer/miniconda3/envs/qubekit/lib/python3.10/site-packages/click/core.py", line 1055, in main
rv = self.invoke(ctx)
File "/work/scratch/b_mayer/miniconda3/envs/qubekit/lib/python3.10/site-packages/click/core.py", line 1657, in invoke
return _process_result(sub_ctx.command.invoke(sub_ctx))
File "/work/scratch/b_mayer/miniconda3/envs/qubekit/lib/python3.10/site-packages/click/core.py", line 1404, in invoke
return ctx.invoke(self.callback, **ctx.params)
File "/work/scratch/b_mayer/miniconda3/envs/qubekit/lib/python3.10/site-packages/click/core.py", line 760, in invoke
return __callback(*args, **kwargs)
File "/work/scratch/b_mayer/miniconda3/envs/qubekit/lib/python3.10/site-packages/qubekit/cli/run.py", line 136, in run
workflow.new_workflow(molecule=molecule, skip_stages=skip_stages, end=end)
File "/work/scratch/b_mayer/miniconda3/envs/qubekit/lib/python3.10/site-packages/qubekit/workflow/workflow.py", line 272, in new_workflow
return self._run_workflow(
File "/work/scratch/b_mayer/miniconda3/envs/qubekit/lib/python3.10/site-packages/qubekit/workflow/workflow.py", line 310, in _run_workflow
molecule = self._run_stage(
File "/work/scratch/b_mayer/miniconda3/envs/qubekit/lib/python3.10/site-packages/qubekit/workflow/workflow.py", line 372, in _run_stage
raise WorkFlowExecutionError(
qubekit.utils.exceptions.WorkFlowExecutionError: The workflow stopped unexpectedly due to the following error at stage: fragmentation
But I guess that doesn't really have to do with the original problem here and is more of a me using the method wrong thing ? (Help required anyways but I guess I will close this here with this comment and wait for your reply in the other issue).
Sorry for opening yet another issue, I would not (yet) have done this if this was only with one molecule but I even get the following error
while running one of the examples:
and as mentioned I also got this with one other molecule. So I have the feeling that something is wrong and this will always be triggered. Hope you have any ideas on how to solve this, I would really like to be able to use the parameters from your very nice toolkit in my MD simulations.