Closed zjing7 closed 2 weeks ago
Be careful the weight should be set based on total QM, not QM-MM.
Thanks, Pengyu
From: zjing7 @.> Sent: Thursday, November 7, 2024 8:31:34 AM To: TinkerTools/poltype2 @.> Cc: Subscribed @.***> Subject: [TinkerTools/poltype2] fix bug in Boltzmann-weighted fitting (PR #59)
There was a typo in the objective function, (model - targetweight). It is corrected as (model-target)weight
You can view, comment on, or merge this pull request online at:
https://github.com/TinkerTools/poltype2/pull/59
Commit Summary
File Changes
(1 filehttps://github.com/TinkerTools/poltype2/pull/59/files)
Patch Links:
— Reply to this email directly, view it on GitHubhttps://github.com/TinkerTools/poltype2/pull/59, or unsubscribehttps://github.com/notifications/unsubscribe-auth/ABNC6XTYSBQZHUCMFU467SLZ7N2UNAVCNFSM6AAAAABRLK4W6SVHI2DSMVQWIX3LMV43ASLTON2WKOZSGY2DCMJZGU2DIMI. You are receiving this because you are subscribed to this thread.Message ID: @.***>
The weight is correct
weightlist=numpy.exp(-numpy.array(qm_energy_list)/poltype.boltzmantemp)
There was a typo in the objective function, (model - targetweight). It is corrected as (model-target)weight