Multiwfn forum

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

You are not logged in.

#1 Re: Multiwfn and wavefunction analysis » Multicenter Bond Index for aromatic macrocycles » Yesterday 11:14:28

Please read Section 3.11.10 of Multiwfn manual, the theoretical background has been carefully introduced.

#2 Re: Multiwfn and wavefunction analysis » Multicenter Bond Index for aromatic macrocycles » Yesterday 00:30:30

For a large ring, please only use multi-center bond order (alternatively, use AV1245), do not use multicenter DI (only supports up to 10 atoms). Your Multiwfn version is new enough.

#3 Re: Multiwfn and wavefunction analysis » Multicenter Bond Index for aromatic macrocycles » 2024-11-22 21:40:39

Please do not use very old version. Due to my special algorithm for Multicenter Bond Index (see "Appendix: The extremely efficient implementation of MCBO in Multiwfn" in Section 3.11.2 of latest Multiwfn manual), current version of Multiwfn is able to calculate Multicenter Bond Index for arbitrarily large rings (even more than 100 atoms!)

#4 Re: Multiwfn and wavefunction analysis » Calculating the contributions of TSCT and TBCT in a molecule » 2024-11-22 21:34:50

Frankly speaking, I don't have any knowledge about TSCT and TBCT. The author Cheng Zhong had developed a script and shared in http://ccc.keinsci.com/thread-18597-1-1.html, which can realize this analysis. Please contact him via E-mail for more information about the script, or use Google translator to try to understand content of his post.

#5 Re: Multiwfn and wavefunction analysis » regarding conceptual dft analysis and "orbital weighted" indicators » 2024-11-21 21:37:06

I don't know why you intend to use a very small delta value, perhaps you only want to consider the contribution from the "most" frontier MO? If yes, you can gradually decrease the delta value from the default value 0.1 a.u., each time you can choose "5 Print current orbital weights used in orbital-weighted (OW) calculation" to check the weights.

Taking examples\benzene.fch as example, if the delta is decreased to 0.01 a.u., you will see the following output after choosing option 5

10 Highest weights in orbital-weighted f+
Orbital    22 (LUMO  )   Weight:  50.00 %   E_diff:     3.410 eV
Orbital    23 (LUMO+1)   Weight:  50.00 %   E_diff:     3.410 eV
Orbital    24 (LUMO+2)   Weight:   0.00 %   E_diff:     5.656 eV
Orbital    25 (LUMO+3)   Weight:   0.00 %   E_diff:     7.147 eV
Orbital    26 (LUMO+4)   Weight:   0.00 %   E_diff:     7.147 eV
Orbital    27 (LUMO+5)   Weight:   0.00 %   E_diff:     7.789 eV
Orbital    28 (LUMO+6)   Weight:   0.00 %   E_diff:     8.202 eV
Orbital    29 (LUMO+7)   Weight:   0.00 %   E_diff:     8.202 eV
Orbital    30 (LUMO+8)   Weight:   0.00 %   E_diff:     8.389 eV
Total weight of above listed orbitals: 100.00 %

10 Highest weights in orbital-weighted f-
Orbital    21 (HOMO  )   Weight:  50.00 %   E_diff:    -3.410 eV
Orbital    20 (HOMO-1)   Weight:  50.00 %   E_diff:    -3.410 eV
Orbital    19 (HOMO-2)   Weight:   0.00 %   E_diff:    -5.924 eV
Orbital    18 (HOMO-3)   Weight:   0.00 %   E_diff:    -5.924 eV
Orbital    17 (HOMO-4)   Weight:   0.00 %   E_diff:    -6.498 eV
Orbital    16 (HOMO-5)   Weight:   0.00 %   E_diff:    -8.011 eV
Orbital    15 (HOMO-6)   Weight:   0.00 %   E_diff:    -8.011 eV
Orbital    14 (HOMO-7)   Weight:   0.00 %   E_diff:    -8.639 eV
Orbital    13 (HOMO-8)   Weight:   0.00 %   E_diff:    -9.084 eV
Orbital    12 (HOMO-9)   Weight:   0.00 %   E_diff:   -10.752 eV
Total weight of above listed orbitals: 100.00 %

Clearly, now only the degenerate two HOMOs and the two LUMOs contribute to orbital weighted Fukui function or dual descriptor.

Best,

Tian

#6 Re: Multiwfn and wavefunction analysis » regarding conceptual dft analysis and "orbital weighted" indicators » 2024-11-20 23:00:21

Hello,

I don't exactly know which data you are talking about, could you provide a screenshot?

Best,

Tian

#7 Re: Multiwfn and wavefunction analysis » Stress tensor analysis now is available in Multiwfn » 2024-11-19 03:23:39

Dear Asdrubal Lozada,

Very nice, I will mentioning your script in Multiwfn manual in the corresponding part of Section 2.7 of manual.

Because the patch may or may not work for other releases of Multiwfn source code, I think it is best to also upload the Multiwfn source code package (version 3.8(dev) Update [2024-Nov-13]) onto https://github.com/aslozada/Stress_tensor/, so that other people can always use your script successfully.

Best regards,

Tian

#8 Re: Multiwfn and wavefunction analysis » on Fukui functions-hirshfeld charge » 2024-11-18 21:21:25

Dear Sedsil,

Just put the three files in current folder, boot up Multiwfn and load any one of them, then enter main function 22, choose "2 Calculate various quantitative indices" to calculate various indices including condensed Fukui function, or choose "3 Calculate grid data of Fukui function, dual descriptor and related functions" to visualize Fukui function.

Best,

Tian

#9 Multiwfn and wavefunction analysis » FOD analysis has been supported in Multiwfn » 2024-11-16 20:02:31

sobereva
Replies: 0

Fractional occupation number weighted electron density (FOD) proposed by Grimme is an important and very convenient method to visually detect and quantify static correlation, this method has been supported by Multiwfn since the version updated on 2024-Nov-13. Please check Section 4.A.7 of new Multiwfn manual for introduction and how to perform the analysis.

Note only for isolated systems, but also for periodic systems this method can be used, molden file produced by CP2K should be used in the latter case.

#10 Re: Multiwfn and wavefunction analysis » Photoinduced Charge Transfer Analysis » 2024-11-16 08:19:32

Mrfahim061 wrote:

Sir, as you previously mentioned, I followed 'Section 3.21.1' and 'Section 4.18.1' for plotting the Photo-Induced Electron Transfer (PET) analysis diagram. However, this procedure only provides a map and does not include the energy state transition figure as illustrated in the attached diagram. I would greatly appreciate it if you could specify which software is used to generate such diagrams.

Additionally, I have another query regarding the Charge Decomposition Analysis (CDA) diagram in Multiwfn. These diagrams often include energy levels indicated in blue. Is it possible to plot the CDA diagram without displaying these energy values?

https://i.postimg.cc/BXfpBBTq/Screenshot-2024-07-30-210852.jpg
https://i.postimg.cc/8sdGHw24/dislin.png


The energy level diagram was plotted manually by the authors (using such as powerpoint), not generated by Multiwfn.

Choose option "6 Disable labelling orbital indices" then replot.

#11 Re: Multiwfn and wavefunction analysis » How to find method/basis for an existing WFN file? » 2024-11-14 22:49:00

Dear Saeed,

Unfortunately, it is not possible to directly derive the method and basis set used for generate the .wfn file.

Best,

Tian

#12 Re: Multiwfn and wavefunction analysis » Forcing cubes with same grid » 2024-11-14 20:06:58

Dear Camps,

When Multiwfn asks you to select a way to set up grid, you can choose option "8 Use grid setting of another cube file" and input path of an existing cube file, then the newly calculated grid data will share exactly the same grid setting with the existing cube file.

Best,

Tian

#13 Re: Multiwfn and wavefunction analysis » Stress tensor analysis now is available in Multiwfn » 2024-11-12 20:30:39

Multiwfn has been updated today. Stress tensor ellipticity has been added as the 118th user-defined function.

#14 Re: Multiwfn and wavefunction analysis » sobEDA analysis » 2024-11-11 14:57:59

Dear Duc,

The atom indices in the two fragments are not contigouous. Please look at atom 198, it is in the large molecule, but it was set to fragment 2 (which should completely correspond to the small molecule).

In addition, implicit solvation model should not be used when performing sobEDA. The solvent contribution should be manually calculated, please check the paragraph "Solvation effect: influence of solvation effect on ..." in page 7025 of sobEDA original paper. In other words, sobEDA.sh can only be directly used to carry out energy decomposition analysis in gas phase.

Best,

Tian

#15 Re: Multiwfn and wavefunction analysis » Calculation of electrostatic potential when using basis with pseudopot » 2024-11-09 10:01:26

Hello,

Multiwfn is able to correctly calculate ESP when fch file was generated with use of ECP. Don't worry about it.

best,

Tian

#16 Re: Multiwfn and wavefunction analysis » Sum of electron density inside an ED isosurface » 2024-11-09 09:58:56

A relatively easy way is:
(1) Using main function 5 to calculate grid data of electron density with grid quality as fine as possible
(2) Then enter main function 13, choose "15 If data value is within certain range, set it to a specified value", input 0,0.025, and then input 0. Now electron density lower than 0.025 a.u. has been set to zero.
(3) Finally, choose "17 Show statistic data of grid points in specific spatial and value ranges", then input "a", now the "Integral of all data" shown on screen is the number of electrons enclosed by the 0.025 a.u. isosurface.

The data obtained in the aforementioned way is calculated by uniform integration grids, note that if there are very heavy atoms, this integration will be inaccurate even if very fine grid is used because of the very sharp electron density around nuclei. In that case the best way of obtaining this quantity is using basin analysis module and use mixed-grid to perform the integration, but modification of source code is needed.

#17 Re: Quantum Chemistry » Calculation of rho, grad(rho) and lap(rho) » 2024-11-09 09:47:21

Dear Prasanta,

You can use Multiwfn to easily calculate these function values and compare between different functionals.
You need to use a quantum chemistry or first-principle code to generate wavefunction files first, such as Gaussian, ORCA, NWChem, GAMESS-US, xtb, CP2K and so on, see Section 2.5 of Multiwfn manual for full list of supported wavefunction file format.
You can perform convergence test of the quantity of interest with respect to number of integration points, and thus to determine the least acceptable quality of integration grid.

Best,

Tian

#18 Re: Multiwfn and wavefunction analysis » Add marks on curve or graphe » 2024-11-07 13:50:08

Dear Jean-Pierre,

Multiwfn doesn't support plotting different spectrum curves in different styles (solid, dash, dot, etc.) to distinguish them. You need to export the curve data and plot via e.g. Origin. Alternatively, manually using e.g. Powerpoint to indicate different curves by arrows and labels.

Best,

Tian

#19 Re: Multiwfn and wavefunction analysis » sobEDA analysis » 2024-11-04 13:45:37

I need your full input files to figure out the reason

#20 Re: Multiwfn and wavefunction analysis » Stress tensor analysis now is available in Multiwfn » 2024-11-02 15:29:12

Dear aslozada,

I will incorporate your code into Multiwfn release several days later and let you know here.

Best,

Tian

#21 Re: Multiwfn and wavefunction analysis » Automatic selection of ELF basins » 2024-11-02 15:23:12

Dear Marcos Verissimo Alves,

Some suggestions:

1 In option "-5 Export basins as cube file", you can use command "a" to export basin.cub, which contains information of all basins.

2 Your quantum region is large, but the reactive region should be much smaller than this size. In the step of "Please select a method for setting up grid", you can consider to choose "8 Set center position, grid spacing and box length", and properly set the center and length of the box so that only reactive region is covered by the grids, the computational cost should be significantly lower than computing the entire quantum region.

3 To avoid calculating grid data twice, you can first use main function 5 to generate grid data for the reactive region, and export the grid data as a cube file, and then return to main menu, enter main function 17 and select 1 to generate basins, then choose "2: Generate the basins by using the grid data stored in memory". This run you can obtain indices of the basins of interest as usual. In the next time of running Multiwfn, you can use the cube file as input file, then still choose "2: Generate the basins by using the grid data stored in memory" when performing basin anaylsis.

Best regards,

Tian

#22 Re: Multiwfn and wavefunction analysis » ETS-NOCV and diffuse basis functions » 2024-11-02 15:01:03

Dear Georg,

Don't worry, you can safely use def2-TZVP for this case.

Best,

Tian

#23 Re: Multiwfn and wavefunction analysis » CHELPG charges for excited EOM-EE-CCSD states (Q-Chem) » 2024-10-26 08:05:47

I don't know, I am not a Q-Chem user. In principle, if Q-Chem is able to store natural orbitals into a file that recognized by Multiwfn (formats like .fch, .molden, .wfn, .mwfn, etc.), then CHELPG charges can be calculated at the corresponding level by Multiwfn.

#25 Re: Quantum Chemistry » Sobtop » 2024-10-22 06:51:38

Dummy atom is irrelvant to sobtop. Only real atoms can occur in the topology file created by sobtop.

#27 Multiwfn and wavefunction analysis » Stress tensor analysis now is available in Multiwfn » 2024-10-20 16:07:43

sobereva
Replies: 6

New feature of Multiwfn updated on 2024-Oct-20:

Stress tensor defined by Bader now can be outputted for any position in main function 1, and can be outputted for critical points by option 7 in topology analysis module. Stress tensor stiffness and stress tensor polarizability are shown together, and they are also available as the 116th and 117th user-defined functions, respectively, see corresponding part of Section 2.7 of updated Multiwfn manual for details. In addition, stiffness of electron density is available as the 115th user-defined function.

#28 Re: Quantum Chemistry » Plotting a Charge Displacement Curve » 2024-10-20 01:30:47

s0fk1 wrote:

Dear Tian, maybe I am missing something important here but I asked for:

5 Output average data of XY planes in a range of Z to output.txt

I inputed to span z from -8 to 8

and the prompt on the screen says:

There are      125 layers within the range
Column 1,2,3 correspond to X,Y,value respectively

So how would I know to what Z each X,Y pair with its corresponding electron density value correspond to, so I can plot the elecytron density vs the z-axis?

This option is not what you need in the present context. You should use "18 Plot (local) integral curve or plane-averaged in X/Y/Z direction"

#29 Re: Quantum Chemistry » Plotting a Charge Displacement Curve » 2024-10-18 15:23:14

Please look at the prompt on screen, Multiwfn always clearly shows meaning of each column when exporting a plain text file.

#30 Re: Quantum Chemistry » Plotting a Charge Displacement Curve » 2024-10-18 10:50:05

Please check Section 4.13.6 of Multiwfn manual, there is an example of plotting charge displacement curve.

If you need to obtain "the point along z′ where the electron densities of the noninteracting fragments become equal.", you can use wavefunction file of fragment 1 to ask Multiwfn to plot plane-averaged curve of electron density of fragment 1, and then repeat the steps for fragment 2, and finally manually find the z' position by plotting the two curves together in e.g. Origin.

Board footer

Powered by FluxBB

Baidu
map