Re: Simulation Using Both NAMD and VMD

From: Axel Kohlmeyer (akohlmey_at_gmail.com)
Date: Fri Apr 20 2012 - 20:28:15 CDT

On Fri, Apr 20, 2012 at 9:19 PM, Ajith Rathnaweera Rajapaksha Mudalige
<arajapak_at_purdue.edu> wrote:
> Hi Axel,
>
> measure sasa has other capabilities...
>
> It can compute the points on the solvent accessible surface and this can be done for each individual atoms.
> Forces are calculated using the coordinates of these points. So when the forces are in action surface changes and therefore forces are to be updated.

ok. i wasn't looking a the argument list of the code only at the return value.

> Everything I need can be done in NAMD, except computing measure sasa.
>
> Any suggestions?

if this is the case, passing the data from NAMD to VMD
and back is EVEN MORE STUPID, if that is at all possible.

why on earth do you insist on this insanity?
the SASA algorithm is not very difficult.
running two executables in step and
communicating data between them is,
not to mention the insane overhead of
formatting numbers into text, communicating
them and then parsing it and converting
to Tcl and back.

axel.

>
> Ajith
>
> ----- Original Message -----
> From: "Axel Kohlmeyer" <akohlmey_at_gmail.com>
> To: "Ajith Rathnaweera Rajapaksha Mudalige" <arajapak_at_purdue.edu>
> Cc: "Aron Broom" <broomsday_at_gmail.com>, namd-l_at_ks.uiuc.edu
> Sent: Friday, April 20, 2012 8:42:49 PM
> Subject: Re: namd-l: Simulation Using Both NAMD and VMD
>
> On Fri, Apr 20, 2012 at 7:14 PM, Ajith Rathnaweera Rajapaksha Mudalige
> <arajapak_at_purdue.edu> wrote:
>> Hi
>>
>> Thanks Alex, Aron,
>>
>> Yes, practically I will want to update the forces after every 500 time steps or so.
>>
>> I know this will slow down the simulation but it is an essential step to run the simulation correctly.
>>
>> The forces I am applying to my system is depending on Solvent Accessible Surface (SAS) of the molecule. SAS can be determined using "measure SASA" command in VMD. This is not directly supported in NAMD (please correct me if I am wrong). Therefore, I need to run the simulation in NAMD load the current phase in VMD and determine SAS at that point, determine the direction of the forces, get back to NAMD and run for the next phase.
>
> wait...
>
> there is one more thing that doesn't add up. measure sasa gives
> your a single number: the area. how does this compute into a
> directional force?
>
> axel.
>
>>
>> Is this possible?
>>
>> I will highly appreciate if anyone can give me a little concrete example, if anyone has done this type of a work before.
>>
>> Regards,
>>
>> Ajith
>>
>>
>>
>>
>> ----- Original Message -----
>> From: "Aron Broom" <broomsday_at_gmail.com>
>> To: "Axel Kohlmeyer" <akohlmey_at_gmail.com>
>> Cc: "Ajith Rathnaweera Rajapaksha Mudalige" <arajapak_at_purdue.edu>, namd-l_at_ks.uiuc.edu
>> Sent: Friday, April 20, 2012 6:30:14 PM
>> Subject: Re: namd-l: Simulation Using Both NAMD and VMD
>>
>> Is it possible to use the tcl forces capability of NAMD to run the script you want?  As Axel points out, passing the information into VMD at every timestep would be crippling (unless your simulation takes an extremely long time per timestep already).  If you only need to update the forces every 1000 timesteps or something, that might work.
>>
>> ~Aron
>>
>>
>> On Fri, Apr 20, 2012 at 6:10 PM, Axel Kohlmeyer < akohlmey_at_gmail.com > wrote:
>>
>>
>>
>> On Fri, Apr 20, 2012 at 6:05 PM, Ajith Rathnaweera Rajapaksha Mudalige
>> < arajapak_at_purdue.edu > wrote:
>>> Hi Everybody,
>>>
>>> Is it possible to write a script for a NAMD simulation to access VMD at each time step through Interactive Molecular Dynamics mode (or any other possible way).
>>
>> technically yes. in practice no, since that would be highly inefficient.
>> better to implement what you need into NAMD directly.
>>
>> axel.
>>
>>
>>
>>> In my simulation I need to execute some VMD commands at each time step and define a new force according to the results output from VMD.
>>>
>>> I highly appreciate if someone could shed some light in this direction.
>>
>>
>>
>>
>>> Regards,
>>>
>>> Ajith
>>>
>>
>>
>>
>> --
>> Dr. Axel Kohlmeyer
>> akohlmey_at_gmail.com   http://goo.gl/1wk0
>>
>> College of Science and Technology
>> Temple University, Philadelphia PA, USA.
>>
>>
>>
>>
>> --
>> Aron Broom M.Sc
>> PhD Student
>> Department of Chemistry
>> University of Waterloo
>
>
>
> --
> Dr. Axel Kohlmeyer
> akohlmey_at_gmail.com  http://goo.gl/1wk0
>
> College of Science and Technology
> Temple University, Philadelphia PA, USA.

-- 
Dr. Axel Kohlmeyer
akohlmey_at_gmail.com  http://goo.gl/1wk0
College of Science and Technology
Temple University, Philadelphia PA, USA.

This archive was generated by hypermail 2.1.6 : Mon Dec 31 2012 - 23:21:27 CST