skip to main content


Title: Range-dependence of two-body intermolecular interactions and their energy components in molecular crystals
Routinely assessing the stability of molecular crystals with high accuracy remains an open challenge in the computational sciences. The many-body expansion decomposes computation of the crystal lattice energy into an embarrassingly parallel collection of computations over molecular dimers, trimers, and so forth, making quantum chemistry techniques tractable for many crystals of small organic molecules. By examining the range-dependence of different types of energetic contributions to the crystal lattice energy, we can glean qualitative understanding of solid-state intermolecular interactions as well as practical, exploitable reductions in the number of computations required for accurate energies. Here, we assess the range-dependent character of two-body interactions of 24 small organic molecular crystals by using the physically interpretable components from symmetry-adapted perturbation theory (electrostatics, exchange-repulsion, induction/polarization, and London dispersion). We also examine correlations between the convergence rates of electrostatics and London dispersion terms with molecular dipole moments and polarizabilities, to provide guidance for estimating convergence rates in other molecular crystals.  more » « less
Award ID(s):
1955940 1828187 1852511
NSF-PAR ID:
10353742
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
The Journal of Chemical Physics
Volume:
157
Issue:
8
ISSN:
0021-9606
Page Range / eLocation ID:
084503
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Heng, Jerry (Ed.)
    The morphological evolution of organic crystals during crystallization depends on the face-specific growth rates. Classical growth rate models relate the face-specific growth rates to the crystal lattice, energy of stable facets, growth mechanism, and supersaturation. The complexities of these models have increased over time to account accurately for solution conditions, the structure of growth units, and their attachment rates. Such advanced growth rate models require several layers of computations to obtain attachment energies of facets, nucleation rates, kink density, and attachment rates. Among these, the most intensive and time-consuming computation is for attachment rates, which require molecular dynamic simulations. This substantially increases the overall computation time to predict the absolute growth rate for even one crystallization condition. Since it is nearly impossible to iterate such a growth rate model, optimization schemes cannot be implemented to identify solution conditions that favor specific crystal growth. To reduce the computational time for attachment rate calculations, we implement a group contribution method (GCM) that relates the properties of functional groups in a molecule to their attachment rates to the crystal lattice, thereby rapidly estimating the growth rates of organic crystals. The process of molecular attachment involves partial desolvation of a solvated molecule, referred to as a transition state, followed by total desolvation via spontaneous attachment to a crystal facet. The first step in GCM is to identify the equilibrium states of fully solvated and partially desolvated solute molecules. The degree of supersaturation dictates the extent of this equilibrium and, thereby, the activation barrier for the growth of crystals, according to transition state theory. Identifying this equilibrium phenomenon allows for capturing the functional-group-specific interactions that depend on molecular motion, which could be related to operating conditions such as temperature and pressure. The stochastic optimization technique with Monte-Carlo sampling allows an efficient optimization problem solution to obtain the group interaction parameters. The GCM approach is first validated for the estimation of growth rates of glutamic acid and L-histidine, and then extended to predict growth rates of alanine and glycine rapidly. The optimized parameters and GCM scheme can be used to estimate growth rates in other crystallization systems. 
    more » « less
  2. An impressive collection of accurate two-body interaction energies for small complexes has been assembled into benchmark databases and used to improve the performance of multiple density functional, semiempirical, and machine learning methods. Similar benchmark data on nonadditive three-body energies in molecular trimers are comparatively scarce, and the existing ones are practically limited to homotrimers. In this work, we present a benchmark dataset of 20 equilibrium noncovalent interaction energies for a small but diverse selection of 10 heteromolecular trimers. The new 3BHET dataset presents complexes that combine different interactions including π −π, anion−π, cation−π, and various motifs of hydrogen and halogen bonding in each trimer. A detailed symmetry-adapted perturbation theory (SAPT)-based energy decomposition of the two- and three-body interaction energies shows that 3BHET consists of electrostatics- and dispersion-dominated complexes. The nonadditive three-body contribution is dominated by induction, but its influence on the overall bonding type in the complex (as exemplified by its position on the ternary diagram) is quite small. We also tested the extended SAPT (XSAPT) approach which is capable of including some nonadditive interactions in clusters of any size. The resulting three-body dispersion term (obtained from the many-body dispersion formalism) is mostly in good agreement with the supermolecular CCSD(T)−MP2 values and the nonadditive induction term is similar to the three-body SAPT(DFT) data, but the overall three-body XSAPT energies are not very accurate as they are missing the first-order exchange terms. 
    more » « less
  3. To study the contribution of three-body dispersion to crystal lattice energies, we compute the three-body contributions to the lattice energies for crystalline benzene, carbon dioxide, and triazine using various computational methods. We show that these contributions converge quickly as the intermolecular distances between the monomers grow. In particular, the smallest value among the three pairwise intermonomer closest-contact distances, Rmin, shows a strong correlation with the three-body contribution to the lattice energy, and, here, the largest of the closest-contact distances, Rmax, serves as a cutoff criterion to limit the number of trimers to be considered. We considered all trimers up to Rmax=15Å. The trimers with Rmin<4Å contribute 90.4%, 90.6%, and 93.9% of the total three-body contributions for crystalline benzene, carbon dioxide, and triazine, respectively, for the coupled-cluster singles, doubles, and perturbative triples [CCSD(T)] method. For trimers with Rmin>4Å, the second-order Møller–Plesset perturbation theory (MP2) supplemented with the Axilrod–Teller–Muto (ATM) three-body dispersion correction reproduces the CCSD(T) values for the cumulative three-body contributions with errors of less than 0.1 kJ mol−1. Moreover, three-body contributions are converged within 0.15 kJ mol−1 by Rmax=10Å. From these results, it appears that in molecular crystals where dispersion dominates the three-body contribution to the lattice energy, the trimers with Rmin>4Å can be computed with the MP2+ATM method to reduce the computational cost, and those with Rmax>10Å appear to be basically negligible. 
    more » « less
  4. null (Ed.)
    9-Anthracene carboxylic acid derivatives comprise a family of thermally reversible photomechanical molecular crystals. The photomechanical response relies on a [4 + 4] photodimerization followed by dissociation that occurs on timescales of seconds to minutes. A combined theoretical and experimental investigation is undertaken to better understand how chemical modification of the anthracene core influences energetics of both the isolated molecule and the crystal lattice. We use both density functional theory and dispersion-corrected Moller–Plesset perturbation theory computational methods to establish orbital energies, photodimerization reaction energies, and lattice energies for a set of substituted 9-anthracene carboxylic acid molecules. The calculations reveal that steric interactions play a dominant role in the ability to form photodimers and indicate an energetic threshold of 80–90 kJ per mole for the dimerization reaction. Examination of intermolecular bonding in a subset of fluorinated 9ACs revealed the absence of H⋯F intermolecular bond formation and energy differences that can explain observed trends in the dissociation kinetics and mechanical reset times. Fluorescence recovery after photobleaching experiments shows that the photodimer dissociation kinetics depend on the amount of initial photodimer, preventing a straightforward correlation between halogen atom substitution and dissociation rates using the Bell–Evans–Polanyi principle. The results clarify how molecular structure affects intermolecular interactions and photoreactivity in this family of molecular crystals, but the origin of the complex photodimer dissociation dynamics remains an open question. 
    more » « less
  5. The many-body expansion (MBE) is promising for the efficient, parallel computation of lattice energies in organic crystals. Very high accuracy should be achievable by employing coupled-cluster singles, doubles, and perturbative triples at the complete basis set limit [CCSD(T)/CBS] for the dimers, trimers, and potentially tetramers resulting from the MBE, but such a brute-force approach seems impractical for crystals of all but the smallest molecules. Here, we investigate hybrid or multi-level approaches that employ CCSD(T)/CBS only for the closest dimers and trimers and utilize much faster methods like Møller–Plesset perturbation theory (MP2) for more distant dimers and trimers. For trimers, MP2 is supplemented with the Axilrod–Teller–Muto (ATM) model of three-body dispersion. MP2(+ATM) is shown to be a very effective replacement for CCSD(T)/CBS for all but the closest dimers and trimers. A limited investigation of tetramers using CCSD(T)/CBS suggests that the four-body contribution is entirely negligible. The large set of CCSD(T)/CBS dimer and trimer data should be valuable in benchmarking approximate methods for molecular crystals and allows us to see that a literature estimate of the core-valence contribution of the closest dimers to the lattice energy using just MP2 was overbinding by 0.5 kJ mol−1, and an estimate of the three-body contribution from the closest trimers using the T0 approximation in local CCSD(T) was underbinding by 0.7 kJ mol−1. Our CCSD(T)/CBS best estimate of the 0 K lattice energy is −54.01 kJ mol−1, compared to an estimated experimental value of −55.3 ± 2.2 kJ mol−1. 
    more » « less