This shows you the differences between two versions of the page.
Both sides previous revision Previous revision Next revision | Previous revision Next revision Both sides next revision | ||
calculate_u-function [2012/11/07 13:12] vlada |
calculate_u-function [2012/11/07 13:36] vlada |
||
---|---|---|---|
Line 3: | Line 3: | ||
This page shows how to calculate U - function. | This page shows how to calculate U - function. | ||
There is an simple example of calculation U -function for simple CH5N molecule below. | There is an simple example of calculation U -function for simple CH5N molecule below. | ||
- | Here is an input bas file. | + | Here are input files. |
- | CH5N.bas | + | ---- |
+ | CH5N.bas | ||
<code> | <code> | ||
7 | 7 | ||
Line 17: | Line 18: | ||
1 0.401767 -0.518300 0.926463 | 1 0.401767 -0.518300 0.926463 | ||
</code> | </code> | ||
+ | Calculation of U-function is turned on by iu = 1 in fireball.in. | ||
fireball.in | fireball.in | ||
Line 39: | Line 41: | ||
&END | &END | ||
</code> | </code> | ||
+ | Next you need input file u.inp, where you have to write total number of electron. You can choose molecular orbitals for which you want to calculate U. There are 5 orbitals (6,7,8,10,22) in this case. | ||
u.inp | u.inp | ||
<code> | <code> | ||
Line 49: | Line 51: | ||
10 | 10 | ||
22 | 22 | ||
+ | </code> | ||
+ | As an output file you get u_func.dat file where is writen number of selected orbital in the first column and in the second one are u values. | ||
+ | <code> | ||
+ | 6 2.1590 | ||
+ | 7 2.3341 | ||
+ | 8 1.9742 | ||
+ | 10 1.9721 | ||
+ | 22 1.4045 | ||
</code> | </code> |