Differences

This shows you the differences between two versions of the page.

--- edp [2026/08/25 14:58] – [Step 1: Edit Necessary Files] edit
+++ edp [2026/08/27 18:29] (current) – edit
@@ Line 324: / Line 324: @@
-b) Run ''1_atm_dens.scr''. The first time you run it, it will take a couple of minutes to set up a virtual environment with the MDAnalysis package installed. Then, it will automate all other steps and run for < 1 minute. You should now see the default plots, thickness data, and interdigitation data in the results folder. If you need NSLD data, all of the formatted .sim and .cmp files can be found inside the SimToEXP folder. Read more about SimToEXP in part V of the legacy script instructions.
+b) Run ''1_atm_dens.scr''. The first time you run it, it will take a couple of minutes to set up a virtual environment with the MDAnalysis package installed. Then, it will automate all other steps and run for < 1 minute. You should now see the default plots, thickness data, and interdigitation data in the results folder. If you need NSLD data, all of the formatted .sim and .cmp files can be found inside the SimToEXP folder. Read more about SimToEXP in part VI of the legacy script instructions.
 ===== Additional Information =====
 This section explains the Python workflow in a bit more detail. I recommend you at least skim over it.
 === Overall Pipeline ===
+_setup.py roughly validates the structure of all the fields in config.py before generating necessary files. 1_atm_dens.scr create a virtual Python environment, then runs 1_atm_dens.py to get .sim files with atom densities -> 2_calc_edp.py which applies the electron weights and groups atoms according to the .cmp files -> 3_analysis.py which computes thickness and interdigitation data -> 4_plot_edp.py which creates generic plots based on the edp data -> 5_exp_format.py which formats .cmp and .sim files to be usable in SimToEXP.
-  - **Edit ''config.py''** — resns list, charmm_dir, dyn_dir, dcd_range, bilayer_only, save_csv, EXP_format, charmm_version, job_time_limit
-  - **''0_setup.py''** — validates config, generates ''recenter_vars.str'', ''toppar_temp.str'', ''recenter.csh'', ''system_info.txt''; creates subdirs
-  - **''1_recenter.scr''** — submits ''recenter.csh'' (3-pass CHARMM recentering) -> recentered DCDs
-  - **''2_run_gen_rtf.scr''** — atom labels + group classification (via ''atmlbl.py''/''cmp.py'') -> ''atmlbl_cmp_files/''
-  - **''3_bldsim.py''** — per-atom, per-leaflet density from DCDs -> ''sim_data/''
-  - **''4_calc_edp.py''** — electron weighting + group EDP -> ''edp_data/''
-  - **''5_analysis.py''** — thickness (DHH/DB/2DC) + interdigitation -> ''results/EDP_results.txt''
-  - **''6_plot_edp.py''** — EDP plots -> ''results/*.png''
-  - **''7_exp_format.py''** — SimToEXP formatting, only if ''EXP_format = True'' -> ''SimToEXP/''
-=== Comparison with the Legacy Scripts ===
-Pros:
-The legacy scripts only worked for symmetric bilayers involving some glycero and sphingo lipids. It also used nested loop structures to get z-coordinate data for all atoms in the simulation, whereas this script directly greps the dcds, turning a 1-6 hour job into a 30-second job.
-Cons:
-The plotting scripts are already heavily customized. For those not comfortable with Python and Matplotlib, it may be easier to import *.dat files with edp data into Excel as CSV files and plot there, or use any alternative software. (Removed old generate Excel script because all *.dat files are already in CSV format, which is easily read into Excel).
 === How to Interpret Last 2 Columns in Thickness Results ======
@@ Line 355: / Line 340: @@
 === Formulas for Self and Hetero Interdigitation ===
-rho is the raw EDP, i denotes the lipid
+{{:new_interdig_formulas.png}}
-Self-Interdigitation for Symmetric Bilayers \\
-{{:self-interdigitation-formula.png}}
-Hetero-Interdigitation for Any Bilayer \\
-{{:hetero-interdigitation-formula.png}}