Multiwfn forum

Multiwfn official website: http://sobereva.com/multiwfn. Multiwfn forum in Chinese: http://bbs.keinsci.com/wfn

You are not logged in.

#1 2025-06-07 10:28:09

smile
Member
Registered: 2025-04-15
Posts: 7

Help with Displaying ESP Extrema Values on van der Waals Surface

Dear Prof. Tian Lu,

I hope you're doing well.

I’m visualizing MEP mapped on the van der Waals surface using the suranalysis.pdb file. The colored potential surface, isosurface, extremas sheres, and CPK model all display correctly — but I’m having trouble with the clear visualization of the ESP extrema values, which I need for a publication-quality figure.

Here are the main issues and questions:

I’ve tried displaying all extrema, as well as only the top 3 or 5 positive and negative. But in all cases, the values either appear buried inside the surface or are overlapping heavily, making them unreadable and cluttered.

I’m unsure how extrema are typically selected for display in articles:

Based on threshold values (e.g., >±10, 15, or 20)?

Based on proximity to functional groups, heteroatoms, or hydrogens?

Or purely based on visual clarity?

Is the selection typically done manually (by removing or hiding labels), or is there a script-based way to handle it?

In some publications, I've seen the use of the “backside” label for certain extrema. When I tried this, the label appeared on all values — is it usually manually assigned?

Some articles show 2D arrows pointing to hidden extrema (e.g., on the rear side of the surface). Are these added inside VMD, or with external tools like Illustrator? Would that affect the image resolution?

I also attempted to use cones in VMD to point toward buried extrema but couldn't get it to work well — any advice?

For label coloring, should all values use a single color (e.g., black), or should positive and negative values be colored differently (e.g., red and blue)?

Everything else in the visualization looks good — the ESP zones, isosurface, and molecular model all render as expected. I just need guidance on how to clearly and selectively display the ESP extrema values to finalize the figure for publication.

Thank you very much for your time and support.

Best regards,

Offline

#2 2025-06-07 16:49:56

sobereva
Tian Lu (Multiwfn developer)
From: Beijing
Registered: 2017-09-11
Posts: 2,076
Website

Re: Help with Displaying ESP Extrema Values on van der Waals Surface

Please upload figure(s) to describe the issues you encounted. You raised many questions, but without figures, I can hardly understand your problems and difficulties.

BTW: In the recent Multiwfn introduction paper J. Chem. Phys., 161, 082503 (2024) DOI: 10.1063/5.0216272, Fig. 14 is a very representative ESP colored vdW surface map, I think the graphical effect is good enough.

Offline

#3 2025-06-09 09:16:39

smile
Member
Registered: 2025-04-15
Posts: 7

Re: Help with Displaying ESP Extrema Values on van der Waals Surface

Dear Prof. Tian Lu,

Thank you for your response.

As requested, I’ve attached figures to illustrate the problems I’m encountering:

Figure 1: This is taken from a reference article. The ESP extrema values are very clearly displayed, well separated, and placed on top of the surface, with labels such as “backside” and arrows used to indicate positions that are not directly visible. Global extrema are also marked with an asterisk *.

Figure 2: This is my result using the same molecule as in Figure 1, for comparison. As you can see, the labels are overcrowded, many are stacked on top of each other, and several are buried inside the surface, making them unreadable. Even limiting the display to top 3 or 5 extrema doesn't resolve this — they still appear unclear or inside the surface.

Figure 3: This is from another publication. This author displayed many extrema values, not just 8 like in Figure 1, but the figure remains clean and readable, with good use of spacing, arrows, and asterisks for emphasis.

Based on these comparisons, I would like to ask:

How can I achieve a clear figure like in these references, where extrema are clearly shown on the surface, not buried inside?

Do I need to limit the number of displayed values manually to avoid overcrowding, or is there a more systematic way to select which extrema to show?

How are elements like “backside” labels, arrows, and asterisks typically added — in VMD or via post-processing tools?

Any guidance or recommendations you can provide would be greatly appreciated.

Best regards,

fig 1par1.png
fig 2 fig-2.jpg
fig 3 fig-3.jpg

Offline

#4 2025-06-09 11:10:35

sobereva
Tian Lu (Multiwfn developer)
From: Beijing
Registered: 2017-09-11
Posts: 2,076
Website

Re: Help with Displaying ESP Extrema Values on van der Waals Surface

The labels in Fig. 1 are added manually using post-processing tool, such as Photoshop and powerpoint. I always add labels in this way.

Offline

#5 2025-06-14 13:34:37

smile
Member
Registered: 2025-04-15
Posts: 7

Re: Help with Displaying ESP Extrema Values on van der Waals Surface

Thank you for the clarification about adding labels manually.

I’d like to ask one more thing: in your publications, how do you decide which extrema points to show and which to omit?
Do you apply a threshold (e.g., only showing ESP values >±15), or do you choose a fixed number of top values (e.g., top 3 or 5 positive/negative)? Or is the selection purely based on visual clarity to avoid overcrowding the figure?

Also, do you filter those points manually, or do you use a script to select and reduce them before labeling?

In my case, for example, I noticed that in the previous reference pic (Fig. 1), the author displayed only 8 extrema values clearly and cleanly, while in my attempt (Fig. 2) using the same molecule, many values appeared — often overlapping or placed on top of each other.
So just to confirm: do you mean it’s fully optional for the user to select which extrema to keep and display, purely based on clarity and visual presentation, not necessarily based on a strict scientific rule?

Your advice would really help me finalize my figure. Thank you again.

Offline

#6 2025-06-14 16:27:12

sobereva
Tian Lu (Multiwfn developer)
From: Beijing
Registered: 2017-09-11
Posts: 2,076
Website

Re: Help with Displaying ESP Extrema Values on van der Waals Surface

I prefer to display all surface extrema, but only manually label some of them, namely the relative important ones (i.e. the ones with relative large magnitude of ESP)

Offline

#7 2025-09-18 17:04:05

smile
Member
Registered: 2025-04-15
Posts: 7

Re: Help with Displaying ESP Extrema Values on van der Waals Surface

Subject: Help matching ESP extrema values to positions

Dear Prof. Tian Lu,

I generated the ESP extrema spheres in Multiwfn and viewed them in VMD.
I tried using a Tcl script to print each sphere’s coordinates and values, but I only get errors like “expected integer but got '' ”.

Could you please tell me how to match each sphere with its ESP value, so I can label them manually in Microsoft PowerPoint as you suggested?

Thank you very much.
Best regards,

Offline

#8 2025-09-18 18:27:30

sobereva
Tian Lu (Multiwfn developer)
From: Beijing
Registered: 2017-09-11
Posts: 2,076
Website

Re: Help with Displaying ESP Extrema Values on van der Waals Surface

Hello,

If you find the extreme of interest is index 3 via inquire mode of VMD, then you can input

[atomselect top "index 3"] get "x y z beta"

To print its X,Y,Z coordinate and ESP value (i.e. Beta field in surfanalysis.pdb)

Best regards,

Offline

#9 Yesterday 15:06:59

smile
Member
Registered: 2025-04-15
Posts: 7

Re: Help with Displaying ESP Extrema Values on van der Waals Surface

Thank you for your quick reply
I am now calculating grid data of mapped function (MEP) in Multiwfn 3.8 for a 57-atom molecule, grid spacing 0.2. In my previous try, it reached 50% after 8 hours, but I stopped it because I thought there was a problem.
Can I make it faster?If the calculation takes more than 16 hours, are the results still correct and reliable for analysis?

Best regards,

Offline

#10 Yesterday 15:11:44

sobereva
Tian Lu (Multiwfn developer)
From: Beijing
Registered: 2017-09-11
Posts: 2,076
Website

Re: Help with Displaying ESP Extrema Values on van der Waals Surface

Without more information (all commands inputted in Multiwfn, operation system, basis set, actual number of CPU cores used), I can hardly answer your question.

Offline

#11 Yesterday 16:29:17

smile
Member
Registered: 2025-04-15
Posts: 7

Re: Help with Displaying ESP Extrema Values on van der Waals Surface

Here are the details of my calculation:
Multiwfn version: 3.8
Commands used / steps:
Loaded .fchk file in the Win exe
Chose 12 for quantitative MEP analysis
Set grid spacing 0.2
Pressed 0 to start calculation
After completion, pressed 5 to get the PDB file
Pressed 2 to get surf.cub
Pressed 13 to get the grid data of mapped function (this step is the longest and causing the problem)
Operating system: Windows 10
Basis set / Method: DFT B3LYP/6-31+G(d,p)
CPU: Intel Core i3 (2 physical cores, 4 logical processors), but Multiwfn is running with the default single core
Molecule size: 56 atoms
Grid spacing: 0.2 Å

Thank you very much for your help.
Best regards,

Offline

#12 Yesterday 22:09:27

sobereva
Tian Lu (Multiwfn developer)
From: Beijing
Registered: 2017-09-11
Posts: 2,076
Website

Re: Help with Displaying ESP Extrema Values on van der Waals Surface

By default Multiwfn uses 4 cores (i.e. the default "nthreads" in settings.ini)
You should set "nthreads" in settings.ini to 2 (equal to number of physical cores), if the cost is still too high, please use a computer with better CPU (>=8 cores), or remove diffuse function when generating wavefunction file to save subsequent cost in ESP related analysis (diffuse functions are not quite important for ESP analysis).

Offline

Board footer

Powered by FluxBB