Skip to content

shirtsgroup/alchemical_metadynamics_project

 
 

Folders and files

NameName
Last commit message
Last commit date

Latest commit

 

History

17 Commits
 
 
 
 
 
 
 
 
 
 
 
 

Repository files navigation

alchemical_metadynamics

plumID:23.003 DOI

Description

This is a repository for maintaining the simulation files and analysis codes for the alchemical metadynamics project, which has been published in the Journal of Chemical Theory and Computation, titled Alchemical Metadynamics: Adding Alchemical Variables to Metadynamics to Enhance Sampling in Free energy calculations.

In the paper, we proposed adding alchemical variables to metadynamics to enhance sampling for free energy calculations. The method, termed as alchemical metadynamics, is demonstrated with three different test systems/alchemical processes, including decoupling an argon atom and a 4-site model from water, and the methylation of a nucleoside both in the isolated form and in a duplex. Note that

  • To keep the repository lightweight, we don't store simulation output files, including *.edr, *.log, *.trr, *.xtc, *cpt, *COLVAR*, *HILLS* and most of the *xvg files. We do keep some *txt files storing analysis results and some analysis outputs in the *png format.
  • Some preliminary tests of the project were stored in the archived project repository MetaD_EXE_TestSys.
  • The first version of this repository is the same as the commit 65cd812 of our previous project repository alchemical_MetaD_archived but without the previous git history the files of the old System 3 that have been migrated to MetaD_EXE_TestSys.

This repository includes the folders listed below. For more details about the files in each folder, please refer to the README file in each folder. Note that any COLVAR or HILLS files named in the form of HILLS*fake or COLVAR*fake in this repo are just placeholders used to pass the test on PLUMED-NEST and were not used in our study. (The corresponding changes have also been made to each plumed_sum_bias.dat.) The names of the real files used in our study are without the suffix of _fake.

  • Case_1: In Case 1, we compared 1D alchemical metadynamics with expanded ensemble in decoupling an argon atom in water. This case validates the usage of our algorithm.

  • Case_2: In Case 2, we compared 1D alchemical metadynamics with 2D alchemical metadynamics in decoupling a 4-site model in water. This case demonstrates the advantages of introducing configurational collective variables (CVs) in higher-dimensional alchemical metadynamics.

  • Case_3: In Case 3, we compared Hamiltonian replica exchange with 1D + 2D alchemical metadynamics in methylating adenosine in its isolated form and in a duplex. This case shows that alchemical metadynamics eliminates the need for running multiple Hamiltonian replica exchange simulations to calculate relative free energy differences.

  • project_paper: This folder contains LaTex-relevant files for different versions of the paper submitted to arXiv and the Journal of Chemical Theory and Computation (JCTC). It also includes intermediate files for preparing the final figures in the manuscript and relevant codes for generating the figures.

Contributors of the repo

About

No description, website, or topics provided.

Resources

Stars

Watchers

Forks

Releases

No releases published

Packages

No packages published

Languages

  • Jupyter Notebook 80.9%
  • Python 14.0%
  • TeX 4.9%
  • Shell 0.2%