From NWChem
Viewed 2526 times, With a total of 3 Posts
|
|
|
Clicked A Few Times
Threads 10
Posts 18
|
|
| 6:42:12 PM PST - Thu, Jan 8th 2015 |
|
Hi, I'm new to NWChem and, as a self-taught ORCA user, trying to learn as much as I can while I'm preparing to my masters degree in chemistry.
I'm very sorry if this is a stupid question but, after some research, I wasnt able to find a solution.
I've decided to start from the very basics with NWChem, so, my first experiment was to try to calculate the transition state for H2SO4. Here is my input:
# ================================================================
# NWChem input file made in Gabedit
# ================================================================
charge 0
memory 3000 mb
geometry autosym
S -1.139420 -0.181760 -0.065320
O -0.603930 0.894620 -0.859210
O -0.327300 -1.144590 0.633970
O -2.146750 0.489400 0.987400
O -2.158340 -1.022960 -0.975080
H -2.473210 1.320810 0.584900
H -2.293600 -1.887270 -0.534180
end
basis "ao basis" spherical
* library 6-31G
end
scf
Singlet
maxiter 999
thresh 1e-2
end
task SCF saddle
And here is the error:
------------------------------------------------------------------------[color=red]
scf saddle failed 0[/color]
------------------------------------------------------------------------
------------------------------------------------------------------------
current input line :
27: task SCF saddle
------------------------------------------------------------------------
------------------------------------------------------------------------
This type of error is most commonly associated with calculations not reaching convergence criteria
------------------------------------------------------------------------
For more information see the NWChem manual at http://www.nwchem-sw.org/index.php/NWChem_Documentation
For further details see manual section:
0:0:scf saddle failed:: -1
(rank:0 hostname:localhost.localdomain pid:10924):ARMCI DASSERT fail. src/common/armci.c:ARMCI_Error():208 cond:0
--------------------------------------------------------------------------
MPI_ABORT was invoked on rank 0 in communicator MPI COMMUNICATOR 4 DUP FROM 0
with errorcode -1.
NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
You may or may not see output from other processes, depending on
exactly when Open MPI kills them.
--------------------------------------------------------------------------
It is a convergence issue, right? I was thinking that the problem was going to be solved using more iterations for scf: (lets say, maxiter 999), but, still the convergence problem is not fixed.
What is the correct way to force the scf saddle task to have more iteration in order to reach a convergence?
|
|
Henrique C. S. Junior
Inorganic Chemist - UFRRJ - Rio de Janeiro
|
|
|
|
Gets Around
Threads 18
Posts 133
|
|
| 1:22:55 PM PST - Fri, Jan 9th 2015 |
|
You have increased the allowable iterations for SCF electronic convergence but you need to increase the allowable iterations for geometry optimization. You should also check your transition state solution by running a frequency calculation and looking for a single imaginary frequency. To get a good TS solution you need to make the geometry convergence tighter than default. Your SCF threshold was too loose for this job. Here's an input file that should work better:
charge 0
memory 3000 mb
geometry autosym
S -1.139420 -0.181760 -0.065320
O -0.603930 0.894620 -0.859210
O -0.327300 -1.144590 0.633970
O -2.146750 0.489400 0.987400
O -2.158340 -1.022960 -0.975080
H -2.473210 1.320810 0.584900
H -2.293600 -1.887270 -0.534180
end
basis "ao basis" spherical
* library 6-31G
end
scf
singlet
thresh 1e-4
end
driver
tight
maxiter 999
end
task scf saddle
task scf frequencies
|
|
Edited On 1:27:26 PM PST - Fri, Jan 9th 2015 by Mernst
|
|
|
|
|
|
Clicked A Few Times
Threads 10
Posts 18
|
|
| 3:38:30 PM PST - Fri, Jan 30th 2015 |
|
Thank you for your reply and sorry for my delay (I went on vacation).
So, this block
driver
tight
maxiter 999
end
Is where I can set more iterations for geometry optimization?
|
|
|
|
|
|
Gets Around
Threads 18
Posts 133
|
|
| 1:03:07 PM PST - Sun, Feb 1st 2015 |
|
Yes, you can set more iterations for geometry optimization in the driver block. Or if you are using the non-default stepper geometry optimizer, you can do the same thing in a stepper block.
Also the default scf thresh value is 1e-4, as I made explicit above; I'm not sure why it was loosened in the original input file. Maybe Gabedit put that in for you.
|
|
|
AWC's:
2.5.10 MediaWiki - Stand Alone Forum Extension
Forum theme style by: AWC
