Dear Amber Folks,
I modified the Run.abfqmmm_dmpoh_md to use infe=1 and pmd. The run completed but nothing follows after the following output for the colvars file (i.e., no colvars printed)
# = NFE%PMD ==================================================================
# << anchor(1) : position = -10.000000, 0.000000, 0.000000, 10.000000
# strength = 400.000000, 400.000000 >>
# << anchor(2) : position = 0.000000, 2.792520, 3.141590, 3.141590
# strength = 100.000000, 100.000000 >>
# -----------------------------------------------------------------------------
# MD time (ps), CV(1:2)
# =============================================================================
The energies and constraints at the end of the 10 step test simulation were identical with the test run using nmr restraints (angle restraints were 0), indicating that the simulation propagated properly. I looked at nfe_pmd_hooks.F90 to try to follow the write statements. I think the printout should be happening with following statement but I can't fully follow the bigger picture logic with abfqmmm. I can get full colvars output with adqmmm.
NFE_MASTER_ONLY_BEGIN
if (nfe_real_mdstep) then
if (mod(mdstep, output_freq).eq.0) then
write (unit = pmd_UNIT, fmt = output_fmt) &
sander_mdtime(), cv_inst(1:ncolvars)
call flush_UNIT(pmd_UNIT)
end if
mdstep = mdstep + 1
end if ! nfe_real_mdstep
NFE_MASTER_ONLY_END
Any help is appreciated. I am running Amber18 (and will be in the process of upgrading to Amber20). Attached are the revised Runfile and the disang.cv file used.
Thanks!
Joe
_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
- application/octet-stream attachment: disang.cv
Received on Mon Aug 17 2020 - 19:00:02 PDT