share/top/oplsaa.ff/tip4pew.itp

   1
   2 ; Horn et al. (2004). J. Chem. Phys.120, 9665-9678
   3 ; http://dx.doi.org/10.1063/1.1683075
   4
   5 [ moleculetype ]
   6 ; molname       nrexcl
   7 SOL             2
   8
   9 [ atoms ]
  10 ; id  at type     res nr  res name  at name  cg nr  charge    mass
  11   1   OW_tip4pew  1       SOL       OW       1       0        16.00000
  12   2   HW_tip4pew  1       SOL       HW1      1       0.52422   1.00800
  13   3   HW_tip4pew  1       SOL       HW2      1       0.52422   1.00800
  14   4   MW          1       SOL       MW       1      -1.04844   0.00000
  15
  16 #ifndef FLEXIBLE
  17
  18 [ settles ]
  19 ; i     funct   doh     dhh
  20 1       1       0.09572 0.15139
  21
  22 #else
  23 [ bonds ]
  24 ; i     j       funct   length  force.c.
  25 1       2       1       0.09572 502416.0 0.09572        502416.0
  26 1       3       1       0.09572 502416.0 0.09572        502416.0
  27
  28 [ angles ]
  29 ; i     j       k       funct   angle   force.c.
  30 2       1       3       1       104.52  628.02  104.52  628.02
  31
  32 #endif
  33
  34
  35 [ virtual_sites3 ]
  36 ; Vsite from                    funct   a               b
  37 4       1       2       3       1       0.106676721     0.106676721
  38
  39
  40 [ exclusions ]
  41 1       2       3       4
  42 2       1       3       4
  43 3       1       2       4
  44 4       1       2       3
  45
  46
  47 ; The position of the virtual site is computed as follows:
  48 ;
  49 ;               O
  50 ;
  51 ;               V
  52 ;
  53 ;       H               H
  54 ;
  55 ; Ewald tip4p:
  56 ; const = distance (OV) / [ cos (angle(VOH))    * distance (OH) ]
  57 ;         0.0125 nm     / [ cos (52.26 deg)     * 0.09572 nm    ]
  58 ;       then a = b = 0.5 * const = 0.106676721
  59 ;
  60 ; Vsite pos x4 = x1 + a*(x2-x1) + b*(x3-x1)