Question about ONIOM/QMMM


Click here for full thread
  • Guest -
I'm pretty new to NWChem and I would like to see If it is possible to use ONIOM or QM/MM for my system. For example, in a water dimer ONIOM calculation, I want to do the following:

Water molecule #1
Optimize with RI-MP2

Water dimer #2
Freeze it + use MM or an available Force Field on NWChem

so I know the following so far for an input file

===========================================================

start filename

memory 50MW

geometry noautosym units angstrom
O 0.88 -0.01 0.00
H 1.41 0.78 0.00
H 1.30 -0.87 0.00
H -0.07 0.05 0.00
O -1.03 0.11 0.00
H -1.41 -0.76 0.00
end

basis "ri-mp2 basis" print
BASIS SET: (4s,3p,2d,1f)
H S 8.5128275910 1.0000000
H S 1.8730891166 1.0000000
H S 0.52618426941 1.0000000
H S 0.28973989936 1.0000000
H P 2.3725491635 1.0000000
H P 1.1804084468 1.0000000
H P 0.60382204241 1.0000000
H D 1.8096373189 1.0000000
H D 1.1439726055 1.0000000
H F 1.8063060576 1.0000000
  1. BASIS SET: (8s,6p,5d,3f,1g)
O S 366.54531503 1.0000000
O S 76.693314858 1.0000000
O S 24.573835091 1.0000000
O S 8.4165585072 1.0000000
O S 3.0202080356 1.0000000
O S 1.3331563617 1.0000000
O S 0.81250658727 1.0000000
O S 0.36505253662 1.0000000
O P 52.854386423 1.0000000
O P 15.021795824 1.0000000
O P 3.9137719482 1.0000000
O P 2.1678840611 1.0000000
O P 0.97419243984 1.0000000
O P 0.51845615412 1.0000000
O D 15.907733132 1.0000000
O D 5.3799329215 1.0000000
O D 2.9552623921 1.0000000
O D 1.1624458997 1.0000000
O D 0.50353399421 1.0000000
O F 4.6873846701 1.0000000
O F 2.1512789489 1.0000000
O F 1.0894068169 1.0000000
O G 2.3270964878 1.0000000
end

set geometry:actlist 1:3

oniom
high ri-mp2?? basis ???
low ???
end

task oniom optimize

========================================================

My questions are at the question marks above in the oniom calc.

I did not see a force field available for oniom calculation. is it possible? or can it be done
with QM/MM? if so how?



Thanks a lot in advance!

Sammer