What is the keyword in NWCHEM to calculate reorganization energies

From NWChem

Viewed 198 times, With a total of 4 Posts
Jump to: navigation, search

Forum Regular
Threads 54
Posts 248
Dear All
Since 2016 or perhaps 2014, GAMESS has added the function to calculate reorganization energies of ET and proton-coupled ET reactions, published, respectively, until 2014, where one of the two necessary steps may fail on an old cpu, in the cases of electrode plus solvent, electrode plus self assembled monolayer and solvent, electrode plus double layer and ionic solution, and homogeneous solvent, respectively, which also has an equivalent Python program. Isn't it fascinating that the electrochemical kinetic parameters can be obtained from the free energies in a solvent though quantum chemical methods, e.g., DFT, without the intricate treatment of electrodes and purification of solvents? These original programs were contributed to GAMESS by Profs. and Drs. previously in the world-wide famous Department of Chemistry of University of Illinois at Urbana-Champaign and evaluated aassembled into the existing solvation module of gamess by Profs and Drs. in Prof. Gordon's Group and Ames Laboratory. According to the authors' paper on SAM monolayer modified electrodes, they can produce total reorganization energies in good agreement with experimental data for all 12 metal complexes tested in aqueous solution, useful for molecular electrocatalysts electron transfer rate constants calculation.
I am now enjoying it for biochemical applications, e.g., cofactors, free radicals, etc., in aqueous and aprotic solutions..


It is said in the NWCHEM manual that it also can calculate reorganization energies, but I cannot find any keyword or example about this.
How to do it in NWCHEM?
Thanks a lot!

Very Best Regards!
Edited On 7:36:20 PM PST - Sun, Feb 17th 2019 by Xiongyan21

Forum Regular
Threads 54
Posts 248
First, I need calibrate it using published experimental data or theoretical results. I am very pleased to find the calculated average reorganization energies of hydroquinone oxidation and the corresponding reduction state in the situation of electrode plus double layer and an ionic aprotic solution, the sole environment chosen for now, is only around 0.09 eV deviating from the outer sphere one obtained with a different method in a published paper in 2004, one of whose author was at EMSL then, using Gaussian98, with the same functional, i.e., B3LYP in the VMN3 case in GAMESS without TDA, although the solvent there may be different from what I have tested. Perhaps the new module is more accurate in some cases.


Now I may know the procedure in NWCHEM but I haven‘t tried.
Edited On 7:26:19 PM PST - Sun, Feb 17th 2019 by Xiongyan21

Forum Regular
Threads 54
Posts 248
Thanks a lot for a professor's answer on Github.
Now I can quantum chemically get ET rate constants from MARCUS's theory using GAMESS and NWCHEM perhaps integrately.
Edited On 6:16:33 AM PST - Tue, Jan 29th 2019 by Xiongyan21

Forum Regular
Threads 54
Posts 248
MP2 6-31g(d,p) MP2 optimized hydroquinone has no warning for HESSIAN using GAMESS whose outer reorganization in another solvent is within 0.05 eV deviating from the previously mentioned article.
Edited On 6:25:46 AM PST - Sun, Mar 3rd 2019 by Xiongyan21

Forum Regular
Threads 54
Posts 248
The dftb3 optimized hydroquinoe has no warning for the HESSIAN using GAMESS, and the outer reorganization enegy in another solvent is in good agreement with that in the previously metioned article, i.e., the difference is within 0.06eV.
Edited On 6:26:15 AM PST - Sun, Mar 3rd 2019 by Xiongyan21


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