free energy topology file error

From NWChem

Viewed 2065 times, With a total of 0 Posts
Jump to: navigation, search

Just Got Here
Threads 2
Posts 2
Dear All,

I have been trying to use the free energy calculation capability of QM/MM module wiht NWChem 6.3. I did google a lot before posting here.

I am mostly following the instructions here

and here

First, it seems that we should use "task qmmm dft fep" instead of "task qmmm fep" shown in the example. The latter will results in "task: no task input for theory?" error.

Second, even with "task qmmm dft fep", I got following error

* 0: Error opening topology file 2
0:0::: 2
forrtl: error (76): Abort trap signal

One may suspect something is wrong with my topology file.
However, the same topology file has been used for optimizing the two configurations successfully. Can anyone give me some suggestions?
Where can I find a complete example (with all required files given) for running FEP with NWchem? There is no such test file in the QA directory.

Thanks a lot!

following is my input file. I have top file as, restart file as BB_path.rst

 start BB

system BB_path
noshake solute

print low
nsamples 1000
ncycles 2

set qmmm:fep_geom BB_path-s.xyzi BB_path-e.xyzi
set qmmm:fep_esp BB_path-s.esp BB_path-e.esp
set qmmm:fep_lambda 0.0 0.1
set qmmm:fep_deriv .true.

task qmmm dft fep

Forum >> NWChem's corner >> QMMM

Who's here now Members 0 Guests 1 Bots/Crawler 0

AWC's: 2.5.10 MediaWiki - Stand Alone Forum Extension
Forum theme style by: AWC