Is it correct to convert from a steady state WE to equilibrium WE if no segment reached the target state?

32 views
Skip to first unread message

Ibrahim Mohamed

unread,
Nov 28, 2023, 1:22:32 PM11/28/23
to westpa-users
Greetings,

I was running a steady state WE simulation and then thought that i need to do an equilibrium WE. Is it correct to change from steady state to equilibrium if no segment reached the target state?
if it is correct, the only change will be in the command of w_init by removing the --tstate-file from the command, right?

Thanks

Daniel Zuckerman

unread,
Dec 1, 2023, 8:14:58 PM12/1/23
to westpa-users
Ibrahim, I don't know about the implementation, but you're right that if no segment has reached the pre-set target, then all those existing segments should be perfectly valid for 'equilibrium' simulation.  As you probably know, this just means that the system will relax toward equilibrium, not that the data can be assumed representative of equilibrium without careful testing.

Perhaps someone else can chime in about the correct implementation for you.

--Dan Zuckerman


From: westpa...@googlegroups.com <westpa...@googlegroups.com> on behalf of Ibrahim Mohamed <hollowic...@gmail.com>
Sent: Tuesday, November 28, 2023 10:12 AM
To: westpa-users <westpa...@googlegroups.com>
Subject: [EXTERNAL] [westpa-users] Is it correct to convert from a steady state WE to equilibrium WE if no segment reached the target state?
 
--
You received this message because you are subscribed to the Google Groups "westpa-users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to westpa-users...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/westpa-users/bee058b7-0a04-4b0b-a3a5-f53e4826d23en%40googlegroups.com.

Yang, Darian T

unread,
Dec 3, 2023, 2:41:31 PM12/3/23
to westpa...@googlegroups.com

Ibrahim,

 

You can use w_init if you’re starting a new simulation, thus generating a new h5 file and preparing the first iteration using the specified bstates w/ or w/o tstates. If you want to continue your current simulation with an adjustment to the tstates, this should be possible using w_states. So, you could try replacing the current target state with no target state, which should then update upon the next iteration.

 

Best,

Darian

Reply all
Reply to author
Forward
0 new messages