changing solventDielectric in GBIS

From: Andres Colubri (andres.colubri_at_gmail.com)
Date: Mon Jun 13 2011 - 07:36:41 CDT

Hello,

I'm simulating a heating cycle using the new GBIS model in NAMD,
starting at T=275 K and ending at T=300K, with increments of 5 degrees.

Since the dielectric of water goes from 88 at 0C to 78 at 28C, the
simulation script should update the solventDielectric parameter
accordingly at each temperature step. So I have the following in
my config file :

# Heating
for {set i $temperature0} {$i <= $temperature1} {incr i 5} {
   set temperaturei $i
   langevinTemp $temperaturei
   langevinPistonTemp $temperaturei
   solventDielectric $dielectric($temperaturei)
}

with:

set temperature0 275
set temperature1 300

and

set dielectric(275) 87.31
set dielectric(280) 85.46
set dielectric(285) 83.61
set dielectric(290) 81.76
set dielectric(295) 79.91
set dielectric(300) 78.06

However, I get this error in the first iteration of the loop:

TCL: Setting parameter solventDielectric to 87.31
FATAL ERROR: Setting parameter solventDielectric from script failed!

Any help with this will be greatly appreciated.
Thanks!

This archive was generated by hypermail 2.1.6 : Mon Dec 31 2012 - 23:20:25 CST