Warren Gallin wrote:
You are not specifying the ions to be added using the -pname and -nname options 
with the genion command.

Perhaps that is a problem?


They are not necessary; the default names are used and should be correct for all force fields now that naming has been standardized.

Thus far I have seen no evidence that genion is not doing what it is supposed to. It reports finding a net -50 charge on the system. What happens then? Does the output coordinate file contain ions? The topology will not be modified because genion is not being told to (i.e. through the use of -p).

-Justin

Warren Gallin

On 2011-07-26, at 2:49 PM, Sara baretller wrote:

thanks but i want to neutralize the system , why genion -s file.tpr -conc 0.2 
-neutral -o file.gro -random  does not neutralize the system to 0 .


On Tue, Jul 26, 2011 at 3:48 PM, Justin A. Lemkul <jalem...@vt.edu> wrote:


Sara baretller wrote:
Hi

this a part of the md.log  where the system has -50 charge

 two-body bonded interactions: 0.407 nm, Bond, atoms 514 515
 multi-body bonded interactions: 0.731 nm, G96Angle, atoms 1272 1275
Minimum cell size due to bonded interactions: 0.804 nm
Maximum distance for 9 constraints, at 120 deg. angles, all-trans: 0.810 nm
Estimated maximum distance required for P-LINCS: 0.810 nm
This distance will limit the DD cell size, you can override this with -rcon
Optimizing the DD grid for 4 cells with a minimum initial size of 0.810 nm
The maximum allowed number of cells is: X 11 Y 11 Z 11
Domain decomposition grid 4 x 1 x 1, separate PME nodes 0
Domain decomposition nodeid 0, coordinates 0 0 0

Table routines are used for coulomb: TRUE
Table routines are used for vdw:     TRUE
Using shifted Lennard-Jones, switch between 0.9 and 1.2 nm
Cut-off's:   NS: 1.2   Coulomb: 1.2   LJ: 1.2
System total charge: -50.000
Generated table with 1100 data points for Shift.
Tabscale = 500 points/nm
Generated table with 1100 data points for LJ6Shift.
Tabscale = 500 points/nm
Generated table with 1100 data points for LJ12Shift.
Tabscale = 500 points/nm
Configuring nonbonded kernels...
Configuring standard C nonbonded kernels...
Testing ia32 SSE2 support... present


There's nothing abnormal here.  genion reports that it finds a -50 charge on 
the system.  That's what it's supposed to do.  Based on what it finds in the 
topology, it adds neutralizing ions and any additional ions to reach the 
specified concentration.  You only have a problem if grompp later complains 
that there's still some residual charge (excluding tiny rounding discrepancies).

-Justin

On Tue, Jul 26, 2011 at 2:18 PM, Justin A. Lemkul <jalem...@vt.edu 
<mailto:jalem...@vt.edu>> wrote:



   Sara baretller wrote:

       Hi All

       I used the genion command using like this "genion -s file.tpr
       -conc 0.2 -neutral -o file.gro -random "  again and i checked
       the md.log file and it says that the net charge is negative like
       it was before using the genion command.

       so can anybody tell me what is wrong with the line  genion -s
       file.tpr -conc 0.2 -neutral -o file.gro -random .


   Please copy and paste the exact (pertinent) output from the log file.

   -Justin

       Sara




       On Tue, Jul 26, 2011 at 1:05 PM, Justin A. Lemkul
       <jalem...@vt.edu <mailto:jalem...@vt.edu>
       <mailto:jalem...@vt.edu <mailto:jalem...@vt.edu>>> wrote:



          Sara baretller wrote:

              Hi all

              I used the genion to add a concentration and to neutalize the
              system in the same time by using the

                   *genion -s file.tpr -conc 0.2 -neutral -o file.gro
       -random
                   so it did add the NA and Cl but it did not
       neutralize the
              system, the net charge of the system still the same negative.


          I find this hard to believe.  The application of genion -conc
          -neutral has worked in every instance I've tried it along every
          Gromacs version.  Check your output again.


                   so i tried to use the genion seperate to neutralize the
              charge using the file.tpr and out file as the file.gro,           
          however it reset the file.gro and i lose the NA And
       CL ions
              added for the concentration.


          Without the actual sequence of commands, this is not a useful
          description. genion adds ions based on whatever it finds in
       the .tpr
          file (which is presumably not neutralized).  If you do not
       re-create
          a .tpr file in between different additions of ions, you're
       going to
          be undoing work you thought you did.

          -Justin

          --     ==============================____==========

          Justin A. Lemkul
          Ph.D. Candidate
          ICTAS Doctoral Scholar
          MILES-IGERT Trainee
          Department of Biochemistry
          Virginia Tech
          Blacksburg, VA
          jalemkul[at]vt.edu <http://vt.edu> <http://vt.edu> | (540)

       231-9080 <tel:%28540%29%20231-9080>
          <tel:%28540%29%20231-9080>

          http://www.bevanlab.biochem.____vt.edu/Pages/Personal/justin
       <http://vt.edu/Pages/Personal/justin>

          <http://www.bevanlab.biochem.__vt.edu/Pages/Personal/justin
       <http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin>>

          ==============================____==========
          --     gmx-users mailing list    gmx-users@gromacs.org
       <mailto:gmx-users@gromacs.org>
          <mailto:gmx-users@gromacs.org <mailto:gmx-users@gromacs.org>>


          http://lists.gromacs.org/____mailman/listinfo/gmx-users
       <http://lists.gromacs.org/__mailman/listinfo/gmx-users>
          <http://lists.gromacs.org/__mailman/listinfo/gmx-users
       <http://lists.gromacs.org/mailman/listinfo/gmx-users>>
          Please search the archive at
          http://www.gromacs.org/____Support/Mailing_Lists/Search
       <http://www.gromacs.org/__Support/Mailing_Lists/Search>
          <http://www.gromacs.org/__Support/Mailing_Lists/Search
       <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
       <mailto:gmx-users-requ...@gromacs.org>
          <mailto:gmx-users-request@__gromacs.org
       <mailto:gmx-users-requ...@gromacs.org>>.

          Can't post? Read
       http://www.gromacs.org/____Support/Mailing_Lists
       <http://www.gromacs.org/__Support/Mailing_Lists>
          <http://www.gromacs.org/__Support/Mailing_Lists
       <http://www.gromacs.org/Support/Mailing_Lists>>



   --     ==============================__==========

   Justin A. Lemkul
   Ph.D. Candidate
   ICTAS Doctoral Scholar
   MILES-IGERT Trainee
   Department of Biochemistry
   Virginia Tech
   Blacksburg, VA
   jalemkul[at]vt.edu <http://vt.edu> | (540) 231-9080
   <tel:%28540%29%20231-9080>
   http://www.bevanlab.biochem.__vt.edu/Pages/Personal/justin
   <http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin>

   ==============================__==========
   --     gmx-users mailing list    gmx-users@gromacs.org
   <mailto:gmx-users@gromacs.org>
   http://lists.gromacs.org/__mailman/listinfo/gmx-users
   <http://lists.gromacs.org/mailman/listinfo/gmx-users>
   Please search the archive at
   http://www.gromacs.org/__Support/Mailing_Lists/Search
   <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
   <mailto:gmx-users-requ...@gromacs.org>.
   Can't post? Read http://www.gromacs.org/__Support/Mailing_Lists
   <http://www.gromacs.org/Support/Mailing_Lists>



--
========================================

Justin A. Lemkul
Ph.D. Candidate
ICTAS Doctoral Scholar
MILES-IGERT Trainee
Department of Biochemistry
Virginia Tech
Blacksburg, VA
jalemkul[at]vt.edu | (540) 231-9080
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin

========================================
--
gmx-users mailing list    gmx-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 mailing list    gmx-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


--
========================================

Justin A. Lemkul
Ph.D. Candidate
ICTAS Doctoral Scholar
MILES-IGERT Trainee
Department of Biochemistry
Virginia Tech
Blacksburg, VA
jalemkul[at]vt.edu | (540) 231-9080
http://www.bevanlab.biochem.vt.edu/Pages/Personal/justin

========================================
--
gmx-users mailing list    gmx-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

Reply via email to