Dispersion Control over Molecule Cohesion: Exploiting and Dissecting the Tipping Power of Aromatic Rings

ConspectusWe have learned over the past years how London dispersion forces can be effectively used to influence or even qualitatively tip the structure of aggregates and the conformation of single molecules. This happens despite the fact that single dispersion contacts are much weaker than competing polar forces. It is a classical case of strength by numbers, with the importance of London dispersion forces scaling with the system size. Knowledge about the tipping points, however difficult to attain, is necessary for a rational design of intermolecular forces. One requires a careful assessment of the competing interactions, either by sensitive spectroscopic techniques for the study of the isolated molecules and aggregates or by theoretical approaches. Of particular interest are the systems close to the tipping point, when dispersion interactions barely outweigh or approach the strength of the other interactions. Such subtle cases are important milestones for a scale-up to realistic multi-interaction situations encountered in the fields of life and materials science. In searching for examples that provide ideal competing interactions in complexes and small clusters, aromatic systems can offer a diverse set of molecules with a variation of dispersion and electrostatic forces that control the dominant and peripheral interactions. Our combined spectroscopic and theoretical investigations provide valuable insights into the balance of intermolecular forces because they typically allow us to switch the aromatic substituent on and off. High-resolution rotational spectroscopy serves as a benchmark for molecular structures, as correct calculations should be based on correct geometries. When discussing the competition with other noncovalent interactions, obvious competitors are directional hydrogen bonds. As a second counterweight to aryl interactions, we will discuss aurophilic/metallophilic interactions, which also have a strong stabilization with a small number of atoms involved. Vibrational spectroscopy is most sensitive to interactions of light atoms, and the competition of OH hydrogen bonds with dispersion forces in a molecular aggregate can be judged well by the OH stretching frequency. Experiments in the gas phase are ideal for gauging the accuracy of quantum chemical predictions free of solvent forces. A tight collaboration utilizing these three methods allows experiment vs experiment vs theory benchmarking of the overall influence of dispersion in molecular structures and energetics.

The hydration of an oxy-polycyclic aromatic compound: the case of naphthaldehyde

The study of the intermolecular interactions of polycyclic aromatic compounds, considered as important pollutants of the Earth's atmosphere since they are emitted by the partial combustion of fuels, is essential to understand the formation and aging of their aerosols. In this study, the hydration of α-naphthaldehyde and β-naphthaldehyde isomers was investigated through a combination of Fourier transform microwave spectroscopy and quantum chemical calculations. Monohydrate structures were observed experimentally for both isomers, with two hydrate structures observed for β-naphthaldehyde and only one for α-naphthaldehyde, consistent with computational predictions. Analysis of the monohydrate structures indicated that the β-isomer exhibits higher hydrophilicity compared to the α-isomer, supported by electronic densities, hydration energies, and structural considerations. Further computational calculations were conducted to explore the planarity of the naphthaldehyde hydrates. Different levels of theory were employed, some of these revealing slight deviations from planarity in the hydrate structures. Low-frequency out-of-plane vibrational modes were examined, and the inertial defect was used to assess the planarity of the hydrates. The results suggested that the hydrates possess a predominantly planar structure, in agreement with the highest level of computational calculations and the absence of c-type transitions in the experimental spectra. Additionally, calculations were extended to dihydrate structures by attaching two water molecules to the naphthaldehyde isomers. The most stable dihydrate structures were predicted to be combinations of the observed monohydrate positions. However, experimental observation of the most stable dihydrate structures was challenging due to their very low vapour pressure, calling for complementary experiments using laser ablation nozzles.

Accurate measurement of sequential Ar desorption energies from the dispersion-dominated Ar1-3 complexes of aromatic molecules

We present experimental determination of the energies associated with the gradual desorption of Ar atoms from the aromatic molecular surface. Non-covalently bound 2,2'-pyridylbenzimidazole-Ar_1-3 complexes were produced in the gas phase and characterized using resonant two-photon ionization (R2PI) spectroscopy. The single Ar desorption from the PBI-Ar, PBI-Ar₂ and PBI-Ar₃ complexes were measured as 581 ± 18, 656 ± 30 and 537 ± 31 cm^-1, respectively. The energies were bracketed between the last observed band in the respective R2PI spectra and the disappeared intramolecular modes of PBI. The Ar_n dissociation energies in the S₁ state were measured as 581 ± 18, 1237 ± 48 and 1774 ± 79 cm^-1, respectively, for n = 1, 2 and 3. The calculated dissociation energies of the respective complexes, obtained using three computational methods, show excellent agreement with the experimental data. The ground state dissociation energies were estimated by subtracting the Δν shift of the origin band, and the respective values are 541 ± 18, 1160 ± 48 and 1634 ± 79 cm^-1. Overall, the calculated values resulted in scaling factors ranging from 0.956 to 1.017, which depict the predictive power of the methods to determine dispersion energies. The current investigation describes a unique methodology to accurately determine the dissociation and desorption energies of Ar atoms from the surfaces of bio-relevant aromatic molecules.

A computational inspection of the dissociation energy of mid-sized organic dimers

The gas-phase value of the dissociation energy ( D0) is a key parameter employed in both experimental and theoretical descriptions of noncovalent complexes. The D0 data were obtained for a set of mid-sized organic dimers in their global minima which was located using geometry optimizations that applied ample basis sets together with either the conventional second-order Møller–Plesset (MP2) method or several dispersion-corrected density-functional theory (DFT-D) schemes. The harmonic vibrational zero-point (VZP) and deformation energies from the MP2 calculations were combined with electronic energies from the coupled cluster theory with singles, doubles, and iterative triples [CCSD(T)] extrapolated to the complete basis set (CBS) limit to estimate D0 with the aim of inspecting values that were most recently measured, and an analogous comparison was performed using the DFT-D data. In at least one case (namely, for the aniline⋯methane cluster), the D0 estimate that employed the CCSD(T)/CBS energies differed from experiment in the way that could not be explained by a possible deficiency in the VZP contribution. Curiously, one of the DFT-D schemes (namely, the B3LYP-D3/def2-QZVPPD) was able to reproduce all measured D0 values to within 1.0 kJ/mol from experimental error bars. These findings show the need for further measurements and computations of some of the complexes. In order to facilitate such studies, the physical nature of intermolecular interactions in the investigated dimers was analyzed by means of the DFT-based symmetry-adapted perturbation theory.

A combined spectroscopic and computational investigation on dispersion-controlled docking of Ar atoms on 2-(2'-pyridyl)benzimidazole

The dispersion-controlled docking of inert Ar atoms on the face of polycyclic 2-(2'-pyridyl)-benzimidazole (PBI) was studied experimentally aided by computational findings. The PBI-Ar_n (n = 1-3) complexes were produced in a supersonically jet-cooled molecular beam and probed using resonant two-photon ionization coupled with a time-of-flight mass spectrometric detection scheme and laser-induced fluorescence spectroscopy. The ground state vibrational frequencies were obtained from single vibronic level fluorescence spectroscopy. The formation of multiple isomers was verified using UV-UV hole-burning spectroscopy. The geometries of PBI-Ar_n (n = 1-3) complexes were derived by mutual agreement between experimental findings and computational results such as vibrational frequencies in the ground and excited electronic states, Franck-Condon factors and spectral shift of the S₁← S₀ transitions. All the above analyses provided good agreement between the experimental and simulated spectrum with the most stable PBI-Ar_n (n = 1-3) clusters. The highest intermolecular interaction between PBI and Ar was obtained with the Ar atom positioned above the imidazolyl ring. A second Ar atom was preferably docking on the other side of the imidazolyl ring than the pyridyl ring. The subsequent addition of the third Ar atom preferred the position above the pyridyl ring. The current investigation can be useful to investigate the preferential docking of dispersion-controlled interacting partners in multifunctional aromatic side chains present in biological systems.