surface_area_of_lipid
Differences
This shows you the differences between two versions of the page.
Both sides previous revisionPrevious revisionNext revision | Previous revision | ||
surface_area_of_lipid [2021/05/25 14:13] – admin | surface_area_of_lipid [2023/02/05 10:39] (current) – admin | ||
---|---|---|---|
Line 9: | Line 9: | ||
This simple analysis is used to estimate the overall average surface area (SA) per lipid. If you are analyzing a single component membrane then this is equivalent to the component' | This simple analysis is used to estimate the overall average surface area (SA) per lipid. If you are analyzing a single component membrane then this is equivalent to the component' | ||
- | Those at UMD should use the DT2 path: **/lustre/ | + | Those at UMD should use the ZT1 path: **/afs/ |
Outside of UMD a sample script is given in the gzip file: {{ : | Outside of UMD a sample script is given in the gzip file: {{ : | ||
Line 82: | Line 82: | ||
The example states that you want to do the averages from the last 20 dyn*dcd files. Let's say you have 50 dyn*.dcd files, then the first 30 files are ignored. If each file had 2ns of data, then you are calculating the average from 60-100ns. | The example states that you want to do the averages from the last 20 dyn*dcd files. Let's say you have 50 dyn*.dcd files, then the first 30 files are ignored. If each file had 2ns of data, then you are calculating the average from 60-100ns. | ||
- | 3.2 Calculate the block-average surface areas by command: '' | + | 3.2 Calculate the block-average surface areas by command: '' |
**Note: **For //#define MAX_R 1000000// | **Note: **For //#define MAX_R 1000000// | ||
Line 88: | Line 88: | ||
As you can see, the output average area data are saved in // | As you can see, the output average area data are saved in // | ||
- | 4. Calculate averages, sample standard deviation, and standard error of block averages by command: '' | + | 3.3 Standard Errors from Block Average |
+ | |||
+ | Calculate averages, sample standard deviation, and standard error of block averages by command: '' | ||
The file // | The file // | ||
- | |||
- | ==== Get the equilibrium data (test) ==== | ||
- | |||
- | (author: Yalun Yu, email: alanyu17@terpmail.umd.edu / yalun.research@gmail.com) | ||
- | |||
- | < Email me if anything is wrong > | ||
- | |||
- | This is a brief instruction to obtain the time range of equilibrium (and also the minimum block size can be used for standard error calculation). Based on these, calculation of SA/lipid and Ka is done automatically for your lipid only system. Run Xiao's // | ||
- | |||
- | 1. Install Anaconda (https:// | ||
- | a) Download the python3.7 version for Linux [I used the ' | ||
- | and put the .sh installer in your home directory (or any other place you can access and execute it) on DT2. | ||
- | b) make the .sh file an executable (chmod u+x $NameOfInstaller.sh) | ||
- | c) Install by running " | ||
- | |||
- | 2. Make sure the Anaconda installation path is included (see below) in your ~/ | ||
- | |||
- | In ~/ | ||
- | set path = ($HOME/ | ||
- | or in ~/ | ||
- | | ||
- | |||
- | 3. create a conda environment and activate the env. (see conda' | ||
- | |||
- | 4. Install the pymbar package. | ||
- | conda install -c omnia pymbar | ||
- | |||
- | Download the script: {{ : | ||
- | |||
- | Change the temperature and number of lipids per leaflet in // | ||
- | set temp=323.15 | ||
- | set nlip=36 | ||
- | |||
- | Change the initial guess (or your preference) for block size (in ns) and the steps size (in ns) of your output (usually 1ps, but make sure to match your xst freq): | ||
- | area_handling.py 10 0.001 >& all.dat | ||
- | |||
- | After running // | ||
- | |||
surface_area_of_lipid.1621966406.txt.gz · Last modified: 2021/05/25 14:13 by admin