Multiwfn forum

FLU reference parameter

dummy@example.com (Hamid) — Mon, 13 Jul 2026 18:48:48 +0000

thanks

Large gap between two methods used to measure S0-T1 excitation energy

dummy@example.com (sobereva) — Sat, 11 Jul 2026 21:27:28 +0000

delta-SCF result is more reliable.

1 delta-SCF itself is an iterative process, reaction field is intrinsically updated during the SCF iteration, so external iteration is never needed for consideration of solvent effect.

2 ORCA doesn't have this feature, this is one of evident limitations of ORCA (especially for studying CT excitation in solvent environment as linear response model performs poorly in this case). The corrected linear-response model that available in G16 is also not supported by ORCA.

TDM with nm as the unit

dummy@example.com (sobereva) — Sat, 11 Jul 2026 03:52:54 +0000

Multiwfn doesn't support TDM of this form.

QTAIM with ECP

dummy@example.com (rikaaardoss) — Fri, 10 Jul 2026 23:36:26 +0000

Hi, AIMAll tells you how to resolve this issue, you have to modify the wfx file. You have to follow the suggestions in the "blue window".

Best regards

Main function 12 with relativistic wavefunction

dummy@example.com (sobereva) — Wed, 08 Jul 2026 13:17:54 +0000

Impact of relativistic effects on the ESP analysis result is fully reflected by the relativistic effects on wavefunctions. So, the formulae for calculation of ESP are exactly the same for relativistic and nonrelativistic cases.

Integrating grid data in cube in whole space

dummy@example.com (sobereva) — Tue, 07 Jul 2026 22:05:49 +0000

Hello,

In Multiwfn, all indices, including orbital indices, always starts from 1 (PS: I quite dislike starting from 0)

Multiwfn AdNDP Analysis with ORCA Output File

dummy@example.com (b322qr) — Tue, 07 Jul 2026 01:10:59 +0000

I'm posting my reply here for future user with the same question to see the results.

I used the phenanthrene molecule from section 4.14.3 in the Multiwfn manual and compared what I believe is the same calculation in ORCA. For others to see, here is the input file for the Gaussian calculation:

# b3lyp/3-21g pop=nboread

b3lyp/3-21g opted

0 1
C 0.00000000 3.56061700 -0.29722900
C 0.00000000 2.83932500 0.87979300
C 0.00000000 1.42361400 0.86771500
C 0.00000000 0.72986200 -0.38070000
C 0.00000000 1.49924400 -1.56931800
C 0.00000000 2.88151800 -1.53149000
C 0.00000000 0.67926500 2.09717700
C 0.00000000 -0.72986200 -0.38070000
C 0.00000000 -1.42361400 0.86771500
C 0.00000000 -0.67926500 2.09717700
C 0.00000000 -2.83932500 0.87979300
H 0.00000000 -3.34927800 1.83744900
C 0.00000000 -3.56061700 -0.29722900
C 0.00000000 -2.88151800 -1.53149000
C 0.00000000 -1.49924400 -1.56931800
H 0.00000000 1.23356400 3.02968400
H 0.00000000 4.64400700 -0.27588200
H 0.00000000 3.34927800 1.83744900
H 0.00000000 1.00203800 -2.53051300
H 0.00000000 3.44643700 -2.45642700
H 0.00000000 -1.23356400 3.02968400
H 0.00000000 -4.64400700 -0.27588200
H 0.00000000 -3.44643700 -2.45642700
H 0.00000000 -1.00203800 -2.53051300

$nbo dmnao aonao $end

And here is the input for the ORCA calculation:
! B3LYP 3-21G NBO

%nbo
nbokeylist="$nbo aonao dmnao $end"
end

%MaxCore 7500
%pal nprocs 16
end

* xyz 0 1
C 0.00000000 3.56061700 -0.29722900
C 0.00000000 2.83932500 0.87979300
C 0.00000000 1.42361400 0.86771500
C 0.00000000 0.72986200 -0.38070000
C 0.00000000 1.49924400 -1.56931800
C 0.00000000 2.88151800 -1.53149000
C 0.00000000 0.67926500 2.09717700
C 0.00000000 -0.72986200 -0.38070000
C 0.00000000 -1.42361400 0.86771500
C 0.00000000 -0.67926500 2.09717700
C 0.00000000 -2.83932500 0.87979300
H 0.00000000 -3.34927800 1.83744900
C 0.00000000 -3.56061700 -0.29722900
C 0.00000000 -2.88151800 -1.53149000
C 0.00000000 -1.49924400 -1.56931800
H 0.00000000 1.23356400 3.02968400
H 0.00000000 4.64400700 -0.27588200
H 0.00000000 3.34927800 1.83744900
H 0.00000000 1.00203800 -2.53051300
H 0.00000000 3.44643700 -2.45642700
H 0.00000000 -1.23356400 3.02968400
H 0.00000000 -4.64400700 -0.27588200
H 0.00000000 -3.44643700 -2.45642700
H 0.00000000 -1.00203800 -2.53051300
*

In the Multiwfn manual, the residual density distribution after accepting the first 27 orbitals is:
1 C : 1.0250 2C : 1.0370    3C : 1.0280 4C : 1.0414
5 C : 1.0339    6C : 1.0262    7C : 0.1322    8C : 1.0414
9 C : 1.0280    10C : 0.1322    11C : 1.0370    12H : 0.0117
13C : 1.0250    14C : 1.0262    15C : 1.0339    16H : 0.0121
17H : 0.0113    18H : 0.0117    19H : 0.0126 20H : 0.0111
21H : 0.0121 22H : 0.0113    23H : 0.0111    24H : 0.0126

In my ORCA calculations, I receive the following:
1C : 1.0234 2C : 1.0295 3C : 1.0243 4C : 1.0381
5C : 1.0246 6C : 1.0246 7C : 0.1302 8C : 1.0473
9C : 1.0167 10C : 0.1341 11C : 1.0345 12H : 0.0112
13C : 1.0234 14C : 1.0246 15C : 1.0217 16H : 0.0115
17H : 0.0106 18H : 0.0112 19H : 0.0120 20H : 0.0104
21H : 0.0115 22H : 0.0106 23H : 0.0104 24H : 0.0120

Later, the analysis of orbital 31 shows the following in the Multiwfn manual:
NAO# Center Label Type    Composition
67    8(C ) px    Val( 2p) 2.247%
76    9(C ) px    Val( 2p) 1.929%
94 11(C ) px    Val( 2p) 13.276%
105 13(C ) px    Val( 2p) 32.846%
114 14(C ) px    Val( 2p) 34.484%
123 15(C ) px    Val( 2p) 15.204%

In my ORCA calculations, I receive the following:
NAO# Center Label Type    Composition
67    8(C ) px    Val( 2p) 2.217%
76    9(C ) px    Val( 2p) 1.942%
94 11(C ) px    Val( 2p) 13.387%
105 13(C ) px    Val( 2p) 32.956%
114 14(C ) px    Val( 2p) 34.401%
123 15(C ) px    Val( 2p) 15.082%

As you can see, the differences between the two calculations in not very significant (>0.2 %); however, the visualization of the orbitals is still weird.
These are the orbitals from the Multiwfn manual

And here are the orbitals from my calculation (using a contour value of 0.03)

NICS curve and plane in ORCA

dummy@example.com (sobereva) — Fri, 03 Jul 2026 17:15:57 +0000

Unfortunately it is not possible. The main reason is that in ORCA the point to be calculated must have at least one basis function, therefore when performing NICS scan, the number of basis functions to be calculated may be huge, causing very high extra computational cost, and the basis functions also affect the NICS scan result. If this limitation of ORCA could be solved in the future version, I would like to make NICS scan and ICSS functions of Multiwfn fully compatible with ORCA.

Multiwfn MK RESP calculations

dummy@example.com (Abdirahman) — Mon, 29 Jun 2026 08:28:45 +0000

Okay, thanks for the info.

Wave function analysis at DLPNO-CCSD(T)

dummy@example.com (saeed_E) — Fri, 26 Jun 2026 20:19:40 +0000

Dear Tian,
Thank you very much.

Saeed

ESP Scale Change

dummy@example.com (alexlester1996) — Sun, 21 Jun 2026 19:32:15 +0000

Dr Tian,
Thank you for this helpful comment.
Best wishes,
Alex

Emission for fluorescein-type dyes

dummy@example.com (sobereva) — Sun, 21 Jun 2026 17:19:43 +0000

If in the actual situation the anionic dye is able to tightly bind counterion (strictly speaking, this can be confirmed by molecular dynamics simulation), your treatment is fully physically sound.

I suggest checking what is the nature of the low-energy transitions (>1000 nm) using hole-electron analysis in Multiwfn, which can provide valuable physical insight.

Also, don't forget to confirm that the reference state wavefunction is stable.

IFCT / Electron excitation analysis / ORCA /json / problem

dummy@example.com (sobereva) — Sat, 20 Jun 2026 11:19:40 +0000

Multiwfn can load configurational coefficients from either ORCA output file or json file. As stated in the prompt, because current TDDFT calculation is based on an open-shell reference state (doublet ground state), Multiwfn cannot load coefficients from the json file despite you provided it (because current version of ORCA has a bug in outputting json file) but automatically load coefficients from ORCA output file. In this case, the analysis result may be inaccurate with "TDA false"

It is expected that next release of ORCA revision will solve this bug, at that time, I will update Multiwfn, then this difficulty will be perfectly solved.

Segmentation fault using LIBRETA for RESP calculation

dummy@example.com (sobereva) — Thu, 18 Jun 2026 00:03:20 +0000

Dear Marcos,

I can normally run the calculation using your files on my RockyLinux 10 computer with Multiwfn 2026.6.2. Using 96 cores, the calculation is finished in about half a minute. All outputted information is given below:

 Multiwfn -- A Multifunctional Wavefunction Analyzer
 Version 2026.6.2 (release date is the same as version name)
 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:  96  Current date: 2026-06-17  Time: 23:59:38 )

 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:\Strawberry_Panic\Chikaru_Minamoto.mwfn
 (.mwfn/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
esp_0.molden.input
 Please wait...
 Loading various information of the wavefunction
 This file is recognized to be generated by ORCA because there is "orca" word in title line. Special treatments are applied...
 Loading basis set definition...
 All D basis functions are spherical harmonic type
 Loading orbitals...
 The actual number of orbitals read:       470
 Converting basis function information to GTF information...
 Back converting basis function information from Cartesian to spherical type...
 Generating density matrix...
 Generating overlap matrix...

 Total/Alpha/Beta electrons:    260.0000    130.0000    130.0000
 Net charge:    -4.00000      Expected multiplicity:    1
 Atoms:     43,  Basis functions:    470,  GTFs:    988
 This is a restricted single-determinant wavefunction
 Orbitals from 1 to   130 are occupied

 Loaded esp_0.molden.input successfully!

 Formula: H12 C10 N5 O13 P3      Total atoms:      43
 Molecule weight:       503.14971 Da
 Point group: C1

 "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)
7
 NOTE: There is a review comprehensively introducing various atomic charges:
 Tian Lu, Qinxue Chen, Partial Charges, In Exploring Chemical Concepts Through T
 heory and Computation. WILEY-VCH GmbH: Weinheim (2024); pp. 161-187. DOI: 10.10
 02/9783527843435.ch6

      ============== Population analysis and atomic charges ==============
 -2 Calculate interaction energy between fragments based on atomic charges
 -1 Define fragment
 0 Return
 1 Hirshfeld atomic charge
 2 Voronoi deformation density (VDD) atom population
 5 Mulliken atom & basis function population analysis
 6 Lowdin atom & basis function population analysis
 7 Modified Mulliken atom population defined by Ros & Schuit (SCPA)
 8 Modified Mulliken atom population defined by Stout & Politzer
 9 Modified Mulliken atom population defined by Bickelhaupt
 10 Becke atomic charge with atomic dipole moment correction
 11 Atomic dipole corrected Hirshfeld atomic charge (ADCH) (recommended)
 12 CHELPG ESP fitting atomic charge
 13 Merz-Kollmann (MK) ESP fitting atomic charge
 14 AIM atomic charge
 15 Hirshfeld-I atomic charge
 16 CM5 atomic charge    -16 Generate 1.2*CM5 atomic charge
 17 Electronegativity Equalization Method (EEM) atomic charge
 18 Restrained ElectroStatic Potential (RESP) atomic charge
 19 Gasteiger (PEOE) charge
 20 Minimal Basis Iterative Stockholder (MBIS) charge
18

             ------------ Calculation of RESP charges ------------
 -1 Load list of conformer and weights from external file
 0 Return
 1 Start standard two-stage RESP fitting calculation
 2 Start one-stage ESP fitting calculation with constraints
 3 Set method and parameters for distributing fitting points, current: MK
 4 Set hyperbolic penalty and various other running parameters
 5 Set equivalence constraint in fitting, current: H in CH2 and CH3
 6 Set charge constraint in fitting, current: No constraint
 7 Set the way of determining connectivity, current: Guess from bond length
 8 Toggle if loading fitting points and ESP values from Gaussian output file of pop=MK/CHELPG task with IOp(6/33=2) during the calculation, current: No
 9 Load additional fitting centers, current: None
 10 Choose the atomic radii used in fitting, current: Automatic
 11 Choose ESP type, current: Nuclear + Electronic
-1
 Input path of the file containing conformer list, e.g. C:\conflist.txt
confs_weights
 There are    5 conformers
 Sum of weights:    1.000000

             ------------ Calculation of RESP charges ------------
-1 Reload list of conformers from external file, current:   5 conformers
 0 Return
 1 Start standard two-stage RESP fitting calculation
 2 Start one-stage ESP fitting calculation with constraints
 3 Set method and parameters for distributing fitting points, current: MK
 4 Set hyperbolic penalty and various other running parameters
 5 Set equivalence constraint in fitting, current: H in CH2 and CH3
 6 Set charge constraint in fitting, current: No constraint
 7 Set the way of determining connectivity, current: Guess from bond length
 8 Toggle if loading fitting points and ESP values from Gaussian output file of pop=MK/CHELPG task with IOp(6/33=2) during the calculation, current: No
 9 Load additional fitting centers, current: None
 10 Choose the atomic radii used in fitting, current: Automatic
 11 Choose ESP type, current: Nuclear + Electronic
5

 Please select options 1~3. You can also use options 10 or 11 to generate file containing equivalence constraint, which can then be utilized by option 1
 Note: For standard two-stage RESP fitting, options 0 and 1 only take effect for the first stage

 0 No equivalence constraint will be imposed
 1 Load equivalence constraint setting from external plain text file
 2 Constraint H in each =CH2, -CH2-, CH3 to be equivalent in one-stage fitting
 10 Export equivalence constraint corresponding to "H in each =CH2, -CH2-, CH3" to eqvcons_H.txt in current folder
 11 Generate equivalence constraint according to point group of global or local geometry and write to eqvcons_PG.txt in current folder
1
 Input path of the plain text file, e.g. C:\eqvcons.txt
 If pressing ENTER button directly, eqvcons.txt in current folder will be loaded
eqvcons.txt
 OK, equivalence constraint will be loaded from it during calculation

             ------------ Calculation of RESP charges ------------
-1 Reload list of conformers from external file, current:   5 conformers
 0 Return
 1 Start standard two-stage RESP fitting calculation
 2 Start one-stage ESP fitting calculation with constraints
 3 Set method and parameters for distributing fitting points, current: MK
 4 Set hyperbolic penalty and various other running parameters
 5 Set equivalence constraint in fitting, current: Customized
 6 Set charge constraint in fitting, current: No constraint
 7 Set the way of determining connectivity, current: Guess from bond length
 8 Toggle if loading fitting points and ESP values from Gaussian output file of pop=MK/CHELPG task with IOp(6/33=2) during the calculation, current: No
 9 Load additional fitting centers, current: None
 10 Choose the atomic radii used in fitting, current: Automatic
 11 Choose ESP type, current: Nuclear + Electronic
1
 Atomic radii used:
 Element:H      vdW radius (Angstrom): 1.200
 Element:C      vdW radius (Angstrom): 1.500
 Element:N      vdW radius (Angstrom): 1.500
 Element:O      vdW radius (Angstrom): 1.400
 Element:P      vdW radius (Angstrom): 1.800
 Generating fitting points and calculate ESP for conformer    1
 Number of MK fitting points used:     23208

 Initializing LIBRETA library (fast version) for ESP evaluation ...
 LIBRETA library has been successfully initialized!

 NOTE: The ESP evaluation code based on LIBRETA library is being used. Please cite Multiwfn original papers (J. Comput. Chem., 33, 580-592 (2012) and J. Chem. Phys., 161, 082503 (2024)) and the paper describing the efficient ESP evaluation algorithm adopted by Multiwfn (Phys. Chem. Chem. Phys., 23, 20323 (2021))

 Progress: [##################################################]   100.0 %     \
 Generating fitting points and calculate ESP for conformer    2
 Number of MK fitting points used:     22932

 Initializing LIBRETA library (fast version) for ESP evaluation ...
 LIBRETA library has been successfully initialized!

 NOTE: The ESP evaluation code based on LIBRETA library is being used. Please cite Multiwfn original papers (J. Comput. Chem., 33, 580-592 (2012) and J. Chem. Phys., 161, 082503 (2024)) and the paper describing the efficient ESP evaluation algorithm adopted by Multiwfn (Phys. Chem. Chem. Phys., 23, 20323 (2021))

 Progress: [##################################################]   100.0 %     /
 Generating fitting points and calculate ESP for conformer    3
 Number of MK fitting points used:     24706

 Initializing LIBRETA library (fast version) for ESP evaluation ...
 LIBRETA library has been successfully initialized!

 NOTE: The ESP evaluation code based on LIBRETA library is being used. Please cite Multiwfn original papers (J. Comput. Chem., 33, 580-592 (2012) and J. Chem. Phys., 161, 082503 (2024)) and the paper describing the efficient ESP evaluation algorithm adopted by Multiwfn (Phys. Chem. Chem. Phys., 23, 20323 (2021))

 Progress: [##################################################]   100.0 %     \
 Generating fitting points and calculate ESP for conformer    4
 Number of MK fitting points used:     23186

 Initializing LIBRETA library (fast version) for ESP evaluation ...
 LIBRETA library has been successfully initialized!

 NOTE: The ESP evaluation code based on LIBRETA library is being used. Please cite Multiwfn original papers (J. Comput. Chem., 33, 580-592 (2012) and J. Chem. Phys., 161, 082503 (2024)) and the paper describing the efficient ESP evaluation algorithm adopted by Multiwfn (Phys. Chem. Chem. Phys., 23, 20323 (2021))

 Progress: [##################################################]   100.0 %     /
 Generating fitting points and calculate ESP for conformer    5
 Number of MK fitting points used:     22904

 Initializing LIBRETA library (fast version) for ESP evaluation ...
 LIBRETA library has been successfully initialized!

 NOTE: The ESP evaluation code based on LIBRETA library is being used. Please cite Multiwfn original papers (J. Comput. Chem., 33, 580-592 (2012) and J. Chem. Phys., 161, 082503 (2024)) and the paper describing the efficient ESP evaluation algorithm adopted by Multiwfn (Phys. Chem. Chem. Phys., 23, 20323 (2021))

 Progress: [##################################################]   100.0 %     \
 Reloading the first file when Multiwfn boots up...

 No charge constraint is imposed in this stage
 Loading equivalence constraint setting from eqvcons.txt
 Atom equivalence constraint imposed in this fitting stage:
 Constraint   1:   14(O )   18(O )
 Constraint   2:   15(O )   19(O )
 Constraint   3:   16(O )   20(O )   24(O )
 Constraint   4:   32(H )   33(H )

 **** Stage 1: RESP fitting under weak hyperbolic penalty
 Convergence criterion:  0.0000010000
 Hyperbolic restraint strength (a): 0.000500    Tightness (b): 0.100000
 Iter:   1   Maximum charge variation:    1.5077332405
 Iter:   2   Maximum charge variation:    0.1107636558
 Iter:   3   Maximum charge variation:    0.0039944874
 Iter:   4   Maximum charge variation:    0.0001558212
 Iter:   5   Maximum charge variation:    0.0000063158
 Iter:   6   Maximum charge variation:    0.0000002668
 Successfully converged!

 **** Stage 2: RESP fitting under strong hyperbolic penalty
 Atom equivalence constraint imposed in this fitting stage:
 Constraint   1:   32(H )   33(H )
 Fitting objects: sp3 carbons, methyl carbons and hydrogens attached to them
 Indices of these atoms:
    4C    32H    33H     6C    34H     8C    35H    10C    36H    12C
   37H
 Convergence criterion:  0.0000010000
 Hyperbolic restraint strength (a): 0.001000    Tightness (b): 0.100000
 Iter:   1   Maximum charge variation:    1.5345257854
 Iter:   2   Maximum charge variation:    0.0125232048
 Iter:   3   Maximum charge variation:    0.0001932288
 Iter:   4   Maximum charge variation:    0.0000027221
 Iter:   5   Maximum charge variation:    0.0000000376
 Successfully converged!

   Center       Charge
     1(P )   1.4135630800
     2(P )   1.4241920085
     3(P )   1.3536680793
     4(C )   0.3760715582
     5(O )  -0.6307782017
     6(C )   0.2527058666
     7(O )  -0.7174906128
     8(C )   0.6780984537
     9(O )  -0.9591365413
    10(C )   0.1435696375
    11(O )  -0.9245005823
    12(C )   0.8770431354
    13(N )  -1.0284259720
    14(O )  -0.8952224744
    15(O )  -0.8323241779
    16(O )  -0.9465635221
    17(C )   0.5937863633
    18(O )  -0.8952224744
    19(O )  -0.8323241779
    20(O )  -0.9465635221
    21(N )  -0.7645206376
    22(O )  -0.6580604115
    23(O )  -0.5547605926
    24(O )  -0.9465635221
    25(C )   0.7500835753
    26(C )  -0.4617961113
    27(C )   1.4542593076
    28(N )  -1.5345257854
    29(N )  -0.6146215415
    30(C )   0.6056978690
    31(N )  -0.6829891638
    32(H )  -0.0477356295
    33(H )  -0.0477356295
    34(H )  -0.0545910770
    35(H )  -0.1429844986
    36(H )   0.0232647432
    37(H )  -0.1020603051
    38(H )   0.0253107478
    39(H )   0.0229862788
    40(H )   0.5191095206
    41(H )   0.5865250363
    42(H )   0.5380987896
    43(H )   0.5834631134
 Sum of charges:  -4.0000000000
 Conformer:    1   RMSE:    0.004924   RRMSE:    0.012598
 Conformer:    2   RMSE:    0.005126   RRMSE:    0.013223
 Conformer:    3   RMSE:    0.007300   RRMSE:    0.019247
 Conformer:    4   RMSE:    0.005913   RRMSE:    0.014901
 Conformer:    5   RMSE:    0.008026   RRMSE:    0.020533
 Weighted RMSE:    0.005650   Weighted RRMSE    0.014530

 Note: Because present calculation involves multiple conformers, the result cannot be exported to .chg file

So I still believe Multiwfn was not fully configurated on your system. I also provide relevant settings in my ~/.bashrc file here:

ulimit -s unlimited
export OMP_STACKSIZE=200M
export PATH=$PATH:/sob/Multiwfn_xxx_bin_Linux
export Multiwfnpath=/sob/Multiwfn_xxx_bin_Linux

Best regards,

Tian

Reasonableness of wB97XD for wavefunction analyses

dummy@example.com (sobereva) — Fri, 12 Jun 2026 22:42:00 +0000

A Multiwfn user asked me if wB97XD is suitable for calculation of delocalization index using Multiwfn and its comparison with GGA and hybrid-GGA functionals. My reply is also provided here, which may be also useful for other Multiwfn users.

---------
wB97XD is a very reasonable choice. Its improvement of representation of wavefunction over (meta-)GGA and common hybrid (meta-)GGA depends on the specific functional and the system.

When a (meta-)(hybrid-)GGA shows evident delocalization error (also known as self-interaction error, SIE) for a system, the improvement of wB97XD is significant. Representative examples include [18]annulene and cyclo[18]carbon. In Angew. Chem. Int. Ed. 2004, 43, 4200–4206, it is shown that B3LYP (20% global HF component) cannot give reasonable structure of [18]annulene, while I found the geometry optimized by wB97XD is fully correct. cyclo[18]carbon cannot be reasonably represented by any functional with insufficient HF composition at long-range, e.g. B3LYP (20% HF) and PBE0 (25% HF), while BHandHLYP (50% HF), M06-2X (54% HF) and wB97XD (22.2% to 100% HF from short-range limit to long-range limit) work reasonably, see my study and review about this system: Carbon, 165, 468-475 (2020), Acc. Mater. Res., 6, 1220−1231 (2025).

(hybrid-)(meta-)GGA functionals with insufficient HF composition at long-range tend to severely overestimate electronic delocalization of the aforemention systems, not only the resulting wavefunction is not reasonable, but also the optimized geometry is qualitatively wrong (strong tendency towards planarization and bond length equalization). For more information about the poor performance of these functionals, see review about delocalization error: WIREs Comput Mol Sci. 2022;e1631.