[AMBER] no collective variables output for infe=1 and abfqmmm

From: FLORA, JOSEPH <FLORA.cec.sc.edu>
Date: Tue, 18 Aug 2020 01:52:20 +0000

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

Received on Mon Aug 17 2020 - 19:00:02 PDT
Custom Search