Abstract: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.

Dynamic Properties of Proton Transfer in the Anharmonic-Interaction Hydrogen Bond Systems

Proton soliton in anharmonic interaction approximation

Quantum Molecular Dynamics Study on the Proton Exchange, Ionic Structures, and Transport Properties of Warm Dense Hydrogen-Deuterium Mixtures.

Hydrogen Bond Dynamics in the Solvation Shell on Proton Transfer in Aqueous Solution.

Understanding the Essential Nature of the Hydrated Excess Proton Through Simulation and Interpretation of Recent Spectroscopic Experiments

Perturbative Treatment of Anharmonic Vibrational Effects on Bond Distances: an Extended Langevin Dynamics Method.

Geometric Characteristics and Energy Landscapes of Halogen–water–hydrogen Bridges at Protein–ligand Interfaces

Mapping Hydrogen Positions along the Proton Transfer Pathway in an Organic Crystal by Computational X-ray Spectra

Structural and Dynamic Properties of Solvated Hydroxide and Hydronium Ions in Water from Ab Initio Modeling.

Simulating and investigating various dynamic aspects of $\rm{H}_2\rm{O}$-related hydrogen bond model

Probing the Structure of Stretched Molecular Ions with High-Harmonic Spectroscopy

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

Neural-network-based molecular dynamics simulations reveal that proton transport in water is doubly gated by sequential hydrogen-bond exchange

The Solvation Effect on the Rattling Behaviour of the Hydrated Excess Proton in Water

Experimental quantification of electrostatics in X-H···π hydrogen bonds

Dynamics of Hydrogen Bonds Coupling on the Specific DNA-Protein Interactions

Concerted motion of protons in hydrogen bonds of DNA-type molecules

Theoretical Study of Binding Interactions and Vibrational Raman Spectra of Water in Hydrogen-Bonded Anionic Complexes: (H2o)n- (n = 2 and 3), H2O···X-(X = F, Cl, Br, and I), and H2O···M- (M = Cu, Ag, and Au)

Anisotropy of the proton kinetic energy as a tool for capturing structural transition in nanoconfined H$_2$O

Modeling of the Response of Hydrogen Bond Properties on an External Electric Field: Geometry, NMR Chemical Shift, Spin-Spin Scalar Coupling