Comments (6)
I believe that the default convergence tolerance for the SQM calculation was decreased from AmberTools 14 to AmberTools 15. So while AmberTools 14 decided that calculation was converged, AmberTools 15 did not. You can try modifying common.c
in the antechamber source code and changing the convergence criteria (making it larger) to see if it starts converging again. But a better bet might be to figure out why the molecule is not minimizing in the first place.
Also, the Amber mailing list is a better place for these kinds of questions (and make sure you're using the actual AmberTools 15 release code, not ambermini). See http://lists.ambermd.org/mailman/listinfo/amber for info about the mailing lists.
from ambermini.
@andrrizzi: Note that we should be using the canonical AM1-BCC charging scheme using OpenEye tools for all of our work, though it would be nice to be able to support pure AMBER workflows. This is implemented in get_charges()
in openmoltools
for now.
from ambermini.
@swails is indeed right, though. From here, Ross Walker changed this:
/* RCW Below is the old setting from AMBER 12, 14. It is ridiculously tight.
64 * As of AmberTools 15 remove tight_p_conv and set grms_tol to 0.0005
65 * tests show it makes no major difference to the charges - max diff = 0.001.
66 *
67 * strcpy((*minfo).skeyword, " qm_theory='AM1', grms_tol=0.0002,\n tight_p_conv=1, scfconv=1.d-10, ");
68 *
69 */
70 strcpy((*minfo).skeyword, " qm_theory='AM1', grms_tol=0.0005,\n scfconv=1.d-10, ");
from ambermini.
I see, thanks for the pointer!
@jchodera Yes, the molecule goes through OpenEye's canonical AM1-BCC charging scheme before going to antechamber
(did you mean I should use it afterwards?).
I do use epik
to generate a protonation state after OpenEye though. If antechamber
uses the user charges as a starting point, I could try to repeat the scheme and generate new charges to feed to antechamber
. Maybe this will solve the minimization failure.
from ambermini.
@jchodera Yes, the molecule goes through OpenEye's canonical AM1-BCC charging scheme before going to antechamber (did you mean I should use it afterwards?).
If that's the case, it shouldn't run through sqm
when it hits antechamber
. It should just read the charges from the mol2
if no -c
argument is specified (or maybe it's -c rc
? but I think that is for an external charge file).
I do use epik to generate a protonation state after OpenEye though. If antechamber uses the user charges as a starting point, I could try to repeat the scheme and generate new charges to feed to antechamber. Maybe this will solve the minimization failure.
We should not be using sqm
at all. The OpenEye charging should happen after epik
is used to enumerate protonation states. Then these charges should be used in antechamber
.
from ambermini.
Got it, thanks! I'll close this for now.
from ambermini.
Related Issues (20)
- Update for AmberTools 15? HOT 14
- Ran into issue with binstar build HOT 2
- Conda package tests fails despite succeeding
- Conda recipe tests fails despite succeeding HOT 2
- 'sqm -h' fails HOT 3
- ambermini conda package still failing HOT 3
- Issues with osx ambermini package HOT 14
- Set up appveyor HOT 5
- pdb4amber? HOT 3
- bondtype and atomtype executables complain AMBERHOME is not set HOT 7
- Windows uses? HOT 15
- Add reduce to ambermini? HOT 5
- Update to AmberTools 16 HOT 47
- How to handle AmberMini Development going forward HOT 11
- AmberTools17 update HOT 5
- OPC water and 12_6_4 parameters
- Incorrectly translated equation in gaff
- Compiling problems of arpack with gfortran-10
- How to change sqm.out directory when using antechamber to convert pdb to mol2? HOT 2
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 ambermini.