Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Add DOIs for OPLS atomtypes #158

Merged
merged 11 commits into from
Aug 17, 2018
Merged

Conversation

summeraz
Copy link
Contributor

@summeraz summeraz commented Mar 5, 2018

Our OPLS force field XML (which we created using atomtype definitions and parameters from GROMACS) features several atomtypes with defined SMARTS strings that lack a DOI reference. This PR adds DOIs to these atomtypes so that the source of the relevant parameters is more transparent.

@summeraz
Copy link
Contributor Author

summeraz commented Mar 5, 2018

I'm not finished with this yet, but wanted to go ahead and open the PR so that I can log comments as to issues I run into trying to locate the proper DOI references.

@summeraz
Copy link
Contributor Author

summeraz commented Mar 5, 2018

The sulfoxide OPLS parameters do not appear to have a published source. The GROMACS parameters that we have in our XML appear in a recent OPLS paper, where these are compared with parameters derived for the improved OPLS-AAx force field. It seems most papers in the literature that used these parameters link to this paper, where the parameters are listed in Ref. 18, although the charges used here are different than those in our XML and as listed in the 2017 paper. I went ahead and added both DOIs to this atomtype. I do not believe our XML parser is currently designed to take more than a single DOI, so I will have to update that.

Also, it might be useful to add a note attribute to Types within the AtomTypes section (unless we already have the ability to do this), so that messages (like the above) could be added if citations are not straightforward.

@summeraz
Copy link
Contributor Author

summeraz commented Mar 5, 2018

Parameters for the α-carbon in alkyl pyrroles do not appear to have a published source. It appears with OPLS that for most α-carbons the epsilon and sigma used are the same as for alkane carbons and the charge is simply adjusted to preserve molecular charge neutrality. As per my above comment, this is something we could probably just add to a note attribute and link the DOI to the alkane carbon parameters.

@summeraz
Copy link
Contributor Author

summeraz commented Mar 5, 2018

Parameters for cyclopropanes do not appear to have a published source. However, these parameters are identical to those for alkanes. Bonded parameters are obviously different; however, for now I'm just focusing on nonbonded parameters since we do not yet have the ability to add DOI tags to bonded parameters.

@summeraz
Copy link
Contributor Author

summeraz commented Mar 5, 2018

DOIs are still needed for propylene and ethylene carbonate parameters and for alkyne parameters.

@mattwthompson
Copy link
Member

And to think there was a time I thought adding SMARTS strings was the hard part!

@ctk3b
Copy link
Member

ctk3b commented Mar 7, 2018

Nice! Thank you for this. Now is there an easy way to review these... 🤔

Copy link
Contributor

@justinGilmer justinGilmer left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Phew, verified the DOI links and atom type parameters from each reference to the designated atom type here. They all look good to me.

<Type name="opls_151" class="Cl" element="Cl" mass="35.453" def="[Cl][C;X4]" desc="Cl in alkyl chlorides" />
<Type name="opls_152" class="CT" element="C" mass="12.011" def="[C;X4](Cl)([Cl,C,H])([Cl,C,H])([Cl,C,H])" desc="RCH2Cl in alkyl chlorides" />
<Type name="opls_153" class="HC" element="H" mass="1.008" def="[H][C;X4](Cl)([Cl,C,H])([Cl,C,H])" overrides="opls_140" desc="H in alkyl chlorides" />
<Type name="opls_151" class="Cl" element="Cl" mass="35.453" def="[Cl][C;X4]" desc="Cl in alkyl chlorides" doi="10.1021/jp0484579"/>
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Confirmed

<Type name="opls_152" class="CT" element="C" mass="12.011" def="[C;X4](Cl)([Cl,C,H])([Cl,C,H])([Cl,C,H])" desc="RCH2Cl in alkyl chlorides" />
<Type name="opls_153" class="HC" element="H" mass="1.008" def="[H][C;X4](Cl)([Cl,C,H])([Cl,C,H])" overrides="opls_140" desc="H in alkyl chlorides" />
<Type name="opls_151" class="Cl" element="Cl" mass="35.453" def="[Cl][C;X4]" desc="Cl in alkyl chlorides" doi="10.1021/jp0484579"/>
<Type name="opls_152" class="CT" element="C" mass="12.011" def="[C;X4](Cl)([Cl,C,H])([Cl,C,H])([Cl,C,H])" desc="RCH2Cl in alkyl chlorides" doi="10.1021/jp0484579"/>
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Confirmed

<Type name="opls_153" class="HC" element="H" mass="1.008" def="[H][C;X4](Cl)([Cl,C,H])([Cl,C,H])" overrides="opls_140" desc="H in alkyl chlorides" />
<Type name="opls_151" class="Cl" element="Cl" mass="35.453" def="[Cl][C;X4]" desc="Cl in alkyl chlorides" doi="10.1021/jp0484579"/>
<Type name="opls_152" class="CT" element="C" mass="12.011" def="[C;X4](Cl)([Cl,C,H])([Cl,C,H])([Cl,C,H])" desc="RCH2Cl in alkyl chlorides" doi="10.1021/jp0484579"/>
<Type name="opls_153" class="HC" element="H" mass="1.008" def="[H][C;X4](Cl)([Cl,C,H])([Cl,C,H])" overrides="opls_140" desc="H in alkyl chlorides" doi="10.1021/jp0484579"/>
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Confirmed

@@ -686,8 +686,8 @@
<Type name="opls_715" class="opls_715" element="H" mass="1.008"/>
<Type name="opls_716" class="opls_716" element="C" mass="12.011"/>
<Type name="opls_717" class="opls_717" element="H" mass="1.008"/>
<Type name="opls_718" class="CA" element="C" mass="12.01100" def="[C;X3;r6]1(F)[C;X3;r6][C;X3;r6][C;X3;r6][C;X3;r6][C;X3;r6]1" overrides="opls_145"/>
<Type name="opls_719" class="F" element="F" mass="18.9984" def="F[C;%opls_718]" desc="fluorobenzene F" overrides="opls_965"/>
<Type name="opls_718" class="CA" element="C" mass="12.01100" def="[C;X3;r6]1(F)[C;X3;r6][C;X3;r6][C;X3;r6][C;X3;r6][C;X3;r6]1" overrides="opls_145" doi="10.1021/jp0484579"/>
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Confirmed

<Type name="opls_718" class="CA" element="C" mass="12.01100" def="[C;X3;r6]1(F)[C;X3;r6][C;X3;r6][C;X3;r6][C;X3;r6][C;X3;r6]1" overrides="opls_145"/>
<Type name="opls_719" class="F" element="F" mass="18.9984" def="F[C;%opls_718]" desc="fluorobenzene F" overrides="opls_965"/>
<Type name="opls_718" class="CA" element="C" mass="12.01100" def="[C;X3;r6]1(F)[C;X3;r6][C;X3;r6][C;X3;r6][C;X3;r6][C;X3;r6]1" overrides="opls_145" doi="10.1021/jp0484579"/>
<Type name="opls_719" class="F" element="F" mass="18.9984" def="F[C;%opls_718]" desc="fluorobenzene F" overrides="opls_965" doi="10.1021/jp0484579"/>
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Confirmed

<Type name="opls_906B" class="opls_906B" element="C" mass="12.011"/>
<Type name="opls_907" class="opls_907" element="C" mass="12.011"/>
<Type name="opls_908" class="opls_908" element="C" mass="12.011"/>
<Type name="opls_909" class="H" element="H" mass="1.00800" def="H[N;%opls_900]" desc="H(N) primary amines"/>
<Type name="opls_909" class="H" element="H" mass="1.00800" def="H[N;%opls_900]" desc="H(N) primary amines" doi="10.1021/ja984106u"/>
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Confirmed

<Type name="opls_910" class="opls_910" element="H" mass="1.008"/>
<Type name="opls_911" class="HC" element="H" mass="1.00800" def="HC[N;%opls_900]" desc="H(C) for C bonded to N in amines, diamines (aziridine H2,H3)" overrides="opls_140"/>
<Type name="opls_911" class="HC" element="H" mass="1.00800" def="HC[N;%opls_900]" desc="H(C) for C bonded to N in amines, diamines (aziridine H2,H3)" overrides="opls_140" doi="10.1021/ja984106u"/>
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Confirmed

@@ -801,11 +801,11 @@
<Type name="opls_958" class="opls_958" element="H" mass="1.008"/>
<Type name="opls_959" class="opls_959" element="C" mass="12.011"/>
<Type name="opls_960" class="opls_960" element="C" mass="12.011"/>
<Type name="opls_961" class="CT" element="C" mass="12.01100" def="[C;X4](F)(F)(F)*" desc="perfluoroalkane CF3" overrides="opls_962"/>
<Type name="opls_962" class="CT" element="C" mass="12.01100" def="[C;X4](F)(F)(*)*" desc="perfluoroalkane CF2"/>
<Type name="opls_961" class="CT" element="C" mass="12.01100" def="[C;X4](F)(F)(F)*" desc="perfluoroalkane CF3" overrides="opls_962" doi="10.1021/jp004071w"/>
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Confirmed

<Type name="opls_961" class="CT" element="C" mass="12.01100" def="[C;X4](F)(F)(F)*" desc="perfluoroalkane CF3" overrides="opls_962"/>
<Type name="opls_962" class="CT" element="C" mass="12.01100" def="[C;X4](F)(F)(*)*" desc="perfluoroalkane CF2"/>
<Type name="opls_961" class="CT" element="C" mass="12.01100" def="[C;X4](F)(F)(F)*" desc="perfluoroalkane CF3" overrides="opls_962" doi="10.1021/jp004071w"/>
<Type name="opls_962" class="CT" element="C" mass="12.01100" def="[C;X4](F)(F)(*)*" desc="perfluoroalkane CF2" doi="10.1021/jp004071w"/>
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Confirmed

<Type name="opls_963" class="opls_963" element="C" mass="12.011"/>
<Type name="opls_964" class="opls_964" element="C" mass="12.011"/>
<Type name="opls_965" class="F" element="F" mass="18.9984" def="FC" desc="perfluoroalkane F"/>
<Type name="opls_965" class="F" element="F" mass="18.9984" def="FC" desc="perfluoroalkane F" doi="10.1021/jp004071w"/>
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Confirmed

@summeraz summeraz added this to the 0.5 milestone May 30, 2018
@summeraz summeraz merged commit dfc1288 into mosdef-hub:master Aug 17, 2018
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

None yet

4 participants