Hello Zhicheng,
Thanks for your interest in our work.
In the topology file, the atomic number was set as 1 for the dummy atom, so
the program treat all the dummy atoms as hydrogen atoms, though the vdW
parameters are set ZERO.
You can set the atomic number of the dummy atom as *-1* in the library
file (*.lib). I think it can fix the problem.
Please let me know if it works or not.
All the best,
Qinghua
On 06/21/2018 11:45 PM, Zuo, Zhicheng wrote:
> Hi,
> The recently released AMBER18 adds GPU support for 12-6-4 ion potentials.
> I did some tests, and found that it did works for the case where an ion is treated
> traditionally as a single point charge with the 12-6-4 LJ potential. However, for the simulation
> of the recent multisite ion (dummy atom) model with the 12-6-4 potential (https://pubs.acs.org/doi/10.1021/acs.jpclett.7b02358),
> I got the following error from pmemd.cuda:
> "Too many hydrogens for a hydrogen network, exiting."
>
> Interestingly, there is no any problem with the cpu version of pmemd for the 12-6-4 multisite ion model above.
>
> It seems that this issue is related to implementation of pmemd.cuda toward such a special divalent ion model.
> Can someone help to figure it out? Thanks!
>
> Best,
> Zhicheng Zuo
> Postdoc research associate from College of Pharmacy,
> UNT Health Science Center, Fort Worth, TX, USA.
>
> _______________________________________________
> 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 Jun 21 2018 - 16:00:02 PDT
