Re: [gmx-users] How to suppress the error "X particles communicated to PME node Y are more than a cell length out of the domain decomposition cell of their charge group
On 27/12/2010 9:22 AM, WU Yanbin wrote: Hi, Mark, Thanks for the reply. As suggested, I tried a "water in vacuum" case. Initially the water droplet is a cubic 4nm-by-4nm-by-4nm water box, located in the middle of the simulation box. Everywhere else is just vacuum. The simulation box size is 8nm by 8nm by 8nm. SPC/E model is used to describe interaction between water molecules. Such system is first equilibrated with "steep" option. The "mdrun" simulation goes with no problem with parallel running on 32 cpus, even though there are occasionaly one or two water molecules move very fast (Mean square displacement being 200 times fast than bulk water) in the vacuum. OK, that sounds like behaviour that might be expected... some molecules evaporate. Check the trajectory to be confident this is the situation. When I switch to 64 cpus, the error "X particles communicated to PME node Y are more than a cell length out of the domain decomposition cell of their charge group" appears. The most plausible reason for this is that the above evaporated molecules are moving so fast that they're doing what the error message says - travelling more than width of a box in one integration step. The DD implementation is predicated on that being impossible. You might succeed by arranging for the largest possible internal DD cell diameter, i.e. a 4x4x4 DD. mdrun -npme 0 might choose this, otherwise use the hidden options to mdrun (use mdrun -hidden) to use that DD with -npme 0. Otherwise, you'll need to accept the fact that there are "engineering" constraints on efficient parallelization algorithms, and not every situation can be catered for. For example, a simulation of a bulk LJ fluid would fail if you used so many processors that the cell diameter was too small with respect to the distribution of particle speeds. Mark Below is the parameter file I'm using. integrator = md tinit= 0 dt = 0.002 nsteps = 500 comm-mode= Linear nstcomm = 1 comm-grps= energygrps = nstlist = 1 ns_type = grid pbc = xyz rlist= 1.5 coulombtype = PME rcoulomb-switch = 0 rcoulomb = 1.5 epsilon-r= 1 vdw-type = Cut-off ;Switch rvdw-switch = 1.0 rvdw = 1.5 fourierspacing = 0.12 fourier_nx = 0 fourier_ny = 0 fourier_nz = 0 pme_order= 4 ewald_rtol = 1e-05 ewald_geometry = 3d epsilon_surface = 0 optimize_fft = no Tcoupl = nose-hoover tc-grps = System tau_t= 1.0 ref_t= 300 Pcoupl = no ;Parrinello-Rahman ;no ;berendsen Pcoupltype = isotropic tau_p= 1.0 compressibility = 4.5e-5 ref_p= 1.0 gen_vel = no gen_temp = 300 gen_seed = 2008 constraints = none constraint-algorithm = Lincs Shake-SOR= no shake-tol= 1e-04 lincs-order = 4 lincs-warnangle = 30 energygrp_excl = Is there any parameters configuration wrong with the simulation? Or is there any way to go around this error? Any tip is appreciated. Thank you. Best, Yanbin -- Message: 1 Date: Thu, 02 Dec 2010 17:24:18 +1100 From: Mark Abraham mailto:mark.abra...@anu.edu.au>> Subject: Re: [gmx-users] How to suppress the error "X particles communicated to PME node Y are more than a cell length out of the domain decomposition cell of their charge group" To: Discussion list for GROMACS users mailto:gmx-users@gromacs.org>> Message-ID: <4cf73b92.40...@anu.edu.au <mailto:4cf73b92.40...@anu.edu.au>> Content-Type: text/plain; charset=ISO-8859-1; format=flowed On 2/12/2010 4:16 PM, WU Yanbin wrote: > Dear GMXers, > > I'm running a simulation of water contact angle measurement on top of > graphite surface. > Initially a water cubic box is placed on two-layer graphite surface > with the rest of the box being vacuum. The water droplet is relaxed > during the simulation to develop a spherical shape. > > An error of "X particles communicated to PME node Y are more than a > cell length out of the domain decomposition cell of their charge >
Re: [gmx-users] How to suppress the error "X particles communicated to PME node Y are more than a cell length out of the domain decomposition cell of their charge group
Hi, Mark, Thanks for the reply. As suggested, I tried a "water in vacuum" case. Initially the water droplet is a cubic 4nm-by-4nm-by-4nm water box, located in the middle of the simulation box. Everywhere else is just vacuum. The simulation box size is 8nm by 8nm by 8nm. SPC/E model is used to describe interaction between water molecules. Such system is first equilibrated with "steep" option. The "mdrun" simulation goes with no problem with parallel running on 32 cpus, even though there are occasionaly one or two water molecules move very fast (Mean square displacement being 200 times fast than bulk water) in the vacuum. When I switch to 64 cpus, the error "X particles communicated to PME node Y are more than a cell length out of the domain decomposition cell of their charge group" appears. Below is the parameter file I'm using. integrator = md tinit= 0 dt = 0.002 nsteps = 500 comm-mode= Linear nstcomm = 1 comm-grps= energygrps = nstlist = 1 ns_type = grid pbc = xyz rlist= 1.5 coulombtype = PME rcoulomb-switch = 0 rcoulomb = 1.5 epsilon-r= 1 vdw-type = Cut-off ;Switch rvdw-switch = 1.0 rvdw = 1.5 fourierspacing = 0.12 fourier_nx = 0 fourier_ny = 0 fourier_nz = 0 pme_order= 4 ewald_rtol = 1e-05 ewald_geometry = 3d epsilon_surface = 0 optimize_fft = no Tcoupl = nose-hoover tc-grps = System tau_t= 1.0 ref_t= 300 Pcoupl = no ;Parrinello-Rahman ;no ;berendsen Pcoupltype = isotropic tau_p= 1.0 compressibility = 4.5e-5 ref_p= 1.0 gen_vel = no gen_temp = 300 gen_seed = 2008 constraints = none constraint-algorithm = Lincs Shake-SOR= no shake-tol= 1e-04 lincs-order = 4 lincs-warnangle = 30 energygrp_excl = Is there any parameters configuration wrong with the simulation? Or is there any way to go around this error? Any tip is appreciated. Thank you. Best, Yanbin > > -- > > Message: 1 > Date: Thu, 02 Dec 2010 17:24:18 +1100 > From: Mark Abraham > Subject: Re: [gmx-users] How to suppress the error "X particles > communicatedto PME node Y are more than a cell length out of the >domain decomposition cell of their charge group" > To: Discussion list for GROMACS users > Message-ID: <4cf73b92.40...@anu.edu.au> > Content-Type: text/plain; charset=ISO-8859-1; format=flowed > > On 2/12/2010 4:16 PM, WU Yanbin wrote: > > Dear GMXers, > > > > I'm running a simulation of water contact angle measurement on top of > > graphite surface. > > Initially a water cubic box is placed on two-layer graphite surface > > with the rest of the box being vacuum. The water droplet is relaxed > > during the simulation to develop a spherical shape. > > > > An error of "X particles communicated to PME node Y are more than a > > cell length out of the domain decomposition cell of their charge > > group" was encountered. > > And I have read the suggested solutions at the link below > > > http://www.gromacs.org/Documentation/Errors#X_particles_communicated_to_PME_node_Y_are_more_than_a_cell_length_out_of_the_domain_decomposition_cell_of_their_charge_group > . > > > > I guess the reason for this error in my case is because of the vacuum > > such that the water molecules at the boundary of the droplet can move > > fast. I have check the trajectory and the simulation is OK. > > I doubt the simulation is OK. This error message is one of several that > can happen when the system is not well-enough conditioned for the MD to > be stable. See www.gromacs.org/Documentation/Terminology/Blowing_Up. > Here, you have atoms moving much faster than GROMACS was engineered to > expect. > > You should be confident that a water drop in a vacuum, and your graphite > surface are both stable on their own before you try the wetting simulation. > > > > > For this situation, is there a way of suppressing this error? Or what > > else can I do? > > Work out why it's poorly conditioned. > > Mark > --
Re: [gmx-users] How to suppress the error "X particles communicated to PME node Y are more than a cell length out of the domain decomposition cell of their charge group"
On Dec 2, 2010, at 6:16 AM, WU Yanbin wrote: > Dear GMXers, > > I'm running a simulation of water contact angle measurement on top of > graphite surface. > Initially a water cubic box is placed on two-layer graphite surface with the > rest of the box being vacuum. The water droplet is relaxed during the > simulation to develop a spherical shape. > > An error of "X particles communicated to PME node Y are more than a cell > length out of the domain decomposition cell of their charge group" was > encountered. > And I have read the suggested solutions at the link below > http://www.gromacs.org/Documentation/Errors#X_particles_communicated_to_PME_node_Y_are_more_than_a_cell_length_out_of_the_domain_decomposition_cell_of_their_charge_group. > > I guess the reason for this error in my case is because of the vacuum such > that the water molecules at the boundary of the droplet can move fast. I have > check the trajectory and the simulation is OK. > > For this situation, is there a way of suppressing this error? Or what else > can I do? If the system is small enough, you can run it on a single core and the problem cannot occur. You could also try to use particle decomposition (-pd) instead of domain decomposition. Or use less domains, i.e. less cores in total or at least less PP nodes if you use PME/PP splitting. This will at least reduce the probability for the problem to occur. Carsten > > PS: the GROMACS version I'm using is GROMACS4.5. > > Thank you. > > Best, > Yanbin > -- > gmx-users mailing listgmx-users@gromacs.org > http://lists.gromacs.org/mailman/listinfo/gmx-users > Please search the archive at > http://www.gromacs.org/Support/Mailing_Lists/Search before posting! > Please don't post (un)subscribe requests to the list. Use the > www interface or send it to gmx-users-requ...@gromacs.org. > Can't post? Read http://www.gromacs.org/Support/Mailing_Lists -- 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/home/grubmueller/ihp/ckutzne -- gmx-users mailing listgmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/Search before posting! Please don't post (un)subscribe requests to the list. Use the www interface or send it to gmx-users-requ...@gromacs.org. Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
Antw: [gmx-users] How to suppress the error "X particles communicated to PME node Y are more than a cell length out of the domain decomposition cell of their charge group"
Hello, Have you also tried this simulation with a velocity rescaling thermostat in the equilibration period? Could it be that it then would work better? Bests, Emanuel >>> WU Yanbin 02.12.10 6.16 Uhr >>> Dear GMXers, I'm running a simulation of water contact angle measurement on top of graphite surface. Initially a water cubic box is placed on two-layer graphite surface with the rest of the box being vacuum. The water droplet is relaxed during the simulation to develop a spherical shape. An error of "X particles communicated to PME node Y are more than a cell length out of the domain decomposition cell of their charge group" was encountered. And I have read the suggested solutions at the link below http://www.gromacs.org/Documentation/Errors#X_particles_communicated_to_PME_node_Y_are_more_than_a_cell_length_out_of_the_domain_decomposition_cell_of_their_charge_group. I guess the reason for this error in my case is because of the vacuum such that the water molecules at the boundary of the droplet can move fast. I have check the trajectory and the simulation is OK. For this situation, is there a way of suppressing this error? Or what else can I do? PS: the GROMACS version I'm using is GROMACS4.5. Thank you. Best, Yanbin -- gmx-users mailing listgmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/Search before posting! Please don't post (un)subscribe requests to the list. Use the www interface or send it to gmx-users-requ...@gromacs.org. Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
Re: [gmx-users] How to suppress the error "X particles communicated to PME node Y are more than a cell length out of the domain decomposition cell of their charge group"
On 2/12/2010 4:16 PM, WU Yanbin wrote: Dear GMXers, I'm running a simulation of water contact angle measurement on top of graphite surface. Initially a water cubic box is placed on two-layer graphite surface with the rest of the box being vacuum. The water droplet is relaxed during the simulation to develop a spherical shape. An error of "X particles communicated to PME node Y are more than a cell length out of the domain decomposition cell of their charge group" was encountered. And I have read the suggested solutions at the link below http://www.gromacs.org/Documentation/Errors#X_particles_communicated_to_PME_node_Y_are_more_than_a_cell_length_out_of_the_domain_decomposition_cell_of_their_charge_group. I guess the reason for this error in my case is because of the vacuum such that the water molecules at the boundary of the droplet can move fast. I have check the trajectory and the simulation is OK. I doubt the simulation is OK. This error message is one of several that can happen when the system is not well-enough conditioned for the MD to be stable. See www.gromacs.org/Documentation/Terminology/Blowing_Up. Here, you have atoms moving much faster than GROMACS was engineered to expect. You should be confident that a water drop in a vacuum, and your graphite surface are both stable on their own before you try the wetting simulation. For this situation, is there a way of suppressing this error? Or what else can I do? Work out why it's poorly conditioned. Mark -- gmx-users mailing listgmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/Search before posting! Please don't post (un)subscribe requests to the list. Use the www interface or send it to gmx-users-requ...@gromacs.org. Can't post? Read http://www.gromacs.org/Support/Mailing_Lists
[gmx-users] How to suppress the error "X particles communicated to PME node Y are more than a cell length out of the domain decomposition cell of their charge group"
Dear GMXers, I'm running a simulation of water contact angle measurement on top of graphite surface. Initially a water cubic box is placed on two-layer graphite surface with the rest of the box being vacuum. The water droplet is relaxed during the simulation to develop a spherical shape. An error of "X particles communicated to PME node Y are more than a cell length out of the domain decomposition cell of their charge group" was encountered. And I have read the suggested solutions at the link below http://www.gromacs.org/Documentation/Errors#X_particles_communicated_to_PME_node_Y_are_more_than_a_cell_length_out_of_the_domain_decomposition_cell_of_their_charge_group . I guess the reason for this error in my case is because of the vacuum such that the water molecules at the boundary of the droplet can move fast. I have check the trajectory and the simulation is OK. For this situation, is there a way of suppressing this error? Or what else can I do? PS: the GROMACS version I'm using is GROMACS4.5. Thank you. Best, Yanbin -- gmx-users mailing listgmx-users@gromacs.org http://lists.gromacs.org/mailman/listinfo/gmx-users Please search the archive at http://www.gromacs.org/Support/Mailing_Lists/Search before posting! Please don't post (un)subscribe requests to the list. Use the www interface or send it to gmx-users-requ...@gromacs.org. Can't post? Read http://www.gromacs.org/Support/Mailing_Lists