Conformational polymorphs of organic molecular crystals represent a challenging test for quantum chemistry because they require careful balancing of the intra- and intermolecular interactions. This study examines 54 molecular conformations from 20 sets of conformational polymorphs, along with the relative lattice energies and 173 dimer interactions taken from six of the polymorph sets. These systems are studied with a variety of van der Waals-inclusive density functionals theory models; dispersion-corrected spin-component-scaled second-order Møller–Plesset perturbation theory (SCS-MP2D); and domain local pair natural orbital coupled cluster singles, doubles, and perturbative triples [DLPNO-CCSD(T)]. We investigate how delocalization error in conventional density functionals impacts monomer conformational energies, systematic errors in the intermolecular interactions, and the nature of error cancellation that occurs in the overall crystal. The density functionals B86bPBE-XDM, PBE-D4, PBE-MBD, PBE0-D4, and PBE0-MBD are found to exhibit sizable one-body and two-body errors vs DLPNO-CCSD(T) benchmarks, and the level of success in predicting the relative polymorph energies relies heavily on error cancellation between different types of intermolecular interactions or between intra- and intermolecular interactions. The SCS-MP2D and, to a lesser extent, ωB97M-V models exhibit smaller errors and rely less on error cancellation. Implications for crystal structure prediction of flexible compounds are discussed. Finally, the one-body and two-body DLPNO-CCSD(T) energies taken from these conformational polymorphs establish the CP1b and CP2b benchmark datasets that could be useful for testing quantum chemistry models in challenging real-world systems with complex interplay between intra- and intermolecular interactions, a number of which are significantly impacted by delocalization error.
- Award ID(s):
- 1955554
- NSF-PAR ID:
- 10335048
- Date Published:
- Journal Name:
- Physical Chemistry Chemical Physics
- Volume:
- 24
- Issue:
- 6
- ISSN:
- 1463-9076
- Page Range / eLocation ID:
- 3695 to 3712
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
-
It is demonstrated that a double hybrid density functional approximation, ωB88PTPSS, that incorporates equipartition of density functional theory and the non-local correlation, however with a meta-generalized gradient approximation correlation functional, as well as with the range-separated exchange of ωB2PLYP, provides accurate excitation energies for conventional systems, as well as correct prescription of negative singlet–triplet gaps for non-conventional systems with inverted gaps, without any necessity for parametric scaling of the same-spin and opposite-spin non-local correlation energies. Examined over “safe” excitations of the QUESTDB set, ωB88PTPSS performs quite well for open-shell systems, correctly and fairly accurately [relative to equation-of-motion coupled-cluster singles and doubles (EOM-CCSD) reference] predicts negative gaps for 50 systems with inverted singlet–triplet gaps, and is one of the leading performers for intramolecular charge-transfer excitations and achieves near-second-order approximate coupled cluster (CC2) and second-order algebraic diagrammatic construction quality for the Q1 and Q2 subsets. Subsequently, we tested ωB88PTPSS on two sets of real-life examples from recent computational chemistry literature–the low energy bands of chlorophyll a (Chl a) and a set of thermally activated delayed fluorescence (TADF) systems. For Chl a, ωB88PTPSS qualitatively and quantitatively achieves DLPNO-STEOM-CCSD-level performance and provides excellent agreement with experiment. For TADF systems, ωB88PTPSS agrees quite well with spin-component-scaled CC2 (SCS-CC2) excitation energies, as well as experimental values, for the gaps between the S1 and T1 excited states.more » « less
-
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
-
more » « less
Abstract The interactions between carbon dioxide and cluster models of coordinatively unsaturated metal–organic frameworks (MOFs) were studied using a variety of ab initio methods. Three metal species and three organic linkers in four structures were considered in these models as a representation of the tunable nature of MOFs and the potential multireference character of such systems. Common single‐reference methods, such as MP2 and CCSD(T), were compared with multireference methods based on complete active space self‐consistent field theory, going as far as multireference configuration interaction with single and double excitations (MRCISD). Special consideration is taken to avoid issues of size inconsistency in the CI results, where an alternate reference is used in the interaction energy definition. The benchmark values are used to judge the adequacy of a selection of density functionals for the current systems. Symmetry‐adapted perturbation theory (SAPT) decomposition was performed to elucidate the important effects that comprise the binding interactions. The systems proved to have very limited multireference character, and MP2 values were closer to the CCSD(T) benchmark than the more difficult MRCISD results. Though the SAPT total energies prove to be relatively poor approximations to the benchmark interaction energies, they reveal (in most cases) the correct trends with respect to the choice of the metal. The SAPT energy decompositions indicate that the
CO2 binding is primarily driven by electrostatics, but induction and dispersion also provide sizable, and quite similar, attractive contributions. Importantly, the small diformate model provides a faithful representation of complexes with large aromatic linkers, both in terms of the total interaction energy and the SAPT decomposition. -
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
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1039/D1CP04922D
