“…If we consider the KID (for Koopmans in DFT) procedure presented in our previous works [6][7][8][9][10] together with a finite difference approximation, then the global reactivity descriptors can be calculated starting from the orbital energies of the HOMO and LUMO: electronegativity χ [19,20], global hardness η [19,20], electrophilicty ω [21], electrodonating ω − [22], electroaccepting power ω + [22], and net electrophilicity Δω ± [23]. Applying the same ideas related to the KID procedure, the local reactivity descriptors can be expressed as follows: [24][25][26][27][28][29][30], nucleophilic Parr function: P − (r) � ρ rc s (r) [31,32], and electrophilic Parr function: P + (r) � ρ ra s (r) [31,32], where ρ N+1 (r), ρ N (r), and ρ N−1 (r) are the electronic densities at point r for the system with N + 1, N, and N − 1 electrons, respectively, and ρ rc s (r) and ρ ra s (r) are related to the atomic spin density (ASD) at the r atom of the radical cation or anion of a given molecule, respectively [33].…”