skip to main content


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
NSF-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. null (Ed.)
    To facilitate computational investigation of intermolecular interactions in the solution phase, we report the development of ALMO-EDA(solv), a scheme that allows the application of continuum solvent models within the framework of energy decomposition analysis (EDA) based on absolutely localized molecular orbitals (ALMOs). In this scheme, all the quantum mechanical states involved in the variational EDA procedure are computed with the presence of solvent environment so that solvation effects are incorporated in the evaluation of all its energy components. After validation on several model complexes, we employ ALMO-EDA(solv) to investigate substituent effects on two classes of complexes that are related to molecular CO 2 reduction catalysis. For [FeTPP(CO 2 -κC)] 2− (TPP = tetraphenylporphyrin), we reveal that two ortho substituents which yield most favorable CO 2 binding, –N(CH 3 ) 3 + (TMA) and –OH, stabilize the complex via through-structure and through-space mechanisms, respectively. The coulombic interaction between the positively charged TMA group and activated CO 2 is found to be largely attenuated by the polar solvent. Furthermore, we also provide computational support for the design strategy of utilizing bulky, flexible ligands to stabilize activated CO 2 via long-range Coulomb interactions, which creates biomimetic solvent-inaccessible “pockets” in that electrostatics is unscreened. For the reactant and product complexes associated with the electron transfer from the p -terphenyl radical anion to CO 2 , we demonstrate that the double terminal substitution of p -terphenyl by electron-withdrawing groups considerably strengthens the binding in the product state while moderately weakens that in the reactant state, which are both dominated by the substituent tuning of the electrostatics component. These applications illustrate that this new extension of ALMO-EDA provides a valuable means to unravel the nature of intermolecular interactions and quantify their impacts on chemical reactivity in solution. 
    more » « less
  2. Abstract

    A far‐red absorbing sensitizer, BF2‐chelated azadipyrromethane (azaBODIPY) has been employed as an electron acceptor to synthesize a series of push‐pull systems linked with different nitrogenous electron donors, viz.,N,N‐dimethylaniline (NND), triphenylamine (TPA), and phenothiazine (PTZ) via an acetylene linker. The structural integrity of the newly synthesized push‐pull systems was established by spectroscopic, electrochemical, spectroelectrochemical, and DFT computational methods. Cyclic and differential pulse voltammetry studies revealed different redox states and helped in the estimation of the energies of the charge‐separated states. Further, spectroelectrochemical studies performed in a thin‐layer optical cell revealed diagnostic peaks of azaBODIPY⋅in the visible and near‐IR regions. Free‐energy calculations revealed the charge separation from one of the covalently linked donors to the1azaBODIPY* to yield Donor⋅+‐azaBODIPY⋅to be energetically favorable in a polar solvent, benzonitrile, and the frontier orbitals generated on the optimized structures helped in assessing such a conclusion. Consequently, the steady‐state emission studies revealed quenching of the azaBODIPY fluorescence in all of the investigated push‐pull systems in benzonitrile and to a lesser extent in mildly polar dichlorobenzene, and nonpolar toluene. The femtosecond pump‐probe studies revealed the occurrence of excited charge transfer (CT) in nonpolar toluene while a complete charge separation (CS) for all three push‐pull systems in polar benzonitrile. The CT/CS products populated the low‐lying3azaBODIPY* prior to returning to the ground state. Global target (GloTarAn) analysis of the transient data revealed the lifetime of the final charge‐separated states (CSS) to be 195 ps for NND‐derived, 50 ps for TPA‐derived, and 85 ps for PTZ‐derived push‐pull systems in benzonitrile.

     
    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. Conjugated polymers composed of tricoordinate boron and π-conjugated units possess extended conjugation with relatively low-lying LUMOs arising from p B –π interactions. However, donor–acceptor (D–A) polymers that feature triorganoboranes alternating with highly electron-rich donors remain scarce. We present here a new class of hybrid D–A polymers that combine electron-rich dithienosiloles or dithienogermoles with highly robust tricoordinate borane acceptors. Polymers of modest to high molecular weight are readily prepared by Pd-catalyzed Stille coupling reaction of bis(halothienyl)boranes and distannyldithienosiloles or -germoles. The polymers are obtained as dark red solids that are stable in air and soluble in common organic solvents. Long wavelength UV-vis absorptions at ca. 500–550 nm indicate effective π-conjugation and pronounced D–A interactions along the backbone. The emission maxima occur at wavelengths longer than 600 nm in solution and experience further shifts to lower energy with increasing solvent polarity, indicative of strong intramolecular charge transfer (ICT) character of the excited state. The powerful acceptor character of the borane comonomer units in the polymer structures is also evident from cyclic voltammetry (CV) analyses that reveal relatively low-lying LUMO levels of the polymers, enhancing the D–A interaction. Density functional theory (DFT) calculations on model oligomers further support these experimental observations. 
    more » « less