Multiwfn forum

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

You are not logged in.

#1 Re: Multiwfn and wavefunction analysis » Plot PDOS for s, p AO using diffuse functions in basis set » Yesterday 18:23:31

Dear Alex,

Hirshfeld and Becke methods are two ways of defining atomic spaces, they are not the ways of defining atomic orbitals in molecular environment. So they cannot be used in this context.

I suggest removing diffuse functions while using def2-TZVP instead when performing DOS analysis, which is adequate.

Best,

Tian

#2 Multiwfn and wavefunction analysis » Integrate electron density difference in its isosurfaces » Yesterday 18:20:39

sobereva
Replies: 0

Multiwfn manual is updated today. Section 4.200.14.3 has been added to manual to show how to integrate electron density difference in its isosurfaces using domain analysis module.

#3 Re: Quantum Chemistry » The input Orca file is immediately terminated with error » 2024-07-25 15:43:54

JK and COSX are two ways of accelerating calculation of HF exchange, the former requires auxiliary basis set of JK type. Without RIJCOSX, the HF exchange will be calculated without any approximation.

#4 Re: Multiwfn and wavefunction analysis » MultiWFN crashing while plotting NMR spectra » 2024-07-25 13:12:03

You do not need to rerun, just change the encoding of the output file (using such as Ultraedit, you can also use google to find other ways).

#5 Re: Multiwfn and wavefunction analysis » Multiwfn won't even extract! » 2024-07-25 13:10:12

The information you provided is not adequate for addressing the problem. Try to use Multiwfn 3.8(dev) instead first. Version 3.7 is quite old and no longer maintained.

Tian

#6 Re: Quantum Chemistry » ELF, LOL, laplacian » 2024-07-25 13:08:32

See
(1) "The ELF Topological Analysis Contribution to Conceptual Chemistry and Phenomenological Models" in book "The Quantum Theory of Atoms in Molecules: From Solid State to DNA and Drug Design"
(2) Struct Bond (2016) 170: 213–248
My work "Revealing Molecular Electronic Structure via Analysis of Valence Electron Density" (http://www.whxb.pku.edu.cn/EN/10.3866/PKU.WHXB201709252) also contains application example of topology analysis of ELF.

#8 Re: Quantum Chemistry » How to augment Def2 series with "aug" prefix? » 2024-07-23 13:07:35

Dear Saeed,

I have a blog article to introduce different ways to add diffuse functions to def2 series of basis set, see http://sobereva.com/340

Briefly speaking, such as aug-def2-TZVP borrows the diffuse functions from aug-cc-pVTZ or aug-pc-2 to def2-TZVP.

Best,

Tian

#9 Re: Multiwfn and wavefunction analysis » MultiWFN crashing while plotting NMR spectra » 2024-07-22 18:57:49

I can normally plot NMR spectrum using your output file (simply inputting 11, then 7, then 0).
If your Multiwfn version is not the latest one, please update Multiwfn first.
If you run ORCA under Powershell, you should change encoding of the output file from the default UTF16 to UTF-8 or Unicode, otherwise Multiwfn cannot load any information from it.

#10 Re: Quantum Chemistry » About Orca » 2024-07-22 18:54:24

Dear Saeed,

I am not sure, please ask developer in ORCA forum

Tian

#11 Re: Multiwfn and wavefunction analysis » EDA - example from manual - Error: Unable to find OW » 2024-07-21 16:26:04

Dear Igor,

You inputted wrong file.
Molecular list file (named as mollist.txt in the example in Section 4.21.1.1 of Multiwfn manual) contains paths of molecular type files (e.g. water.txt), which record atom types and atomic charges of each kind of molecule.

Multiwfn first asked you to input molecular list file, however you inputted path of a molecular type file, in which the first word is OW (atom type of oxygen in water), this is why Multiwfn prompted "Unable to find OW". So, you should first input molecular list file.

Best,

Tian

#12 Re: Multiwfn and wavefunction analysis » .wfn file » 2024-07-21 13:34:02

I updated Multiwfn today, now this function can use the file only containing atomic coordinate such as .xyz as input file.

#13 Re: Multiwfn and wavefunction analysis » promolecule » 2024-07-21 13:27:24

Please notice "ispheratm" in settings.ini. If you set it to 0, atomic wavefunction files will not be automatically sphericalized.

#14 Re: Multiwfn and wavefunction analysis » AIM analysis of a metallacyles » 2024-07-20 14:39:00

Dear Paul,

I don't think there should be a CP, as the distance is long. Please check if Poincare-Hopf relationship has been satisfied, if yes, all CPs should have been found.

Best,

Tian

#16 Re: Multiwfn and wavefunction analysis » MultiWFN crashing while plotting NMR spectra » 2024-07-18 13:10:32

Please upload your input file for Multiwfn (or send it to my E-mail), and show me all commands you inputted in Multiwfn, I will try to reproduce the issue.

#17 Re: Multiwfn and wavefunction analysis » Integration inaccuracy for heavy atoms » 2024-07-17 13:38:40

Dear John,

Multiwfn employs uniform grid for integrating ELF basins. The integration accuracy for valence ELF basins should be fine if grid spacing is small, however it is not expected that integration accuracy for core basins is also satisfactory for heavy atoms, because electron density varies quite sharply around their nuclei. Usually only valence ELF basins are of interest, it is not a problem that core basins cannot be well integrated using uniform grids.

Best regards,

Tian

#20 Multiwfn and wavefunction analysis » Two important updates on the function of exporting .47 file » 2024-07-12 10:27:53

sobereva
Replies: 0

Multiwfn updated today has two important changes about exporting NBO .47 file:

IMPROVEMENT: f-type spherical-harmonic basis functions have been supported in exporting NBO .47 file.

BUG FIXED: The contraction coefficients of spherical-harmonic D shells in the .47 file exported by Multiwfn are all zero.

#21 Re: Multiwfn and wavefunction analysis » Crash while converting molden into .47 » 2024-07-12 09:45:49

Kerelsso wrote:

I'm sorry for the delay. I performed a singlepoint of my molecule on ORCA with def2-SV(P), converted .gbw into .molden, change Pt number of electrons to 18 (for the ECP) and tried to export it to .47 with Multifwn. The results have been the same, but this time the program crashes wheteher I set "iloadasCart" to 0 or 1. It also crashes when I tried to export it to .wfx, or if I export it to .fchk and then try to use that .fchk to export it to .47.

I'll send you my input to your e-mail. I'm very sorry for the incoveniences.

It seems that the Multiwfn version you used is too old. I can normally convert your .molden file to .47 via latest version of Multiwfn (updated on 2024-Jul-12), see below output:

 Multiwfn -- A Multifunctional Wavefunction Analyzer
 Version 3.8(dev), update date: 2024-Jul-12
 Developer: Tian Lu (Beijing Kein Research Center for Natural Sciences)
 Below paper ***MUST BE CITED IN MAIN TEXT*** if Multiwfn is used in your work:
          Tian Lu, Feiwu Chen, J. Comput. Chem., 33, 580-592 (2012)
 See "How to cite Multiwfn.pdf" in Multiwfn binary package for more information
 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:  24  Current date: 2024-07-12  Time: 17:35:08 )
 OpenMP stacksize for each thread:     953.67 MB

 Input file path, for example E:\Sakurako-san_no_Ashimoto_ni_wa_Shitai_ga_Umatteiru\Sakurako.wfn
 (.mwfn/wfn/wfx/fch/molden/pdb/xyz/mol2/cif/cub... see Section 2.5 of manual)
 Hint: Pressing ENTER button directly can select file in a GUI window. To reload the file last time used, simply input the letter "o". Input such as ?miku.fch can open the miku.fch in the same folder as the file last time used.
C:\Users\sober\Desktop\sp_bw_BPt_benz_H2act_HH.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:       985
 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:    496.0000    248.0000    248.0000
 Net charge:     0.00000      Expected multiplicity:    1
 Atoms:    117,  Basis functions:    985,  GTFs:   1803
 This is a restricted single-determinant wavefunction
 Orbitals from 1 to   248 are occupied
 Loading electron density functions (EDF) information from built-in EDF library
 (The library is available at https://github.com/zorkzou/Molden2AIM)
 Pt(   20)      Core electrons: 60     EDF primitive GTFs: 36
 The number of total inner-core electrons:    60
 The number of total EDF primitive GTFs:    36
 Loading EDF library finished!

 Loaded C:\Users\sober\Desktop\sp_bw_BPt_benz_H2act_HH.molden.input successfully!

 Formula: H56 B1 C42 F15 P2 Pt1      Total atoms:     117
 Molecule weight:      1113.71465 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)
100

               ============ Other functions (Part 1) ============
 0 Return
 1 Draw scatter graph between two functions and generate their cube files
 2 Export various files (mwfn/pdb/xyz/wfn/wfx/molden/fch/47/mkl...) or generate input file of quantum chemistry programs
 3 Calculate molecular van der Waals volume
 4 Integrate a function in whole space
 5 Show overlap integral between alpha and beta orbitals
 6 Monitor SCF convergence process of Gaussian
 7 Auxiliary tools for CP2K (CP2Kmate)
 8 Generate Gaussian input file with initial guess from fragment wavefunctions
 9 Evaluate interatomic connectivity and atomic coordination number
 11 Calculate overlap and centroid distance between two orbitals
 12 Biorthogonalization between alpha and beta orbitals
 14 Calculate LOLIPOP (LOL Integrated Pi Over Plane)
 15 Calculate intermolecular orbital overlap
 17 Generate Fock/KS matrix based on orbital energies and coefficients
 18 Yoshizawa's electron transport route analysis
 19 Generate new wavefunction by combining fragment wavefunctions
 20 Calculate Hellmann-Feynman forces
 21 Calculate properties based on geometry information for specific atoms
 22 Detect pi orbitals, set occupation numbers and calculate pi composition
 23 Fit function distribution to atomic value
2
 0 Return
                 Export system to various formats of files:
 1 Output current structure to .pdb file
 2 Output current structure to .xyz file
 3 Output current structure and atomic charges to .chg file
 4 Output current wavefunction as .wfx file
 5 Output current wavefunction as .wfn file
 6 Output current wavefunction as Molden input file (.molden)
 7 Output current wavefunction as .fch file
 8 Output current wavefunction as .47 file
 9 Output current wavefunction as old Molekel input file (.mkl)
 31 Output current structure to .cml file
 32 Output current wavefunction as .mwfn file
 33 Output current structure and cell information as .cif file
 34 Output current structure and cell information as .gro file
              Generate input file of quantum chemistry codes:
 10 Gaussian
 11 GAMESS-US
 12 ORCA         13 NWChem       14 MOPAC
 15 PSI4         16 MRCC         17 CFOUR
 18 Molpro       19 Dalton       20 Molcas
 21 Q-Chem
 25 CP2K         26 Quantum ESPRESSO      27 VASP (POSCAR)
8
 Input path for exporting file, e.g. C:\ltwd.47
t.47
 Exporting, please wait...

 If writing Fock matrix to .47 file?
 0 Do not write
 1 Generate Fock matrix based on MO energies and coefficients via F=SCE(C^-1) and write it to .47 file
1
 Generating Fock/KS matrix based on orbital energies, expansion coefficients and overlap matrix...
 Fock/KS matrix has been successfully generated!
 Generating electric dipole moment integral matrix...

 Exporting .47 file finished!

I have tested, the .47 file exported in this way can be normally analyzed by GENNBO 7.0.

By the way, the version updated on 2024-Jul-12 fixed a bug, namely the contraction coefficients of pure type D shells in the .47 file exported by Multiwfn are all zero.

#22 Re: Quantum Chemistry » jun-Cheap (jChS) composite method » 2024-07-11 15:16:52

Dear Saeed,

It is not a built-in method in G16. There is no evident necessarity to use this method in most cases.

Tian

#23 Re: Multiwfn and wavefunction analysis » sobEDA combined with Distortion Energy » 2024-07-11 01:41:35

Dear Mudit Jain,

I do not comment the result of ADF, I never use it. According to the original paper of sobEDA, this method works fairly well for wide variety of systems.

ETS-NOCV is supported by Multiwfn, which is able to further decompose orbital interaction energy given by sobEDA to provide deeper insight.

Best,

Tian

#24 Re: Multiwfn and wavefunction analysis » Crash while converting molden into .47 » 2024-07-11 01:38:23

Note that def2-SVP contains f basis functions for Pt, while def2-SV(P) doesn't. If def2-SV(P) also doesn't work, you can send me your ORCA input file and molden file to my E-mail, I will check.

#25 Re: Multiwfn and wavefunction analysis » sobEDA combined with Distortion Energy » 2024-07-10 06:09:13

Dear Mudit Jain,

If you want to fully decompose energy variation, you can calculate deformation energies and perform sobEDA analysis for every point of IRC to check how various physical components change along the reaction coordinate.

Best,

Tian

#27 Re: Multiwfn and wavefunction analysis » Very strange AIM results! » 2024-07-10 05:54:46

Dear Saeed,

I can confirm that O-B is much more covalent than O-Se. If you calculate fuzzy bond order for O-B, you will find the bond order is close to 0.9, clearly indicating its strong covalency. In contrast, the bond order of O-Se is not so large. As you mentioned, this is also supported by sobEDA analysis.

Regarding why some BCP properties such as Laplacian of rho give contradictory conclusion, you can understand it easily by plotting contour line map of Laplacian of rho, see below image for the case of O-B. It can be seen that there is an evident electron density concentration area, as highlighted by red circle, clearly indicating its prominent covalent character. However, since the BCP occurs in the region that far away from the actual bonding region, the Laplacian of rho at the BCP is positive. It obviously doesn't mean that the O-B interaction doesn't show evident covalent character.

Clipboard01.jpg

Best,

Tian

#28 Re: Multiwfn and wavefunction analysis » Installing Multiwfn on Linux » 2024-07-10 02:14:50

Dear Mudit,

I have a blog article to describe what is the electronic energy frequently referred in computational chemistry literatures: http://sobereva.com/488 (in Chinese, you can use google translator).
For DFT calculation (except for double-hybrid functionals) in Gaussian, electronic energy simply refers to the value after "SCF Done".

Best,

Tian

#29 Re: Multiwfn and wavefunction analysis » Crash while converting molden into .47 » 2024-07-10 02:10:42

The Multiwfn version you used is too old. If you use latest version, you will see much richer prompt.
Briefly speaking, for ORCA user, it is best to use a basis set without angular moment >=f. Note that most NBO analyses are quite insensitive to basis set, high angular moment basis functions are generally not needed.

#30 Re: Multiwfn and wavefunction analysis » modify function » 2024-07-09 05:28:45

Derivative at nuclear position cannot be obtained. In the case of good promolecular density, for each element, Multiwfn has an array recording radial density of its isolated state, and Lagrangian interpolation is used to obtain the value and derivatives at any radial distance. At nuclear position, the interpolation is not feasible, because it is the starting position of radial density.

In fact, it is well-known that electron density at nuclear position is ill-defined, because there is a cusp of electron density at this point (see https://en.wikipedia.org/wiki/Kato_theorem).

Board footer

Powered by FluxBB