Multiwfn forum

Multiwfn official website: http://sobereva.com/multiwfn. Multiwfn forum in Chinese: http://bbs.keinsci.com/wfn. E-mail of admin: sobereva[at]sina.com

You are not logged in.

#1 2020-03-19 10:20:59

Mehdi
Member
Registered: 2020-03-19
Posts: 11

energy decomposition analysis using combined fragment wavefunctions

Dear Admin
I am doing energy decomposition analysis using the section”4.100.8 Perform simple energy decomposition by using combined fragment wavefunctions” in Multiwfn package.
I run a single point calculation with Gaussian09, revision A.02. My input file contains fragment-combined wavefunction as initial guess. Part of my input file is as follows:
When I run this input file, I receive an error (Error termination via Lnk1e in C:\G09W\l401.exe
). Could you please kindly guide me in this problem? I appreciate in advance for the time you devote to answer my question.
Best regards,
Mehdi Shakourian


Chk=new.chk
%mem=700mb
%nprocshared=2
#p rcam-b3lyp/gen pseudo=read nosymm guess=cards

Please check this file to ensure validity

  2  1
C     1.90992500    3.16721400   -0.50060600
C     0.54789600    2.74660600   -0.44522800
C     0.19538600    1.36319800   -0.49336500
C     1.25707900    0.46851800   -0.75010000
C     2.59624300    0.86019300   -0.71804900
C     2.93638500    2.21919900   -0.59908800
C    -1.20085300    0.98099500   -0.29788200
C    -1.66491700   -0.33704700   -0.09830500
C    -1.05101200   -2.62529500   -0.15795600
C     0.02904800   -3.53159000   -0.41070300
C     1.40536500   -3.15465800   -0.60944800
C     1.87884200   -1.81648000   -0.80490800
C     3.22699900   -1.49410400   -0.76687600
C     3.60293700   -0.12763500   -0.78042900
C    -0.47035900    3.75364400   -0.35159700
C    -1.83057600    3.38527900   -0.31726600
C    -2.20234800    1.99958600   -0.28581900
C    -3.59026600    1.67436500   -0.23011600
C    -3.99860200    0.33507800   -0.19172200
C    -3.01935600   -0.67180300   -0.13135900
C    -2.37847600   -3.02318500   -0.20008200
C    -3.39102800   -2.03386800   -0.12896800
C     2.26924200    4.54648900   -0.43784800
C     4.30754300    2.59176200   -0.60208700
C     4.65641100    3.95365500   -0.53107300
C     3.63124900    4.93130100   -0.45735500
C     5.31259000    1.60285900   -0.67171400
C     4.95884500    0.23362700   -0.73382300
C     5.95028900   -0.77580700   -0.73003900
C     5.57480400   -2.14614800   -0.72175100
C     4.21039800   -2.50612800   -0.71503000
C     2.40959500   -4.17137600   -0.63363800
C     3.79869800   -3.84685100   -0.68822600
C    -0.32664800   -4.91577500   -0.44955900
C     2.03733700   -5.54933900   -0.59725100
C     0.68804700   -5.91503900   -0.54820300
C    -2.71940600   -4.38906800   -0.31263300
C    -1.68935800   -5.33710800   -0.39764900
C    -4.74500400   -2.40228400   -0.17770400
C    -5.09348300   -3.77203100   -0.24082600
C    -4.07947200   -4.76513900   -0.30740600
C    -5.37102000   -0.03357400   -0.19102100
C    -5.74064800   -1.39578600   -0.18152000
C    -4.59566600    2.68650100   -0.23268800
C    -5.96609200    2.33228500   -0.21534100
C    -6.35873000    0.96926400   -0.20312000
C    -2.83779700    4.39378800   -0.27848200
C    -4.21587000    4.04367300   -0.24931000
C    -0.10823800    5.13257200   -0.32917700
C    -1.11058900    6.12683400   -0.28094300
C    -2.47288900    5.75843800   -0.26504900
C     1.25817400    5.52530700   -0.35923900
C     6.56651700   -3.16191800   -0.71505300
C     7.94916700   -2.76323900   -0.71782100
H     8.70976500   -3.53694200   -0.71794400
C     8.30788100   -1.45534400   -0.71437400
H     9.35728200   -1.17873700   -0.71009300
C     7.32198900   -0.40538000   -0.71254400
C     4.80199800   -4.85975200   -0.69108300
C     6.15412100   -4.49824900   -0.70493000
H     6.90486200   -5.28263700   -0.70200600
C     3.06419000   -6.57130700   -0.63971400
C     4.38550800   -6.22209900   -0.67683600
H     5.14305000   -6.99936600   -0.69842100
C     0.31584100   -7.31593300   -0.56888500
C     1.36468100   -8.31080200   -0.62344100
H     1.07229700   -9.35552200   -0.63703600
C     2.65986900   -7.96004300   -0.64758400
H     3.43860500   -8.71483200   -0.68094400
C    -1.00071700   -7.68129000   -0.53153300
H    -1.26397200   -8.73440200   -0.55736400
C    -2.04620700   -6.71645900   -0.45667500
C    -3.39654000   -7.08540400   -0.44291300
H    -3.65037900   -8.13986400   -0.49304900
C    -4.42507300   -6.14048100   -0.37478400
C    -6.46534800   -4.14190100   -0.25970400
C    -6.78869400   -5.54402200   -0.31898100
H    -7.83467500   -5.83339600   -0.32417400
C    -5.82048600   -6.49220000   -0.37134400
H    -6.08838700   -7.54257000   -0.41703700
C    -7.43580900   -3.13846000   -0.23061400
H    -8.48408600   -3.42530100   -0.24374300
C    -7.11206700   -1.76557700   -0.20159100
C    -8.09085800   -0.73904900   -0.20340900
H    -9.14122700   -1.01579900   -0.21063200
C    -7.74400100    0.60004400   -0.20093900
C    -8.72656900    1.65923800   -0.20155500
H    -9.77676100    1.38705800   -0.19781000
C    -6.96625200    3.35315000   -0.21314300
C    -8.35622700    2.95686200   -0.20528600
H    -9.10796200    3.73934700   -0.20404100
C    -5.21662200    5.05187500   -0.23009200
C    -6.57593200    4.67911500   -0.21595600
H    -7.33176300    5.45855900   -0.20706100
C    -3.49214700    6.77391200   -0.23072100
C    -4.81679700    6.41093400   -0.22111500
H    -5.58212500    7.18089900   -0.20122400
C    -0.74158600    7.51769500   -0.25329700
C    -1.78972800    8.51265900   -0.21459400
H    -1.49917200    9.55794300   -0.19419200
C    -3.08799600    8.16171000   -0.20688400
H    -3.86542700    8.91823200   -0.18039300
C     0.58543200    7.87152200   -0.26483100
H     0.85898700    8.92201700   -0.23810900
C     1.61512600    6.89992900   -0.31782700
C     2.97697100    7.26295300   -0.33523700
H     3.23770400    8.31639700   -0.29758700
C     3.98147100    6.31576100   -0.39879400
C     6.03911000    4.33229200   -0.53555700
C     5.38205700    6.67350800   -0.41085100
H     5.63769300    7.72726600   -0.36594200
C     6.35467200    5.74082400   -0.47333100
H     7.39916700    6.03490200   -0.47877300
C     6.68388400    1.97311700   -0.65352300
C     7.01286700    3.35143200   -0.59266500
H     8.05965900    3.64127600   -0.58675200
C     7.65507000    0.95005100   -0.68438900
H     8.70529500    1.22994300   -0.67268100
O     1.00595500   -0.82777600   -1.05517500
O    -0.79162200   -1.34393700    0.14890500
Hg    1.57958200   -0.64661800    3.20888800

C H O 0
  6-31G**
  ****
Hg     0
S   3   1.00
     48.013786990            0.58008164912E-02     
     21.239875095           -0.17328165235   
     15.876100879            0.36416685034   
S   1   1.00
      5.5180531820           1.0000000       
S   1   1.00
      1.5099145479           1.0000000       
S   1   1.00
      0.70683482159          1.0000000       
S   1   1.00
      0.16266265576          1.0000000       
S   1   1.00
      0.56960589112E-01            1.0000000       
P   4   1.00
     17.500000000           -0.86457187331E-01     
     15.252594855            0.15327318234   
      6.4404715169          -0.30274955249   
      1.8180159920           0.54361285026   
P   1   1.00
      0.90067981801          1.0000000       
P   1   1.00
      0.41304540122          1.0000000       
P   1   1.00
      0.11845702655          1.0000000       
P   1   1.00
      0.36087619349E-01            1.0000000       
D   4   1.00
     10.028197701            0.39922353165E-01     
      7.5920661493          -0.94774230522E-01     
      1.9144256118           0.31164157463   
      0.88641552102          0.45567078406   
D   1   1.00
      0.38154235320          1.0000000       
D   1   1.00
      0.14880740086          1.0000000       
F   1   1.00
      1.5521700              1.0000000       
F   1   1.00
      0.5243100              1.0000000       
G   1   1.00
      1.2895000              1.0000000       
****

HG     0
HG-ECP     3     60
f-ul potential
  1
2      3.88579112            30.36499643       
s-ul potential
  3
2     12.98154870           275.73721174       
2      6.49077440            49.08921249       
2      3.88579112           -30.36499643       
p-ul potential
  3
2     10.53809580           241.54007398       
2      5.26904790            27.39659081       
2      3.88579112           -30.36499643       
d-ul potential
  3
2      8.10172051           127.86700761       
2      4.05086026            16.60831151       
2      3.88579112           -30.36499643       

5(E16.9)
-1
! Orbital:     1 Occ:  2.000000 from fragment   1
0.100000000E-04 0.500000000E-04-0.100000000E-04-0.300000000E-04 0.000000000E+00
-0.110000000E-02 0.150000000E-03 0.301000000E-02 0.150000000E-03-0.200000000E-04
0.000000000E+00 0.100000000E-04 0.300000000E-04 0.000000000E+00 0.000000000E+00
0.000000000E+00-0.300000000E-04 0.000000000E+00-0.200000000E-04 0.000000000E+00
0.157000000E-02 0.117000000E-02-0.162000000E-02-0.110000000E-03 0.200000000E-04
0.300000000E-04 0.100000000E-04 0.000000000E+00 0.000000000E+00-0.100000000E-04
0.000000000E+00 0.100000000E-03 0.500000000E-04 0.000000000E+00-0.200000000E-04
-0.101000000E-02-0.283000000E-02-0.128000000E-02 0.130000000E-03-0.400000000E-04
-0.800000000E-04 0.000000000E+00 0.100000000E-04 0.000000000E+00-0.100000000E-04
-0.500000000E-04 0.300000000E-04-0.130000000E-03-0.480000000E-03-0.170000000E-03
0.234000000E-02 0.250000000E-02 0.406000000E-02 0.850000000E-03 0.300000000E-04
0.160000000E-03-0.900000000E-04 0.900000000E-04 0.300000000E-04 0.110000000E-03
0.200000000E-04 0.230000000E-03-0.300000000E-04 0.400000000E-04 0.100000000E-04
-0.340000000E-02 0.36

Offline

#2 2020-03-19 11:04:06

sobereva
Tian Lu (Multiwfn developer)
From: Beijing
Registered: 2017-09-11
Posts: 1,611
Website

Re: energy decomposition analysis using combined fragment wavefunctions

Please use latest version of Multiwfn on Multiwfn official site (i.e. 3.7(dev) currently). If it still doesn't work, please directly upload entire Gaussian input file as attachment, I will check the reason.

Offline

#3 2020-03-21 20:34:32

Mehdi
Member
Registered: 2020-03-19
Posts: 11

Re: energy decomposition analysis using combined fragment wavefunctions

Dear Admin,
Thank you so much for your quick reply. I used the last version of Multiwfn 3.7. Unfortunately, I see the error again. The input file has been attached. I appreciate if you see the error.
Best regards,
Mehdi Shakourian

Offline

#4 2020-03-22 05:47:07

sobereva
Tian Lu (Multiwfn developer)
From: Beijing
Registered: 2017-09-11
Posts: 1,611
Website

Re: energy decomposition analysis using combined fragment wavefunctions

I don't find your attachment...

Offline

#5 2020-03-22 09:11:36

Mehdi
Member
Registered: 2020-03-19
Posts: 11

Re: energy decomposition analysis using combined fragment wavefunctions

Dear Admin,
Thank you so much for your quick reply. I used the last version of Multiwfn 3.7. Unfortunately, I see the error again. The input file was too large. Thus, the file was divided into seven sections. Please use the seven links below, download the files, put them in a folder and extract. I appreciate if you see the error.
Best regards,
Mehdi Shakourian
File.part01.rar
File.part02.rar
File.part03.rar
File.part04.rar
File.part05.rar
File.part06.rar
File.part07.rar

Offline

#6 2020-03-22 14:21:26

sobereva
Tian Lu (Multiwfn developer)
From: Beijing
Registered: 2017-09-11
Posts: 1,611
Website

Re: energy decomposition analysis using combined fragment wavefunctions

You should replace 5(E16.9) in the gjf file by (5E16.9).

I think the Multiwfn you are using is not the latest version on Multiwfn website, because this problem has been identified about one month ago and fixed.

In addition, when you calculating the Hg, plate and then the complex, please always add "5d" keyword. (Gaussian uses 6d type of basis function when 6-31G* is employed, while uses 5d type of basis function when "gen" keyword is employed, this causes evident inconsistency, while using "5d" keyword can fully eliminiate this problem).

If after doing above modifications, the calculation of complex is still unsuccessful, please upload the gjf file for Hg and that for plate, I will check.

Offline

#7 2020-03-31 10:29:24

Mehdi
Member
Registered: 2020-03-19
Posts: 11

Re: energy decomposition analysis using combined fragment wavefunctions

Dear Admin
Thank you so much for your help. I used the modifications in the input files but I couldn't solve the problem. I appreciate if you see the input files (attached files) and check it again.

Best regards
Mehdi Shakourian-Fard

inputs.part04.rar
inputs.part03.rar
inputs.part02.rar
inputs.part01.rar
inputs.part05.rar
inputs.part06.rar

Offline

#8 2020-03-31 17:33:20

sobereva
Tian Lu (Multiwfn developer)
From: Beijing
Registered: 2017-09-11
Posts: 1,611
Website

Re: energy decomposition analysis using combined fragment wavefunctions

In the new.gjf, you should use genecp instead of gen, otherwise the pseudopotential information will not be loaded.

Offline

#9 2020-04-01 08:02:25

Mehdi
Member
Registered: 2020-03-19
Posts: 11

Re: energy decomposition analysis using combined fragment wavefunctions

Dear Admin
The problem was solved with "genecp" keyword. I really appreciate for your help.

Best regards,
Mehdi Shakourian

Offline

#10 2020-04-01 09:58:18

Mehdi
Member
Registered: 2020-03-19
Posts: 11

Re: energy decomposition analysis using combined fragment wavefunctions

Dear Admin,
I calculated the ΔEtot, ΔEorb and ΔEsteric as follows:
ΔEtot= -387.98
ΔEorb= -373.24
ΔEsteric = -14.74
I think ΔEsteric should be positive, while i have calculated it with negative value. Am i right? How is the sign of ΔEdis in energy decomposition analysis?
Best regards,
Mehdi Shakourian

Offline

#11 2020-04-01 22:46:21

sobereva
Tian Lu (Multiwfn developer)
From: Beijing
Registered: 2017-09-11
Posts: 1,611
Website

Re: energy decomposition analysis using combined fragment wavefunctions

The -14.74 corresponds to the sum of ΔE_electrostatic and ΔE_exchange (see Section 3.100.8 of Multiwfn manual), in principle it can be positive or negative, only in some papers the authors labelled it as ΔE_steric for discussion convenience; note that the word "steric" in this context is not consistent with the steric effect in common sense (in fact the steric effect is represented by the ΔE_exchange). So, I suggest you report -14.74 in your work as ΔE_electrostatic + ΔE_exchange instead of ΔE_steric, which is somewhat misleading.

The ΔE_dis must be negative. You can obtain it in terms of DFT-D dispersion correction, as mentioned in Section 3.100.8 of Multiwfn manual.

Offline

#12 2020-04-02 08:07:06

Mehdi
Member
Registered: 2020-03-19
Posts: 11

Re: energy decomposition analysis using combined fragment wavefunctions

Dear Admin
Thank you so much for all your help.

Best regards,
Mehdi Shakourian-Fard

Offline

#13 2021-01-05 17:24:02

Mehdi
Member
Registered: 2020-03-19
Posts: 11

Re: energy decomposition analysis using combined fragment wavefunctions

Dear Admin,
I have plotted a RDG surface for my structure. We are considering the interaction of Choline chloride (ChCl):Urea complex with Ag, Au and Cu clusters. I know that the interaction between the Cl- anion in the ChCl:Urea complex with nanoclusters is attractive in nature (blue color or electrostatic interaction, whereas we can see some red colors (repulsive force). I appreciate if you give us a reasonable explanation for seeing the red observation?
Picture2.zip


Bset regards,
Mehdi shakourian-Fard

Offline

#14 2021-01-05 18:23:28

sobereva
Tian Lu (Multiwfn developer)
From: Beijing
Registered: 2017-09-11
Posts: 1,611
Website

Re: energy decomposition analysis using combined fragment wavefunctions

Dear Mehdi shakourian-Fard,

Please simply ignore the red color, it never means that the interaction is repulsive, it just comes from the drawback of the definition of sign(lambda2)rho function. In the circular red color region, the sign of lambda2 happens to change to positive, this observation doesn't have any noteworthy chemical significance.

If you slightly adjust isovalue of RDG to make the isosurface slightly shrink, the red edge may disappear.

By the way, the grid spacing of your RDG analysis is too large, therefore the isosurfaces look unsmooth. I suggest properly decreasing the grid spacing to make the graph finer.

Best regards,

Tian

Offline

#15 2021-01-05 18:49:52

Mehdi
Member
Registered: 2020-03-19
Posts: 11

Re: energy decomposition analysis using combined fragment wavefunctions

Dear Tian
Thank you so much for your help.

Best regards,
Mehdi Shakourian-Fard

Offline

#16 2021-01-10 18:09:53

Mehdi
Member
Registered: 2020-03-19
Posts: 11

Re: energy decomposition analysis using combined fragment wavefunctions

Dear Tian,
We have investigated the interaction of Mn (M = Cu, Ag, and Au; n =1-4) nanoparticles with ChCl:Urea molecule. We have done the EDA using Multiwfn program. In the energy decomposition analysis of the complex (attached file), why the electrostatic+exchange repulsion terms are completely positive for the Au nanoparticles, while it is mostly negative for the Cu and Ag systems. What is structurally different among these systems that contribute to this trend.
Thank you in advance for the time you devote to answer my question.

Picture1.zip

Best regards,
Mehdi Shakourian-Fard

Offline

#17 2021-01-10 19:29:27

sobereva
Tian Lu (Multiwfn developer)
From: Beijing
Registered: 2017-09-11
Posts: 1,611
Website

Re: energy decomposition analysis using combined fragment wavefunctions

Dear Mehdi,

It is really difficult to give you a definitive answer without more information, there are many possible reasons. For example, incorrect setting of basis set and pseudopotential in input file, the SCF was not converged to stable wavefunction, incorrect calculation step, etc.

Best regards,

Tian

Offline

#18 2021-01-10 21:29:21

Mehdi
Member
Registered: 2020-03-19
Posts: 11

Re: energy decomposition analysis using combined fragment wavefunctions

Dear Tian
Thank you for your quick response. We are sure about the basis sets and the calculations. We have used the 6-311++G(d,p) for C, H, N, O and Cl- as well as Def2-TZVP for metal atoms. These basis sets have used for these interactions in literature. You have already said me that this term can be positive and negative in the EDA analysis. In fact, I want to know when the Eelectrostatic+exchange term can be positive or negative.  If there if a reason for that, please let me know.
Best Regards,
mehdi

Offline

#19 2021-01-11 11:42:51

sobereva
Tian Lu (Multiwfn developer)
From: Beijing
Registered: 2017-09-11
Posts: 1,611
Website

Re: energy decomposition analysis using combined fragment wavefunctions

It is really difficult to give an explanation of your observation without additional information. It is best to send me your optimized structures via E-mail, also providing your Gaussian input files will be useful.

Offline

Board footer

Powered by FluxBB