Improved Description of Intra- and Intermolecular Interactions through Dispersion-Corrected Second-Order Møller-Plesset Perturbation Theory.

ConspectusThe quantum chemical modeling of organic crystals and other molecular condensed-phase problems requires computationally affordable electronic structure methods which can simultaneously describe intramolecular conformational energies and intermolecular interactions accurately. To achieve this, we have developed a spin-component-scaled, dispersion-corrected second-order Møller-Plesset perturbation theory (SCS-MP2D) model. SCS-MP2D augments canonical MP2 with a dispersion correction which removes the uncoupled Hartree-Fock dispersion energy present in canonical MP2 and replaces it with a more reliable coupled Kohn-Sham treatment, all evaluated within the framework of Grimme's D3 dispersion model. The spin-component scaling is then used to improve the description of the residual (nondispersion) portion of the correlation energy.The SCS-MP2D model improves upon earlier corrected MP2 models in a few ways. Compared to the highly successful dispersion-corrected MP2C model, which is based solely on intermolecular perturbation theory, the SCS-MP2D dispersion correction improves the description of both inter- and intramolecular interactions. The dispersion correction can also be evaluated with trivial computational cost, and nuclear analytic gradients are computed readily to enable geometry optimizations. In contrast to earlier spin-component scaling MP2 models, the optimal spin-component scaling coefficients are only mildly sensitive to the choice of training data, and a single global parametrization of the model can describe both thermochemistry and noncovalent interactions.The resulting dispersion-corrected, spin-component-scaled MP2 (SCS-MP2D) model predicts conformational energies and intermolecular interactions with accuracy comparable to or better than those of many range-separated and double-hybrid density functionals, as is demonstrated on a variety of benchmark tests. Among the functionals considered here, only the revDSD-PBEP86-D3(BJ) functional gives consistently smaller errors in benchmark tests. The results presented also hint that further improvements of SCS-MP2D may be possible through a more robust fitting procedure for the seven empirical parameters.To demonstrate the performance of SCS-MP2D further, several applications to molecular crystal problems are presented. The three chosen examples all represent cases where density-driven delocalization error causes GGA or hybrid density functionals to artificially stabilize crystals exhibiting more extended π-conjugation. Our pragmatic strategy addresses the delocalization error by combining a periodic density functional theory (DFT) treatment of the infinite lattice with intramolecular/conformational energy corrections computed with SCS-MP2D. For the anticancer drug axitinib, applying the SCS-MP2D conformational energy correction produces crystal polymorph stabilities that are consistent with experiment, in contrast to earlier studies. For the crystal structure prediction of the ROY molecule, so named for its colorful red, orange, and yellow crystals, this approach leads to the first plausible crystal energy landscape, and it reveals that the lowest-energy polymorphs have already been found experimentally. Finally, in the context of photomechanical crystals, which transform light into mechanical work, these techniques are used to predict the structural transformations and extract design principles for maximizing the work performed.

Measuring and Modeling Highly Accurate 15 N Chemical Shift Tensors in a Peptide.

NMR studies measuring chemical shift tensors are increasingly being employed to assign structure in difficult-to-crystallize solids. For small organic molecules, such studies usually focus on 13 C sites, but proteins and peptides are more commonly described using 15 N amide sites. An important and often neglected consideration when measuring shift tensors is the evaluation of their accuracy against benchmark standards, where available. Here we measure 15 N tensors in the dipeptide glycylglycine at natural abundance using the slow-spinning FIREMAT method with SPINAL-64 decoupling. The accuracy of these 15 N tensors is evaluated by comparing to benchmark single crystal NMR 15 N measurements and found to be statistically indistinguishable. These FIREMAT experimental results are further used to evaluate the accuracy of theoretical predictions of tensors from four different density functional theory (DFT) methods that include lattice effects. The best theoretical approach provides a root mean square (rms) difference of ±3.9 ppm and is obtained from a fragment-based method and the PBE0 density functional.

Symmetry breaking in benzene and larger aromatic molecules within generalized valence bond coupled cluster methods.

The origin of symmetry breaking (SB) in benzene in generalized valence bond methods is investigated within a coupled cluster formalism that correlates all valence electrons. Retention of a limited number of pair correlation amplitudes (as in the perfect- and imperfect-pairing models) that incompletely describes interpair correlations leads to symmetry breaking as the orbitals and amplitudes are optimized. Local correlation models that are exact for one, two, and three interacting pairs at the doubles excitation level are compared against the exact pair correlation treatment, which correlates four interacting pairs at once in the connected double substitution operator. For simplicity, this comparison is performed with a second-order model of electron correlation, which is reasonably faithful to the infinite-order result. The significant SB known for the one-pair model (perfect pairing) is not eliminated at the two-pair level, but is virtually eliminated at the three-pair level. Therefore, a tractable hybrid model is proposed, which combines three-pair correlations at the second-order level and infinite-order treatment for the strong imperfect-pairing correlations involving one and two-pair correlations. This model greatly reduces SB in benzene and larger delocalized pi systems such as naphthalene and the phenalenyl cation and anion. The resulting optimized orbitals are localized in the sigma space but exhibit significant delocalization in the pi space. This means that correlation effects associated with different resonance structures are treated in a more balanced way than if the pi orbitals localize, leading to reduced SB.