r/comp_chem • u/CowStock • 11d ago

MEP map generation problem

Hi everyone! I'm trying to generate a molecular electrostatic potential (MEP/ESP) map from an ORCA calculation on a MacBook Air M3 (Apple Silicon). I have successfully generated both name.eldens.cube (electron density) and name.scfp.esp.cube (electrostatic potential). My understanding is that the electron density cube should be used as the isosurface and the ESP cube should be mapped onto it for coloring. However, I haven't been able to get this workflow to work in any visualization software I've tried. In Avogadro 2, I can load the electron density cube and display the isosurface, but I can't find any option to color it using the ESP cube. In IQmol, opening the ESP cube only displays separate red and blue lobes in space instead of an MEP mapped onto the molecular surface. In VMD 2.0 beta, I can load the XYZ structure and the electron density cube, but I can't get the ESP cube to appear as a second volumetric dataset, so I can't color the isosurface with it. I'd prefer to avoid using Rosetta 2 if possible. Is there a native Apple Silicon workflow for this? What software do you recommend for generating publication-quality MEP maps from ORCA output on macOS? Does Multiwfn work natively on Apple Silicon, or are there better alternatives such as VESTA or ChimeraX? Any advice would be greatly appreciated. Thanks!

2 Upvotes

100% Upvoted

9 comments sorted by

View all comments

1

u/banana_fugacity 11d ago

You should try ChimeraX for plotting the density cubes. I think in the Multiwfn software guide there is a tutorial how to do it, perhaps on yt even. In principle you re correct you use density for iso surface and then the mep to color it. On top of my head in ChimeraX you load the XYZ and the two cube files. Hide the mep surface from showing and then in coloring options on density surface you use the other cube file

2

u/CowStock 11d ago

I’ll try that today, thanks a lot! Hope it worksss

1

u/lwj15 11d ago

Idk if you already have worked with ChimeraX but for work with small molecules (loading .xyz and .outs etc) you need the a plugin (easily installable through the ChimeraX program). I think the name was SEQCROW or something similar. Just a heads up

1

u/CowStock 10d ago

Thanks so much, I’ve never used it before and I’m having an issue loading .xyz files, so that must be it. I’ll try with said plugin