Thank you very much! Everything is working now!

Best regards,
Timofey

This problem has been fixed, please download the latest version of Multiwfn.

Best  regards,

Tian

Thanks for the implementation of Local Hartree-Fock exchange energy. I have performed integration over the whole space for the open-shell He atom. Unfortunately, the result is 2 times more than the HF exchange energy given by Gaussian. At the same time, for the open-shell Be atom, the energy in Mulfiwfn is the same as in Gaussian.

Here, the example UHF-He.wfn where the problem is:
UHF-He.wfn

That is RHF-He.wfn results:
RHF-He.wfn

$ echo "1000 2 999 100 4 100" | tr " " "\n" | Multiwfn UHF-He.wfn 2>/dev/null | grep "Final result"
Final result:       -2.0409120547
$ echo "1000 2 999 100 4 100" | tr " " "\n" | Multiwfn RHF-He.wfn 2>/dev/null | grep "Final result"
Final result:       -1.0204560273

And two more examples for Be atom:
UHF-Be.wfn
RHF-Be.wfn

$ echo "1000 2 999 100 4 100" | tr " " "\n" | Multiwfn UHF-Be.wfn 2>/dev/null | grep "Final result"
Final result:       -2.6674331900
$ echo "1000 2 999 100 4 100" | tr " " "\n" | Multiwfn RHF-Be.wfn 2>/dev/null | grep "Final result"
Final result:       -2.6674331900

Best regards,
Timofey

Thank you so much! Now, it is working!

Best regards,
Igor

I just updated again, this time it should work.

Best regards,

Tian

I downloaded most recent //www.umsyar.com/multiwfn/misc/Multi … _noGUI.zip. Unfortunately, it continues to produce NaN local Hartree-Fock energy with

$ echo "1000 2 999 100 4 100" | tr " " "\n" | ./Multiwfn Al_large.wfn 2>/dev/null | grep "Final result"
 Final result:                 NaN

Additionally, running that version of MultiWFN with attached to archive settings.ini leads to error:

forrtl: severe (59): list-directed I/O syntax error, unit -5, file Internal List-Directed Read
Image              PC                Routine            Line        Source             
Multiwfn           0000000001F13C4B  Unknown               Unknown  Unknown
Multiwfn           0000000001F4037E  Unknown               Unknown  Unknown
Multiwfn           0000000001F3F0A5  Unknown               Unknown  Unknown
Multiwfn           000000000079DFBF  Unknown               Unknown  Unknown
Multiwfn           0000000000802012  Unknown               Unknown  Unknown
Multiwfn           000000000042B022  Unknown               Unknown  Unknown
libc-2.31.so       00007FC51DAFD0B3  __libc_start_main     Unknown  Unknown
Multiwfn           000000000042AF29  Unknown               Unknown  Unknown

Best regards,
Igor

I have fixed this problem, please download latest version.

Best regards,

Tian

While we tested a new code, we found that some basis sets lead to unreliable results while local Hartree-Fock energy is integrated over the whole space.

Here, the example .wfn where the problem is:
Al_large.wfn

At the same time, that .wfn have reliable energy:
Al_small.wfn

With binaries, provided in www.umsyar.com, we got the following energies:

$ echo "1000 2 999 100 4 100" | tr " " "\n" | ./Multiwfn Al_large.wfn 2>/dev/null | grep "Final result"
 Final result:                 NaN
$ echo "1000 2 999 100 4 100" | tr " " "\n" | ./Multiwfn Al_small.wfn 2>/dev/null | grep "Final result"
 Final result:      -18.0695836639

Best regards,
Igor

I forgot to mention, if this function is utilized in research, please not only cite Multiwfn but also cite original paper of Libreta liraray: J. Chem. Theory Comput., 14, 572−587 (2018). The developer of this library provided help about modifying Libreta interface during my implementation of this function.

Best regards,

Tian

Thank you so much!

Now, that is working!

Best regards,
Igor

I just fixed this issue and updated Multiwfn. Now you can also use 1000 to choose local HF exchange energy.

Best regards,

Tian

I found the problem. I do not use settings.ini.

I run Multiwfn as:

./Multiwfn H2O.wfn
1000
2
999
100
4
100

In that case, libreta is not initialized and genPprim routine is not called, so, Pprim_tot is not initialized.

I can write a code in 1000 section, but I'm not sure that it is a good idea since at the next time someone will need to fix several places.
What will be a better place for the correct initialization of needed libraries and matrixes?

Best regards,
Igor

At least for my tested system (examples\dopamine.wfn), both official Windows and Linux version (in RHEL6 U1) of Multiwfn work normally.
In order to implement this function, the interface of libreta library has been largely reorganized, make pay attention to proper compliation, the Makefile and "COMPLIATION_METHOD.txt" in source code package have been modified.

Best regards,

Tian

Thank you so much!

Unfortunately, I got SIGSEGV: Segmentation fault - invalid memory reference.

I will look what is the problem with my input file or with the implementation or with the compiler.

Best regards,
Igor

I am glad to inform you that this function has been supported as the 999th user-defined function in the Multiwfn updated today, see corresponding part of Section 2.7 of Multiwfn manual for detail.

Note that the -0.5 in the expression should be -1/4 for closed-shell case, this can be easily demonstrated. Also, the code segment you mentioned only works for closed-shell system, while the local HF exchange energy in Multiwfn also support open-shell system.

PS: To validate my implementation, I have performed integration for this real space function over the whole space using subfunction 4 of main function 100, the result is exactly identical to the HF exchange energy given by Gaussian.

Best regards,

Tian