Energy decomposition analysis based on broken symmetry unrestricted density functional theory

Theoretical Study on the Correlation between Open-Shell Electronic Structures and Third-Order Nonlinear Optical Properties in One-Dimensional Chains of π-Radicals

This paper theoretically investigated the correlation between the open-shell electronic structure and third-order nonlinear optical (NLO) properties of one-dimensional (1D) stacked chains of π-radicals. By employing the finite N-mer models consisting of methyl or phenalenyl radicals with different stacking distances, we evaluated the average and standard deviation of diradical characters yi for N-mer models of π-radicals (yav and ySD). Then, we estimated these diradical characters at the limit of N → ∞. These y-based indices were helpful in discussing the correlation between the open-shell electronic structures and the second hyperpolarizability per dimer at the limit N → ∞, γ∞ for the 1D chains with stacking distance alternation (SDA). The calculated γ∞ values and the polymer/dimer ratio γ∞/γ(N = 2) were enhanced significantly when both the stacking distance and SDA are small. We also found that the spin-unrestricted long-range-corrected (LC-)UBLYP method with the range-separating parameter μ = 0.47 bohr-1 well reproduced the trend of γ∞ of this type of 1D chain estimated at the spin-unrestricted coupled-cluster levels. The present study is expected to contribute to establishing the design guidelines for future high-performance open-shell molecular NLO materials.

Solvent Effect on Dative and Ionic Bond Strengths: A Unified Theory from Potential Analysis and Valence-Bond Computations

Solvent molecules interact with a solute through various intermolecular forces. Here we employed a potential energy surface (PES) analysis to interpret the solvent-induced variations in the strengths of dative (Me3NBH3) and ionic (LiCl) bonds, which possess both ionic and covalent (neutral) characteristics. The change of a bond is driven by the gradient (force) of the solvent-solute interaction energy with respect to the focused bond length. Positive force shortens the bond length and increases the bond force constant, leading to a blue-shift of the bond stretching vibrational frequency upon solvation. Conversely, negative force elongates the bond, resulting in a reduced bond force constant and red-shift of the stretching vibrational frequency. The different responses of Me3NBH3 and LiCl to solvation are studied with valence bond (VB) theory, as Me3NBH3 and LiCl are dominated by the neutral covalent VB structure and the ionic VB structure, respectively. The dipole moment of an ionic VB structure increases along the increasing bond distance, while the dipole moment of a neutral covalent VB structure increases with the decreasing bond distance. The roles of the dominating VB structures are further examined by the geometry optimizations and frequency calculations with the block-localized wavefunction (BLW) method.

A many-body energy decomposition analysis (MB-EDA) scheme based on a target state optimization self-consistent field (TSO-SCF) method

In this paper, we combine an energy decomposition analysis (EDA) scheme with many-body expansion (MBE) to develop a MB-EDA method to study the cooperative and anti-cooperative effects in molecular cluster systems. Based on the target state optimization self-consistent field (TSO-SCF) method, the intermolecular interaction energy can be decomposed into five chemically meaningful terms, i.e., electrostatic, exchange, polarization, charge transfer and dispersion interaction energies. MB-EDA can decompose each of these terms in MBE. This MB-EDA has been applied to 3 different cluster systems: water clusters, ionic liquid clusters, and acetonitrile-methane clusters. This reveals that electrostatic, exchange, and dispersion interactions are highly pairwise additive in all systems. In water and ionic liquid clusters, the many-body effects are significant in both polarization and charge transfer interactions, but are cooperative and anti-cooperative, respectively. For acetonitrile-methane clusters, which do not involve hydrogen bonds or charge-charge Coulombic interactions, the many-body effects are quite small. The chemical origins of different many-body effects are deeply analyzed. The MB-EDA method has been implemented in Qbics (https://qbics.info) and can be a useful tool for understanding the many-body behavior in molecular aggregates at the quantum chemical level of theory.

A theoretical study on π-stacking and ferromagnetism of the perylene diimide radical anion dimer and tetramer

Ferromagnetism is rare in pure organic materials. Recently, the perylene diimide radical anion (PDI-) salt prepared through solvothermal reduction by hydrazine hydrate has shown room-temperature ferromagnetism in our work [Jiang et al., Adv. Mater., 2022, 34, 2108103]. Based on this, herein we conduct a theoretical study based on density functional theory (DFT) to reveal the stacked geometries between two NH4PDI monomers for low-spin (LS) and high-spin (HS) states and their magnetic exchange interactions (JAB) using Yamaguchi's approximate spin projection. It is observed that the pancake-bonded dimer of NH4PDI is the most stable pimer compared to others on both LS and HS potential energy surfaces. A transition of magnetic properties from strong antiferromagnetic (-1333.9 cm-1) to moderate ferromagnetic (67.0 cm-1) appears after increasing the interplanar distance between monomers and their relative rotation angle to access the HS state. According to energy decomposition analysis, the enhanced hydrogen bond formation and decrease of Pauli repulsion is able to counteract the decrease of attraction induced by electron correlation after accessing the HS state. Stacking patterns of exchange-coupled chain consisting of the NH4PDI tetramer are obtained for the HS state after geometry optimization of the structure constructed by two most stable HS pimers. The exchange interactions (51.8 cm-1, 381.2 cm-1 and 53.2 cm-1) between adjacent NH4PDI monomers are ferromagnetic in the HS state, which is in accordance with the experimentally observed room-temperature ferromagnetism.

Theoretical analysis of expanded porphyrins: Aromaticity, stability, and optoelectronic properties

Expanded porphyrin systems are capable of binding a variety of substrates due to their increased cavity size and aromatic nature, holding important applications as magnetic resonance imaging contrast agents and as sensitizers for photodynamic therapy. It is there of fundamental interest to know the photoelectrical properties of expanded porphyrins using quantum chemistry calculations. In this work, we theoretically designed and screened a series of expanded porphyrins by incorporating terthiophene (TTH) and dithienothiophene (DTT) moieties. Our calculations showed that all the designed molecules exhibit excellent optoelectronic performance than the reference molecule. It is suggested that the porphyrin molecule with TTH moiety has better stability than the one with DTT moiety. Finally, we demonstrated that molecule 2 features with TTH moiety and the inverted selenophene ring outperform other molecules because it exhibits increased HOMO-LUMO gap, planar geometry, and strengthened aromaticity. We expect that this work can provide theoretical guidelines for the design of novel porphyrin materials.

Energy Decomposition Analysis of the Nature of Coordination Bonding at the Heme Iron Center in Cytochrome P450 Inhibition

Drug compounds or their metabolic intermediates (MIs) sometimes inhibit the function of cytochrome P450 enzymes (P450s) by forming a coordination bond with the Fe(III) heme or Fe(II) heme of P450s. Such inhibition is one of the major causes of drug-drug interactions (DDIs), a subject of longstanding academic and practical interest. However, such coordination bonding is not fully understood at the quantum mechanical level, thus hampering rational improvement of the accuracy of DDI-related predictions. In this work, we employed density functional theory (DFT) and the generalized Kohn-Sham energy decomposition analysis (GKS-EDA) scheme to investigate the nature of the coordination bonding formed in the reversible and quasi-irreversible inhibition of P450s. The GKS-EDA results highlighted a previously unrecognized role of the electron correlation effect in P450 inhibition. The correlation effect tends to be larger in Fe(II) complexes of MI-type inhibitors and is particularly prominent for the nitrosoalkane ligand. An additional natural bond orbital (NBO) analysis provided insight into the relative significance of the σ donation and π backdonation effects in various heme-inhibitor complexes.

XEDA, a fast and multipurpose energy decomposition analysis program

A fast and multipurpose energy decomposition analysis (EDA) program, called XEDA, is introduced for quantitative analysis of intermolecular interactions. This program contains a series of variational EDA methods, including LMO-EDA, GKS-EDA and their extensions, to analyze non-covalent interactions and strong chemical bonds in various environments. XEDA is highly efficient with a similar computational scaling of single point energy calculations. Its efficiency and universality are validated by a series of test examples including van der Waals interactions, hydrogen bonds, radical-radical interactions and strong covalent bonds.

Stability, Aromaticity, and Photophysical Behaviors of Macrocyclic Molecules: A Theoretical Analysis

The macrocyclic molecules with terthiophene (TTH) isomers unit exhibit intriguing properties in terms of aromaticity, stability, and absorption. In this work, we theoretically designed a series of macrocyclic molecules featured with TTH and dithienothiophene (DTT) π-conjugated building units, which are used to permute pyrrole unit in porphyrin skeleton. Density functional theory and time-dependent DFT methods are used to evaluate the performance of the designed molecules. Our simulations show that molecules 1–3 exhibit excellent optoelectronic performance. Specifically, the molecule with the DTT unit is more stable than the one with TTH unit in terms of aromaticity and aromatic stabilization energy. This is because DTT unit enhances the coplanarity of the molecular material, facilitating electronic communication. Calculation of vertical electronic excitations suggests the absorption feature of these molecules is mainly contributed by the electronic excitations of higher occupied molecular orbital (HOMO) → lowest unoccupied molecular orbital (LUMO)+1 and HOMO-1 → LUMO. Judging from the key parameters determining the overall performance, 3 stands out because of its good planarity, large HOMO–LUMO gap, and strong aromaticity among all molecules. Interestingly, molecule 1 has the current density flow distributes around the outer section of TTH unit; in contrast, molecule 3 with DTT unit has the current density flow located at the inner section of DTT, which is beneficial for stability and aromaticity. Second-order perturbation energies are calculated to rationalize this observation. We expect that these research results can provide valuable insights into the rational design of novel molecular materials for a variety of applications.

Probing radical-molecule interactions with a second generation energy decomposition analysis of DFT calculations using absolutely localized molecular orbitals

Intermolecular interactions between radicals and closed-shell molecules are ubiquitous in chemical processes, ranging from the benchtop to the atmosphere and extraterrestrial space. While energy decomposition analysis (EDA) schemes for closed-shell molecules can be generalized for studying radical-molecule interactions, they face challenges arising from the unique characteristics of the electronic structure of open-shell species. In this work, we introduce additional steps that are necessary for the proper treatment of radical-molecule interactions to our previously developed unrestricted Absolutely Localized Molecular Orbital (uALMO)-EDA based on density functional theory calculations. A "polarize-then-depolarize" (PtD) scheme is used to remove arbitrariness in the definition of the frozen wavefunction, rendering the ALMO-EDA results independent of the orientation of the unpaired electron obtained from isolated fragment calculations. The contribution of radical rehybridization to polarization energies is evaluated. It is also valuable to monitor the wavefunction stability of each intermediate state, as well as their associated spin density profiles, to ensure the EDA results correspond to a desired electronic state. These radical extensions are incorporated into the "vertical" and "adiabatic" variants of uALMO-EDA for studies of energy changes and property shifts upon complexation. The EDA is validated on two model complexes, H₂O˙F and FH˙OH. It is then applied to several chemically interesting radical-molecule complexes, including the sandwiched and T-shaped benzene dimer radical cation, complexes of pyridine with benzene and naphthalene radical cations, binary and ternary complexes of the hydroxyl radical with water (˙OH(H₂O) and ˙OH(H₂O)₂), and the pre-reactive complexes and transition states in the ˙OH + HCHO and ˙OH + CH₃CHO reactions. These examples suggest that this second generation uALMO-EDA is a useful tool for furthering one's understanding of both energetic and property changes associated with radical-molecule interactions.

Fragment Exchange Potential for Realizing Pauli Deformation of Interfragment Interactions

In fragment-based methods, the lack of explicit short-range exchange interactions between monomers can result in unphysical deformation in charge density. In this study, we describe a fragment exchange potential (XFP) to explicitly account for interfragmental Pauli deformation. In our implementation, a Kohn-Sham exchange potential is adopted along with the Yukawa potential. The method has been validated by comparison of the computed exchange energies using the XFP potential with results obtained from antisymmetrized fragmental orbitals on the S66×8 data set containing 528 bimolecular interactions of equilibrium and arbitrary geometries. It was also found that it is only necessary to deploy numerical grids on atoms within their van der Waals contacts, significantly reducing the small, albeit extra, computational cost. We anticipate that the XFP presented here may be applied to molecular dynamics simulations of macromolecules using a fragment-based quantum mechanical potential with improved SCF convergence and computational accuracy.