Re: [AMBER] sqm convergence failure

From: Hannes Loeffler <Hannes.Loeffler.stfc.ac.uk>
Date: Thu, 8 May 2014 13:49:09 +0100

On Thu, 8 May 2014 08:22:41 -0400
Jason Swails <jason.swails.gmail.com> wrote:

> On Thu, 2014-05-08 at 11:16 +0100, Hannes Loeffler wrote:
> > Hi,
> >
> > I can't get the attached structure (sqm.in) converged. What
> > happens is that the wave function won't converge somewhere along
> > geometry optimisation, at some arbitrary step (not the first step)
> > depending on name list parameters.
> >
> > I have no particular interest in this one. The problem just came up
> > in a larger test set. The structure is ZINC03814828. Same result
> > with ZINC03814832 which is almost the same.
> >
> > The current name list is an outcome of playing around with various
> > parameters but without luck. What I wonder is if this could still
> > be made to converge or if there is a set of known structures which
> > just cannot be converged? If I replace the carboxy-ethyl group
> > with propyl the structure converges at scfconv=1.0D-9.
>
> I think I may be missing something obvious or perhaps I postprocessed
> your sqm input file wrong. If I delete the sqm input variables and
> put '55' at the top of the file (to try and make the file look like
> an XYZ), I get an image which appears to have 3 hydrogen atoms
> floating around somewhere in space (image attached).
>
> If I didn't make a mistake in my hacky sqm-to-xyz file conversion (it
> is _very_ possible that I did), this could certainly explain why the
> structure fails to converge (most of the structure looks fine, but not
> the 3 labeled H atoms), and more reflects a problem with the starting
> structure than with sqm.

Many thanks Jason but I believe the structure is sane. Not sure how
you got this as looks good in my vmd. The C-C distance is 1.54 A and
the C-H distances are all 1.09.

The structure has been generated as part of a convergence protocol
which includes MM pre-minimisation steps with the MMFF94 force field.
The original SDF file contains connectivities and bond order
information so that should be fine, see attached mol2 file (ignore
Openbabel's feeble attempt at recognizing residues :-).


-- 
Scanned by iCritical.


Content-Transfer-Encoding: base64
Content-Disposition: attachment; filename=ZINC03814828.mol2

QDxUUklQT1M+TU9MRUNVTEUKWklOQzAzODE0ODI4CiA1NiA1OCAwIDAgMApTTUFMTApHQVNURUlH
RVIKCkA8VFJJUE9TPkFUT00KICAgICAgMSBDICAgICAgICAgIC0xLjE3MjkgICAtMS4zMzA0ICAg
LTAuMjU4MSBDLjMgICAgIDAgIFVOSzAgICAgICAgLTAuMDM4NQogICAgICAyIEMgICAgICAgICAg
LTIuNDQxMCAgIC0wLjQ1ODggICAtMC4yMDU1IEMuYXIgICAgMCAgVU5LMCAgICAgICAtMC4wMjkw
CiAgICAgIDMgQyAgICAgICAgICAtMy42NzYzICAgLTEuMDM1NSAgIC0wLjU3MDEgQy5hciAgICAw
ICBVTkswICAgICAgIC0wLjAyNjQKICAgICAgNCBDICAgICAgICAgIC00Ljg1MzIgICAtMC4yOTU2
ICAgLTAuNTQyOCBDLmFyICAgIDAgIFVOSzAgICAgICAgLTAuMDUwNgogICAgICA1IEMgICAgICAg
ICAgLTQuODM1MyAgICAxLjAzNTcgICAtMC4xNzY5IEMuYXIgICAgMCAgVU5LMCAgICAgICAgMC4w
MTA4CiAgICAgIDYgQyAgICAgICAgICAtMy42MjkwICAgIDEuNjMzMSAgICAwLjE4MTUgQy5hciAg
ICAwICBVTkswICAgICAgICAwLjIxOTUKICAgICAgNyBDICAgICAgICAgIC0yLjQxOTYgICAgMC45
MDIwICAgIDAuMTk3MyBDLmFyICAgIDAgIFVOSzAgICAgICAgIDAuMDg5MQogICAgICA4IEMgICAg
ICAgICAgLTEuMjYyNCAgICAxLjYxNTEgICAgMC42MDkxIEMuYXIgICAgMCAgVU5LMCAgICAgICAg
MC4xMzc5CiAgICAgIDkgTiAgICAgICAgICAtMS4yOTE2ICAgIDIuOTI4MSAgICAwLjkyMDMgTi5h
ciAgICAwICBVTkswICAgICAgIC0wLjE1NjYKICAgICAxMCBDICAgICAgICAgIC0yLjQ4MDYgICAg
My41NDk4ICAgIDAuODY1MyBDLmFyICAgIDAgIFVOSzAgICAgICAgIDAuMzcxOQogICAgIDExIE4g
ICAgICAgICAgLTIuNDcwMSAgICA0Ljg3MDggICAgMS4yMDgyIE4ucGwzICAgMCAgVU5LMCAgICAg
ICAtMC4yODgyCiAgICAgMTIgTiAgICAgICAgICAtMC4wMjg1ICAgIDEuMDY1OSAgICAwLjc5OTgg
Ti5wbDMgICAwICBVTkswICAgICAgIC0wLjM0MjMKICAgICAxMyBDICAgICAgICAgIC0zLjgxMTQg
ICAtMi40ODQ5ICAgLTAuOTgzOCBDLjMgICAgIDAgIFVOSzAgICAgICAgIDAuMDM0NgogICAgIDE0
IE4gICAgICAgICAgLTMuNzE0NiAgIC0zLjM1MDYgICAgMC4xODc1IE4ucGwzICAgMCAgVU5LMCAg
ICAgICAtMC4zMzk3CiAgICAgMTUgQyAgICAgICAgICAtMy43MTE2ICAgLTQuNzE0OSAgICAwLjA4
NDMgQy5hciAgICAwICBVTkswICAgICAgICAwLjAyNzQKICAgICAxNiBDICAgICAgICAgIC0zLjc0
OTAgICAtNS4zNTcwICAgLTEuMTU3OCBDLmFyICAgIDAgIFVOSzAgICAgICAgLTAuMDQwMQogICAg
IDE3IEMgICAgICAgICAgLTMuNzUzNCAgIC02Ljc0MzMgICAtMS4yMzQxIEMuYXIgICAgMCAgVU5L
MCAgICAgICAtMC4wNDc5CiAgICAgMTggQyAgICAgICAgICAtMy43MjY3ICAgLTcuNTE1NyAgIC0w
LjA2NzggQy5hciAgICAwICBVTkswICAgICAgICAwLjAzNjgKICAgICAxOSBDICAgICAgICAgIC0z
LjY4MTcgICAtNi44Njg5ICAgIDEuMTcwMSBDLmFyICAgIDAgIFVOSzAgICAgICAgLTAuMDQ3OQog
ICAgIDIwIEMgICAgICAgICAgLTMuNjczNCAgIC01LjQ3OTggICAgMS4yNDgzIEMuYXIgICAgMCAg
VU5LMCAgICAgICAtMC4wNDAxCiAgICAgMjEgQyAgICAgICAgICAtMy43NjQzICAgLTkuMDE1MCAg
IC0wLjEwMjQgQy4yICAgICAwICBVTkswICAgICAgICAwLjI0NDYKICAgICAyMiBPICAgICAgICAg
IC0zLjc4NjAgICAtOS42NzY4ICAgIDAuOTE4NyBPLjIgICAgIDAgIFVOSzAgICAgICAgLTAuMjcw
MgogICAgIDIzIE4gICAgICAgICAgLTMuODU4NyAgIC05LjU3NzEgICAtMS4zMjUwIE4uYW0gICAg
MSAgR0xVMSAgICAgICAtMC4zMDM0CiAgICAgMjQgQ0EgICAgICAgICAtMy43Njk2ICAtMTEuMDE3
MCAgIC0xLjU4NzIgQy4zICAgICAxICBHTFUxICAgICAgICAwLjA2MzIKICAgICAyNSBIQSAgICAg
ICAgIC00LjI4OTYgIC0xMS41NDc3ICAgLTAuNzg1MyBIICAgICAgIDEgIEdMVTEgICAgICAgIDAu
MDU2MAogICAgIDI2IENCICAgICAgICAgLTQuNDQ2MyAgLTExLjM4MTMgICAtMi45MzQ4IEMuMyAg
ICAgMSAgR0xVMSAgICAgICAtMC4wMjQxCiAgICAgMjcgQ0cgICAgICAgICAtNS45NDQ1ICAtMTAu
OTc2MCAgIC0yLjk3MzQgQy4zICAgICAxICBHTFUxICAgICAgIC0wLjAxMTIKICAgICAyOCBDRCAg
ICAgICAgIC02LjY0NTUgIC0xMS4zNTEyICAgLTQuMjcwOCBDLjIgICAgIDEgIEdMVTEgICAgICAg
IDAuMDQxOAogICAgIDI5IE9FMSAgICAgICAgLTYuMDk5MyAgLTExLjgzNjcgICAtNS4yNDE3IE8u
Y28yICAgMSAgR0xVMSAgICAgICAtMC41NDk4CiAgICAgMzAgQyAgICAgICAgICAtMi4zMTMzICAt
MTEuNDc0OSAgIC0xLjU0ODEgQy4yICAgICAxICBHTFUxICAgICAgICAwLjA2MjAKICAgICAzMSBP
ICAgICAgICAgIC0xLjM1ODkgIC0xMC43MjIxICAgLTEuNDY0NCBPLmNvMiAgIDEgIEdMVTEgICAg
ICAgLTAuNTQ3OQogICAgIDMyIEggICAgICAgICAgLTAuODYwMiAgIC0xLjYwMTUgICAgMC43NTMy
IEggICAgICAgMCAgVU5LMCAgICAgICAgMC4wMjc4CiAgICAgMzMgSCAgICAgICAgICAtMS4yODM2
ICAgLTIuMjYwNCAgIC0wLjgwOTkgSCAgICAgICAwICBVTkswICAgICAgICAwLjAyNzgKICAgICAz
NCBIICAgICAgICAgIC0wLjM1NjggICAtMC44MTY5ICAgLTAuNzY3NSBIICAgICAgIDAgIFVOSzAg
ICAgICAgIDAuMDI3OAogICAgIDM1IEggICAgICAgICAgLTUuNzk4OSAgIC0wLjc2NDUgICAtMC44
MTIyIEggICAgICAgMCAgVU5LMCAgICAgICAgMC4wNjIzCiAgICAgMzYgSCAgICAgICAgICAtNS43
NjkzICAgIDEuNTkzMiAgIC0wLjE3MTYgSCAgICAgICAwICBVTkswICAgICAgICAwLjA2ODQKICAg
ICAzNyBIICAgICAgICAgIC0zLjMyNTYgICAgNS40MTA2ICAgIDEuMjcxOSBIICAgICAgIDAgIFVO
SzAgICAgICAgIDAuMTQ5MwogICAgIDM4IEggICAgICAgICAgLTEuNTk4NiAgICA1LjI5NDcgICAg
MS41MTM2IEggICAgICAgMCAgVU5LMCAgICAgICAgMC4xNDkzCiAgICAgMzkgSCAgICAgICAgICAg
MC43NTMzICAgIDEuNTI0NyAgICAwLjM0MzQgSCAgICAgICAwICBVTkswICAgICAgICAwLjE0MzcK
ICAgICA0MCBIICAgICAgICAgICAwLjAyNjEgICAgMC4wNTMyICAgIDAuODMwMyBIICAgICAgIDAg
IFVOSzAgICAgICAgIDAuMTQzNwogICAgIDQxIEggICAgICAgICAgLTQuNzc1MiAgIC0yLjY2MjEg
ICAtMS40NjkxIEggICAgICAgMCAgVU5LMCAgICAgICAgMC4wNTExCiAgICAgNDIgSCAgICAgICAg
ICAtMy4wNDM3ICAgLTIuNzUyMyAgIC0xLjcxMzQgSCAgICAgICAwICBVTkswICAgICAgICAwLjA1
MTEKICAgICA0MyBIICAgICAgICAgIC0zLjY2ODQgICAtMi45MjcxICAgIDEuMTA3NiBIICAgICAg
IDAgIFVOSzAgICAgICAgIDAuMTQ2NAogICAgIDQ0IEggICAgICAgICAgLTMuNzc5MyAgIC00Ljc5
NzUgICAtMi4wODc5IEggICAgICAgMCAgVU5LMCAgICAgICAgMC4wNjM2CiAgICAgNDUgSCAgICAg
ICAgICAtMy43ODUxICAgLTcuMjA1MiAgIC0yLjIxNzcgSCAgICAgICAwICBVTkswICAgICAgICAw
LjA2MjYKICAgICA0NiBIICAgICAgICAgIC0zLjY1ODkgICAtNy40NTI4ICAgIDIuMDg2NSBIICAg
ICAgIDAgIFVOSzAgICAgICAgIDAuMDYyNgogICAgIDQ3IEggICAgICAgICAgLTMuNjQ0MSAgIC01
LjAwMjYgICAgMi4yMjQ5IEggICAgICAgMCAgVU5LMCAgICAgICAgMC4wNjM2CiAgICAgNDggSCAg
ICAgICAgICAtMy43OTUwICAgLTguOTkwMCAgIC0yLjEzOTggSCAgICAgICAxICBHTFUxICAgICAg
ICAwLjE0OTYKICAgICA0OSBIQjEgICAgICAgIC0zLjkxMzkgIC0xMC44OTYxICAgLTMuNzU1NiBI
ICAgICAgIDEgIEdMVTEgICAgICAgIDAuMDI5MAogICAgIDUwIEhCMiAgICAgICAgLTQuMzcyMyAg
LTEyLjQ1OTkgICAtMy4wODkxIEggICAgICAgMSAgR0xVMSAgICAgICAgMC4wMjkwCiAgICAgNTEg
SEcxICAgICAgICAtNi40Njk0ICAtMTEuNDUzMyAgIC0yLjE0MjQgSCAgICAgICAxICBHTFUxICAg
ICAgICAwLjAzMjUKICAgICA1MiBIRzIgICAgICAgIC02LjA0MzQgICAtOS44OTczICAgLTIuODMx
OCBIICAgICAgIDEgIEdMVTEgICAgICAgIDAuMDMyNQogICAgIDUzIE9FMiAgICAgICAgLTcuODg2
OCAgLTEwLjkzNzMgICAtNC4zMDk5IE8uY28yICAgMSAgR0xVMSAgICAgICAtMC41NDk4CiAgICAg
NTQgT1hUICAgICAgICAtMi4xNzczICAtMTIuNzcyNyAgIC0xLjQ2MzIgTy5jbzIgICAxICBHTFUx
ICAgICAgIC0wLjU0NzkKICAgICA1NSBOICAgICAgICAgIC0zLjY0NzYgICAgMi45NDg3ICAgIDAu
NTMyOCBOLmFyICAgIDAgIFVOSzAgICAgICAgLTAuMTc0MQogICAgIDU2IEggICAgICAgICAgLTQu
NTg5OCAgICAzLjUxNjEgICAgMC41NTIwIEggICAgICAgMCAgVU5LMCAgICAgICAgMC40NTY1CkA8
VFJJUE9TPkJPTkQKICAgICAxICAgICAxICAgICAyICAgIDEKICAgICAyICAgICAxICAgIDMyICAg
IDEKICAgICAzICAgICAxICAgIDMzICAgIDEKICAgICA0ICAgICAxICAgIDM0ICAgIDEKICAgICA1
ICAgICAyICAgICA3ICAgYXIKICAgICA2ICAgICAyICAgICAzICAgYXIKICAgICA3ICAgICAzICAg
ICA0ICAgYXIKICAgICA4ICAgICAzICAgIDEzICAgIDEKICAgICA5ICAgICA0ICAgICA1ICAgYXIK
ICAgIDEwICAgICA0ICAgIDM1ICAgIDEKICAgIDExICAgICA1ICAgICA2ICAgYXIKICAgIDEyICAg
ICA1ICAgIDM2ICAgIDEKICAgIDEzICAgICA2ICAgICA3ICAgYXIKICAgIDE0ICAgICA2ICAgIDU1
ICAgYXIKICAgIDE1ICAgICA3ICAgICA4ICAgYXIKICAgIDE2ICAgICA4ICAgICA5ICAgYXIKICAg
IDE3ICAgICA4ICAgIDEyICAgIDEKICAgIDE4ICAgICA5ICAgIDEwICAgYXIKICAgIDE5ICAgIDEw
ICAgIDExICAgIDEKICAgIDIwICAgIDEwICAgIDU1ICAgYXIKICAgIDIxICAgIDExICAgIDM3ICAg
IDEKICAgIDIyICAgIDExICAgIDM4ICAgIDEKICAgIDIzICAgIDEyICAgIDM5ICAgIDEKICAgIDI0
ICAgIDEyICAgIDQwICAgIDEKICAgIDI1ICAgIDEzICAgIDE0ICAgIDEKICAgIDI2ICAgIDEzICAg
IDQxICAgIDEKICAgIDI3ICAgIDEzICAgIDQyICAgIDEKICAgIDI4ICAgIDE0ICAgIDE1ICAgIDEK
ICAgIDI5ICAgIDE0ICAgIDQzICAgIDEKICAgIDMwICAgIDE1ICAgIDIwICAgYXIKICAgIDMxICAg
IDE1ICAgIDE2ICAgYXIKICAgIDMyICAgIDE2ICAgIDE3ICAgYXIKICAgIDMzICAgIDE2ICAgIDQ0
ICAgIDEKICAgIDM0ICAgIDE3ICAgIDE4ICAgYXIKICAgIDM1ICAgIDE3ICAgIDQ1ICAgIDEKICAg
IDM2ICAgIDE4ICAgIDE5ICAgYXIKICAgIDM3ICAgIDE4ICAgIDIxICAgIDEKICAgIDM4ICAgIDE5
ICAgIDIwICAgYXIKICAgIDM5ICAgIDE5ICAgIDQ2ICAgIDEKICAgIDQwICAgIDIwICAgIDQ3ICAg
IDEKICAgIDQxICAgIDIxICAgIDIyICAgIDIKICAgIDQyICAgIDIxICAgIDIzICAgYW0KICAgIDQz
ICAgIDIzICAgIDI0ICAgIDEKICAgIDQ0ICAgIDIzICAgIDQ4ICAgIDEKICAgIDQ1ICAgIDI0ICAg
IDI1ICAgIDEKICAgIDQ2ICAgIDI0ICAgIDI2ICAgIDEKICAgIDQ3ICAgIDI0ICAgIDMwICAgIDEK
ICAgIDQ4ICAgIDI2ICAgIDI3ICAgIDEKICAgIDQ5ICAgIDI2ICAgIDQ5ICAgIDEKICAgIDUwICAg
IDI2ICAgIDUwICAgIDEKICAgIDUxICAgIDI3ICAgIDI4ICAgIDEKICAgIDUyICAgIDI3ICAgIDUx
ICAgIDEKICAgIDUzICAgIDI3ICAgIDUyICAgIDEKICAgIDU0ICAgIDI4ICAgIDI5ICAgYXIKICAg
IDU1ICAgIDI4ICAgIDUzICAgYXIKICAgIDU2ICAgIDMwICAgIDMxICAgYXIKICAgIDU3ICAgIDMw
ICAgIDU0ICAgYXIKICAgIDU4ICAgIDU1ICAgIDU2ICAgIDEK


_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Thu May 08 2014 - 06:00:04 PDT
Custom Search