Frustrated Lewis pair-mediated hydro-dehalogenation: crucial role of non-covalent interactions

Organic halides stand as invaluable reagents with diverse applications in synthetic chemistry and various industrial processes. Despite their utility, concerns arise due to their inherent toxicity. Addressing these apprehensions, hydro-dehalogenation has emerged as a promising strategy involving the replacement of halogen atoms with hydrogen atoms to transform toxic organic halides into hydrocarbons. This study delves into the computational exploration of hydro-dehalogenation reactions of benzyl halide, mediated by frustrated Lewis pairs (FLPs), using density functional theory (DFT). The reactions entail the formation of FLP1 or FLP2 in the presence of TMP or lutidine with B(C₆F₅)₃, respectively. This is followed by heterolytic cleavage of dihydrogen and subsequent reaction with benzyl halides. Non-covalent interaction analysis underscores the significance of π-π stacking and CH-π interactions in stabilizing transition states. Additionally, the activation strain model (ASM) dissects activation energies, revealing the substantial impact of strain energy on reaction barriers. Energy decomposition analysis (EDA) offers insights into the contributions of electrostatic, orbital, and dispersion energies to the overall attractive interaction energy. The investigation extends to hydro-dehalogenation reactions of ethyl halides, uncovering distinct mechanisms and activation barriers. This comprehensive analysis illuminates the intricacies of hydro-dehalogenation reactions, providing valuable insights into their mechanisms and paving the way for future studies in this field.

Geometry optimizations were carried out at the M06-2X/def2-SVP level of theory, which was performed using the Gaussian 16 program. Solvent-corrected single-point energies were also calculated using the polarizable continuum model (PCM) at the PCM(chloroform)-M06-2X/def2-TZVP//M06-2X/def2-SVP level of theory. The Gibbs free energy correction was determined from computations performed at the M06-2X/def2-SVP level of theory. Principal interacting orbital (PIO) analysis was conducted using the NBO 6.0 software. The nature of bonding in the respective transition state (TS) structures was analyzed using atoms-in-molecules (AIM) analyses. Additionally, the presence of non-covalent interactions (NCI) was exemplified using Multiwfn software.