[gmx-users] set charge as 0 in gromacs code
Hello everyone, My system is 512 waters and 1 ion CL-. Now I would like to change the CL charge as 0 in the gromacs code for the short non-bonded interaction. While for the long-range electrostatics (pme) the CL is -1. My gromacs version is 4.6.2, which automatically starts 8 threads to compute. In src/kernel/md.c, at line 1191 just before function do_force() , I added codes as following which could find out which atom is CL and change its charge as 0 in a certain thread. After this, when I use g_energy, the Coul.SR of SOL and CL is -88 kJ, which should be 0 kJ. Does anyone could help me out? int stone_a; for (stone_a=mdatoms-start; stone_amdatoms-start+mdatoms-homenr; stone_a++) { if (mdatoms-chargeA[stone_a] == -1.0) { fprintf(stderr,the atom %d found with charge %f\n, stone_a,mdatoms-chargeA[stone_a]); mdatoms-chargeA[stone_a]=0.0; break; } } Thanks, Stony -- 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] set charge as 0 in gromacs code
Hello Mark, Thanks for your reply!!! Not self energy term. For my research reason, I would like to turn off the long-ranged electrostatics between SOL and CL. That is to say, SOL and CL only interacts with erfc(eta*r)/r potential, without reciprocal electrostatics parts. I use PME in the mdp file for coulomb, and in topol file the CL is -1 charged. Hence I have to set CL charge as 0 in the Gromacs code to turn off the reciprocal calculation of the pme scheme. E.g before gmx_pme_do(), I set Cl charge as 0, but this fails to work. So do you some idea to do this? Stony From: gmx-users-boun...@gromacs.org [gmx-users-boun...@gromacs.org] on behalf of Mark Abraham [mark.j.abra...@gmail.com] Sent: Tuesday, August 06, 2013 8:20 PM To: Discussion list for GROMACS users Subject: Re: [gmx-users] set charge as 0 in gromacs code Ewald self energy, perhaps? Mark On Wed, Aug 7, 2013 at 1:16 AM, Shi, Yu (shiy4) sh...@mail.uc.edu wrote: Hello everyone, My system is 512 waters and 1 ion CL-. Now I would like to change the CL charge as 0 in the gromacs code for the short non-bonded interaction. While for the long-range electrostatics (pme) the CL is -1. My gromacs version is 4.6.2, which automatically starts 8 threads to compute. In src/kernel/md.c, at line 1191 just before function do_force() , I added codes as following which could find out which atom is CL and change its charge as 0 in a certain thread. After this, when I use g_energy, the Coul.SR of SOL and CL is -88 kJ, which should be 0 kJ. Does anyone could help me out? int stone_a; for (stone_a=mdatoms-start; stone_amdatoms-start+mdatoms-homenr; stone_a++) { if (mdatoms-chargeA[stone_a] == -1.0) { fprintf(stderr,the atom %d found with charge %f\n, stone_a,mdatoms-chargeA[stone_a]); mdatoms-chargeA[stone_a]=0.0; break; } } Thanks, Stony -- 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 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 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 set charge 0 in the source code for PME
Dear gmx-users, My system is composed of an ion(CL-) and waters. In the .top file the charge of CL is -1. PME-USER is used to deal with electrostatics. Now I would like set CL charge as 0 in the source code, in order to turn off the ELECTROSTATIC FORCE of Reciprocal part between CL and waters. Anyone has the idea where to change the code? Or, I could set CL charge as 0 in .top file, and PME-USER is still used. Then in the code, how to set CL as -1 in the code for ELECTROSTATIC FORCE of Real part? Yu -- 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] mdrun no error, but hangs no results
Dear gmx-users, My problem is weird. My mdrun worked well using the old serial version 4.5.5 (about two years ago). And I have these top, ndx, mdp, and gro files. Basing on those old files, for the serial 4.6.2, the grompp works through, resulting the .tpr file successfully. After that when I make the mdrun, mdrun -v -s em-nv.tpr -deffnm ss it only shows: Reading file em-nv.tpr, VERSION 4.6.2 (double precision) Using 8 MPI threads Killed and there is no further processing. Later, it is killed. Also, using cmake my installation process works well, so does anyone meet this problem before? And part of the logfile is: Log file opened on Wed Jul 17 12:17:30 2013 Host: opt-login03.osc.edu pid: 32177 nodeid: 0 nnodes: 1 Gromacs version:VERSION 4.6.2 Precision: double Memory model: 64 bit MPI library:thread_mpi OpenMP support: disabled GPU support:disabled invsqrt routine:gmx_software_invsqrt(x) CPU acceleration: SSE2 FFT library:fftw-3.3-sse2 Large file support: enabled RDTSCP usage: enabled Built on: Wed Jul 17 10:51:22 EDT 2013 Built by: ucn1...@opt-login03.osc.edu [CMAKE] Build OS/arch: Linux 2.6.18-308.11.1.el5 x86_64 Build CPU vendor: AuthenticAMD Build CPU brand:Dual-Core AMD Opteron(tm) Processor 8218 Build CPU family: 15 Model: 65 Stepping: 2 Build CPU features: apic clfsh cmov cx8 cx16 htt lahf_lm mmx msr pse rdtscp sse2 sse3 C compiler: /usr/bin/cc GNU cc (GCC) 4.1.2 20080704 (Red Hat 4.1.2-48) C compiler flags: -msse2-Wextra -Wno-missing-field-initializers -Wno-sign-compare -Wall -Wno-unused -Wunused-value -Wno-unknown-pragmas -fomit-frame-pointer -funroll-all-loops -O3 -DNDEBUG . . . . Initializing Domain Decomposition on 8 nodes Dynamic load balancing: no Will sort the charge groups at every domain (re)decomposition Using 0 separate PME nodes, as there are too few total nodes for efficient splitting Optimizing the DD grid for 8 cells with a minimum initial size of 0.000 nm Domain decomposition grid 8 x 1 x 1, separate PME nodes 0 PME domain decomposition: 8 x 1 x 1 Domain decomposition nodeid 0, coordinates 0 0 0 Using 8 MPI threads Detecting CPU-specific acceleration. Present hardware specification: Vendor: AuthenticAMD Brand: Dual-Core AMD Opteron(tm) Processor 8218 Family: 15 Model: 65 Stepping: 2 Features: apic clfsh cmov cx8 cx16 htt lahf_lm mmx msr pse rdtscp sse2 sse3 Acceleration most likely to fit this hardware: SSE2 Acceleration selected at GROMACS compile time: SSE2 Table routines are used for coulomb: FALSE Table routines are used for vdw: FALSE Will do PME sum in reciprocal space. PLEASE READ AND CITE THE FOLLOWING REFERENCE U. Essmann, L. Perera, M. L. Berkowitz, T. Darden, H. Lee and L. G. Pedersen A smooth particle mesh Ewald method J. Chem. Phys. 103 (1995) pp. 8577-8592 --- Thank You --- -- 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] Default Electrostatic Algorithm for User Defined Potential
Hello everyone, For a system of water (SOL) and ion (NA), if a speical potential needs to be applied only to SOL and NA, we have to make a user-define tabulation potental table_SOL_NA.xvg. And then the table.xvg is for the other groups of SOL SOL and NA NA. Here is the question, for groups of SOL SOL and NA NA, when table.xvg is adopted, the default electrostatic interaction (1/r) algorithm is the cut-off scheme. How could we use other schemes such as PME for that under periodic boundary condition? Stone -- 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