Re: [AMBER] ATP charge convergence

From: Eric Pettersen <pett.cgl.ucsf.edu>
Date: Thu, 18 Mar 2010 17:17:58 -0700

I'll just cut and paste the sqm.in file, since it's short:

Run semi-empirical minimization
  &qmmm
     qm_theory='AM1', grms_tol=0.02,
   tight_p_conv=0, scfconv=1.d-7, qmcharge=-4,
  /
   15 PG 7.7270 -127.7270 -164.2590
    8 O1G 7.1390 -129.0740 -164.0840
    8 O2G 7.1960 -126.6730 -163.1630
    8 O3G 7.3760 -127.1060 -165.6990
   15 PB 10.1260 -126.3960 -164.6610
    8 O1B 9.4570 -125.1840 -164.1420
    8 O2B 10.2690 -126.5480 -166.1250
    8 O3B 9.3310 -127.6920 -164.1370
   15 PA 11.7330 -126.4560 -162.3570
    8 O1A 12.3780 -125.1570 -162.0590
    8 O2A 10.4420 -126.7610 -161.7010
    8 O3A 11.5660 -126.5470 -163.9550
    8 O5' 12.7770 -127.6400 -162.0360
    6 C5' 13.0970 -128.6550 -162.9930
    6 C4' 13.0520 -130.0350 -162.3350
    8 O4' 13.7310 -131.0510 -163.0970
    6 C3' 13.7730 -129.9820 -160.9880
    8 O3' 12.8470 -129.5630 -159.9830
    6 C2' 14.0860 -131.4660 -160.8180
    8 O2' 12.8750 -132.1970 -160.6110
    6 C1' 14.6320 -131.7500 -162.2180
    7 N9 16.0030 -131.1990 -162.3840
    6 C8 16.3250 -129.9360 -162.6610
    7 N7 17.6490 -129.8280 -162.7800
    6 C5 18.1790 -131.0310 -162.5780
    6 C6 19.4730 -131.5300 -162.6070
    7 N6 20.5090 -130.7300 -162.8610
    7 N1 19.6790 -132.8970 -162.3780
    6 C2 18.5740 -133.7230 -162.1230
    7 N3 17.3560 -133.2130 -162.1020
    6 C4 17.1230 -131.9070 -162.3210
    1 H5' 14.0980 -128.4760 -163.3860
    1 H5'' 12.3770 -128.6210 -163.8100
    1 H4' 12.0140 -130.3300 -162.1800
    1 H3' 14.6710 -129.3640 -161.0160
    1 HO3' 13.2910 -129.5270 -159.1330
    1 H2' 14.8240 -131.6500 -160.0370
    1 HO2' 13.0310 -132.9030 -159.9790
    1 H1' 14.6200 -132.8200 -162.4260
    1 H8 15.6210 -129.1240 -162.7720
    1 H61 20.3580 -129.7460 -163.0310
    1 H62 21.4450 -131.1080 -162.8840
    1 H2 18.7200 -134.7780 -161.9440

On Mar 18, 2010, at 5:11 PM, Eric Pettersen wrote:

> Hmm, I don't see the attachments that were in the mail I sent. Is
> there some magic to avoid having them stripped?
>
> --Eric
>
> On Mar 18, 2010, at 5:07 PM, Eric Pettersen wrote:
>
>> Hi,
>> I cannot get sqm to produce SCF convergence for an ATP derived
>> from PDB 1xfa, even with the convergence criteria cranked way
>> down. Is it because of the high charge (-4)? Is there anything
>> that can be done to get it to converge? I've attached the sqm.in
>> file (with the lowered criteria), sqm.out file, and the mol2 file
>> that was fed to Antechamber originally. The sqm command line I
>> used was the same as antechamber uses: sqm -O -i sqm.in -o sqm.out.
>>
>> --Eric
>>
>> Eric Pettersen
>> UCSF Computer Graphics Lab
>> http://www.cgl.ucsf.edu
>>
>>
>>
>> _______________________________________________
>> AMBER mailing list
>> AMBER.ambermd.org
>> http://lists.ambermd.org/mailman/listinfo/amber
>
>
> _______________________________________________
> AMBER mailing list
> AMBER.ambermd.org
> http://lists.ambermd.org/mailman/listinfo/amber


_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Thu Mar 18 2010 - 17:30:04 PDT
Custom Search