## Definitions

Rather than use the phrase "dimerization energy", which is the energy difference specific to two chemically identical species coming together, such as $\ce{2 BH3 -> B2H6}$, it is better to use the phrase _[interaction energy](https://en.wikipedia.org/wiki/Interaction_energy)_. See [here](http://chemistry.stackexchange.com/questions/31/do-terminal-amides-participate-in-%cf%80-stacking/38#38) and [here](http://chemistry.stackexchange.com/questions/24326/what-tools-can-be-used-to-estimate-binding-energies/24376#24376) for a detailed distinction between interaction energies and binding energies. I will simplify a bit some of the text there and in [this](http://chemistry.stackexchange.com/questions/32560/significance-of-single-point-energy-when-calculating-interaction-energies) answer.

The _single point energy_ (SP energy, SPE) is the total energy of a system at one and only one configuration of the nuclear positions. In the case of two monomers (sometimes called fragments), the interaction energy will be the energy difference between the complex (sometimes called the supersystem) and the sum of the individual fragments or monomers _at the same nuclear positions as in the complex_. This gives rise to the equation

$$
\Delta E_{\text{int}}(\text{AB}) = E(\text{AB}) - \left[ E(\text{A}) + E(\text{B}) \right]
$$

where there are two molecules, A and B, leading to three calculations, one for each of the monomers and one for the AB complex. In the dimerization of borane from above, or your specific example of the $\ce{HF}$ dimer, assuming that the _internal coordinates_ (bond lengths, bond angles) of each monomer are identical, then the equation above can be simplified as A and B are identical:

$$
\Delta E_{\text{int}} = E_{\text{dimer}} - 2 E_{\text{monomer}}
$$

To take your specific example of the $\ce{HF}$ dimer, you mention the role of the "orientation of atoms" in both the single point and dimerization energies. For a single monomer, assuming again that internal coordinates are not changing, the relative positions of the nuclei are not changing, and the total energy is _invariant_ to translations and rotations. However, in the dimer or complex geometry, the relative orientation of one monomer to the other will change the total energy of the dimer, leading to a change in the interaction energy. For example, in addition to the linear form $\ce{F-H}\cdots\ce{F-H}$, you might be interested in the "flipped" versions $\ce{F-H}\cdots\ce{H-F}$ or $\ce{H-F}\cdots\ce{F-H}$, or forming a T-shape, etc. The interaction energy of each will tell you about the relative stability of that nuclear configuration relative to the other configurations.

The most important variables then in an interaction energy calculation are the relative orientation of the interaction species and the _distance between each of the fragments_. How do you decide this? Usually some intermolecular coordinate is scanned, such as the distance between closest atoms on the monomers, still while leaving the monomer internal coordinates fixed to give an _unrelaxed potential energy surface scan_.

## Calculation Examples

For example, here is the interaction energy for the $\ce{F-H}\cdots\ce{F-H}$ configuration with a variety of methods using the cc-pVTZ basis set and an intraatomic distance of 0.927 angstroms, scanning the $\ce{H}\cdots\ce{F}$ distance:

[![interaction energy plot: different methods][1]][1]

From a plot like this, we see that MP2 overestimates the interaction and CCSD underestimates it compared to CCSD(T); the shape of the B97-D potential is probably incorrect. HF is missing some repulsive interaction at long-range that the other methods include. It isn't dispersion, as the total dispersion interaction is always attractive and the minimum of the potential is qualitatively correct. You could certainly use any density functional or other method of your choice here (SCS-MP2, SOS-MP2, CISD, ...).

You could similar plots for other configurations, or choose a single method and place multiple configurations on the same plot (here, CCSD(T)/cc-pVTZ):

[![interaction energy plot: different conformers][2]][2]

From a plot like this, it's pretty clear that only the $\ce{F-H}\cdots\ce{F-H}$ configuration is ever bound; it would be interesting to compare the T-shaped configurations to the linear ones.

## Interaction Energy Decomposition

You mention "electrostatic monomeric interactions", and we also see some interesting long-range behavior in plot [1][1]. One of the most powerful techniques for understanding intermolecular interactions is [symmetry-adapted perturbation theory](http://www.psicode.org/psi4manual/master/sapt.html) (SAPT), which gives a decomposition of interaction energies in terms of electrostatics, Pauli exchange, polarization (called induction), dispersion, and charge transfer (though the definition of the latter is often ambiguous).

Here are each of the interaction energy components (SAPT0/jun-cc-pVTZ, monomer-centered basis):

[![SAPT0 energy decomposition plot: separate][3]][3]

Combining all the electrostatic (1), exchange (1), induction(3), and dispersion (2) terms gives

[![SAPT0 energy decomposition plot: together][4]][4]

This makes it pretty clear that all points, $\ce{F-H}\cdots\ce{F-H}$ is bound mostly by electrostatic attraction, and dispersion is negligible. Maybe you knew that, but I didn't. My guess is that the presence of dynamic correlation in all the wavefunction methods damps the long-range electrostatics compared to Hartree-Fock.

## References

### Interaction energies

* [J. Řezáč and P. Hobza, _Chem. Rev._, **2016**, _116_ (9), pp 5038–5071](http://dx.doi.org/10.1021/acs.chemrev.5b00526)

### SAPT

* [B. Jeziorski, R. Moszynski, and K. Szalewicz, _Chem. Rev._ **94**, 1887 (1994)](http://dx.doi.org/10.1021/cr00031a008)
* [E. G. Hohenstein and C. D. Sherrill, J. Chem. Phys. 132, 184111 (2010)](http://dx.doi.org/10.1063/1.3426316)
* [E. G. Hohenstein, R. M. Parrish, C. D. Sherrill, J. M. Turney, and H. F. Schaefer, _J. Chem. Phys._ **135**, 174107 (2011)](http://dx.doi.org/10.1063/1.3656681)
* [T. M. Parker, L. A. Burns, R. M. Parrish, A. G. Ryno, and C. D. Sherrill, _J. Chem. Phys._ **140**, 094106 (2014)](http://dx.doi.org/10.1063/1.4867135)

### Technical notes

All of these calculations were run with Psi4. None of the interaction or SAPT energies are CP-corrected.


  [1]: https://i.sstatic.net/tu00i.png
  [2]: https://i.sstatic.net/FhbWM.png
  [3]: https://i.sstatic.net/Osntu.png
  [4]: https://i.sstatic.net/6Ow7r.png