Hello Carsten,

Thanks for your reply! The system I'm testing is a 128x128x128 nm box 
containing only water. I tried using double precision, but the error is the 
same. There are lots of errors like this:

Program:     gmx mdrun, version 2019.3
Source file: src/gromacs/mdlib/nbnxn_grid.cpp (line 1164)
MPI rank:    9 (out of 16)

Fatal error:
grid cell cx -207 cy 0 out of range (max 23 372)
atom 0.000000 0.000000 0.000000, grid->c0 72.000000 0.000000


There are many atoms with coordinates 0.000000 0.000000 0.000000, which I think 
should be an error.

Zhang

> -----原始邮件-----
> 发件人: "Kutzner, Carsten" <ckut...@gwdg.de>
> 发送时间: 2019-07-09 22:03:06 (星期二)
> 收件人: "<gmx-us...@gromacs.org> GROMACS users" <gmx-us...@gromacs.org>
> 抄送: 
> 主题: Re: [gmx-users] what's the maximum number of atoms that GROMACS can       
> simulate?
> 
> Hi Zhang,
> 
> > On 9. Jul 2019, at 15:16, 张驭洲 <zhangyuzho...@mails.ucas.edu.cn> wrote:
> > 
> > Hello,
> > 
> > I want to know if there is a maximum number of atoms that GROMACS can 
> > simulate, or if there are any bugs in the code that cause the error in 
> > computing the coordinate of atoms when the number of atoms is very large? 
> > I'm trying to simulate some very large systems using GROMACS-2019.3. 
> > Unfortunately, the energy minimization crashed during the run. I tried 
> > using more than one compute nodes, for example, 4, 8 or 32, and I'm sure 
> > the memory is not exhausted. The system I tried contains more than 100 
> > million atoms.
> There is no limit on the number of atoms that you can simulate in GROMACS.
> 
> However, if the simulation systems become very large in at least one 
> dimension,
> you might need double precision, as single precision could be insufficient 
> to accurately store the positions of atoms that are far off the origin.
> This can then lead to a positive energy drift which could become so large
> that the system gets instable.
> 
> A while time ago we have observed such problems with large systems of mostly
> vacuum, which contained several water droplets. The droplets that were 
> farthest
> away from the origin would start to heat up and eventually evaporate after 
> some time, unless a double precision version of mdrun was used.
> 
> Carsten
> 
> 
> > 
> > Zhang
> > -- 
> > Gromacs Users mailing list
> > 
> > * Please search the archive at 
> > http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!
> > 
> > * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
> > 
> > * For (un)subscribe requests visit
> > https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send 
> > a mail to gmx-users-requ...@gromacs.org.
> 
> --
> Dr. Carsten Kutzner
> Max Planck Institute for Biophysical Chemistry
> Theoretical and Computational Biophysics
> Am Fassberg 11, 37077 Goettingen, Germany
> Tel. +49-551-2012313, Fax: +49-551-2012302
> http://www.mpibpc.mpg.de/grubmueller/kutzner
> http://www.mpibpc.mpg.de/grubmueller/sppexa
> 
> -- 
> Gromacs Users mailing list
> 
> * Please search the archive at 
> http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!
> 
> * Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
> 
> * For (un)subscribe requests visit
> https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
> mail to gmx-users-requ...@gromacs.org.
-- 
Gromacs Users mailing list

* Please search the archive at 
http://www.gromacs.org/Support/Mailing_Lists/GMX-Users_List before posting!

* Can't post? Read http://www.gromacs.org/Support/Mailing_Lists

* For (un)subscribe requests visit
https://maillist.sys.kth.se/mailman/listinfo/gromacs.org_gmx-users or send a 
mail to gmx-users-requ...@gromacs.org.

Reply via email to