Comments (11)
Yes, I agree. I'll check how things are implemented in SOMD and see if @jmichel80 has any input. It's an easy fix to make once we agree on the correct approach. (It's almost like you want to do the partitioning before the merge, i.e. just add the dummies afterwards.)
from biosimspace.
Thanks for this. I hadn't thought about dummies when writing the code and agree with your suggested solution. I'll check what SOMD does too, since I believe I largely followed the implementation there, or at least tested mine against that one when first written.
from biosimspace.
I think the issue is that the repartitioning code generates a connectivity object to working the bonding when performing the HMR. This doesn't work properly if an atom is a dummy, since it's based on equilibrium bond radii for specific elements. I think I'll need to adjust this to use the elements from the non-dummy end state instead. Hopefully this should be easy enough to fix.
Thanks for reporting.
from biosimspace.
Thanks @lohedges, would this work if both endstates have dummies? I think that's the common use case.
from biosimspace.
Yes, assuming the velocity conversion is okay I think capping with a warning would be an okay solution.
from biosimspace.
I think there are two approaches:
- Include dummies in the repartitioning by using the element property from the opposite end state. In this case dummies are part of the repartitioning so the mass of the end state is consistent before and after repartitioning, but the mass of a dummy atom isn't necessarily the same between end states as it was before.
- Exclude dummies from the repartitioning, then set their mass to that of the opposite end state after repartitioning. In this case, the overall molecule mass of the end state may differ before / after repartitioning (if the sum includes dummies) but the dummy mass is consistent between end states.
from biosimspace.
Yes I was thinking about that as well, it seems to me that the second approach makes more sense, as one would not expect conservation of mass in an alchemical perturbation setting anyway as you are adding dummies regardless.
from biosimspace.
SOMD uses the second approach above, i.e. excluding dummy atoms and replacing with the repartitioned mass at the other end state. I'll implement this and check that things are working as expected.
from biosimspace.
Okay, this is now fixed, i.e. I get the same behaviour as SOMD, which agrees with the approach agreed upon above. This requires a new build of Sire, so I'll submit a PR when that's done. (Likely Monday.)
Thanks again for catching this.
from biosimspace.
Many thanks for that @lohedges, I am off until Jan 3rd, so I will only be able to test it then. Happy holidays!
from biosimspace.
Thanks, you too. Have a good break!
from biosimspace.
Related Issues (20)
- Storage quota exceeded on the Anaconda cloud HOT 1
- HILLS file not updating during restarting simulation HOT 4
- Is it possible to add ions to an already solvated system, without adding any more water molecules? HOT 4
- Problem when converting Gro/Top to Rst7 HOT 12
- Make _removeDummies() recognise dummy bonds HOT 10
- Would it be possible to include HiMap as an alternative to LoMap in the generateNetwork method? HOT 2
- Unable to import BioSimSpace after mamba installation HOT 12
- BioSimSpace.Align.viewMapping: py3Dmol.view not rendering in Jupyter Notebook HOT 11
- `_toRegularMolecule` failing with `convert_amber_dummies=True` HOT 15
- Absolute Binding Free Energy Calculations with GROMACS HOT 7
- Pickling a parametrised molecule from an SDF file with GAFF2 HOT 4
- Are atomtypes meant to be updated when using repartitionHydrogenMass? HOT 16
- [BUG] OSError "It looks like you failed to include a topology file." when reading SDF which doesn't contain redundant bond order information HOT 10
- [BUG] Cannot import BioSimSpace after installation of biosimspace 2023.1.2 using mamba HOT 3
- AnalysisError: SOMD free-energy analysis failed! HOT 5
- [BUG] HOT 10
- Reporting a vulnerability
- [BUG] Water molecule coordinates are setting to (0,0,0) HOT 4
- Issue in reproducing the funnel methadynamics tutorial HOT 2
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from biosimspace.