share/top/amber03.ff/tip4pew.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_tip4pew  1       SOL       OW       1       0        16.00000
   8   2   HW_tip4pew  1       SOL       HW1      1       0.52422   1.00800
   9   3   HW_tip4pew  1       SOL       HW2      1       0.52422   1.00800
  10   4   MW          1       SOL       MW       1      -1.04844   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 [ bonds ]
  20 ; i     j       funct   length  force.c.
  21 1       2       1       0.09572 502416.0 0.09572        502416.0
  22 1       3       1       0.09572 502416.0 0.09572        502416.0
  23
  24 [ angles ]
  25 ; i     j       k       funct   angle   force.c.
  26 2       1       3       1       104.52  628.02  104.52  628.02
  27
  28 #endif
  29
  30
  31 [ virtual_sites3 ]
  32 ; Vsite from                    funct   a               b
  33 4       1       2       3       1       0.106676721     0.106676721
  34
  35
  36 [ exclusions ]
  37 1       2       3       4
  38 2       1       3       4
  39 3       1       2       4
  40 4       1       2       3
  41
  42
  43 ; The position of the virtual site is computed as follows:
  44 ;
  45 ;               O
  46 ;
  47 ;               V
  48 ;
  49 ;       H               H
  50 ;
  51 ; Ewald tip4p:
  52 ; const = distance (OV) / [ cos (angle(VOH))    * distance (OH) ]
  53 ;         0.0125 nm     / [ cos (52.26 deg)     * 0.09572 nm    ]
  54 ;       then a = b = 0.5 * const = 0.106676721
  55 ;
  56 ; Vsite pos x4 = x1 + a*(x2-x1) + b*(x3-x1)