Dear Professor Tian,
Thank you for the quick answer! I also wonder, is it possible to plot the second map (or even the first) using the main function 13 (Process grid data), so I would be grateful for your professional feedback
Best wishes,
Alex

Dear Alex,

Main function 13 (grid data processing function) is completely irrelevant to this purpose.

In main function 3 of Multiwfn, you can scan any real space function (e.g. ESP) along a given line, and you can manually define coordinate of the two ends of the line. Please check Section 4.3 of Multiwfn manual for example. According to the first map, one end should be set to nuclear position of Se, and another end should be set to the given point "d".

To reproduce the second map, you need to write your own script to generate a batch of points in the scanning path, then you can use the special feature of main function 5 (see "Special case: Calculate data for a set of arbitrarily distributed points" in Section 3.6 of Multiwfn manual) to calculate ESP values at all points at once.