An accurate determination of the electron correlation energy is an essential prerequisite for describing the structure, stability, and function in a wide variety of systems. Therefore, the development of efficient approaches for the calculation of the correlation energy (and hence the dispersion energy as well) is essential and such methods can be coupled with many density-functional approximations, local methods for the electron correlation energy, and even interatomic force fields. In this work, we build upon the previously developed many-body dispersion (MBD) framework, which is intimately linked to the random-phase approximation for the correlation energy. We separate the correlation energy into short-range contributions that are modeled by semi-local functionals and long-range contributions that are calculated by mapping the complex all-electron problem onto a set of atomic response functions coupled in the dipole approximation. We propose an effective range-separation of the coupling between the atomic response functions that extends the already broad applicability of the MBD method to non-metallic materials with highly anisotropic responses, such as layered nanostructures. Application to a variety of high-quality benchmark datasets illustrates the accuracy and applicability of the improved MBD approach, which offers the prospect of first-principles modeling of large structurally complex systems with an accurate description of the long-range correlation energy.

The sixth blind test of organic crystal structure prediction (CSP) methods has been held, with five target systems: a small nearly rigid molecule, a polymorphic former drug candidate, a chloride salt hydrate, a co-crystal and a bulky flexible molecule. This blind test has seen substantial growth in the number of participants, with the broad range of prediction methods giving a unique insight into the state of the art in the field. Significant progress has been seen in treating flexible molecules, usage of hierarchical approaches to ranking structures, the application of density-functional approximations, and the establishment of new workflows and `best practices' for performing CSP calculations. All of the targets, apart from a single potentially disorderedZ′ = 2 polymorph of the drug candidate, were predicted by at least one submission. Despite many remaining challenges, it is clear that CSP methods are becoming more applicable to a wider range of real systems, including salts, hydrates and larger flexible molecules. The results also highlight the potential for CSP calculations to complement and augment experimental studies of organic solid forms.

The development and application of computational methods for studying molecular crystals, particularly density-functional theory (DFT), is a large and ever-growing field, driven by their numerous applications. Here we expand on our recent study of the importance of many-body van der Waals interactions in molecular crystals [A. M. Reilly and A. Tkatchenko, J. Phys. Chem. Lett. 4, 1028 (2013)10.1021/jz400226x], with a larger database of 23 molecular crystals. Particular attention has been paid to the role of the vibrational contributions that are required to compare experiment sublimation enthalpies with calculated lattice energies, employing both phonon calculations and experimental heat-capacity data to provide harmonic and anharmonic estimates of the vibrational contributions. Exact exchange, which is rarely considered in DFT studies of molecular crystals, is shown to have a significant contribution to lattice energies, systematically improving agreement between theory and experiment. When the vibrational and exact-exchange contributions are coupled with a many-body approach to dispersion, DFT yields a mean absolute error (3.92 kJ/mol) within the coveted “chemical accuracy” target (4.2 kJ/mol). The role of many-body dispersion for structures has also been investigated for a subset of the database, showing good performance compared to X-ray and neutron diffraction crystal structures. The results show that the approach employed here can reach the demanding accuracy of crystal-structure prediction and organic material design with minimal empiricism.

In this perspective we discuss recent advances in the understanding of collective and many-body van der Waals interactions and their role and impact for molecular materials.

Aspirin has been used and studied for over a century but has only recently been shown to have an additional polymorphic form, known as form II. Since the two observed solid forms of aspirin are degenerate in terms of lattice energy, kinetic effects have been suggested to determine the metastability of the less abundant form II. Here, first-principles calculations provide an alternative explanation based on free-energy differences at room temperature. The explicit consideration of many-body van der Waals interactions in the free energy demonstrates that the stability of the most abundant form of aspirin is due to a subtle coupling between collective electronic fluctuations and quantized lattice vibrations. In addition, a systematic analysis of the elastic properties of the two forms of aspirin rules out mechanical instability of form II as making it metastable.