Comments (6)
@eneas77 Hi, could you try the v1.4.2 development branch, and see if the problem still persists?
from autode.
Hi, I suspect this is the same issue as #261 , however, now autodE detects the problem much before it happens and then throws the appropriate error. I am currently working on it, but it will be some time before the coordinate problems get completely fixed.
from autode.
Hi,
Did you try to run the DA example when you released v1.4.0 & v1.4.1?
I suppose the answer is you didn't.
Perhaps you should include that basic QC: every case ending in bug/enhancement should be tested previously to the deployment, as well as the basic QC examples, like the DA. This is the second time it happens...
Two more shilings on this matter.
Thanks a lot.
E77
from autode.
Hi @shoubhikraj,
I´ve tried but unfrotunately some issue is triggered...
--------------------
Traceback (most recent call last):
File "/home/data/amiduros/142/test_thf-78C/test_thf_-78C.py", line 11, in
rxn.calculate_reaction_profile(free_energy=True)
File "/apps/anaconda3/envs/autode_142/lib/python3.10/site-packages/autode/utils.py", line 433, in wrapped_function
return func(*args, **kwargs)
File "/apps/anaconda3/envs/autode_142/lib/python3.10/site-packages/autode/reactions/reaction.py", line 168, in calculate_reaction_profile
calculate(self)
File "/apps/anaconda3/envs/autode_142/lib/python3.10/site-packages/autode/utils.py", line 239, in wrapped_function
result = func(*args, **kwargs)
File "/apps/anaconda3/envs/autode_142/lib/python3.10/site-packages/autode/reactions/reaction.py", line 156, in calculate
reaction.find_lowest_energy_conformers()
File "/apps/anaconda3/envs/autode_142/lib/python3.10/site-packages/autode/utils.py", line 668, in wrapped_function
result = func(reaction)
File "/apps/anaconda3/envs/autode_142/lib/python3.10/site-packages/autode/reactions/reaction.py", line 622, in find_lowest_energy_conformers
mol.find_lowest_energy_conformer(hmethod=h_method)
File "/apps/anaconda3/envs/autode_142/lib/python3.10/site-packages/autode/utils.py", line 239, in wrapped_function
result = func(*args, **kwargs)
File "/apps/anaconda3/envs/autode_142/lib/python3.10/site-packages/autode/species/species.py", line 1443, in find_lowest_energy_conformer
self.conformers.optimise(hmethod)
File "/apps/anaconda3/envs/autode_142/lib/python3.10/site-packages/autode/conformers/conformers.py", line 306, in optimise
return self._parallel_calc("optimise", method, keywords)
File "/apps/anaconda3/envs/autode_142/lib/python3.10/site-packages/autode/conformers/conformers.py", line 288, in _parallel_calc
self[idx] = res.result()
File "/apps/anaconda3/envs/autode_142/lib/python3.10/concurrent/futures/_base.py", line 458, in result
return self.__get_result()
File "/apps/anaconda3/envs/autode_142/lib/python3.10/concurrent/futures/_base.py", line 403, in __get_result
raise self.exception
autode.exceptions.CoordinateTransformFailed: DIC->Cart iterative back-transform did not converge
--------------------
This is my input file in case you want to try:
import autode as ade
ade.Config.n_cores = 4
ade.Config.max_core = 4000
ade.Config.hcode = 'NWChem'
ade.Config.lcode = 'XTB'
r1 = ade.Reactant(smiles='CC(C)=O')
r2 = ade.Reactant(smiles='CN(C)[Li-][O+]1CCCC1')
p1 = ade.Product(smiles='CC(=C)O[Li-][O+]1CCCC1')
p2 = ade.Product(smiles='[H]N(C)C')
rxn = ade.Reaction(r1, r2, p1, p2, name='test_thf_-78C', solvent_name='thf', temp=195.15)
rxn.calculate_reaction_profile(free_energy=True)
quit()
Thanks
from autode.
@eneas77 Hi, does the DA example work? The optimiser and the coordinate systems need some more improvements before it can work on every possible system. I suspect in your case above, one or more angles are becoming linear(180) during optimisation, or becoming <90 when starting from linear. Both of these cases will fail currently.
from autode.
@shoubhikraj Yes, the DA example now works.
from autode.
Related Issues (20)
- ImportError: cannot import name 'cached_property' from 'functools' HOT 2
- Temperature should be a value
- Better optimisers for NEB
- Performing NEB given exact atom-mapping HOT 5
- Diffrence in frequencies from ORCA and autode HOT 5
- Optimiser bugs
- Implement image-pair regeneration with DHS (-GS) and i-EIP
- Increase default SCF cycles in Q-Chem HOT 2
- SMILES input -> molecule bug HOT 3
- Turn of truncation of structure HOT 2
- Failed to create a complete set of delocalised internal coordinates HOT 3
- ORCA error ? CouldNotGetProperty: Could not get energy HOT 2
- Codecov pytest runs are not using xTB on windows CI runs
- Debugging TS failures via file type generated
- Bump Python version in CI
- Problem with decarboxylation reactions HOT 4
- Solvent not being added with G09 for TS quick-reaction coordinate HOT 1
- Autode error "Could not get energy" for the XTB calculation HOT 2
- autode deals with free radicals? HOT 3
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from autode.