Skip to content

Commit

Permalink
Merge pull request #29 from WelbornGroup/usage_procedure
Browse files Browse the repository at this point in the history
adds a paragraph about usage procedure to the paper
  • Loading branch information
valeriewelborn authored Oct 12, 2020
2 parents b645005 + d85d560 commit bfee1f8
Showing 1 changed file with 2 additions and 0 deletions.
2 changes: 2 additions & 0 deletions paper.md
Original file line number Diff line number Diff line change
Expand Up @@ -33,6 +33,8 @@ Improved electrostatics with AMOEBA give a unique opportunity to accurately comp

`ELECTRIC` [@electric] is a MolSSI Driver Interface (MDI) [@mdi_repo; @barnes_taylor_arnold_2020_3659285] driver that utilizes Tinker [@doi:10.1021/acs.jctc.8b00529] to analyze specific components of electric fields that are modeled using the AMOEBA force field. `ELECTRIC` parses Tinker trajectories and orchestrates additional Tinker calculations in order to project components of the electric fields onto user-defined bonds (specified by two atoms). It outputs the field in MV/cm, which is the sum of the direct field (from permanent electrostatics) and the induced field (from the induce dipole term), projected onto the bond unit vector (i.e. normalized by the bond length). `ELECTRIC` enables splitting of the total field into contributions from different components of the system, by molecules or by residues as specified in a reference PDB file. In summary, `ELECTRIC` was designed to expand quantitative system characterization via the computation of electric fields with user-friendly processing tools of Tinker-AMOEBA simulations.

In practice, the user needs a Tinker trajectory file (i.e. `filename.arc`), a Tinker input file (i.e. `tinker.key`) stripped from any keywords corresponding to periodic boundary conditions and a Tinker reference snapshot (`filename.xyz` - this can be the first frame in `filename.arc`) where the dimensions of the box, usually printed on line 2, has been deleted. Tinker is then launched as an MDI engine before the `ELECTRIC` driver. Options regarding the electric fields calculations, such as probe index number, number of frames at the beginning of the trajectory file to skip (equilibration procedure), reference `PDB` file, etc. are specified as command line arguments when launching the driver. A complete usage procedure is provided in the `README` file.

# Mathematics
The electric field at atom $i$, $\vec{E}^i$, has components defined as:
\begin{equation}
Expand Down

0 comments on commit bfee1f8

Please sign in to comment.