Multiwfn official website: http://sobereva.com/multiwfn. Multiwfn forum in Chinese: http://bbs.keinsci.com/wfn
You are not logged in.
Hello,
If you choose two atomic nuclei to define the line to be plot in Multiwfn, then the value of X-axis will corresponds to the distance from starting point, which is irrelevant to Cartesian coordinate of the atoms.
Hello,
The color transition methods are hard-coded, you cannot change it without modifying source code. Please look at plot.f90 in source code package, in which "subroutine setcolortable" contains the code for defining the coloring transitions.
Change the last line to: $NBO AONAO DMNAO $END
Dear Prof. Tian Lu,
Thank you very much for your clear answer. So, for a set of different molecules that have the same nucleophilic reactive group, for example -SH (as in my case), can I also use Hirshfeld atomic charges or the Fukui function to evaluate the nucleophilicity of sulfur S between different molecules and make a rank between the molecules from the most nucleophile to the less one right ?
my best
Alessio
Dear Alessio,
Fukui function cannot be used to evaluate relative nucleophilicity between different molecules, but local softness can, which includes global softness on the top of Fukui function. Hirshfeld charges can also be compared between different molecules.
Tian
This depends on your practical purpose. I suggest reading original paper of Seminario (International Journal of Quantum Chemistry: Quantum Chemistry Symposium 30, 1271 -1277 (1996)), and review of relaxed force constant (Chem. Soc. Rev., 2008, 37, 1558–1567).
Bond force constant can be calculated by my code Sobtop (http://sobereva.com/soft/Sobtop) based on Gaussian .fch file produced by freq task (namely this .fch file contains Hessian), Seminario and modified Seminario methods are supported in Sobtop for this purpose.
If you need relaxed force constant, you can use http://www.oc.tu-bs.de/Grunenberg/compliance.html
This is no other post about AIMD.
Without detailed information, including ORCA input and output files, input file for Multiwfn, and all commands you inputted in Multiwfn, I cannot easily answer this question.
No. The ESP code in current version of Multiwfn is much faster than orca_plot, see comparison of speed (Table 1): Phys. Chem. Chem. Phys., 23, 20323–20328 (2021)
Dear Alessio,
1 In principle, considering implicit solvent should be more rigorous, but it doesn't necessarily improve fitting quality. So you can try both cases, with and without implicit solvation model.
2 Nucleophilicity index should only be calculated for reactant molecule, protonated amine NH3 should be ignored. Nucleophilic index is not necessarily limited to neutral species, though you can simultaneously consider other descriptors, such as Hirshfeld atomic charge of key atom(s), average local ionization energy (ALIE), etc.
Best regards,
Tian
Perhaps you loaded a wrong file. You should use Gaussian output file rather than fch file as input file in this analysis, as illustrated in Section 4.9.4 of Multiwfn manual.
Today a Multiwfn user asked me how to generate wavefunction file of excited state calculated by ORCA, my reply is given below, which may also be useful for other people:
Run the following task, assume that input file is H2CO.inp and output file is H2CO.out
! PBE0 def2-SVP
%tddft
tprint 1E-8
end
*xyz 0 1
C 0.00000000 0.00000000 -0.52887900
H 0.00000000 0.93777000 -1.12367000
O 0.00000000 0.00000000 0.67757700
H 0.00000000 -0.93777000 -1.12367000
*Then run this command to obtain H2CO.molden.input:
orca_2mkl H2CO -moldenLoad H2CO.molden.input into Multiwfn, and input
18
13
H2CO.out
1 //The excited state of interestNow you have NO_0001.mwfn in current folder, which records natural orbitals of S1 state, and can be used as input file of Multiwfn for wavefunction analysis on the S1 state.
If you want to convert it to .wfx format, boot up Multiwfn and load NO_0001.mwfn, then input
100
2
4
S1.wfxNow you have S1.wfx in current folder, which may be used in other code to perform wavefunction analysis for S1 state.
Note that using SDD for B is not good, 6-31+G(d,p) should also be applied to B. In addition, as you already specified genecp, pseudo=read is redundant.
Using revised input file, I performed calculation using G16, the output file is attached: a.out
and I didn't encounter the same problem when using the FiPC-NICS analysis in latest version of Multiwfn. Full output information is shown below for your check:
Multiwfn -- A Multifunctional Wavefunction Analyzer
Version 3.8(dev), update date: 2025-Aug-14
Developer: Tian Lu (Beijing Kein Research Center for Natural Sciences)
Multiwfn official website: http://sobereva.com/multiwfn
Multiwfn English forum: http://sobereva.com/wfnbbs
Multiwfn Chinese forum: http://bbs.keinsci.com/wfn
( Number of parallel threads: 16 Current date: 2025-09-01 Time: 10:25:19 )
Both following papers ***MUST BE CITED IN MAIN TEXT*** if Multiwfn is used:
Tian Lu, Feiwu Chen, J. Comput. Chem., 33, 580 (2012) DOI: 10.1002/jcc.22885
Tian Lu, J. Chem. Phys., 161, 082503 (2024) DOI: 10.1063/5.0216272
See "How to cite Multiwfn.pdf" in Multiwfn binary package for more information
Now input file path, for example, E:\Shakugan_no_Shana\Shana.mwfn
(.wfn/wfn/wfx/fch/molden/pdb/xyz/mol2/cif/cub... see Section 2.5 of manual)
Hint: Pressing ENTER button directly can select a file in a GUI window. To reload the past file, inputting "o". Input such as ?miku.fch can open the miku.fch in the same folder as the past file
C:\Users\sober\Desktop\a.out
Please wait...
Trying to load geometry from this file...
Geometry (final, input orientation) has been loaded from this file
Number of alpha/beta/total electrons: 181 181 362
Loaded C:\Users\sober\Desktop\a.out successfully!
Formula: Bq200 B1 C20 O20 Mn5 Total atoms: 246
Molecule weight: 845.70406 Da
"q": Exit program gracefully "r": Load a new file
************ Main function menu ************
0 Show molecular structure and view orbitals
1 Output all properties at a point 2 Topology analysis
3 Output and plot specific property in a line
4 Output and plot specific property in a plane
5 Output and plot specific property within a spatial region (calc. grid data)
6 Check & modify wavefunction
7 Population analysis and calculation of atomic charges
8 Orbital composition analysis 9 Bond order analysis
10 Plot total DOS, PDOS, OPDOS, local DOS, COHP and photoelectron spectrum
11 Plot IR/Raman/UV-Vis/ECD/VCD/ROA/NMR spectrum
12 Quantitative analysis of molecular surface
13 Process grid data (No grid data is presented currently)
14 Adaptive natural density partitioning (AdNDP) analysis
15 Fuzzy atomic space analysis
16 Charge decomposition analysis (CDA) and plot orbital interaction diagram
17 Basin analysis 18 Electron excitation analysis
19 Orbital localization analysis 20 Visual study of weak interaction
21 Energy decomposition analysis 22 Conceptual DFT (CDFT) analysis
23 ETS-NOCV analysis 24 (Hyper)polarizability analysis
25 Electron delocalization and aromaticity analyses
26 Structure and geometry related analyses
100 Other functions (Part 1) 200 Other functions (Part 2)
300 Other functions (Part 3)
25
================ Delocalization and aromaticity analyses ===============
0 Return to main menu
1 Multicenter bond order
-1 Multicenter bond order in NAO basis
2 AV1245 index
3 Iso-chemical shielding surface (ICSS)
4 NICS_ZZ for non-planar or tilted system
5 ELF-sigma/pi and LOL-sigma/pi
6 Harmonic oscillator measure of aromaticity (HOMA) and Bird indices
6a HOMAc (reparameterized HOMA) 6b HOMER (HOMA for excited states)
7 Shannon aromaticity index
8 Para-delocalization index (PDI)
9 Aromatic fluctuation index (FLU) and FLU-pi
10 Para linear response index (PLR)
11 Information-theoretic (ITA) aromaticity index
12 Properties of ring critical point
13 NICS-1D scan curve map, integral NICS (INICS) and FiPC-NICS
14 NICS-2D scan plane map
13
---------- NICS-1D scan and integral ----------
0 Return
Choose the way of defining the two end points of the line for scanning
1 Directly input Cartesian coordinates of two end points
2 The two end points are above and below the center of a plane fitted for specific atoms, and the line perpendicularly passes through their center
1
Input X,Y,Z of the first end point in Angstrom, e.g. 3.2,0.1,-9.5
-0.799429 0.96579 0.115633
Input X,Y,Z of the second end point in Angstrom, e.g. 3.2,0.1,-9.5
-0.799429 0.96579 -9.776147
How many points evenly distributing in the scanning line? e.g. 50
If press ENTER button directly, 99 points will be used, which is fine enough
200
0 Exit
1 Generate Gaussian input file for NICS-1D scanning
2 Load Gaussian output file of NICS-1D scanning
2
Input path of Gaussian output file of NICS-1D scanning task
e.g. D:\Aqours\Mari\shiny.out
C:\Users\sober\Desktop\a.out
Loading Bq 1
Loading Bq 2
Loading Bq 3
Loading Bq 4
Loading Bq 5
Loading Bq 6
Loading Bq 7
Loading Bq 8
Loading Bq 9
Loading Bq 10
Loading Bq 11
Loading Bq 12
Loading Bq 13
Loading Bq 14
Loading Bq 15
Loading Bq 16
Loading Bq 17
Loading Bq 18
Loading Bq 19
Loading Bq 20
Loading Bq 21
Loading Bq 22
Loading Bq 23
Loading Bq 24
Loading Bq 25
Loading Bq 26
Loading Bq 27
Loading Bq 28
Loading Bq 29
Loading Bq 30
Loading Bq 31
Loading Bq 32
Loading Bq 33
Loading Bq 34
Loading Bq 35
Loading Bq 36
Loading Bq 37
Loading Bq 38
Loading Bq 39
Loading Bq 40
Loading Bq 41
Loading Bq 42
Loading Bq 43
Loading Bq 44
Loading Bq 45
Loading Bq 46
Loading Bq 47
Loading Bq 48
Loading Bq 49
Loading Bq 50
Loading Bq 51
Loading Bq 52
Loading Bq 53
Loading Bq 54
Loading Bq 55
Loading Bq 56
Loading Bq 57
Loading Bq 58
Loading Bq 59
Loading Bq 60
Loading Bq 61
Loading Bq 62
Loading Bq 63
Loading Bq 64
Loading Bq 65
Loading Bq 66
Loading Bq 67
Loading Bq 68
Loading Bq 69
Loading Bq 70
Loading Bq 71
Loading Bq 72
Loading Bq 73
Loading Bq 74
Loading Bq 75
Loading Bq 76
Loading Bq 77
Loading Bq 78
Loading Bq 79
Loading Bq 80
Loading Bq 81
Loading Bq 82
Loading Bq 83
Loading Bq 84
Loading Bq 85
Loading Bq 86
Loading Bq 87
Loading Bq 88
Loading Bq 89
Loading Bq 90
Loading Bq 91
Loading Bq 92
Loading Bq 93
Loading Bq 94
Loading Bq 95
Loading Bq 96
Loading Bq 97
Loading Bq 98
Loading Bq 99
Loading Bq 100
Loading Bq 101
Loading Bq 102
Loading Bq 103
Loading Bq 104
Loading Bq 105
Loading Bq 106
Loading Bq 107
Loading Bq 108
Loading Bq 109
Loading Bq 110
Loading Bq 111
Loading Bq 112
Loading Bq 113
Loading Bq 114
Loading Bq 115
Loading Bq 116
Loading Bq 117
Loading Bq 118
Loading Bq 119
Loading Bq 120
Loading Bq 121
Loading Bq 122
Loading Bq 123
Loading Bq 124
Loading Bq 125
Loading Bq 126
Loading Bq 127
Loading Bq 128
Loading Bq 129
Loading Bq 130
Loading Bq 131
Loading Bq 132
Loading Bq 133
Loading Bq 134
Loading Bq 135
Loading Bq 136
Loading Bq 137
Loading Bq 138
Loading Bq 139
Loading Bq 140
Loading Bq 141
Loading Bq 142
Loading Bq 143
Loading Bq 144
Loading Bq 145
Loading Bq 146
Loading Bq 147
Loading Bq 148
Loading Bq 149
Loading Bq 150
Loading Bq 151
Loading Bq 152
Loading Bq 153
Loading Bq 154
Loading Bq 155
Loading Bq 156
Loading Bq 157
Loading Bq 158
Loading Bq 159
Loading Bq 160
Loading Bq 161
Loading Bq 162
Loading Bq 163
Loading Bq 164
Loading Bq 165
Loading Bq 166
Loading Bq 167
Loading Bq 168
Loading Bq 169
Loading Bq 170
Loading Bq 171
Loading Bq 172
Loading Bq 173
Loading Bq 174
Loading Bq 175
Loading Bq 176
Loading Bq 177
Loading Bq 178
Loading Bq 179
Loading Bq 180
Loading Bq 181
Loading Bq 182
Loading Bq 183
Loading Bq 184
Loading Bq 185
Loading Bq 186
Loading Bq 187
Loading Bq 188
Loading Bq 189
Loading Bq 190
Loading Bq 191
Loading Bq 192
Loading Bq 193
Loading Bq 194
Loading Bq 195
Loading Bq 196
Loading Bq 197
Loading Bq 198
Loading Bq 199
Loading Bq 200
Loading finished!
---------- Post-processing menu ----------
-3 Invert direction of X-axis
-2 Multiply NICS data by a factor
-1 Select component of NICS, current: Along the scanning line
0 Exit to main menu
1 Plot NICS curve along the line
2 Save image file of NICS curve along the line
3 Export NICS curve data along the line
4 Calculate integral of NICS along the line
5 Find minima and maxima along NICS curve
6 Calculate FiPC-NICS
1
Integral of NICS component: -76.59 ppm*Angstrom
---------- Post-processing menu ----------
-3 Invert direction of X-axis
-2 Multiply NICS data by a factor
-1 Select component of NICS, current: Along the scanning line
0 Exit to main menu
1 Plot NICS curve along the line
2 Save image file of NICS curve along the line
3 Export NICS curve data along the line
4 Calculate integral of NICS along the line
5 Find minima and maxima along NICS curve
6 Calculate FiPC-NICS
6
The scanning direction is found to be Z
FiPC-NICS.txt has been generated in current folder, meaning of each column:
Column 1: Point index
Column 2: Scanning distance in Angstrom
Column 3: In-plane component of NICS
Column 4: Out-of-plane component of NICS
FiPC-NICS is -1.116487 ppm, at 4.182 Angstrom
Perhaps the commands you inputted in Multiwfn were incorrect.
In TDDFT calculation, M06-2X represents triplet excited state better than most other functionals, you can have a try.
In addition, I strongly suggest using UKS instead of TDDFT to optimize T1, the former is not only significantly faster, but also the result is more robust.
Hello,
I have no experience in using the EFP feature in GAMESS. I suggest:
(1) Check if the wavefunction generated with EFP can be properly loaded into Multiwfn. Please check section "6.5 Check sanity of wavefunction" of Multiwfn manual.
(2) Plot electron density difference (EDD) between the wavefunction with and without the presence of EFP, and check if the EDD makes sense (i.e. properly reflecting the polarization effect caused by EFP). Please check Section 5.5 of Multiwfn manual on how to plot isosurface map of EDD.
If the above tests can be successfully passed, then you can normally perform other analyses on the wavefunction generated with EFP, including calculating electronic spatial extent.
Increasing number of grid point will not evidently affect the observed result.
Please take care of the calculation equation of Pauli potential used in Multiwfn (part 60 of Section 2.7 of Multiwfn manual). As you can see there are three component terms in the case of ispecial=0, you can individually plot them via Multiwfn to understand the detail. In which, Weizsacker potential (Vw) corresponds to iuserfunc=4, exchange-correlation potential (V_XC) corresponds to iuserfunc=1100.
BTW: Note that the position of nucleus is a singular point of the nuclear attractive potential in ESP, Multiwfn arbitrarily sets nuclear contribution to ESP at this point to 1000 a.u.
Because of approximations, such as approximated V_XC and finite-sized basis set, the resulting curve may not be well consistent with analyical solution.
In this case, your steps are correct.
In this way it is impossible to check if the Pauli potential is actually spherically symmetric. You should view isosurface map (via main function 5. Using isovalue of e.g. 0.7). You will find the isosurface is perfectly spherical.
Using iuserfunc=60 is correct. But as I mentioned, I cannot reproduce your observation, at least for Ar atom... More detailed information, such as input file, all commands you inputted in Multiwfn, the result, should be given.
Without description of your calculation steps in Multiwfn, it is difficult for me to make comment. I just tried, using Pauli potential (iuserfunc=53 with ispecial=0 or 1), to calculate grid data for Ar atom, the isosurface of Pauli potential is indeed fully spherically symmetric.
Hello,
You can use .molden.input files exported by orca_2mkl instead of the .fch files in this post.
.fch and .molden files provide essentially the same information for Multiwfn, see Section 2.5 of Multiwfn manual for more information.
There is no necessary contradiction. For example, MOLOPT basis set in CP2K is fully generally contracted, while it is used in combination with GTH pseudopotential. I suggest checking original paper of MOLOPT: THE JOURNAL OF CHEMICAL PHYSICS, 127, 114105 (2007).
I suggest changing to def2 basis set using "Extrapolate(3/4,def2) def2-QZVPP/C", the ECP should be automatically assigned for Rh.
Multiwfn doesn't directly support .dat file. For GAMESS-US users, one can obtain .wfn file from .dat file and then use it as Multiwfn input file, please check beginning part of Chapter 4 of Multiwfn manual on how to do.
After that, you can use Multiwfn to perform a wide variety of analyses, including calculating grid data of electron density around the interested region using main function 5, please check Section 4.5 of manual for examples. The region to be calculated can be flexibly defined, see Section 3.6 of manual about the different modes for setting up grids.
If your purpose is generating path from BCP, you should select option "8 Generating the paths connecting (3,-3) and (3,-1) CPs" rather than "9 Generating the paths connecting (3,+1) and (3,+3) CPs", the former doesn't need (3,+3).
It is not possible to use wB97X-D3(BJ) in G16, but ORCA supports it.
Hello,
This is no corresponding option, however, you can delete other type of CPs (option "-4 Modify or export CPs" and then "2 Delete some CPs") before printing out properties of BCP.
Also, note that in option -3 of topology analysis module there are suboptions "4 Save CPs to CPs.txt in current folder" and "5 Load CPs from a file", via them you can restore all CPs after deleting some of them.
Best regards,
Tian
As I noted earlier, this depends on how the surface is defined. If it is defined as isosurface of electron density, then "a" will be a constant for all atoms and all molecules, making the discussion completely meaningless.
Detailed description about these modes can be found in Section 3.6 of Multiwfn manual. The most convenient way is using option 10, by which you can visually setting up the box.
HELP and HELV, in their original paper, were used to characterize lone pair, it seems that you hope to use them to study bonding ELF basin, this might also make sense but I am not fully sure.
You may visualize other basins surrounding the An to try to identify their meaning, I think they should not be significant and thus can be ignored if your purpose is simply studying bonding.
By the way, to greatly save computational cost of basin analysis of ELF for large systems like this, you can properly define the box when setting up grid so that the box only cover the An atom and coordination atoms around it.