Closed Alexsp32 closed 2 months ago
NQCModels.mobileatoms
can be implemented to return the indices of the atoms that are unconstrained. This function can then be used with the integrators to ensure only the mobileatoms
move. This technique is already used in some places such as IESH but only supports atoms that are fixed in place.
Implemented in #38
This would be useful to carry over constraints such as fixed atomic positions to dynamics simulations in NQCD and avoid unintended behaviour with stochastic dynamics methods.
Add-on to NQCD/NQCDynamics.jl#295.