CHarMM c32b2
CHARMM Development and Information
CHARMM (Chemistry at HARvard Macromolecular Mechanics) is a program for macromolecular dynamics and mechanics. It performs standard molecular dynamics in many different ensembles (e.g., NVE, NVT, NPT) using state-of-the-art algorithms for timestepping, long range force calculation and periodic images. CHARMM can be used for energy minimization, normal modes and crystal optimizations as well. The potential energy functions available for use with CHARMM have been extensively parameterized for simulations of proteins, nucleic acids and lipids. Free energy methods for chemical and conformational free energy calculations are also fully developed and available in CHARMM. Many other novel tools have been developed and are available in CHARMM, these include: replicas (multiple copies); many types of restraints and constraints, including fixed atoms, atomic, NOE, dihedral and internal coordinate restraints and generalized SHAKE for bonds and arbitrary internal coordinates; minimium energy path following and transition state optimization; etc.
CHARMM development began in the research group of Professor Martin Karplus at Harvard University. Development continues throughout the world with contributing developers in over 20 universities, research institutes and companies. CHARMM is available to academic users in full source form, plus extensive documentation and test input scripts, through Harvard University. A license for CHARMM can be obtained by inquiring to:
Professor Martin Karplus
Department of Chemistry & Chemical Biology
Harvard University
12 Oxford Street
Cambridge, MA 02138
(Contact Martin Karplus)
The appropriate references for CHARMM in publications are:
CHARMM: A Program for Macromolecular Energy, Minimization, and Dynamics Calculations, J. Comp. Chem. 4, 187-217 (1983), by B. R. Brooks, R. E. Bruccoleri, B. D. Olafson, D. J. States, S. Swaminathan, and M. Karplus.
CHARMM: The Energy Function and Its Parameterization with an Overview of the Program, in The Encyclopedia of Computational Chemistry, 1, 271-277, P. v. R. Schleyer et al., editors (John Wiley & Sons: Chichester, 1998), by A. D. MacKerell, Jr., B. Brooks, C. L. Brooks, III, L. Nilsson, B. Roux, Y. Won, and M. Karplus.
from http://brooks.chem.lsa.umich.edu/charmm_docs/charmm.htm
Usage
If you wish to test your input file before submitting production jobs to PBS you should start by:
module load charmm
charmm < filename.inp > & filename.out [[param:value] ...] [ & ]
where:
- charmm
- The actual name of the program (or script which runs the program) on the computer system being used
- filename.inp
- A text file containing CHARMM input commands, using the command syntax described in usage.doc. The filename is up to the user.
- filename.out
- The log file for the CHARMM run, containing echoed commands, and various amounts of command output. The output print level may be increased or decreased in general, and procedures such as minimization and dynamics have printout frequency specifications.
- param:value
- CHARMM script parameters (essentially user defined variables) may optionally be defined on Unix shell command line, using the form param:value; the value is substituted for occurrences of @param in the input script.
- &
- The optional ampersand will place the program in the background
For foreground usage, e.g. graphics, start the program by typing it's name alone, wait for the program "header" to be printed, and enter a title. Then use a "stream" file, e.g.
stream psf-crd.str
will read CHARMM commands from the file named psf-crd.str as if they were being typed at the keyboard; as implied by the name, the stream file reads the PSF, PARAms, COORdinates, etc., the basic starting point for all work with CHARMM.
-from http://brooks.chem.lsa.umich.edu/charmm_docs/c32docs/c32a1/
Remember that a run on the head node must complete within 5 minutes; so you'll need to run your production codes from PBS.
As before start with:
module load charmm
Then submit a script, such as the one below,
qsub ex.scr
#PBS -N charmm_test
#PBS -q route
#PBS -l nodes=2:ppn=2,walltime=24:00:00
#PBS -S /bin/csh
#PBS -M your-email-address
#PBS -m abe
#PBS -V
#
echo "I ran on:"
cat $PBS_NODEFILE
#
#cd to your execution directory first
cd /home/your-user-name
#
mpirun -np 4 $chmroot/exec/gnu/charmm < $chmroot/test/c28test/gridtest.inp
#
however if your input file to charm contains a goto you should use the following form of PBS script.
#PBS -N charmm_test
#PBS -q route
#PBS -l nodes=2:ppn=2,walltime=24:00:00
#PBS -S /bin/csh
#PBS -M your-email-address
#PBS -m abe
#PBS -V
#
echo "I ran on:"
cat $PBS_NODEFILE
#
#cd to your execution directory first
cd /home/your-user-name
#
echo "* CHARMM test stream file for parallel runs" >stream.inp
echo "* where rewinding stdin is not possible" >>stream.inp
echo "stream @stream" >>stream.inp
echo "stop" >>stream.inp
echo $chmroot
mpirun -np 4 $chmroot/exec/gnu/charmm stream=$chmroot/test/c28test/gridtest.inp
