drops below some threshold value τ BE (where N cons is the total number of constraints). One can then obtain the BE energy value using the 1 and 2-PDMs of the fragments. The accuracy of the calculation will ultimately depend on the choice of the embedding fragments and the high-level solver.