This shows you the differences between two versions of the page.
Next revision | Previous revision | ||
scissor_operator [2011/12/05 10:10] prokop created |
scissor_operator [2011/12/06 13:19] (current) |
||
---|---|---|---|
Line 9: | Line 9: | ||
===== Free molecule computation ===== | ===== Free molecule computation ===== | ||
- | As example you can use fireball.in like this | + | As example you can use fireball.in like this ( importaint is to se iwrtcdcoefs = 1 ) |
<code> | <code> | ||
&OPTION | &OPTION | ||
basisfile = freemol.bas | basisfile = freemol.bas | ||
+ | lvsfile = cel.lvs | ||
+ | icluster = 0 | ||
+ | nstepf = 1 | ||
+ | sigmatol = 0.000001 | ||
+ | max_scf_iterations = 100 | ||
+ | iqout = 0 | ||
+ | ifixcharge = 0 | ||
+ | &END | ||
+ | |||
+ | &OUTPUT | ||
+ | iwrtcdcoefs = 1 | ||
+ | &END | ||
+ | </code> | ||
+ | |||
+ | <code> | ||
+ | 12 | ||
+ | 8 -0.000000 0.006964 -0.307366 | ||
+ | 6 0.000000 -0.001848 1.561491 | ||
+ | 1 0.000000 -2.253742 1.713401 | ||
+ | 1 0.000000 2.250737 1.711678 | ||
+ | 6 0.000000 -1.242741 2.296730 | ||
+ | 6 0.000000 1.240398 2.295749 | ||
+ | 6 -0.000000 -1.240740 3.721606 | ||
+ | 6 -0.000000 1.239884 3.720292 | ||
+ | 1 -0.000000 -2.253492 4.301542 | ||
+ | 1 -0.000001 2.253787 4.299402 | ||
+ | 6 -0.000000 0.000210 4.460464 | ||
+ | 8 0.000000 0.000069 6.330528 | ||
+ | </code> | ||
+ | |||
+ | ===== SYSTEM computation ===== | ||
+ | |||
+ | - copy cdcoeffs.dat from molecule computation to the SYSTEM directory | ||
+ | - Set koopmans shifts for orbitals of molecule. To do this you should specify start and end atom of molecule inside the system geometry file. Ten set number of molecular orbitals (the same number as first number of molecular cdcoefs.da) to be shifted and specify shifts for each of the orbitals. | ||
+ | |||
+ | koopman_shift | ||
+ | <code> | ||
+ | 19 first_atom | ||
+ | 30 last_atom | ||
+ | 76 orbitals | ||
+ | -1.00000 | ||
+ | -1.00000 | ||
+ | -1.00000 | ||
+ | -1.00000 | ||
+ | -1.00000 | ||
+ | -1.00000 | ||
+ | -1.00000 | ||
+ | -1.00000 | ||
+ | -1.00000 | ||
+ | -1.00000 | ||
+ | -1.00000 | ||
+ | -1.00000 | ||
+ | -1.00000 | ||
+ | -1.00000 | ||
+ | -1.00000 | ||
+ | -1.00000 | ||
+ | -1.00000 | ||
+ | -1.00000 | ||
+ | -1.00000 | ||
+ | -1.00000 | ||
+ | 1.00000 | ||
+ | 1.00000 | ||
+ | 1.00000 | ||
+ | 1.00000 | ||
+ | 1.00000 | ||
+ | 1.00000 | ||
+ | 1.00000 | ||
+ | 1.00000 | ||
+ | 1.00000 | ||
+ | 1.00000 | ||
+ | 1.00000 | ||
+ | 1.00000 | ||
+ | 1.00000 | ||
+ | 1.00000 | ||
+ | 1.00000 | ||
+ | 1.00000 | ||
+ | 1.00000 | ||
+ | 1.00000 | ||
+ | 1.00000 | ||
+ | 1.00000 | ||
+ | 1.00000 | ||
+ | 1.00000 | ||
+ | 1.00000 | ||
+ | 1.00000 | ||
+ | 1.00000 | ||
+ | 1.00000 | ||
+ | 1.00000 | ||
+ | 1.00000 | ||
+ | 1.00000 | ||
+ | 1.00000 | ||
+ | 1.00000 | ||
+ | 1.00000 | ||
+ | 1.00000 | ||
+ | 1.00000 | ||
+ | 1.00000 | ||
+ | 1.00000 | ||
+ | 1.00000 | ||
+ | 1.00000 | ||
+ | 1.00000 | ||
+ | 1.00000 | ||
+ | 1.00000 | ||
+ | 1.00000 | ||
+ | 1.00000 | ||
+ | 1.00000 | ||
+ | 1.00000 | ||
+ | 1.00000 | ||
+ | 1.00000 | ||
+ | 1.00000 | ||
+ | 1.00000 | ||
+ | 1.00000 | ||
+ | 1.00000 | ||
+ | 1.00000 | ||
+ | 1.00000 | ||
+ | 1.00000 | ||
+ | 1.00000 | ||
+ | 1.00000 | ||
+ | </code> | ||
+ | |||
+ | then you can run computation with (igap=3), see this fireball.in as a reference | ||
+ | |||
+ | fireball.in | ||
+ | <code> | ||
+ | &OPTION | ||
+ | basisfile = answer_small.bas | ||
lvsfile = cel.lvs | lvsfile = cel.lvs | ||
icluster = 0 | icluster = 0 | ||
Line 23: | Line 147: | ||
ismeagol=0 | ismeagol=0 | ||
ifixcharge = 0 | ifixcharge = 0 | ||
+ | igap = 3 | ||
&END | &END | ||
&OUTPUT | &OUTPUT | ||
- | iwrtcdcoefs = 1 | ||
iwrtHSrho = 0 | iwrtHSrho = 0 | ||
iwrteigen = 0 | iwrteigen = 0 | ||
Line 32: | Line 156: | ||
&END | &END | ||
</code> | </code> | ||
- | |||