Skip to content

Commit

Permalink
Update paper.md
Browse files Browse the repository at this point in the history
  • Loading branch information
Sulstice committed Dec 21, 2021
1 parent d68a1c2 commit eb400f0
Showing 1 changed file with 2 additions and 3 deletions.
5 changes: 2 additions & 3 deletions paper.md
Original file line number Diff line number Diff line change
Expand Up @@ -191,15 +191,14 @@ bond, angle, dihedral are associated with a charge increment value subtracted fr
This value is associated with how these two atoms behave, electronically, in the environment space it is in. It is noted
that the dihedral charge increment is set to a limit to of 50.

The lower the distribution is to 0 the more performant the forcefield is. The distributions are reported in accordance with bonds, angles, dihedrals, charge classifications of the charmm potential energy
equation.
The lower the distribution is to 0 the more performant the forcefield is. The distributions are reported in accordance with bonds, angles, dihedrals, charge classifications of the charmm potential energy equation.

<p align="center">
<img width="1400" height="400" src="images/figures/figure_3.png">
<i>Figure 4: Charmm Potential Energy Equation Picture Depiction: Top row and bond row is bonded and non-bonded terms respectively. </i>
</p>

We passed each object individually into `CGenFF` recorded the results. For anything `CGenFF` determines it cannot find an adequate substitute within reason
We passed each object individually into `CGenFF` recorded the results. <b>Note</b> that due to our initial input being SMILES we had to offer up some of the original functionality of `CGenFF` that was used to process `mol2` for bond assignment. Whereas now we rely on `SDF` and go through our own in-house decision tree. For anything `CGenFF` determines it cannot find an adequate substitute within reason
returns a failed response. These failed responses are captured. It is noted, that the IUPAC blue book is a list of radicals (Object 14 in the tests) which is the reason for it's repeated compound failures which you can see in the log file. Other noteworthy failures that would be of particular interest are listed in `Figure 5`. For example, cyclobutadiene is a non-traditional ring system with a lot of ring strain where the carbon atom types are common. `CGenFF` might determine that this particular ring system with it's existing atom type network is not allowed or detrimental to the network if added and needs to be handled with care. Silicon-based systems seems to be more ubuiqituious in these data sets and a lot of the failures are attributed to silicon-based compounds. An interesting functional group to handle would be the allene-based compounds and perhaps warrant a new carbon atom type.

<p align="center">
Expand Down

0 comments on commit eb400f0

Please sign in to comment.