Pressure coupling and system heating

From: DimitryASuplatov (genesup_at_gmail.com)
Date: Mon Apr 20 2009 - 03:38:36 CDT

Hello,

I want to heat the system while preserving cell dimensions constant
(to avoid bubbles in the water).
I use Langevin temperature and pressure coupling algorithms.

1/ To my understanding to achieve my goal I need to adjust
LangevinPistonTemp parameter
to meet the current heating temperature.

# Minimization
minimize 5000

#Heat
for { set TEMP $start_temperature } { $TEMP < $hold_temperature } {
incr TEMP 25 } {
        reassignTemp $TEMP
        langevinPistonTemp $TEMP
        langevinTemp $TEMP
        run 5000
}

reassignTemp $hold_temperature
langevinTemp $hold_temperature
langevinPistonTemp $TEMP
run 5000000

Is this correct?

2/ langevinPistonTemp may not be <1. That means that I can not start
from 0K, I should start from 1K instead. Am I correct?

Thank you for your time.
SDA

This archive was generated by hypermail 2.1.6 : Wed Feb 29 2012 - 15:50:47 CST