From NWChem
Viewed 1552 times, With a total of 2 Posts
|
|
|
Just Got Here
Threads 1
Posts 1
|
|
| 11:05:25 AM PST - Sat, Feb 4th 2012 |
|
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
|
|
|
-
Marat Forum:Admin, Forum:Mod, NWChemDeveloper, bureaucrat, sysop
|
|
|
Clicked A Few Times
Threads 2
Posts 43
|
|
| 10:29:51 PM PST - Thu, Feb 23rd 2012 |
|
You just add
charge -1
somewhere in the beginning of your input file
|
|
|
|
|
|
Clicked A Few Times
Threads 12
Posts 43
|
|
| 6:15:44 AM PST - Mon, Dec 24th 2012 |
|
Thanks, this suggestion also works for me.
Jingbo
|
|
Edited On 6:16:05 AM PST - Mon, Dec 24th 2012 by Wjb0920
|
|
|
AWC's:
2.5.10 MediaWiki - Stand Alone Forum Extension
Forum theme style by: AWC
