share/top/amber96.ff/tip4p.itp

   1 [ moleculetype ]
   2 ; molname       nrexcl
   3 SOL             2
   4
   5 [ atoms ]
   6 ; id  at type     res nr  res name  at name  cg nr  charge    mass
   7   1   OW_tip4p    1       SOL       OW       1       0        16.00000
   8   2   HW_tip4p    1       SOL       HW1      1       0.52      1.00800
   9   3   HW_tip4p    1       SOL       HW2      1       0.52      1.00800
  10   4   MW          1       SOL       MW       1      -1.04      0.00000
  11
  12 #ifndef FLEXIBLE
  13
  14 [ settles ]
  15 ; i     funct   doh     dhh
  16 1       1       0.09572 0.15139
  17
  18 #else
  19
  20 [ bonds ]
  21 ; i     j       funct   length  force.c.
  22 1       2       1       0.09572 502416.0 0.09572        502416.0
  23 1       3       1       0.09572 502416.0 0.09572        502416.0
  24
  25 [ angles ]
  26 ; i     j       k       funct   angle   force.c.
  27 2       1       3       1       104.52  628.02  104.52  628.02
  28
  29 #endif
  30
  31
  32 [ virtual_sites3 ]
  33 ; Vsite from                    funct   a               b
  34 4       1       2       3       1       0.128012065     0.128012065
  35
  36
  37 [ exclusions ]
  38 1       2       3       4
  39 2       1       3       4
  40 3       1       2       4
  41 4       1       2       3
  42
  43
  44 ; The position of the virtual site is computed as follows:
  45 ;
  46 ;               O
  47 ;
  48 ;               V
  49 ;
  50 ;       H               H
  51 ;
  52 ; const = distance (OV) / [ cos (angle(VOH))    * distance (OH) ]
  53 ;         0.015 nm      / [ cos (52.26 deg)     * 0.09572 nm    ]
  54 ;
  55 ; Vsite pos x4 = x1 + a*(x2-x1) + b*(x3-x1)