fep
Differences
This shows you the differences between two versions of the page.
Next revision | Previous revision | ||
fep [2017/11/26 20:23] – created admin | fep [2020/12/10 08:16] (current) – admin | ||
---|---|---|---|
Line 1: | Line 1: | ||
- | test | + | **Free Energy Perturbation Calculation in CHARMM** |
+ | |||
+ | Link to FEP files in DT2: | ||
+ | / | ||
+ | |||
+ | To run an FEP calculation, | ||
+ | - A coordinate (.crd) file for your molecule in a vacuum, equilibrated. Call it “moleculename_gas.crd, | ||
+ | - A coordinate file of your molecule in a waterbox that is an appropriate size for the molecule. Call it “moleculename_wbox.crd.” It’s important that this system is NOT equilibrated, | ||
+ | |||
+ | Copy the files and follow these steps: | ||
+ | - Check for proper permission to all .run files. Navigate to the appropriate folder in DT2. Then, you can use the " | ||
+ | - Update the path to the CHARMM executable file in the following files: | ||
+ | - / | ||
+ | - Edit the files / | ||
+ | - Move your moleculename_gas.crd file into the /vac folder and the /water folder, and your moleculename_wbox.crd file into the /water folder. | ||
+ | - Submit the jobs by running “fep.run” in two separate folders, once in the /vac folder and once in the /water folder. Let all the jobs complete. THEN tally up the totals using “xtract.run” separately in each folder. It will create a “job.1” folder where it will store all the outputs. BE SURE TO BACK UP THE FILES before doing this, as they will be deleted. (Alternately, | ||
+ | - To compute FE, you use the “total” under “final results.” Simply subtract initial from final. That is, subtract the total in the gas phase FROM the total in the water phase. Units of energy are kcal/mol. | ||
+ | - Lastly, for multiple runs, be sure to change “seed” in the file “fep.run, | ||
+ | |||
+ | The basic idea is that you use fep.run to submit a series of simulations of your gas and water systems separately. The file xtract.run does a WHAM analysis (see CHARMM documentation) and a thermodynamic integration analysis of the results. You subtract to find the change in FE from the gas state to the solvated state (water result – vacuum result). | ||
+ | |||
+ | On building an appropriate water box: | ||
+ | - You can use step4.2_waterbox.inp in the /water folder. Edit parameters “A,” “B,” “C,” and “watboxZ” in the file step3_size.str. Make all of these the same length, the length of your box in angstroms. Make sure your box is big enough to allow for a 7-or-so-angstrom buffer on all sides of your molecule, which should be placed centered at zero. | ||
+ | - In the command line of DT2, you can just use CHARMM to read step4.2_waterbox.inp. Assuming the path to the CHARMM executable is " | ||
+ | - Once you have the resulting files, equilibrate the water box using step6.1_equilibration.inp: | ||
+ | - Place your molecule into the water box using step5_assembly.inp. You’ll need to change the name of your gas molecule. Mine was DMOE, so you can just use edit/ |
fep.1511745794.txt.gz · Last modified: 2017/11/26 20:23 by admin