fep
Differences
This shows you the differences between two versions of the page.
Both sides previous revisionPrevious revisionNext revision | Previous revision | ||
fep [2018/01/16 12:20] – edit | fep [2020/12/10 08:16] (current) – admin | ||
---|---|---|---|
Line 1: | Line 1: | ||
**Free Energy Perturbation Calculation in CHARMM** | **Free Energy Perturbation Calculation in CHARMM** | ||
- | Link to FEP files in DT2: /a/fs-3/ | + | |
+ | Link to FEP files in DT2: | ||
+ | /lustre/jbklauda/FEP/FEP_example.zip | ||
To run an FEP calculation, | To run an FEP calculation, | ||
- A coordinate (.crd) file for your molecule in a vacuum, equilibrated. Call it “moleculename_gas.crd, | - 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, | - 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: | 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 " | - Check for proper permission to all .run files. Navigate to the appropriate folder in DT2. Then, you can use the " | ||
Line 13: | Line 17: | ||
- 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. | - 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, | - 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). | 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: | 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. | - 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: | + | - 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: | + | - 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/ | - 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.1516123218.txt.gz · Last modified: 2018/01/16 12:20 by edit