skip to main content


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 matter are discussed.  more » « less
Award ID(s):
2102568
NSF-PAR ID:
10331860
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Electronic Structure
Volume:
4
Issue:
1
ISSN:
2516-1075
Page Range / eLocation ID:
014003
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    ipso‐Arylative ring‐opening polymerization of 2‐bromo‐8‐aryl‐8H‐indeno[2,1‐b]thiophen‐8‐ol monomers proceeds to Mnup to 9 kg mol−1with conversion of the monomer diarylcarbinol groups to pendent conjugated aroylphenyl side chains (2‐benzoylphenyl or 2‐(4‐hexylbenzoyl)phenyl), which influence the optical and electronic properties of the resulting polythiophenes. Poly(3‐(2‐(4‐hexylbenzoyl)phenyl)thiophene) was found to have lower frontier orbital energy levels (HOMO/LUMO=−5.9/−4.0 eV) than poly(3‐hexylthiophene) owing to the electron‐withdrawing ability of the aryl ketone side chains. The electron mobility (ca. 2×10−3 cm2 V−1 s−1) for poly(3‐(2‐(4‐hexylbenzoyl)phenyl)thiophene) was found to be significantly higher than the hole mobility (ca. 8×10−6 cm2 V−1 s−1), which suggests such polymers are candidates for n‐type organic semiconductors. Density functional theory calculations suggest that backbone distortion resulting from side‐chain steric interactions could be a key factor influencing charge mobilities.

     
    more » « less
  2. Abstract

    ipso‐Arylative ring‐opening polymerization of 2‐bromo‐8‐aryl‐8H‐indeno[2,1‐b]thiophen‐8‐ol monomers proceeds to Mnup to 9 kg mol−1with conversion of the monomer diarylcarbinol groups to pendent conjugated aroylphenyl side chains (2‐benzoylphenyl or 2‐(4‐hexylbenzoyl)phenyl), which influence the optical and electronic properties of the resulting polythiophenes. Poly(3‐(2‐(4‐hexylbenzoyl)phenyl)thiophene) was found to have lower frontier orbital energy levels (HOMO/LUMO=−5.9/−4.0 eV) than poly(3‐hexylthiophene) owing to the electron‐withdrawing ability of the aryl ketone side chains. The electron mobility (ca. 2×10−3 cm2 V−1 s−1) for poly(3‐(2‐(4‐hexylbenzoyl)phenyl)thiophene) was found to be significantly higher than the hole mobility (ca. 8×10−6 cm2 V−1 s−1), which suggests such polymers are candidates for n‐type organic semiconductors. Density functional theory calculations suggest that backbone distortion resulting from side‐chain steric interactions could be a key factor influencing charge mobilities.

     
    more » « less
  3. 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. 
    more » « less
  4. Three implicit solvation models, the conductor‐like polarizable continuum model (C‐PCM), the conductor‐like screening model (COSMO), and universal implicit solvent model (SMD), combined with a hybrid two layer QM/QM approach (ONIOM), were utilized to calculate the pKavalues, using a direct thermodynamic scheme, of a set of Group 10 transition metal (TM) hydrides in acetonitrile. To obtain the optimal combination of quantum methods for ONIOM calculations with implicit solvation models, the influence of factors, such as the choice of density functional and basis set, the atomic radii used to build a cavity in the solvent, and the size of the model system in an ONIOM scheme, was examined. Additionally, the impact of Grimme's empirical dispersion correction and exact exchange was also investigated. The results were calibrated by experimental data. This investigation provides insight about effective models for the prediction of thermodynamic properties of TM‐containing complexes with bulky ligands. © 2019 Wiley Periodicals, Inc.

     
    more » « less
  5. We derive and investigate numerically a reduced model for wave–vortex interactions involving non-dispersive waves, which we study in a two-dimensional shallow water system with an eye towards applications in atmosphere–ocean fluid dynamics. The model consists of a coupled set of nonlinear partial differential equations for the Lagrangian-mean velocity and the wave-related pseudomomentum vector field defined in generalized Lagrangian-mean theory. It allows for two-way interactions between the waves and the balanced flow that is controlled by the distribution of Lagrangian-mean potential vorticity, and for strong solutions it features a desirable exact energy conservation law for the sum of wave energy and mean flow energy. Our model goes beyond standard ray tracing as we can derive weak solutions that contain discontinuities in the pseudomomentum field, using the theory of weakly hyperbolic systems. This allows caustics to form without predicting infinite wave amplitudes, as would be the case in the standard ray-tracing theory. Suitable wave forcing and dissipation terms are added to the model and a numerical scheme for the model is implemented as a coupled set of pseudo-spectral and finite-volume integrators. Idealized examples of interactions between wavepackets and simple vortex structures are presented to illustrate the model dynamics. The unforced and non-dissipative simulations suggest a heuristic rule of ‘greedy’ waves, i.e. in the long run the wave field always extracts energy from the mean flow. 
    more » « less