An in depth discussion of the theory underlying the LED scheme can be found elsewhere. In this section, we describe how this approach can be used to obtain an exact decomposition of the DLPNO-CCSD(T) binding energy of subsystems (with labels A, B, C, …) in a supramolecular complex into contributions originating from the interaction of pairs of fragments (with labels X, Y, Z, …), as shown schematically in Fig. . Each subsystem is divided into a number of fragments (denoted , , etc. for subsystem A, B, C, … respectively), resulting in a total of fragments. The DLPNO-CCSD(T) binding energy between M interacting subsystems can be computed using a supramolecular approach as: where is the total electronic energy of the system, is the electronic energy of subsystem A, and is the sum of the energy of the isolated noninteracting subsystems. The geometric preparation energy term (also called -strain energy‖) and the thermodynamic corrections to the Gibbs free energy can also be computed separately and added to when needed.