share/top/oplsaa.ff/tip5pe.itp

   1 [ moleculetype ]
   2 ; molname       nrexcl
   3 SOL             2
   4
   5 [ atoms ]
   6 ; id    at type res nr  residu name     at name         cg nr   charge
   7 1       opls_118e    1       SOL              OW             1       0
   8 2       opls_119     1       SOL             HW1             1       0.241
   9 3       opls_119     1       SOL             HW2             1       0.241
  10 4       opls_120     1       SOL             LP1             1      -0.241
  11 5       opls_120     1       SOL             LP2             1      -0.241
  12
  13 [ settles ]
  14 ; i     funct   doh     dhh
  15 1       1       0.09572 0.15139
  16
  17 [ virtual_sites3 ]
  18 ; The position of the virtual site is computed as follows:
  19 ;
  20 ; The distance from OW to L is 0.07 nm, the geometry is tetrahedral
  21 ; (109.47 deg)
  22 ; Therefore, a = b = 0.07 * cos (109.47/2) / | xOH1 + xOH2 |
  23 ;            c = 0.07 * sin (109.47/2) / | xOH1 X xOH2 |
  24 ;
  25 ; Using | xOH1 X xOH2 | = | xOH1 | | xOH2 | sin (H1-O-H2)
  26 ;       | xOH1 + xOH2 | = 2 | xOH1 | cos (H1-O-H2)
  27 ; Vsite pos x4 = x1 + a*x21 + b*x31 + c*(x21 X x31)
  28
  29 ; Vsite from                    funct   a       b               c
  30 4       1       2       3       4       -0.344908  -0.344908  -6.4437903493
  31 5       1       2       3       4       -0.344908  -0.344908   6.4437903493
  32
  33 [ exclusions ]
  34 1       2       3       4       5
  35 2       1       3       4       5
  36 3       1       2       4       5
  37 4       1       2       3       5
  38 5       1       2       3       4