The energy decomposition consists of estimating the interaction energy either per-residue or per-wise. In the first case, the DELTA energy corresponds to the sum of the interaction energy of that residue with the rest of the selected residues and itself. In the second case, it is computed in pairs, i.e., the interaction energy between two residues. This is an approximation, which is not equivalent to the relative binding free energy value since it is computed only between the selected residues. The decomposition aims to estimate which residues are apparently more important for the interaction, useful when generating mutants (e.g. peptide design, in silico affinity maturation, etc.) In your case, the energy value reported for the ligand corresponds to the sum of the interaction energies between the ligand and the selected residues and itself.
According to one of the figures, the interaction energy is approximately 40 kcal/mol. If you add up the interaction energy of all the residues shown in the decomposition plot, you get -41.28 kcal/mol, which is expected.
There is no way to remove the ligand from the graph in gmx_MMPBSA_ana, mainly because it would not be correct for the analysis. If your goal is to show only the receptor residues, edit the csv and use your favorite graphing program/library. If you have any doubts, I can help you set up the graph in the same format as shown in gmx_MMPBSA_ana.
To perform the decomposition analysis at least one residue of the ligand and the receptor is required, otherwise, there is no way to calculate the interaction energy, as the name suggests.
Let me know any doubt
Mario S.