In general you can differentiate between binding energy and interaction energy.

Binding energy is the difference between the energy of the complex minus the energy of the isolated monomers in their minima configuration.

Interaction energy is the difference between the energy of the complex minus the energy of the isolated monomers in the geometry of the complex.

The difference between binding and interaction energy is the distortion energy of the monomers.

The binding and interaction energies can be corrected or not from the Basis Set Superposition Error (BSSE).

In general you can differentiate between binding energy and interaction energy.

Binding energy is the difference between the energy of the complex minus the energy of the isolated monomers in their minima configuration.

Interaction energy is the difference between the energy of the complex minus the energy of the isolated monomers in the geometry of the complex.

The difference between binding and interaction energy is the distortion energy of the monomers.

The binding and interaction energies can be corrected or not from the Basis Set Superposition Error (BSSE).

the interaction energy ∆E(R) = EAB(R) _ EA _ EB, in which EA and EB represent the energy of the optimized cation and anion. The interaction energy, Eint (r), for two A and B systems can simply be given as:

Eint (r) =EAB(A…B) – EAB(A + B) -------------(1)

EAB(A…B) AB is the energy of the A...B system (total energy of ion pairs) at the (r) distance while EAB(A + B) is the energy of the two isolated components (cation and anion), at infinity (r =∞) .

But this equation neglect the energy of cation and anion when they exists in ionic liquid and takes into account the energy of cation and anion alone at (r=∞).Therefore, corrected this equation to eliminated this error by using counterpoise method (CP) and become:

Eint =EAB(A…B) – EAB(A + B) + ∆ECP -------------(2)

Where

∆ECP =[EA(A+B) – EA(A.B)] + [EB(A+B) –EB(A.B)] -------------(3)

EA(A+B), EB(A+B):energy of isolated A and B both alone (r=∞).

EA(A.B), EB(A.B):energy of both A and B alone but in ion pair.

okay, you can calculate BSSE in easy way by put separation between each two complex 10 angestrom and get the energy or do it in normal way, I do not know what is the problem to get BSSE since you have large system, you should correct the interaction energy

Many thanks to all for  the answers; actually, I want to compare the corrected energies using the BSSE-CP correction to the uncorrected ones, and my question was how can calculate the uncorrected energies (without using the counterpoise technique). I have previously optimized my systems, now to obtain the energy of each monomer I perform a SP calculation at the same level of theory or I should optimize each monomer alone.

Thank you

If you want to obtain the interaction energy, you should used the monomers in the geometry of the complex (no optimization of the monomers).

If you want to report the binding energy, you should optimized monomers.

Many thanks Pr Ibon Alkorta

first can calculation the minimization the geometry of the complex (ex: A---B) , then delete the A monomer from complex and used single point calculation for B monomer in complex and thus result the energy of B monomer in complex. the same steps for A monomer in complex. or used the ghost atom option    

Tanks a lot Ali Hamed.

Dear Prof. Ibon Alkorta can you please give us some important and useful references on this topic

Thank you

Regards

The steps for predicting a binding energy is relatively straightforward in computational chemistry but the process itself can become extremely difficult depending on

· The types of molecules involved (size, complexity, etc.)

· The desired accuracy of your final answer

· The environment that the molecules reside in (gas phase, aqueous phase, etc.)

For very small molecules in the gas phase, the process is easy and routine assuming that electron correlation is trivial. Aqueous phase environments, complex/large molecules and high accuracy quickly makes the problem more difficult if not impossible.

I have written up a previous answer which addresses how to predict an interaction energy (Eint) which is completely different from a binding energy (Ebind). However, much of the concepts will remain the same. I will try to outline here what in fact is the difference between an interaction energy and a binding energy. Unfortunately, 'binding energy' is tossed around a lot in the literature and consequently grossly misused. I hope to clarify the subtle differences between these two terms in my answer to your question.

I will focus my discussion around an elementary system of two water molecules to make this conceptually easier. In the gas phase, an isolated water molecule will adopt a particular geometry (at 0K and in a vacuum). We can determine what this minimum energy configuration will be at a particular level of theory by approximating the Schrodinger equation. The geometry that gives rise to the lowest possible energy for an isolated water molecule will be the geometry that you will observe under those conditions. Remember that the internal coordinates of water include things like the O-H bonds, and the H-O-H angle. We don't need to guess at what these parameters will be in the end, solving the Schrodinger equation will determine that.

Now, we know exactly what an isolated water molecule looks like and what the corresponding energy of that water molecule will be. Both components are equally important when it comes to determining a binding energy because (as the name suggests) we are trying to determine an energy. So, the energies that we start off with are crucially important.

We can now consider what happens when two water molecules interact (form a hydrogen bond). What will be the binding energy of these two water molecules (i.e. what will be the strength of this bond that forms between then?). It is common to approximate the strength of this bond by simply computing an interaction energy (its computationally easier to determine Eint rather than Ebind). How would we compute an interaction energy? In this case, you wouldn't need to figure out what geometry or energy an isolated water molecule in its lowest energy configuration. You would simply need to determine what the geometry of two interacting water molecules (i.e. dimer) would be (and corresponding energy) as well as the energy of each water molecule isolated from each other in the geometry that it adopts in the dimer.

I'll explain. We can determine the lowest energy configuration of a single water molecule in the gas phase. This geometry is inevitably different than a water that is not isolated... a water that is interacting with something else such as another water. If you want to get an interaction energy, you only need...

· Energy of the lowest energy configuration of the water dimer

· Energy of each water (in the water dimer) isolated from each other but only at the geometry that the water takes on within the dimer.

In this case, you are NOT allowing each isolated water molecule to 'relax' into its lowest energy configuration. You are simply going to neglect the contributions of this energy relaxation effect to the overall bond strength. Therefore,

Eint=Eopt-dimer−(Emon1+Emon2)Eint=Eopt-dimer−(Emon1+Emon2)

where mon1 is the first water molecule in your dimer and mon2 is the second water molecule in your dimer. (opt-dimer means optimized dimer)

However, in the case of a binding energy, you MUST account for the relaxation effects that occur when you break the dimer apart and allow the waters to become isolated from each other (and everything else). To do this, you simply carry out what I've stated above... you simply optimize each water and then take its corresponding energy. (Determining an interaction energy means you don't have to do this extra optimization(s)... you would simply do a single point energy on each monomer). Clearly for a water dimer, you would only do one optimization as each water will become identical when isolated. For a heterodimer, you would have to optimize each piece separately.

To determine the binding energy, you would use the following equation:

Ebind=Eopt-dimer−(Eopt-mon1+Eopt-mon2)Ebind=Eopt-dimer−(Eopt-mon1+Eopt-mon2)

Of course the problem becomes more complex when taking into account basis-set superposition error (you will need to perform counterpoise corrections though there exists a heated debate in the literature about the usefulness of this type of correction). Also, zero-point vibrational energy corrected binding energies will require that the zero-point energy of the dimer and each monomer will need to be determined (which comes from taking the 2nd derivative of the wavefunction which is computationally demanding).

As far as tools go, you can choose from several available quantum software packages. Gaussian 09 is a popular choice but very expensive. Q-Chem (paid) has the benefit of having the most beautiful, gorgeous, well-laid out documentation of anything I've ever seen (makes for a good supplemental textbook to boot). NWChem is opensource. Google can help you find plenty of options.

A quick note: Basis set superposition error (BSSE) should NOT be referred to as an error as 'error' suggests deviation from the RIGHT answer. BSSE has been shown to actually make non-corrected answers even WORSE.

One last important note: When determining differences of energies (such as you do to get an interaction or binding energy), ALL energies being compared must be evaluated at the SAME LEVEL OF THEORY (method/basis set, etc.).