Multiwfn official website: //www.umsyar.com/multiwfn. Multiwfn forum in Chinese: http://bbs.keinsci.com/wfn
You are not logged in.
There is no internal limitation in Multiwfn. But if the file is large, while your physical memory is limited, or you are using Linux but you didn't properly configuration your system according to Section 2.1.2 of Multiwfn manual, loading the file may be crash. Please take care of the information shown on screen when Multiwfn crashes. In addition, if the file is corrupted, loading will also be failed.
First, you should change drawing method to "Points", rather than using the default "Line" style. Then you need to fine-tune coloring method, material, visible status and so on.
You didn't properly set up the drawing method in VMD. I strongly suggest using the method described in the following video tutorial to plot the AIM topology map via Multiwfn+VMD, it is very convenient and the graphical quality is ideal.
Using Multiwfn with VMD to quickly and easily plot high quality atoms-in-molcules (AIM) topology map
https://youtu.be/mgsnhvWH5SI
If you refer to sobEDA, there is no special point for host:guest system, just treat it as a molecular dimer.
Dear Saeed,
The ETS-NOCV in Multiwfn is also able to study Pauli repulsion interaction, but it doesn't represent the interaction as fragment orbital interactions.
In the CDA interface, you should also need to first identify the fragment orbitals that heavily participlate in the repulsion term. I mention "4 Export overlap matrix between fragment orbitals" because you do not need to use subfunction 15 of main function 100.
I don't know how did you represent fragments in IRC.log. If something like C(fragment=2) appears in the log file, GauIRC2xyz doesn't formally support this very special case, and I don't have intention to modify the code.
Best,
Tian
Dear Saeed,
Your description is correct. However, I would like to mention that in fact you can directly obtain overlap integral between fragment orbitals via option "4 Export overlap matrix between fragment orbitals" in Multiwfn CDA module, it is more convenient. Also it should be noted that magnitude of overlap integral is just a very crude estimation of contribution to Pauli repulsion, the value should not be overly discussed.
Best,
Tian
There is an important update of Multiwfn (Windows version only) on 2025-Apr-28, which makes changing viewpoint significantly easier!
Now in the GUI window of Windows version of Multiwfn, one can use the left mouse button to drag the system to freely rotate it. In addition, with holding the Ctrl key and dragging the system horizontally and vertically using left mouse button, one can rotate the system along the screen and zoom in/out the system, respectively. With holding the Shift key and dragging the system using left mouse button, the system can be translated in the drawing region. These improvements make changing viewpoint significantly easier!
Dear Saeed,
The CDA analysis in Multiwfn is able to provide overlap integral S(i,j) between fragment orbitals of the two fragments.
Best,
Tian
Hello,
Please take care of the following points:
(1) As you can see from prompt, Multiwfn command was not found, that means Multiwfn has not been properly installed. Please check Section 2.1.2 on how to install Multiwfn under Linux. You need to guarantee that Multiwfn can be normally invoked using "Multiwfn" command (or "Multiwfn_noGUI" command if you use noGUI version of Multiwfn).
(2) G09 is not formally supported by sobEDA.sh, using G16 is strongly recommended.
(3) I didn't test sobEDA.sh in WSL, if it doesn't work finally, please use Linux machine, or install a Linux guest system using VMware virtual machine in Windows.
You need to modify "subroutine outwfx" in source code file fileIO.f90 of Multiwfn and recompile.
In this subroutine, please note lines including "if (MOocc(imo)/=0D0) " and "if (MOocc(imo)==0D0)", you can easily realize your aim by properly modifying corresponding parts.
%Chk=DPB_gr.chk
# B3LYP/cc-pvtz Freq(ReadFC,FC,ReadFCHT) NoSymm SCRF(Solvent=Cyclohexane) Geom=Checkpoint Guess=Read
FC simulation of UV absorption spectrum for trans,trans-1,4-diphenyl-l,3-butadiene
0 1
SpecHwHm=400 SpecRes=20 InpDEner=0.133 Print=Matrix=JK
DPB_ex.chkDon't worry, you can input e.g. 0.002 in the text box of color scale of VMD and press ENTER button, it will take effect. Although then you will see 0.00, it doesn't matter, because the text box only show two decimal places.
Alternatively, you can use ESPiso_eV.bat along with ESPiso_eV.txt in "examples\drawESP\" instead of ESPiso.bat along with ESPiso.txt. In this case, the grid data of ESP exported by Multiwfn will be in eV, so you can input the color scale in VMD in a much larger magnitude.
Alternatively, you can use VMD command to directly set the color scale. In ESPiso.vmd you can find the following two lines, you can change the values and run them in VMD console window manually.
set colorlow -0.03
set colorhigh 0.03The new update of Multiwfn today brings an important new feature:
The electron density polarization analysis based on electron excitations proposed in J. Phys. Chem. A, 124, 633 (2020) has been implemented as subfunction 17 of main function 18. This method is able to provide very valuable insight into the nature of electron density polarization under an external perturbation (e.g. point charge), and can be used to study substitution effect, mechanism of electrophilic/nucleophilic reactions, atomic polarizability, and so on. See Section 3.21.17 of Multiwfn manual for introduction and Section 4.18.17 for example.
To force ORCA to optimize T1 in the case 2, you should add "irootmult triplet", then the optimization will only consider the triplet excited states.
I am not a NCIplot user. The NCI data exported by Multiwfn are exactly in line with the standard definition, see equations in Multiwfn manual (Section 3.23.1 of Multiwfn manual) and my reviews:
Tian Lu, Qinxue Chen, Visualization Analysis of Weak Interactions in Chemical Systems (2023), Comprehensive Computational Chemistry, Vol. 2 pp. 240-264. Oxford: Elsevier. DOI: 10.1016/B978-0-12-821978-2.00076-3 (preprint: //www.umsyar.com/attach/Visualizatio … ystems.pdf)
Tian Lu, Visualization Analysis of Covalent and Noncovalent Interactions in Real Space, ChemRxiv (2025) DOI: 10.26434/chemrxiv-2025-9t442
I also strongly suggest you reproducing the NCI examples given in Multiwfn manual, and then apply the analysis to your own system.
There is no any special prefactor like "100" in the outputted data of Multiwfn.
When you ask Multiwfn to export a text file, including the file exported by option "2 Output scatter points to output.txt in current folder", meaning of each column in the exported file is always very clearly shown on screen. Please carefully check these information on screen.
Transition density is very different to charge density difference, there is no direct relationship among them. Please check Section 3.21.1 of Multiwfn manual to understand definition of transition density.
To calculate electron density difference between S1 and T1 via Multiwfn, you should do:
(1) Using hole-electron analysis module (subfunction 1 of main function 18) to calculate electron density difference between S1 and S0, then export the grid data as CDD.cub and manually rename it as S1CDD.cub. See Section 4.18.1 of Multiwfn manual for example.
(2) Similarly, generate electron density difference between T1 and S0, and store as T1CDD.cub.
(3) Using Multiwfn to generate difference grid data between S1CDD.cub and T1CDD.cub (see the grid data mathematical operations exemplified in Section 4.13.2 of Multiwfn manual), then the grid data in memory will corresponds to electron density difference between S1 and T1. Then you can directly visualize it as isosurface map in Multiwfn.
As I mentioned in past posts, you may also perform NTO analysis for T1 and S1 respectively, if their hole-NTO happen to be similar, then you can simply visually compare electron-NTO of T1 and that of S1 to study electron transition in terms of orbital wavefunctions.
Please describe which quantum chemistry program and computational level you are using.
Multiwfn doesn't support wavefunction of VASP. Currently, for VASP users, Multiwfn can only do two kinds of analyses:
(1) Analyses solely based on geometries: Such as mIGM, IGM, promolecular NCI, van der Waals potentials, etc. See Section 2.9.3 of Multiwfn manual for full list. In this case, you can use POSCAR of VASP as input file.
(2) Analyses based on grid data, such as basin analysis, grid data manipulations (including calculating plane averaged curve, etc.), plotting plane map based on interpolation of 3D grid data, etc. The grid data can be loaded from such as CONTCAR, ELFCAR.
If you use a large-core ECP, then there will only be valence ELF basins. If a small-core ECP is used, core basins will also be present, which correspond to subvalence shell electrons. You can use option 0 in basin analysis module to visualize the basins to better understand their characters.
It is found to be a bug, I have fixed it, please download the latest version from Multiwfn website. Thank you for bringing this bug to my attention.
You Google driver only contains Be2Br6_vacuum_high.log, I don't find .wfn file.
.wfn file produced under implicit solvation model can be normally analyzed by Multiwfn.
ELF gradient is calculated analytically in Multiwfn.
Integration of ELF basins in Multiwfn is carried out based on uniform grids. However, for very heavy atoms, because electron density around their nuclei varies very sharply, it is impossible to accurate integrate core ELF basins based on the uniform grids, even if lunatic quality grid is used. Since core basin is usually not of chemical interest (and if they are accurately integrated, then their populations must be very close to integer, that means you can easily predict their populations without any calculation), you can only focues on valence basins. Alternatively, using pseudopotential for these heavy atoms, then at least innermost ELF basins will not be presented.
It is fully possible to use implicit solvation model when generating wavefunction files.
The unit of ESP directly generated by Multiwfn is independent of charged state of the system. Only when you ask Multiwfn to export surface extrema or vertices as pdb file (in which B-factor field records ESP), the unit will be different for neutral and charged systems, because for charged systems, the ESP value on vdW surface is significantly larger than neutral systems, while the number of columns for recording B-factor in pdb format is very limited, so eV should be used instead of kcal/mol in this situation.
You can easily manually convert the unit. In addition, if you choose option "8 Export all surface vertices and surface extrema as vtx.pqr and extrema.pqr", then ESP will always be recorded in a.u. in the exported files, irrespective of charged state.
For a charged system, you should manually modify upper and lower limits of color scale in VMD, so that different colors can distinguish ESP in various surface areas. It is irrelevant to setting of Multiwfn.
Unfortunately, MR-SF- and SF- variants of TDDFT have not been explicitly supported yet.
IOp(3/33=3) extralinks=L316 scf=conventional noraff
You understanding is correct, an IOp is used to suppress orthogonalization. I didn't notice any freely available program that can suppress the orthogonalization...but I believe there should exist.
Dear Zander,
In fact TDA case is supported. Previously I didn't explicitly test compatibility with TDA, I just made a test, Multiwfn load the information correctly.
Best,
Tian
I just updated sobEDA tutorial package (//www.umsyar.com/soft/sobEDA_tutorial.zip), Section 6 is added to sobEDA_tutorial.pdf, which discussed how to solve SCF non-convergence problem.