Re: restart file - langevin options removed?

From: Maxim Belkin (mbelkin_at_ks.uiuc.edu)
Date: Thu Apr 28 2016 - 15:17:35 CDT

Both errors explicitly say that variable temperature (used as $temperature) is not set. Add the following line at the beginning of the file:

set temperature 310

Also, keep in mind that "temperature $temperature" has nothing to do with thermostats. It only assigns random velocities to all the atoms so that the total kinetic energy corresponds to desired temperature.

> On Apr 28, 2016, at 15:13, Olya Kravchenko <ovkrav_at_gmail.com> wrote:
>
> Thank you!
>
> I don't quite understand then - if I comment out only the temperature
> option (i.e. the line that sets certain temperature, f.e. 310) then
> namd complains about missing variable:
>
> Processor 0 Exiting: Called CmiAbort ------------
> Reason: FATAL ERROR: can't read "temperature": no such variable
> while executing
> "temperature $temperature"
> (file "restart.conf" line 40)
>
> And when I comment out that line (i.e. temperature $temperature) it
> complains about
> langevin temperature:
>
> Processor 0 Exiting: Called CmiAbort ------------
> Reason: FATAL ERROR: can't read "temperature": no such variable
> while executing
> "langevinTemp $temperature"
> (file "restart.conf" line 68)
>
> And when I comment out that one, it says that langevin requires
> temperature to be specified. And when I comment out the whole routine
> it runs with no complaints...
>
>
> On Thu, Apr 28, 2016 at 4:00 PM, Maxim Belkin <mbelkin_at_ks.uiuc.edu> wrote:
>> Hi Olga,
>>
>> You should have the same temperature control mechanism throughout the production simulations. So, if you maintain constant temperature using Langevin thermostat, you have to have the langevin thermostat lines in all your NAMD configure files.
>>
>> langevinPiston part should be there as long as you want to run NPT (constant P - pressure) rather than NVP (constant V - volume).
>>
>> Also, make sure you use "commotion yes", otherwise NAMD will subtract the average velocity of the system when you restart your simulations.
>>
>> Maxim
>>
>>
>>
>>
>>> On Apr 28, 2016, at 14:49, Olya Kravchenko <ovkrav_at_gmail.com> wrote:
>>>
>>> Hi all,
>>>
>>> just a quick question, do I understand it correctly that langevin
>>> section should not be in a restart configuration file at all? I.e. to
>>> restart my simulation I should remove/comment out this whole section
>>> (as opposed to only removing the langevinTemp $temperature line):
>>>
>>> # Constant Temperature Control
>>> langevin on
>>> langevinDamping 1
>>> langevinTemp $temperature
>>> langevinHydrogen no
>>>
>>> What about langevinPistonTemp $temperature, should that part be
>>> commented out as well?
>>>
>>> Thanks in advance!
>>>
>>> Olga
>>>
>>> P.S. Does this restart script look right?
>>>
>>> structure ionized.psf
>>> coordinates ionized.pdb
>>> outputName test_1
>>>
>>> if {1} {
>>> set inputname test
>>> binCoordinates $inputname.restart.coor
>>> binVelocities $inputname.restart.vel ;# remove the "temperature"
>>> entry if you use this!
>>> extendedSystem $inputname.xsc
>>> }
>>>
>>> firsttimestep 26500
>>>
>>> paraTypeCharmm on
>>> parameters my_parameters.prm
>>>
>>> wrapWater on #what about this option, is it needed?
>>> wrapAll on
>>>
>>> exclude scaled1-4
>>> 1-4scaling 1.0
>>> cutoff 12.0
>>> switching on
>>> switchdist 10.0
>>> pairlistdist 14.0
>>>
>>>
>>>
>>> timestep 2.0 ;# 2fs/step
>>> rigidBonds all
>>> nonbondedFreq 1
>>> fullElectFrequency 2
>>> stepspercycle 10
>>>
>>>
>>> #PME
>>> if {0} {
>>> PME yes
>>> PMEGridSpacing 1.0
>>> }
>>>
>>> if {0} {
>>> useGroupPressure yes ;
>>> useFlexibleCell no ;
>>> useConstantArea no
>>>
>>> langevinPiston on
>>> langevinPistonTarget 1.01325 ;
>>> langevinPistonPeriod 100.0
>>> langevinPistonDecay 50.0
>>> langevinPistonTemp $temperature
>>> }
>>>
>>> restartfreq 500 ;# 500steps = every 1ps
>>> dcdfreq 250
>>> xstFreq 250
>>> outputEnergies 100
>>> outputPressure 100
>>>
>>>
>>> run 10000000
>>

This archive was generated by hypermail 2.1.6 : Tue Dec 27 2016 - 23:22:05 CST