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-14 12:54:58

may01dz
Member
Registered: 2018-07-17
Posts: 54

binding energies

Hello;

I have read the following sentences:

(A further interesting comparison is that of the binding energies, which we have calculated with respect to the fragments, both charged and neutral. We find that the binding energy of carboplatin is higher in both cases, by AEc = 2.9 eV and AE=1.8 eV, respectively. These results reflect the higher instability of the dicarboxylate ligand in the gas phase, both in the charged and the neutral state.)

DOI: https://doi.org/10.1016/0009-2614(95)01099-2

My question is:

How I figured out the binding energies, and whether it is achievable with Multiwfn.

thanks you very much

Offline

#2 2020-03-14 15:36:13

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

Re: binding energies

Commonly, you need to do following things:
1 Optimize complex structure using medium calculation level
2 Calculate single point energy of complex using high calculation level (EAB)
3 Extract the two fragments from the optimized complex structure, and then respectively calculate their single point energies using high calculation level (EA and EB)
4 Then the binding energy can be calculated as: E_bind=EAB-EA-EB.

For H-bond dimer system, you can also use electron density at bond critical point to evaluate the binding energy due to H-bond interaction, thus only step 1 is needed (used to generate wavefunction of the complex). This calculation is quite easy by using Multiwfn, see Section 4.2.1 of latest version of Multiwfn manual for introduction and example, as well as my work J. Comput. Chem., 40, 2868–2881 (2019).

Offline

#3 2020-03-14 19:05:33

may01dz
Member
Registered: 2018-07-17
Posts: 54

Re: binding energies

and, the most stable complex has the highest negative E_bind value?

Offline

#4 2020-03-16 03:02:01

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

Re: binding energies

To compare stability, commonly you should compare free energy, the one with lowest free energy is the most (thermodynamic) stable.

Offline

#5 2020-03-16 17:04:49

may01dz
Member
Registered: 2018-07-17
Posts: 54

Re: binding energies

Gibbs free energy change?

sobereva wrote:

To compare stability, commonly you should compare free energy, the one with lowest free energy is the most (thermodynamic) stable.

Offline

#6 2020-03-17 04:31:16

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

Re: binding energies

Not "change", you just need to compute Gibbs free energy for each configuation, and then directly compare the values.

Offline

#7 2020-03-18 00:12:08

may01dz
Member
Registered: 2018-07-17
Posts: 54

Re: binding energies

good, thank you very much

sobereva wrote:

Not "change", you just need to compute Gibbs free energy for each configuation, and then directly compare the value.

Offline

#8 2021-01-24 06:29:22

fyzan
Member
Registered: 2020-01-25
Posts: 21

Re: binding energies

sobereva wrote:

Not "change", you just need to compute Gibbs free energy for each configuation, and then directly compare the values.

To compute accurate Gibbs free energy, is this right approach to compute frequency at medium level and add thermal/entropy corrections from this medium level frequency calculation to the single point energy calculated at higher level? For example opt and freq job at pbe0/def2-tzvp and then add thermal/entropy corrections to single point energy computed at PWPB95 D3 def2-QZVPP level will yield better Gibbs free energy.?

Offline

#9 2021-01-24 06:57:35

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

Re: binding energies

This is a very common and standard way of calculating accurate Gibbs free energy.

However, if the system contains many low frequency modes, it is also important to apply quasi-RRHO thermodynamic model, which is supported by Shermo code (http://sobereva.com/soft/shermo/)

Offline

#10 2021-01-24 07:45:42

fyzan
Member
Registered: 2020-01-25
Posts: 21

Re: binding energies

thanks professor, can the BSSE error energy from gaussian counterpoise method be also added to zero point and Gibbs free energies to get BSSE corrected Gibbs free energy?

Offline

#11 2021-01-24 09:10:11

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

Re: binding energies

“BSSE corrected Gibbs free energy” makes no sense. BSSE is defined for interfragment interaction, rather than for a single system. You can say e.g. "BSSE corrected binding Gibbs free energy"

Offline

#12 2021-01-26 08:48:07

fyzan
Member
Registered: 2020-01-25
Posts: 21

Re: binding energies

sobereva wrote:

This is a very common and standard way of calculating accurate Gibbs free energy.

However, if the system contains many low frequency modes, it is also important to apply quasi-RRHO thermodynamic model, which is supported by Shermo code (http://sobereva.com/soft/shermo/)

sir the result of G from text file is different from shm file generated by shermo.
I have given text file as input with content gaussian-freq-job.log;-2812.272184125518 and I got value of G as -2811.969688 au. then I gave shm file exported by shermo as input and I get value of G as -2812.4253519 a.u.

Offline

#13 2021-01-26 13:36:03

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

Re: binding energies

I don't understand your meaning. You do not need to export shm file by Shermo, you can directly read thermal correction to G from output of Shermo.

Note that if you request Shermo to employ RRHO model (ilowfreq=0), then the thermal corrections printed by Shermo are identical to those by Gaussian.

Offline

#14 2021-01-27 03:22:13

fyzan
Member
Registered: 2020-01-25
Posts: 21

Re: binding energies

sobereva wrote:

I don't understand your meaning. You do not need to export shm file by Shermo, you can directly read thermal correction to G from output of Shermo.

Note that if you request Shermo to employ RRHO model (ilowfreq=0), then the thermal corrections printed by Shermo are identical to those by Gaussian.

Thank you sir, I was exporting shm file to read ZPE corrections as from text file of multiple entries Shermo don't print ZPE corrections. The exported shm file takes electronic energy value from Gausssian log file and not the value provided in text file input. Manually changing electronic energy value in shm file printed same results.

Offline

Board footer

Powered by FluxBB