Torsion Meeting Summary (Feb 19, 2019)

Summary from torsion subgroup meeting on February 19, 2019

Atom Order: There can be multiple ways to tag SMILES involving self-symmetry and C. Stern raised the question of how to ensure the same geometry always gets assigned the same atom ordering. D. Smith said that a string comparison is currently possible, where string is contained in a QCEngine job description. J. Chodera thinks that having identical geometries with different atom ordering is highly improbable and should not be a major concern.

Fragmentation: The underlying assumption of torsion fitting is that the MM energy can be decomposed and decoupled from other terms, leaving the clean torsion energy profile as a single function to fit. J. Chodera is interested in how to validate this assumption. D. Smith suggested looking into AMOEBA SAPT, which identifies specific orbital pairs to calculate energies between, and label them as required to correspond to certain MM terms. David Sherrill might be a good person to reach out to for more details about SAPT calculations.

Due to licensing issues with using OpenEye tools, it’s desirable that fragmentation is performed in OE-independent way in the future. D. Smith suggested fragmenting based on Wiberg bond orders and potentially training a neural network based on that data.

QCFractal: MolSSI expects to launch the QCFractal closed beta soon, although later than planned, but that gives more time to integrate QCFractal into everyone’s CI and for local testing. This period will likely introduce many API breaking changes.

TorsionDrive / geomeTRIC: geomeTRIC now achieves faster convergence and redundant files are removed. TorsionDrive can now run with limited ranges for torsion scans. It can break the symmetry of the torsion map, but the user should know what they’re doing. For 2D scans, it is possible to set ranges independently.

Published by in science and tagged atom order, fitting, fragmentation, geomeTRIC, QCFractal, SAPT, SMILES, torsion and TorsionDrive using 285 words.