Previous Thread
Next Thread
Print Thread
Bilayer Center
#18597 07/22/08 05:25 PM
Joined: Jan 2005
Posts: 86
M
mzzt Offline OP
Forum Member
OP Offline
Forum Member
M
Joined: Jan 2005
Posts: 86
Dear Charmm users,

I did a simulation of a lipid bilayer by keeping the bilayer center at 0,0,0 while system setup. But while production run, since there was no constraint on the bilayer at all, so is the center of bilayer still 0,0,0?

I wish to calculate the average distance between a methyl carbon atom of DMPC (i.e. say C13) and bilayer center. Simply evaluating the Z coordinate of C13 does not sound appropriate!! Kindly help.

MzZt.

Re: Bilayer Center
mzzt #18598 07/22/08 05:38 PM
Joined: Sep 2003
Posts: 8,580
Likes: 11
rmv Online Content
Forum Member
Online Content
Forum Member
Joined: Sep 2003
Posts: 8,580
Likes: 11
I often use the symmetry of the time series for the Z coords of the P atoms for each leaflet as a quick gauge of drift through the course of a simulation. It's a quick and easy analysis with CORREL, indicated by the commands below--

:
define leafa sele type P .and. prop Z .gt. 0. end
define leafb sele type P .and. prop Z .lt. 0. end
correl ...
enter lfa atom Z sele leafa end
enter lfb atom Z sele leafb end
:


Moderated by  BRBrooks, lennart, rmv 

Link Copied to Clipboard
Powered by UBB.threads™ PHP Forum Software 7.7.4
(Release build 20200307)
Responsive Width:

PHP: 7.3.31-1~deb10u1 Page Time: 0.004s Queries: 18 (0.002s) Memory: 0.7178 MB (Peak: 0.7535 MB) Data Comp: Off Server Time: 2021-11-30 23:39:35 UTC
Valid HTML 5 and Valid CSS