Unraveling chemical bonding mechanisms through dipole moment variations under external electric fields

In a recent study [Angew. Chem. Int. Ed. 62, e202312078 (2023)], it was suggested that the maxima of the dipole polarizability along the bond dissociation coordinate can be used to define chemical bond cleavage. In this work, we systematically analyze this hypothesis by performing a spatial partition of the local derivatives of the electronic densities and dipole moments with respect to the external field along the bond axis, considering a set of chemically bonded and non-covalently bonded systems. Through this approach, we advance the characterization of three interaction regions: (I) the bonded region, where the system responds to the field as a single molecular entity; (II) the bond-breaking region, characterized by the additional formation of a dipole moment, in between the atoms participating in the bond, that opposes the total molecular dipole; and (III) the dissociated state, corresponding to the sum of independent dipoles. Furthermore, our decomposition allows us to identify one critical point that might be used to roughly mark the transition between region (II) and region (III). In particular, our decomposition provides a more exhaustive understanding of the second region within a first principles framework, unravelling the nature of a bond by its electronic response properties, rather than by analyzing the hybridization of the electronic structure.