Static and Dynamical Quantum Studies of CX3-AlX2 and CSiX3-BX2 (X = F, Cl, Br) Complexes with Hydrocyanic Acid: Unusual Behavior of Strong π-Hole at Triel Center

• Published in: International journal of molecular sciences 2023-04, Vol.24 (9), p.7881
• Main Authors: Michalczyk, Mariusz, Wojtkowiak, Kamil, Panek, Jarosław J., Jezierska, Aneta, Zierkiewicz, Wiktor
• Format: Article
• Language: English
• Subjects: Acids; BOMD; Dimerization; Electrostatic properties; Geometry; Germanium; Hydrocyanic acid; Molecular dynamics; molecular electrostatic potential; NCI; QTAIM; SAPT; triel bond
• ISSN: 1422-0067; 1661-6596; 1422-0067
• DOI: 10.3390/ijms24097881

The set of TX3-TrX2 (T = C, Si, Ge; Tr = B, Al, Ga; X = F, Cl, Br) molecules offers a rather unique opportunity to study both σ-hole and π-hole dimerization on the tetrel and triel ends, respectively. According to the molecular electrostatic potential (MEP) distribution, the π-hole extrema (acidic sites) were more intense than their σ-hole counterparts. The molecules owning the most (CX3-AlX2) and least (SiX3-BX2) intense π-holes were chosen to evaluate their capacities to attract one and two HCN molecules (Lewis bases). We discovered that the energetic characteristics of π-hole dimers severely conflict with the monomers MEP pattern since the weakest π-hole monomer forms a dimer characterized by interaction energy compared to those created by the monomers with noticeably greater power in the π-hole region. This outcome is due to the deformation of the weakest π-hole donor. Furthermore, the MEP analysis for monomers in the geometry of respective dimers revealed a “residual π-hole” site that was able to drive second ligand attachment, giving rise to the two “unusual trimers” examined further by the NCI and QTAIM analyses. Apart from them, the π-hole/π-hole and σ-hole/π-hole trimers have also been obtained throughout this study and described using energetic and geometric parameters. The SAPT approach revealed details of the bonding in one of the “unusual trimers”. Finally, Born-Oppenheimer Molecular Dynamics (BOMD) simulations were carried out to investigate the time evolution of the interatomic distances of the studied complexes as well as their stability.