From NWChem
You are viewing a single post from the thread title above
-
Huub Forum:Admin, Forum:Mod, NWChemDeveloper, bureaucrat, sysop
|
|
|
Forum Regular
Threads 1
Posts 185
|
|
| 6:53:11 PM PDT - Tue, Jul 1st 2014 |
|
Hi,
Indeed the code does not print the HOMO-LUMO gap explicitly. However, at the end of the SCF or DFT calculation it will print the orbital analysis. For H2 (I have used SCF instead of DFT so I could simply use "task ccsd energy", and I used the cc-pvdz basis set) I got:
Vector 1 Occ=2.000000D+00 E=-5.924110D-01 Symmetry=a1g
MO Center= 0.0D+00, 0.0D+00, -4.0D-17, r^2= 7.1D-01
Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
----- ------------ --------------- ----- ------------ ---------------
1 0.403665 1 H s 6 0.403665 2 H s
2 0.173763 1 H s 7 0.173763 2 H s
Vector 2 Occ=0.000000D+00 E= 1.974400D-01 Symmetry=a2u
MO Center= 3.3D-17, 6.2D-17, -9.0D-16, r^2= 2.7D+00
Bfn. Coefficient Atom+Function Bfn. Coefficient Atom+Function
----- ------------ --------------- ----- ------------ ---------------
2 1.958656 1 H s 7 -1.958656 2 H s
1 0.150244 1 H s 6 -0.150244 2 H s
In this table "Occ" refers to the orbital occupation and "E" to the corresponding orbital energy. So it is relatively easy to identify vector 1 as the HOMO and vector 2 as the LUMO, giving a HOMO-LUMO gap of 0.19744 - -0.59241 = 0.78985 Hartree. Of course if you want you can write a little script to extract that information. I am more than happy to put that in our contrib directory, so you have it in future.
Best wishes, Huub
|
|
|
AWC's:
2.5.10 MediaWiki - Stand Alone Forum Extension
Forum theme style by: AWC
