QM/MM calculation with a charge qm region

From NWChem

You are viewing a single post from the thread title above
Jump to: navigation, search

Click here for full thread
Just Got Here
Threads 1
Posts 1
Hi,

I want to perform a qmmm optimization with nwchem 6.1. First I prepare my system that come from amber with this script.

The system is completely uncharged.

start test
prepare
source complex.pdb
new_top new_seq
new_rst
solvent name WAT model spce
modify segment 1 quantum
update lists
ignore
write test_ref.pdb
write test_ref.rst
end
task prepare


I have put the resid 1 as quantum, this resid is charged (-1)

How I put the charge for the qm region in the following script qmmm.nw and is it mandatory to give the charge for the qm region or not.



permanent_dir ./perm
scratch_dir ./data
start test

md
msa 200000
system test_ref
end
basis
* library 6-31++g  
end

dft
xc b3lyp
end

qmmm
region all
maxiter 100 3000
ncycles 5
density espfit
xyz foo
end

task qmmm dft optimize

Thx very much


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


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