Quasi-relativistic two-component computations of intermolecular dispersion energies

Advances and Prospects in Understanding London Dispersion Interactions in Molecular Chemistry

London dispersion (LD) interactions are the main contribution of the attractive part of the van der Waals potential. Even though LD effects are the driving force for molecular aggregation and recognition, the role of these omnipresent interactions in structure and reactivity had been largely underappreciated over decades. However, in the recent years considerable efforts have been made to thoroughly study LD interactions and their potential as a chemical design element for structures and catalysis. This was made possible through a fruitful interplay of theory and experiment. This review highlights recent results and advances in utilizing LD interactions as a structural motif to understand and utilize intra- and intermolecularly LD-stabilized systems. Additionally, we focus on the quantification of LD interactions and their fundamental role in chemical reactions.

TURBOMOLE: Modular program suite for ab initio quantum-chemical and condensed-matter simulations

TURBOMOLE is a collaborative, multi-national software development project aiming to provide highly efficient and stable computational tools for quantum chemical simulations of molecules, clusters, periodic systems, and solutions. The TURBOMOLE software suite is optimized for widely available, inexpensive, and resource-efficient hardware such as multi-core workstations and small computer clusters. TURBOMOLE specializes in electronic structure methods with outstanding accuracy-cost ratio, such as density functional theory including local hybrids and the random phase approximation (RPA), GW-Bethe-Salpeter methods, second-order Møller-Plesset theory, and explicitly correlated coupled-cluster methods. TURBOMOLE is based on Gaussian basis sets and has been pivotal for the development of many fast and low-scaling algorithms in the past three decades, such as integral-direct methods, fast multipole methods, the resolution-of-the-identity approximation, imaginary frequency integration, Laplace transform, and pair natural orbital methods. This review focuses on recent additions to TURBOMOLE's functionality, including excited-state methods, RPA and Green's function methods, relativistic approaches, high-order molecular properties, solvation effects, and periodic systems. A variety of illustrative applications along with accuracy and timing data are discussed. Moreover, available interfaces to users as well as other software are summarized. TURBOMOLE's current licensing, distribution, and support model are discussed, and an overview of TURBOMOLE's development workflow is provided. Challenges such as communication and outreach, software infrastructure, and funding are highlighted.

Benchmark Experimental Gas-Phase Intermolecular Dissociation Energies by the SEP-R2PI Method

The gas-phase ground-state dissociation energy D0( S0) of an isolated and cold bimolecular complex is a fundamental measure of the intermolecular interaction strength between its constituents. Accurate D0 values are important for the understanding of intermolecular bonding, for benchmarking high-level theoretical calculations, and for the parameterization of dispersion-corrected density functionals or force-field models that are used in fields ranging from crystallography to biochemistry. We review experimental measurements of the gas-phase D0( S0) and D0( S1) values of 55 different M⋅S complexes, where M is a (hetero)aromatic molecule and S is a closed-shell solvent atom or molecule. The experiments employ the triply resonant SEP-R2PI laser method, which involves M-centered ( S0 → S1) electronic excitation, followed by S1 → S0 stimulated emission spanning a range of S0 state vibrational levels. At sufficiently high vibrational energy, vibrational predissociation of the M⋅S complex occurs. A total of 49 dissociation energies were bracketed to within ≤1.0 kJ/mol, providing a large experimental database of accurate noncovalent interactions.

Can Combined Electrostatic and Polarization Effects Alone Explain the F···F Negative-Negative Bonding in Simple Fluoro-Substituted Benzene Derivatives? A First-Principles Perspective

The divergence of fluorine-based systems and significance of their nascent non-covalent chemistry in molecular assemblies are presented in a brief review of the field. Emphasis has been placed to show that type-I and -II halogen-centered F···F long-ranged intermolecular distances viable between the entirely negative fluorine atoms in some fluoro-substituted dimers of C6H6 can be regarded as the consequence of significant non-covalent attractive interactions. Such attractive interactions observed in the solid-state structures of C6F6 and other similar fluorine-substituted aromatic compounds have frequently been underappreciated. While these are often ascribed to crystal packing effects, we show using first-principles level calculations that these are much more fundamental in nature. The stability and reliability of these interactions are supported by their negative binding energies that emerge from a supermolecular procedure using MP2 (second-order Møller-Plesset perturbation theory), and from the Symmetry Adapted Perturbation Theory, in which the latter does not determine the interaction energy by computing the total energy of the monomers or dimer. Quantum Theory of Atoms in Molecules and Reduced Density Gradient Non-Covalent Index charge-density-based approaches confirm the F···F contacts are a consequence of attraction by their unified bond path (and bond critical point) and isosurface charge density topologies, respectively. These interactions can be explained neither by the so-called molecular electrostatic surface potential (MESP) model approach that often demonstrates attraction between sites of opposite electrostatic surface potential by means of Coulomb’s law of electrostatics, nor purely by the effect of electrostatic polarization. We provide evidence against the standalone use of this approach and the overlooking of other approaches, as the former does not allow for the calculation of the electrostatic potential on the surfaces of the overlapping atoms on the monomers as in the equilibrium geometry of a complex. This study thus provides unequivocal evidence of the limitation of the MESP approach for its use in gaining insight into the nature of reactivity of overlapped interacting atoms and the intermolecular interactions involved.