skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: A non-perturbative pairwise-additive analysis of charge transfer contributions to intermolecular interaction energies
Energy decomposition analysis (EDA) based on absolutely localized molecular orbitals (ALMOs) decomposes the interaction energy between molecules into physically interpretable components like geometry distortion, frozen interactions, polarization, and charge transfer (CT, also sometimes called charge delocalization) interactions. In this work, a numerically exact scheme to decompose the CT interaction energy into pairwise additive terms is introduced for the ALMO-EDA using density functional theory. Unlike perturbative pairwise charge-decomposition analysis, the new approach does not break down for strongly interacting systems, or show significant exchange–correlation functional dependence in the decomposed energy components. Both the energy lowering and the charge flow associated with CT can be decomposed. Complementary occupied–virtual orbital pairs (COVPs) that capture the dominant donor and acceptor CT orbitals are obtained for the new decomposition. It is applied to systems with different types of interactions including DNA base-pairs, borane-ammonia adducts, and transition metal hexacarbonyls. While consistent with most existing understanding of the nature of CT in these systems, the results also reveal some new insights into the origin of trends in donor–acceptor interactions.  more » « less
Award ID(s):
1955643
PAR ID:
10231407
Author(s) / Creator(s):
; ; ; ;
Date Published:
Journal Name:
Physical Chemistry Chemical Physics
Volume:
23
Issue:
2
ISSN:
1463-9076
Page Range / eLocation ID:
928 to 943
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Energy decomposition analysis (EDA) is a useful method to unravel an intermolecular interaction energy into chemically meaningful components such as geometric distortion, frozen interactions, polarization, and charge transfer. A further decomposition of the polarization (POL) and charge transfer (CT) energy into fragment-wise contributions would be useful to understand the significance of each fragment during these two processes. To complement the existing exact pairwise decomposition of the CT term, this work describes formulation and implementation of a non-perturbative polarization analysis that decomposes the POL energy into an exactly fragment-wise additive sum based on the absolutely localized molecular orbital energy decomposition analysis (ALMO-EDA). These fragment-wise contributions can be further decomposed into chemically intuitive molecular orbital pairs using complementary occupied-virtual pairs (COVP) analysis. A very useful phase convention is established for each COVP such that constructive interference of occupied and virtual corresponds to electron flow into that region, whilst destructive interference corresponds to electron outflow. A range of model problems are used to demonstrate that the polarization process is typically a collective behavior of the electrons that is quite different from the charge transfer process. This provides another reason in addition to their different distance-dependence on fragment separation for separating these two processes in EDA. 
    more » « less
  2. In this article, we introduce the occupied-virtual orbitals for chemical valence (OVOCV). The OVOCVs can replace or complement the closely related idea of the natural orbitals for chemical valence (NOCV). The input is a difference density matrix connecting any initial single determinant to any final determinant, at a given molecular geometry, and a given one-particle basis. This arises in problems such as orbital rearrangement or charge-transfer in energy decomposition analysis. The OVOCVs block-diagonalize the density difference operator into 2 × 2 blocks which are spanned by one level that is filled in the initial state (the occupied OVOCV) and one which is empty (the virtual OVOCV). By contrast, the NOCVs fully diagonalize the density difference matrix, and therefore are orbitals with mixed occupied-virtual character. Use of the OVOCVs makes it much easier to identify the donor and acceptor orbitals. We also introduce two different types of energy decomposition analysis (EDA) methods with the OVOCVs, and most importantly, a charge decomposition analysis (CDA) method that fixes the unreasonably large charge transfer amount obtained directly from NOCV analysis. 
    more » « less
  3. As known, small HCl–water nanoclusters display a particular dissociation behaviour, whereby at least four water molecules are required for the ionic dissociation of HCl. In this work, we examine how intermolecular interactions promote the ionic dissociation of such nanoclusters. To this end, a set of 45 HCl–water nanoclusters with up to four water molecules is introduced. Energy decomposition analysis based on absolutely localized molecular orbitals (ALMO-EDA) is employed in order to study the importance of frozen interaction, dispersion, polarization, and charge-transfer for the dissociation. The vertical ALMO-EDA scheme is applied to HCl–water clusters along a proton-transfer coordinate varying the amount of spectator water molecules. The corresponding ALMO-EDA results show a clear preference for the dissociated cluster only in the presence of four water molecules. Our analysis of adiabatic ALMO-EDA results reveals a push–pull mechanism for the destabilization of the HCl bond based on the synergy between forward and backward charge-transfer. 
    more » « less
  4. null (Ed.)
    Quantum chemistry in the form of density functional theory (DFT) calculations is a powerful numerical experiment for predicting intermolecular interaction energies. However, no chemical insight is gained in this way beyond predictions of observables. Energy decomposition analysis (EDA) can quantitatively bridge this gap by providing values for the chemical drivers of the interactions, such as permanent electrostatics, Pauli repulsion, dispersion, and charge transfer. These energetic contributions are identified by performing DFT calculations with constraints that disable components of the interaction. This review describes the second-generation version of the absolutely localized molecular orbital EDA (ALMO-EDA-II). The effects of different physical contributions on changes in observables such as structure and vibrational frequencies upon complex formation are characterized via the adiabatic EDA. Example applications include red- versus blue-shifting hydrogen bonds; the bonding and frequency shifts of CO, N 2 , and BF bound to a [Ru(II)(NH 3 ) 5 ] 2 + moiety; and the nature of the strongly bound complexes between pyridine and the benzene and naphthalene radical cations. Additionally, the use of ALMO-EDA-II to benchmark and guide the development of advanced force fields for molecular simulation is illustrated with the recent, very promising, MB-UCB potential. 
    more » « less
  5. p-Donor/Acceptor charge-transfer (CT) interactions between redox-complementary p-systems often give rise to non-native optical and electronic properties that are beneficial for modern electronics and energy technologies. However, the formation of extended supramolecular p-donor/acceptor stacks capable of long-range charge transport requires ingenious design strategies that can help reinforce otherwise weak p-donor/acceptor noncovalent interactions. Herein, we demonstrate that a large tetragonal prismatic metal–organic cage (MOC28+) having two parallel p-donor tetrakis(4- carboxyphenyl)-Zn-porphyrin (ZnTCPP) faces located B14 Å apart can accommodate up to three redox-complementary planar aromatic guests (either three p-acceptor guests or two p-acceptors surrounding one p-donor guest) between the ZnTCPP faces, forming extended p-donor/acceptor stacks. While empty MOC28+ behaves as an insulator due to the lack of charge delocalization across its large cavity, its inclusion complexes saturated with p-acidic hexaazatriphenylene hexacarbonitrile (HATHCN) and hexacyanotriphenylene (HCTP) displayed noticeably higher electrical conductivity (8.7   10 6 and 1.3   10 6 S m 1, respectively) owing to more facile charge transport through the p-donor/ acceptor stacks composed of the p-acidic guests intercalated between the ZnTCPP faces. Thus, this work demonstrates that tetragonal prismatic metallacages with two parallel electroactive faces can facilitate the creation of extended p-donor/acceptor stacks by encapsulating redox-complementary planar guests, which in turn facilitates through-space charge delocalization, generating non-native electrical conductivity. 
    more » « less