On 16/03/11, "C.Y. Chang" <chiayun.ch...@gmail.com> wrote: > Hi, > > I try to calculate the Hbond in the ligand and remove the duplicate atoms in > the .ndx file dumped from g_hbond. > This is the contents in the .ndx file. > > [ HBD ] > 17385 17382 17422 17419 17459 17456 17496 17493 17533 17530 > > [ HBA ] > 17367 17375 17381 17384 17387 17388 17400 17404 17412 17418 17421 > 17424 17425 17437 17441 > 17449 17455 17458 17461 17462 17474 17478 17486 17492 17495 17498 > 17499 17511 17515 17523 > > 17529 17532 17535 17536 17548 > > And then, I perfrom the grompp. > It shows the error msg. > > Fatal error: > atoms 17366 and 17367 in charge group 1 of molecule type 'LIG' are in > different energy groups > > > How should I deal with the problem? >
Redefine charge or energy groups as appropriate for whatever you're trying to do. Non-bonded interactions are calculated by looping over atom-atom interactions grouped by charge group and then by energy group. Every charge group must be a subset of an energy group or fully disjoint from it. > On the other side, I read the information in the web page > (http://www.mail-archive.com/gmx-users@gromacs.org/msg34617.html) > > I remove the lines in the [bonds] region in .top files. > But I still get these terms > > 1 Angle 2 G96Angle 3 Proper-Dih. 4 > Ryckaert-Bell. > 5 Improper-Dih. 6 LJ-14 7 Coulomb-14 8 LJ-(SR) > > 9 Disper.-corr. 10 Coulomb-(SR) 11 Coul.-recip. 12 > Position-Rest. > 13 Potential 14 Kinetic-En. 15 Total-Energy 16 Temperature > 17 Pres.-DC 18 Pressure 19 Constr.-rmsd 20 Box-X > > 21 Box-Y 22 Box-Z 23 Volume 24 Density > 25 pV 26 Enthalpy 27 Vir-XX 28 Vir-XY > 29 Vir-XZ 30 Vir-YX 31 Vir-YY 32 Vir-YZ > > 33 Vir-ZX 34 Vir-ZY 35 Vir-ZZ 36 Pres-XX > 37 Pres-XY 38 Pres-XZ 39 Pres-YX 40 Pres-YY > 41 Pres-YZ 42 Pres-ZX 43 Pres-ZY 44 Pres-ZZ > > 45 #Surf*SurfTen 46 Mu-X 47 Mu-Y 48 Mu-Z > 49 Coul-SR:LIG-LIG 50 LJ-SR:LIG-LIG > 51 Coul-14:LIG-LIG 52 LJ-14:LIG-LIG > 53 Coul-SR:LIG-DPPC_SOL 54 LJ-SR:LIG-DPPC_SOL > > 55 Coul-14:LIG-DPPC_SOL 56 LJ-14:LIG-DPPC_SOL > 57 Coul-SR:DPPC_SOL-DPPC_SOL 58 LJ-SR:DPPC_SOL-DPPC_SOL > 59 Coul-14:DPPC_SOL-DPPC_SOL 60 LJ-14:DPPC_SOL-DPPC_SOL > 61 T-DPPC 62 T-SOL 63 T-LIG 64 Lamb-DPPC > > 65 Lamb-SOL 66 Lamb-LIG > > And I try to remove lines in [bonds], [pairs], [angles] and [ dihedrals ], > but it is the same. > Where is the wrong process? And how could I get the angle and dihedral energy > in the intramolecules? > > Thanks for your any comments. > You're decomposing PME electrostatics. That produces garbage values for quantities that have little meaning, like it says in the thread I suggested you read the whole of. I'm going to stop helping :-) Mark > > Best, > > Chia-yun > > > 2011/3/14 Mark Abraham <mark.abra...@anu.edu.au> > > > > > > > > On 14/03/11, "C.Y. Chang" <chiayun.ch...@gmail.com> wrote: > > > > > > Hi, > > > > > > I try to calculate hydrogen bond (HB) energy. > > > The g_energy does not have this term. > > > And I find the g_hbond function in Gromacs. > > > But the HB energy calculation is not in g_hbond. > > > > > > > > > > > There's good reasons for this. How would you define the "HB energy" in > > terms of the kind of information accessible to MD simulations? > > > > > > > > > > Therfore, I also try to dump the .ndx file including the HB_donor, > > > HB_acceptor and HB_system from g_hbond, and perfrom the grompp > > > > > > But there is a error msg, > > > > > > Atom 17380 in multiple Energy Mon. groups > > > > > > > > > Look up energy groups in the manual - start of section 3.3. Your .mdp file > > is defining an illegal combination of atoms and energy monitor groups. > > > > > > > > > Another problem is about calculating the intramolecular energy e.g. > > > 1,4-nonbonded, van der waal, electrostatic etc. of the ligand in lipid > > > layer-ligand complex system. > > > > > > > > > I could set up the energy_grp and calculate energy between the ligand > > > group and the lipid layer group. > > > But I need the intramolecular energy in the groups. > > > How should I deal with these problems? > > > > > > > > > > > An inter-group energy doesn't mean anything much, so don't bother. Please > > read the whole of this thread > > http://lists.gromacs.org/pipermail/gmx-users/2010-November/055687.html > > > > > > Mark > > > > -- > > > > 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