skip to main content

This content will become publicly available on January 25, 2023

Title: Dispersion size-consistency
Abstract A multivariate adiabatic connection (MAC) framework for describing dispersion interactions in a system consisting of N non-overlapping monomers is presented. By constraining the density to the physical ground-state density of the supersystem, the MAC enables a rigorous separation of induction and dispersion effects. The exact dispersion energy is obtained from the zero-temperature fluctuation–dissipation theorem and partitioned into increments corresponding to the interaction energy gained when an additional monomer is added to a K -monomer system. The total dispersion energy of an N -monomer system is independent of any partitioning into subsystems. This statement of dispersion size consistency is shown to be an exact constraint. The resulting additive separability of the dispersion energy results from multiplicative separability of the generalized screening factor defined as the inverse generalized dielectric function. Many-body perturbation theory (MBPT) is found to violate dispersion size-consistency because perturbative approximations to the generalized screening factor are nonseparable; on the other hand, random phase approximation-type methods produce separable generalized screening factors and therefore preserve dispersion size-consistency. This result further explains the previously observed increase in relative errors of MBPT for dispersion interactions as the system size increases. Implications for electronic structure theory and applications to supramolecular materials and condensed more » matter are discussed. « less
; ; ;
Award ID(s):
Publication Date:
Journal Name:
Electronic Structure
Page Range or eLocation-ID:
Sponsoring Org:
National Science Foundation
More Like this
  1. Molecular crystal engineering seeks to tune the material properties by controlling the crystal packing. However, the range of achievable properties is constrained by the limited energy range of polymorphs which can be crystallized. Here, computational modeling highlights that a solid-state crystal-to-crystal chemical reaction in 9- tert -butyl anthracene ester (9TBAE) nanorods [Al-Kaysi et al. , J. Am. Chem. Soc. , 2006, 128 , 15938] imparts “synthetic memory” into the crystal structure that allows reproducible formation of a highly metastable, yet long-lived polymorph. Specifically, whereas the vast majority of known polymorphs exhibit lattice energy differences below 10 kJ mol −1 , the conformational polymorph formed via solid state reaction chemistry lies 14 kJ mol −1 higher in energy than the form grown from solution, according to calculations that combine a dispersion-corrected second-order Møller–Plesset perturbation theory (MP2D) treatment of the monomer and photodimer with a density functional theory treatment (B86bPBE-XDM) of the intermolecular interactions in the crystal. Moreover, the solid-state reaction environment traps a highly unstable intramolecular photodimer conformation which defies the conventional wisdom surrounding conformational polymorphs. These observations suggest that solid-state reaction chemistry represents an under-appreciated strategy for producing polymorphs that would likely be unobtainable otherwise.
  2. 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-bodymore »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.

    « less
  3. The defining characteristic of hole-doped cuprates is d -wave high temperature superconductivity. However, intense theoretical interest is now focused on whether a pair density wave state (PDW) could coexist with cuprate superconductivity [D. F. Agterberg et al., Annu. Rev. Condens. Matter Phys. 11, 231 (2020)]. Here, we use a strong-coupling mean-field theory of cuprates, to model the atomic-scale electronic structure of an eight-unit-cell periodic, d -symmetry form factor, pair density wave (PDW) state coexisting with d -wave superconductivity (DSC). From this PDW + DSC model, the atomically resolved density of Bogoliubov quasiparticle states N r , E is predicted at the terminal BiO surface of Bi 2 Sr 2 CaCu 2 O 8 and compared with high-precision electronic visualization experiments using spectroscopic imaging scanning tunneling microscopy (STM). The PDW + DSC model predictions include the intraunit-cell structure and periodic modulations of N r , E , the modulations of the coherence peak energy Δ p r , and the characteristics of Bogoliubov quasiparticle interference in scattering-wavevector space q - space . Consistency between all these predictions and the corresponding experiments indicates that lightly hole-doped Bi 2 Sr 2 CaCu 2 O 8 does contain a PDW + DSC state. Moreover,more »in the model the PDW + DSC state becomes unstable to a pure DSC state at a critical hole density p *, with empirically equivalent phenomena occurring in the experiments. All these results are consistent with a picture in which the cuprate translational symmetry-breaking state is a PDW, the observed charge modulations are its consequence, the antinodal pseudogap is that of the PDW state, and the cuprate critical point at p * ≈ 19% occurs due to disappearance of this PDW.« less
  4. We study statistical inference and distributionally robust solution methods for stochastic optimization problems, focusing on confidence intervals for optimal values and solutions that achieve exact coverage asymptotically. We develop a generalized empirical likelihood framework—based on distributional uncertainty sets constructed from nonparametric f-divergence balls—for Hadamard differentiable functionals, and in particular, stochastic optimization problems. As consequences of this theory, we provide a principled method for choosing the size of distributional uncertainty regions to provide one- and two-sided confidence intervals that achieve exact coverage. We also give an asymptotic expansion for our distributionally robust formulation, showing how robustification regularizes problems by their variance. Finally, we show that optimizers of the distributionally robust formulations we study enjoy (essentially) the same consistency properties as those in classical sample average approximations. Our general approach applies to quickly mixing stationary sequences, including geometrically ergodic Harris recurrent Markov chains.
  5. From local mode stretching force constants and topological electron density analysis, computed at either the UM06/6-311G(d,p), UM06/SDD, or UM05-2X/6– 31++G(d,p) level of theory, we elucidate on the nature/strength of the parallel π- stacking interactions (i.e. pancake bonding) of the 1,2-dithia-3,5-diazolyl dimer, 1,2-diselena-3,5-diazolyl dimer, 1,2-tellura-3,5-diazolyl dimer, phenalenyl dimer, 2,5,8-tri-methylphenalenyl dimer, and the 2,5,8-tri-t-butylphenalenyl dimer. We use local mode stretching force constants to derive an aromaticity delocalization index (AI) for the phenalenyl-based dimers and their monomers as to determine the effect of substitution and dimerization on aromaticity, as well as determining what bond property governs alterations in aromaticity. Our results reveal the strength of the C⋯C contacts and of the rings of the di-chalcodiazoyl dimers investigated decrease in parallel with decreasing chalcogen⋯chalcogen bond strength. Energy density values Hb suggest the S⋯S and Se⋯Se pancake bonds of 1,2-dithia-3,5- diazolyl dimer and the 1,2-diselena-3,5-diazolyl dimer are covalent in nature. We observe the pancake bonds, of all phenalenyl-based dimers investigated, to be electrostatic in nature. In contrast to their monomer counterparts, phenalenyl- based dimers increase in aromaticity primarily due to CC bond strengthening. For phenalenyl-based dimers we observed that the addition of bulky substituents steadily decreased the system aromaticity predominately due to CC bond weakening.