We show that electrodynamic dipolar interactions, responsible for long-range fluctuations in matter, play a significant role in the stability of molecular crystals. Density functional theory calculations with van der Waals interactions determined from a semilocal "atom-in-a-molecule" model result in a large overestimation of the dielectric constants and sublimation enthalpies for polyacene crystals from naphthalene to pentacene, whereas an accurate treatment of nonlocal electrodynamic response leads to an agreement with the measured values for both quantities. Our findings suggest that collective response effects play a substantial role not only for optical excitations, but also for cohesive properties of noncovalently bound molecular crystals. Polyacene molecular crystals form a fundamental class of aromatic solids, and have been extensively studied as potential materials for organic electronics. [1][2][3] It is understood that the optical properties of polyacenes are very sensitive to longrange intra-and intermolecular electrodynamic interactions. This is reflected by shifts in the optical absorption frequencies upon increasing the molecule size or upon solid formation, 4 and is further exhibited by the visible color of oligoacene crystals, which changes from transparent in naphthalene and anthracene, to bright orange in tetracene, and deep blue in pentacene. 4,5 The optical absorption spectrum is directly related to the polarizability through the Kramers-Kronig transformation.6 Therefore, the observed changes in the optical spectrum upon crystallization of polyacenes are accompanied by a change in the molecular polarizability. In addition, these changes in polarization should directly impact the crystal lattice energy. However, the effect of electrodynamic intermolecular interactions on the cohesive properties of molecular crystals remains poorly understood. In this Rapid Communication, we show that the dipolar electrodynamic coupling between polyacene molecules reduces the solid dielectric constant by 15%, and has an impact of up to 0.5 eV per molecule on the computed van der Waals (vdW) energies and sublimation enthalpies of these molecular crystals. Our results imply that electrodynamic response is crucial for describing both the cohesive energy and the optical properties of molecular crystals, also providing strong quantitative support to empirical relations between stability and refractive index of molecular crystals.
7Polyacene crystals are extended aromatic networks characterized by polarizable π clouds. Therefore, an appreciable part of the crystal lattice energy stems from ubiquitous attractive vdW dispersion interactions. When studying the cohesion of molecular systems, for example, using densityfunctional theory (DFT) 8,9 or classical potentials, 10 the vdW energy is typically computed using effective polarizabilities for hybridized "atoms" inside a molecule. It is common to approximate the frequency-dependent polarizability of every atom using a single effective excitation frequency (also called t...