Test Cases & Examples¶
This section contains brief descriptions of the sample calculations found in the TEST subdirectory of the Tinker distribution. These test cases exercise several of the current Tinker programs and are intended to provide an overview of the capabilities of the package. In addition, a number of input files for various molecular systems are provided in the EXAMPLE subdirectory.
Computes an estimation of the hydration free energy difference between Cl- anion vs. Br- anion via a 2 picosecond simulation on a “hybrid” anion in a box of water followed by a free energy perturbation calculation.
Performs an initial energy minimization on a periodic box containing 150 argon atoms followed by 6 picoseconds of a molecular dynamics using a modified Beeman integration algorithm and a Bersedsen thermostat.
Computes an estimation of the hydration free energy difference between Rb+ cation vs. Cs+ cation via a 2 picosecond simulation of each cation in a box of water followed by a BAR free energy calculation.
Performs a set of 10 Gaussian density annealing (GDA) trials on a cluster of 13 argon atoms in an attempt to locate the global minimum energy structure.
Generates a Tinker file from a PDB file, followed by a single point energy computation and determination of the molecular volume and surface area.
First approximately locates the transition state between chair and boat cyclohexane, followed by subsequent refinement of the transition state and a final vibrational analysis to show that a single negative frequency is associated with the saddle point.
Performs 10 steps of molecular dynamics on a pre-equilibrated system of DHFR protein in a box or water using the AMOEBA force field. Note this test case is the so-called Joint Amber-CHARMM “JAC” benchmark containing 23558 total atoms.
Finds all the local minima of alanine dipeptide via a potential energy surface scan using torsional modes to jump between the minima.
Produces coordinates from the met-enkephalin amino acid sequence and phi/psi angles, followed by truncated Newton energy minimization and determination of the lowest frequency normal mode.
Performs fitting of torsional parameter values for the ethanol C-C-O-H bond based on relative quantum mechanical (G09) energies for rotating the C-O bond.
Generates a unit cell from fractional coordinates, followed by full crystal energy minimization and determination of optimal carbonyl oxygen energy parameters from a fit to lattice energy and structure.
Finds the lowest-frequency normal mode of bacteriorhodopsin using vibrational analysis via a sliding block iterative matrix diagonalization. Alter the gpcr.run script to save the file gpcr.001 for later viewing of the mode.
Performs a rigid-body optimization of the packing of two idealized polyalanine helices using only van der Waals interactions.
Performs a short MD simulation of the monoclinic ice V crystal form using the iAMOEBA water model, pairwise neighbor lists and PME electrostatics.
Generates three distance geometry structures for intestinal fatty acid binding protein from a set of NOE distance restraints and torsional restraints.
Processes distributed multipole analysis (DMA) output to extract coordinates and permanent multipoles, set local frames and polarization groups, remove intramolecular polarization, detect and average equivalent atomic sites.
Calculates the self-diffusion constant and the N-N radial distribution function for liquid nitrogen via analysis of a 50ps MD trajectory.
Converts a sodium chloride assymetric unit to the corresponding unit cell, then runs a crystal minimization starting from the initial diffraction structure using Ewald summation to model the long-range electrostatic interactions.
Generates capped alanine tetrapeptide in an extended conformation, then use Monte Carlo Minimization with random torsional moves to find the global minimum energy structure.
Compares analytical and finite difference numerical gradients determined over Cartesian and internal coordinates for vasopressin using the AMOEBA force field model.
Fits the electrostatic potential for the TIP3P, AMOEBA and HIPPO water models to a QM-derived potential (MP2/aug-cc-pVTZ) on a grid of points outside the molecular surface.