“…Many efforts have been made to model accute toxicity of NACs applying a set of molecular descriptors, such as the energies of the highest occupied (E HO-MO ) and the lowest unoccupied (E LUMO ) molecular orbitals, the LUMO-HOMO energy gap (∆E H-L ) (or chemical hardness, η = E LUMO -E HOMO , or softness, S = 1/η), vertical and/or adiabatic ionization potentials (IPs) and electron affinities (EAs), electro-static densities, charges and potentials, molecular area and/or volume, (topological) polar surface area ((T)PSA), the numbers of hydrogen-bond donor (HBD) or hydrogen-bond acceptor (HBA), lipophilicity (log P and/or log D), a set of topological indices, etc. [13][14][15][16][17]. In addition, the efforts have been made to explore the DFT-based reactivity indices of nitroaromatics for prediction of their toxicity, viz., the chemical hardness (η), the chemical potential (µ) (or electronegativity, χ = -µ) and the electrophilicity (ω) ( [18], and refs.…”