Please, see the following manuscript: Understanding the local reactivity in polar organic reactions through electrophilic and nucleophilic Parr functions. RSC Adv. 3, 1486-1494 (2013)
Hi, Multiwfn program (http://sobereva.com/multiwfn) can easily do this. A practial example is given in Section 4.5.4 of the manual. You only need to provide .wfn files (or .fch, .wfx, .molden files) for the related two states to Multiwfn.
The program DMol3 from Accelrys, Inc can do this, as described in https://www.researchgate.net/publication/200612981_On_the_use_of_fractional_charges_for_computing_Fukui_functions?ev=prf_pub
Article On the use of fractional charges for computing Fukui functions
Gabedit is capable to display fukui functions from molecular orbitals save as cube files or other formats. Then, you can obtain the square of homo or lumo to obtain the respective fukui funtions, as indicated in literature. In addition, Gabedit is capable to obtain fukui function in a extended squeme.
UCA-FUKUI software is a good supporting software. Nice to share an outstanding video about the combination of UCA-FUKUI, Gaussian, and Gaussview to calculate Fukui functions as follow:
You can use Gaussian 16 W to calculate FUKUI functions. Either Mulliken charges or Hirshfeld charges can be used to calculate the FUKUI function. But the Hirshfeld charges based FUKUI functions would give more accurate results.
You can compute the electrophilic and nucleophilic Parr functions, which are more accurate than the Fukui functions and more easy to obtained (RSC Adv.2013,3,1486).
You need only to do single point calculations of the radical anion and radical cation of the molecular optimized geometry, and analyse the Mulliken atomic spin densities.
For organic and small molecules the procedure might be like this -compute the single points in the neutral, radical cation and radical anion species and consider the relative density differences (form the cube files) . Just look for the populations ( say NPA, Hirschfeld or simply Mulliken) in the respective atoms in all the cases and derive fukui functions( electrophilic or nucleophilic) from the formulae set available. Once the electrophilic and nucleophilic functions are computed one can calculate other functions based on them (just school level mathematics). The use of specific software to generate the fukui functions is really not important it you have a electronic structure calculations program like orca, gausian etc.
Diganta Raychaudhuri Can you provide a literature reference for the fukui functions (and their derivative functions)? It would be of great help for a newcomer like me
Any ab initio software you can use depending on the systems of interest. From the calculation results you can calculate the respective fukui functions from the HOMO-LUMO gap itself. The attached paper could be of help.