How Substitution Combines with Non-Covalent Interactions to Modulate 1,4-Naphthoquinone and Its Derivatives Molecular Features-Multifactor Studies

Unraveling the Nature of Hydrogen Bonds of "Proton Sponges" Based on Car-Parrinello and Metadynamics Approaches

The nature of intra- and intermolecular non-covalent interactions was studied in four naphthalene derivatives commonly referred to as "proton sponges". Special attention was paid to an intramolecular hydrogen bond present in the protonated form of the compounds. The unsubstituted "proton sponge" served as a reference structure to study the substituent influence on the hydrogen bond (HB) properties. We selected three compounds substituted by methoxy, amino, and nitro groups. The presence of the substituents either retained the parent symmetry or rendered the compounds asymmetric. In order to reveal the non-covalent interaction properties, the Hirshfeld surface (HS) was computed for the crystal structures of the studied compounds. Next, quantum-chemical simulations were performed in vacuo and in the crystalline phase. Car-Parrinello molecular dynamics (CPMD), Path Integral Molecular Dynamics (PIMD), and metadynamics were employed to investigate the time-evolution changes of metric parameters and free energy profile in both phases. Additionally, for selected snapshots obtained from the CPMD trajectories, non-covalent interactions and electronic structure were studied. Quantum theory of atoms in molecules (QTAIM) and the Density Overlap Regions Indicator (DORI) were applied for this purpose. It was found based on Hirshfeld surfaces that, besides intramolecular hydrogen bonds, other non-covalent interactions are present and have a strong impact on the crystal structure organization. The CPMD results obtained in both phases showed frequent proton transfer phenomena. The proton was strongly delocalized in the applied time-scale and temperature, especially in the PIMD framework. The use of metadynamics allowed for tracing the free energy profiles and confirming that the hydrogen bonds present in "proton sponges" are Low-Barrier Hydrogen Bonds (LBHBs). The electronic and topological analysis quantitatively described the temperature dependence and time-evolution changes of the electronic structure. The covalency of the hydrogen bonds was estimated based on QTAIM analysis. It was found that strong hydrogen bonds show greater covalency, which is additionally determined by the proton position in the hydrogen bridge.

Comprehensive Empirical Model of Substitution—Influence on Hydrogen Bonding in Aromatic Schiff Bases

In this work, over 500 structures of tri-ring aromatic Schiff bases with different substitution patterns were investigated to develop a unified description of the substituent effect on the intramolecular hydrogen bridge. Both proximal and distal effects were examined using Density Functional Theory (DFT) in the gas phase and with solvent reaction field (Polarizable Continuum Model (PCM) and water as the solvent). In order to investigate and characterize the non-covalent interactions, a topological analysis was performed using the Quantum Theory of Atoms In Molecules (QTAIM) theory and Non-Covalent Interactions (NCI) index. The obtained results were summarized as the generalized, empirical model of the composite substituent effect, assessed using an additional group of simple ring-based Schiff bases. The composite substituent effect has been divided into separate increments describing the different interactions of the hydrogen bridge and the substituent: the classical substituent effect, involving resonance and induction mediated through the ring, steric increment based on substituent proximity to the bridge elements, and distal increment, derived from substitution on the distal ring.