Many-Body Dispersion

Computational predictions of cocrystal formation: A benchmark study of 28 assemblies comparing five methods from high-throughput to advanced models

Cocrystals are assemblies of more than one type of molecule stabilized through noncovalent interactions. They are promising materials for improved drug formulation in which the stability, solubility, or biocompatibility of the active pharmaceutical ingredient (API) is improved by including a coformer. In this work, a range of density functional theory (DFT) and density functional tight binding (DFTB) models are systematically compared for their ability to predict the lattice enthalpy of a broad range of existing pharmaceutically relevant cocrystals. These range from cocrystals containing model compounds 4,4'-bipyridine and oxalic acid to those with the well benchmarked APIs of aspirin and paracetamol, all tested with a large set of alternative coformers. For simple cocrystals, there is a general consensus in lattice enthalpy calculated by the different DFT models. For the cocrystals with API coformers the cocrystals, enthalpy predictions depend strongly on the DFT model. The significantly lighter DFTB models predict unrealistic values of lattice enthalpy even for simple cocrystals.

Extension of the D3 and D4 London dispersion corrections to the full actinides series

Efficient dispersion corrections are an indispensable component of modern density functional theory, semi-empirical quantum mechanical, and even force field methods. In this work, we extend the well established D3 and D4 London dispersion corrections to the full actinides series, francium, and radium. To keep consistency with the existing versions, the original parameterization strategy of the D4 model was only slightly modified. This includes improved reference Hirshfeld atomic partial charges at the ωB97M-V/ma-def-TZVP level to fit the required electronegativity equilibration charge (EEQ) model. In this context, we developed a new actinide data set called AcQM, which covers the most common molecular actinide compound space. Furthermore, the efficient calculation of dynamic polarizabilities that are needed to construct CAB6 dispersion coefficients was implemented into the ORCA program package. The extended models are assessed for the computation of dissociation curves of actinide atoms and ions, geometry optimizations of crystal structure cutouts, gas-phase structures of small uranium compounds, and an example extracted from a small actinide complex protein assembly. We found that the novel parameterizations perform on par with the computationally more demanding density-dependent VV10 dispersion correction. With the presented extension, the excellent cost-accuracy ratio of the D3 and D4 models can now be utilized in various fields of computational actinide chemistry and, e.g., in efficient composite DFT methods such as r2SCAN-3c. They are implemented in our freely available standalone codes (dftd4, s-dftd3) and the D4 version will be also available in the upcoming ORCA 6.0 program package.

Distinction and Quantification of Noncovalent Dispersive and Hydrophobic Effects

The possibilities of comparing computational results of noncovalent interactions with experimental data are discussed, first with respect to intramolecular interactions. For these a variety of experimental data such as heats of formation, crystal sublimation heats, comparison with energy minimized structures, and spectroscopic data are available, but until now largely have not found widespread application. Early force field and QM/MP2 calculations have already shown that the sublimation heats of hydrocarbons can be predicted with an accuracy of ±1%. Intermolecular interactions in solution or the gas phase are always accompanied by difficult to compute entropic contributions, like all associations between molecules. Experimentally observed T∆S values contribute 10% to 80% of the total ∆G, depending on interaction mechanisms within the complexes, such as, e.g., hydrogen bonding and ion pairing. Free energies ∆G derived from equilibrium measurements in solution allow us to define binding increments ∆∆G, which are additive and transferable to a variety of supramolecular complexes. Data from more than 90 equilibrium measurements of porphyrin receptors in water indicate that small alkanes do not bind to the hydrophobic flat surfaces within a measuring limit of ∆G = ±0.5 kJ/mol, and that 20 functions bearing heteroatoms show associations by dispersive interactions with up to ∆G = 8 kJ/mol, roughly as a function of their polarizability. Aromatic systems display size-dependent affinities ∆G as a linear function of the number of π-electrons.

On the role of non-additive interactions in three-body recombination

Non-additive forces are a cornerstone of molecular spectroscopy and reaction dynamics. However, the relevance of non-additive forces in three-body recombination remains largely unexplored. In this work, we present a global study on the impact of non-additive interactions in three-body recombination: atom-atom-atom and ion-atom-atom. Our study explores these reactions in a wide range of energies, from the cold to the hyperthermal regime, finding no effect of non-additive interactions. Therefore, pair-wise interactions are enough to describe the three-body recombination dynamics adequately.

libMBD: A general-purpose package for scalable quantum many-body dispersion calculations

Many-body dispersion (MBD) is a powerful framework to treat van der Waals (vdW) dispersion interactions in density-functional theory and related atomistic modeling methods. Several independent implementations of MBD with varying degree of functionality exist across a number of electronic structure codes, which both limits the current users of those codes and complicates dissemination of new variants of MBD. Here, we develop and document libMBD, a library implementation of MBD that is functionally complete, efficient, easy to integrate with any electronic structure code, and already integrated in FHI-aims, DFTB+, VASP, Q-Chem, CASTEP, and Quantum ESPRESSO. libMBD is written in modern Fortran with bindings to C and Python, uses MPI/ScaLAPACK for parallelization, and implements MBD for both finite and periodic systems, with analytical gradients with respect to all input parameters. The computational cost has asymptotic cubic scaling with system size, and evaluation of gradients only changes the prefactor of the scaling law, with libMBD exhibiting strong scaling up to 256 processor cores. Other MBD properties beyond energy and gradients can be calculated with libMBD, such as the charge-density polarization, first-order Coulomb correction, the dielectric function, or the order-by-order expansion of the energy in the dipole interaction. Calculations on supramolecular complexes with MBD-corrected electronic structure methods and a meta-review of previous applications of MBD demonstrate the broad applicability of the libMBD package to treat vdW interactions.

Review of Applications of Density Functional Theory (DFT) Quantum Mechanical Calculations to Study the High-Pressure Polymorphs of Organic Crystalline Materials

Since its inception, chemistry has been predominated by the use of temperature to generate or change materials, but applications of pressure of more than a few tens of atmospheres for such purposes have been rarely observed. However, pressure is a very effective thermodynamic variable that is increasingly used to generate new materials or alter the properties of existing ones. As computational approaches designed to simulate the solid state are normally tuned using structural data at ambient pressure, applying them to high-pressure issues is a highly challenging test of their validity from a computational standpoint. However, the use of quantum chemical calculations, typically at the level of density functional theory (DFT), has repeatedly been shown to be a great tool that can be used to both predict properties that can be later confirmed by experimenters and to explain, at the molecular level, the observations of high-pressure experiments. This article's main goal is to compile, analyze, and synthesize the findings of works addressing the use of DFT in the context of molecular crystals subjected to high-pressure conditions in order to give a general overview of the possibilities offered by these state-of-the-art calculations.

Three-in-one: exploration of co-encapsulation of cabazitaxel, bicalutamide and chlorin e6 in new mixed cyclodextrin-crosslinked polymers

Three-in-one: a single bCyD polymer easily prepared in water is used to co-encapsulate cabazitaxel and bicalutamide with chlorin e6 affording a nanoplatform to implement multimodal cancer therapy.

The chemical bond in solids-revisited

This article complements an earlier topical review of the chemical bond (Jones 2018J. Phys.: Condens. Matter30153001), starting in the mid-19th century and seen from the perspective of a condensed matter physicist. The discussion of applications focused on the structure and properties of phase change materials. We review here additional aspects of chemistry, particularly some that have raised interest recently in this context. Concepts such as 'electron-rich', 'electron-deficient (excess orbital)', 'hypervalent', 'three-centre', and 'metavalent' bonds, and 'multicentre hyperbonding' are now found in the condensed matter literature. They are surveyed here, as well as the bond in metals and the 'Peierls' distortion. What are these concepts, are they related, and are they sometimes new labels for established, but unfamiliar ideas? 'Half bonds' and 'fractional valencies' play a central role in this discussion. It is remarkable that they were introduced 100 years ago, but ignored or forgotten, and have needed to be rediscovered more than once.

Current Status of Quantum Chemical Studies of Cyclodextrin Host-Guest Complexes

This article aims to review the application of various quantum chemical methods (semi-empirical, density functional theory (DFT), second order Møller-Plesset perturbation theory (MP2)) in the studies of cyclodextrin host-guest complexes. The details of applied approaches such as functionals, basis sets, dispersion corrections or solvent treatment methods are analyzed, pointing to the best possible options for such theoretical studies. Apart from reviewing the ways that the computations are usually performed, the reasons for such studies are presented and discussed. The successful applications of theoretical calculations are not limited to the determination of stable conformations but also include the prediction of thermodynamic properties as well as UV-Vis, IR, and NMR spectra. It has been shown that quantum chemical calculations, when applied to the studies of CD complexes, can provide results unobtainable by any other methods, both experimental and computational.

O(N) Stochastic Evaluation of Many-Body van der Waals Energies in Large Complex Systems

We propose a new strategy to solve the key equations of the many-body dispersion (MBD) model by Tkatchenko, DiStasio Jr., and Ambrosetti. Our approach overcomes the original O(N3) computational complexity that limits its applicability to large molecular systems within the context of O(N) density functional theory. First, to generate the required frequency-dependent screened polarizabilities, we introduce an efficient solution to the Dyson-like self-consistent screening equations. The scheme reduces the number of variables and, coupled to a direct inversion of the iterative subspace extrapolation, exhibits linear-scaling performances. Second, we apply a stochastic Lanczos trace estimator resolution to the equations evaluating the many-body interaction energy of coupled quantum harmonic oscillators. While scaling linearly, it also enables communication-free pleasingly parallel implementations. As the resulting O(N) stochastic massively parallel MBD approach is found to exhibit minimal memory requirements, it opens up the possibility of computing accurate many-body van der Waals interactions of millions-atoms' complex materials and solvated biosystems with computational times in the range of minutes.

Direct orbital selection within the domain-based local pair natural orbital coupled-cluster method

Domain-based local pair natural orbital coupled cluster (DLPNO-CC) has become increasingly popular to calculate relative energies (e.g., reaction energies and reaction barriers). It can be applied within a multi-level DLPNO-CC-in-DLPNO-CC ansatz to reduce the computational cost and focus the available computational resources on a specific subset of the occupied orbitals. We demonstrate how this multi-level DLPNO-CC ansatz can be combined with our direct orbital selection (DOS) approach [M. Bensberg and J. Neugebauer, J. Chem. Phys. 150, 214106 (2019)] to automatically select orbital sets for any multi-level calculation. We find that the parameters for the DOS procedure can be chosen conservatively such that they are transferable between reactions. The resulting automatic multi-level DLPNO-CC method requires no user input and is extremely robust and accurate. The computational cost is easily reduced by a factor of 3 without sacrificing accuracy. We demonstrate the accuracy of the method for a total of 61 reactions containing up to 174 atoms and use it to predict the relative stability of conformers of a Ru-based catalyst.

Can We Predict the Isosymmetric Phase Transition? Application of DFT Calculations to Study the Pressure Induced Transformation of Chlorothiazide

Isosymmetric structural phase transition (IPT, type 0), in which there are no changes in the occupation of Wyckoff positions, the number of atoms in the unit cell, and the space group symmetry, is relatively uncommon. Chlorothiazide, a diuretic agent with a secondary function as an antihypertensive, has been proven to undergo pressure-induced IPT of Form I to Form II at 4.2 GPa. For that reason, it has been chosen as a model compound in this study to determine if IPT can be predicted in silico using periodic DFT calculations. The transformation of Form II into Form I, occurring under decompression, was observed in geometry optimization calculations. However, the reverse transition was not detected, although the calculated differences in the DFT energies and thermodynamic parameters indicated that Form II should be more stable at increased pressure. Finally, the IPT was successfully simulated using ab initio molecular dynamics calculations.

Data-efficient machine learning for molecular crystal structure prediction

The combination of modern machine learning (ML) approaches with high-quality data from quantum mechanical (QM) calculations can yield models with an unrivalled accuracy/cost ratio. However, such methods are ultimately limited by the computational effort required to produce the reference data. In particular, reference calculations for periodic systems with many atoms can become prohibitively expensive for higher levels of theory. This trade-off is critical in the context of organic crystal structure prediction (CSP). Here, a data-efficient ML approach would be highly desirable, since screening a huge space of possible polymorphs in a narrow energy range requires the assessment of a large number of trial structures with high accuracy. In this contribution, we present tailored Δ-ML models that allow screening a wide range of crystal candidates while adequately describing the subtle interplay between intermolecular interactions such as H-bonding and many-body dispersion effects. This is achieved by enhancing a physics-based description of long-range interactions at the density functional tight binding (DFTB) level-for which an efficient implementation is available-with a short-range ML model trained on high-quality first-principles reference data. The presented workflow is broadly applicable to different molecular materials, without the need for a single periodic calculation at the reference level of theory. We show that this even allows the use of wavefunction methods in CSP.

Density functional theory predictions of the mechanical properties of crystalline materials

The DFT-predicted mechanical properties of crystalline materials are crucial knowledge for their screening, design, and exploitation.