Roundtable updates | CC Got people to commit to getting outstanding LiveCOMS sections by the end of the week. The last section left undone is the biopolymer benchmarking one, so I’ll have more clarity on those details soon Submitted dipeptide torsiondrives. Updates TSCC workers environment, so those should be good to go shortly. Â
MT Short week last week. Tried and made some progress on the reference energies package. Quickly ran into cases/details where the toolkit doesn’t handle things correctly, though it’s some question as to whether the toolkit should. Also found some issues with the SMIRNOFF spec. So I opened up some OFF-EPs, but this took a lot of work since it’s a lot of writing. The first ones were fairly straightforward and I expect easy approval. But there are some hard ones coming. Got a support request from Roitberg group member. The issue was fairly unclear/somewhat low quality. Also the email went to my personal email which I don’t like. Generally I’ve seen a lot more issue/question traffic from external users, which is good, but we should make sure to be disciplined about it so we still get our own work done. Worked on toolkit refactor to make OpenMM optional. Big lift here is to remove openMM units and replace with openff units. Except for some possible edge cases with GBSA and VSites, this should be in good shape. Most of the test suite looks good. For the topology and units refactor, we could use feedback and user testing – Please get in touch if you’re able to help here. M1 builds are almost supported. One remaining big sticking point netcdf-fortran Wrote+published interchange blog post Interchange project update LW – Re: API breaking change – Should we stage other API breaking changes with the Feb release as well? JW – We don’t really have a plan/policy, now that I think about it. For now, the best thing to do would be to open API-breaking PRs to the topology-refactor branch instead of master so that we can consolidate all the API breaks in the Feb release.
DD LW Short week last week Refactored PsiRESP to use QCFractal (a while ago) and made docs (last fortnight), now asking for feedback if anyone uses it (last week):
Past two weeks: debugging with Simon (unsuccessfully so far) Yank writing to PDB error: File "/data/homezvol0/lilyw7/miniconda3/envs/psiresp-3.8/lib/python3.8/site-packages/mdtraj/formats/pdb/pdbfile.py", line 353, in write
line = "ATOM %5d %-4s %3s %1s%4d %s%s%s 1.00 %5s %-4s%2s " % ( # Right-justify atom symbol
ValueError: cannot convert float NaN to integerCannot reproduce on local machine (following what we think is the logic flow – loops etc make this a little unobvious) I’ve run into this sort of error before – It’s probably not actually MDTraj’s fault, more likely YANK. Previous solution seemed to be downgrading to openmm 7.5.1 but this may not be solving it for us. Now clear why, and error is deeply buried in Evaluator workfow so it’s hard to pin down. JW – I wonder what we can extract from this investigation - Some sort of principle when we engineer workflows to always have “access hatches” to intermediate data? DD – Does this only happen on distributed DASK tasks? LW – Looking into it. It’s so far happened on remote hosts, but I’m trying to reproduce on my laptop. Somewhat tricky since it took a lot of simulation to get there. DD – My general solution in this kind of situation is to run things serially so I can use python -m pdb to get right to the error and look around. LW – I’m not sure if setting it to one process will run it in serial. I’d need to look at this. JW – I’d love to know if it’s possible to get a debugger itno there. We have lots of distribution architectures in our stack and each one has different tricks forinspecting the underlying runs.
Successfully used Evaluator + ForceBalance to automatically train a force field for properties that are not SolvationFreeEnergys on Slurm UCI infrastructure (hpc3) Catch for hpc3: need to explicitly request gpus Tiny dataset of properties for C/H molecules (alkanes, alkenes, and aromatics) did not show a lot of improvement / difference between am1bcc/resp2 charges. Trying again on a larger dataset of just alkane densities, as there are systematic difference between am1bcc/resp2 charges here and lots of density data
JW – Worked on ring stereo/pseudosymmetric representation problem - Looks like we should be able to make it work
Some molecule sanitization/canonicalization work for benchmarking - OPLS “sanitizes” inputs to have different graphs in some cases, like R[S+2]([O-])([O-])R → RS(=O)(=O)R, which breaks our analysis workflows. Used MolStandardize.rdMolStandardize.Normalize(rdmol) Biopolymer topology refactor – Implementing quick molecule-identity checks, incorporating into parameterization loops. Worked with AMBER dev (Kellon) to standardize on CCD. Will meet again this coming Fri. Â
|