mkccsd geometry optimization

From NWChem

Viewed 823 times, With a total of 2 Posts
Jump to: navigation, search

Just Got Here
Threads 1
Posts 1
Dear NWChem developers,

I've compiled 6.5 version but unfirtunatelly cannot run optimization at mkccsd level.
The same problem I've met with version 6.3 where after http://www.nwchem-sw.org/images/Mrcctmp.patch.gz
patching situation even situation became worse - it could not reach convergence with mkccsd for hte startin geometry


input:

echo
start tce_mrcc_mkcc
memory stack 1000 mb heap 100 mb global 500 mb verify
geometry units au
H 0.00000000 -2.27289450 -1.58834700
O 0.00000000 0.00000000 -0.01350000
H 0.00000000 2.27289450 -1.58834700
end
basis spherical
O library cc-pvtz file /work/02018/sokovyty/Nwchem-6.5/src/basis/libraries/cc-pvtz
H library cc-pvtz file /work/02018/sokovyty/Nwchem-6.5/src/basis/libraries/cc-pvtz
end
charge 0
scf
rohf
singlet
thresh 1.0e-10
tol2e 1.0e-10
end
tce
mkccsd
thresh 1.0e-10
targetsym a1
maxiter 500
io ga
tilesize 18
end
mrccdata
se4t
root 1
cas 2 2 # Please make sure the references generated are correct.
end

task tce OPTIMIZE


error massage:


tce_hash: key not found                   0
------------------------------------------------------------------------
------------------------------------------------------------------------
current input line :
34: task tce OPTIMIZE
------------------------------------------------------------------------
------------------------------------------------------------------------
This error has not yet been assigned to a category
------------------------------------------------------------------------

  • Karol Forum:Admin, Forum:Mod, NWChemDeveloper, bureaucrat, sysop
    Profile
    Send PM
Clicked A Few Times
Threads 1
Posts 30
MRCC formalisms have been designed for single point calculations only.

Best,
Karol

Gets Around
Threads 15
Posts 51
Quote:Karol Feb 13th 3:30 am
MRCC formalisms have been designed for single point calculations only.

Best,
Karol

It is not true. One can optimize geometry, saddle point and calculate frequencies with MRCC formalisms by using finite differences.
In this case I think MRCC iteration is unstable "Eigenvectors swaps", so you shell use another set of references orbitals (not cas 2 2).
Edited On 7:04:36 PM PST - Thu, Feb 12th 2015 by Vladimir


Forum >> NWChem's corner >> NWChem functionality



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