Local Electron Attachment Energy and Its Use for Predicting Nucleophilic Reactions and Halogen Bonding

Perylene Bisimide-Functionalized Triphenylmethyl Radicals Showing High Stability and Reversible Electrochemical Redox Properties

This study presents a series of triphenylmethyl monoradicals incorporating varying numbers of peripheral perylene bisimide (PBI) substituents (1PBI-TTM⋅, 2PBI-TTM⋅ and 3PBI-TTM⋅). The incorporation of electron-withdrawing PBI substituents significantly enhances the stability of these carbon radicals, enabling them to display exceptional electrochemical redox reversibility. Notably, the electronic interplay between the PBI substituents and the central triphenylmethyl core facilitates unique and reversible multi-step redox reactions. Among the reported radicals, the tris-PBI-functionalized radical (3PBI-TTM⋅) demonstrates the remarkable ability to accommodate up to seven electrons under negative potentials, forming high valence anions. This research promotes the development of highly stable carbon radicals with superior electrochemical oxidation-reduction processes, presenting promising avenues for the advancement of electric energy storage technologies.

Intermolecular interaction potential maps from energy decomposition for interpreting reactivity and intermolecular interactions

The electrostatic potential (ESP) has been widely used to visualize electrostatic interactions about a molecule. However, electrostatic effects are often insufficient for capturing the entirety of an interaction or a reaction of interest. In this investigation, intermolecular interaction potential maps (IMIPs), constructed from the potentials derived from energy decomposition analysis (EDA) using density functional theory, were developed and applied to provide unique insight into molecular interactions and reactivity. To this end, rather than constructing a potential map from probe point charge interactions, IMIPs were constructed from probe interactions with small molecular fragments, including CH3+, CH3-, benzene, and atomic probes including alkali metals, transition metals, and halides. The interaction potentials are further decomposed producing IMIPs for each interaction component using EDA (electrostatic, orbital, steric, etc.). The IMIPs are applied to the study of various interactions including cation-π and anion-π interactions, electrophilic and nucleophilic aromatic substitution, Lewis acid activation, π-stacking, endohedral fullerenes, and select organometallics which reveal fundamental insight into the positional preferences and physical origins of the interactions that otherwise would be difficult to uncover through other surface analyses.

Lignin-Based Visible Light-Triggered Nitric Oxide Nanogenerator for Antibacterial Applications

Nitric oxide (NO) has received growing attention as an effective antibacterial agent with broad-spectrum activity and a low risk of resistance. However, it remains challenging to develop effective, controllable, and biocompatible NO-releasing materials. Here, we report a novel NO nanogenerator (AL-BNN6-PEG) self-assembled by lignin, a UV-absorbing and hydrophobic NO donor (N,N'-disec-butyl-N,N'-dinitroso-1,4-phenylenediamine, BNN6), and PEG-DSPE2000. It was discovered that upon visible light irradiation (450-460 nm), BNN6 can be decomposed by lignin within micellar nanoparticles via a photoinduced electron transfer mechanism in the aqueous medium. Lignin not only served as a sustainable carrier, enhancing the water dispersity of BNN6, but also acted as a biocompatible photosensitizer, triggering BNN6 decomposition with the concomitant release of NO. As a result, the micellar nanoparticles displayed superior antibacterial effects against Gram-negative and Gram-positive bacteria upon visible light illumination. Moreover, MTT assay revealed the negligible cytotoxic effect of the micellar nanoparticles to the mouse fibroblast cells (L929). This research provides more insight into the BNN6 decomposition mechanism and demonstrates a straightforward, effective, and biocompatible strategy for controlled NO-mediated antibacterial applications.

A comprehensive electron wavefunction analysis toolbox for chemists, Multiwfn

Analysis of electron wavefunction is a key component of quantum chemistry investigations and is indispensable for the practical research of many chemical problems. After more than ten years of active development, the wavefunction analysis program Multiwfn has accumulated very rich functions, and its application scope has covered numerous aspects of theoretical chemical research, including charge distribution, chemical bond, electron localization and delocalization, aromaticity, intramolecular and intermolecular interactions, electronic excitation, and response property. This article systematically introduces the features and functions of the latest version of Multiwfn and provides many representative examples. Through this article, readers will be able to fully understand the characteristics and recognize the unique value of Multiwfn. The source code and precompiled executable files of Multiwfn, as well as the manual containing a detailed introduction to theoretical backgrounds and very rich tutorials, can all be downloaded for free from the Multiwfn website (http://sobereva.com/multiwfn).

Fully Bio-Based 2,5-Furandicarboxylic Acid Polyester toward Plastics with Mechanically Robust, Excellent Gas Barrier and Fast Degradation

Aliphatic-aromatic copolyesters offer a promising solution to mitigate plastic pollution, but high content of aliphatic units (>40 %) often suffer from diminished comprehensive performances. Poly(butylene oxalate-co-furandicarboxylate) (PBOF) copolyesters were synthesized by precisely controlling the oxalic acid content from 10 % to 60 %. Compared with commercial PBAT, the barrier properties of PBOF for H2O and O2 increased by more than 6 and 26 times, respectively. The introduction of the oxalic acid units allowed the water contact angle to be reduced from 82.5° to 62.9°. Superior hydrophilicity gave PBOF an excellent degradation performance within a 35-day hydrolysis. Interestingly, PBO20F and PBO30F also displayed obvious decrease of molecular weight during hydrolysis, with elastic modulus >1 GPa and tensile strength between 35-54 MPa. PBOF achieved the highest hydrolysis rates among the reported PBF-based copolyesters. The hydrolytic mechanism was further explored based on Fukui function analysis and density functional theory (DFT) calculation. Noncovalent analysis indicated that the water molecules formed hydrogen bonding interaction with adjacent ester groups and thus improved the reactivity of carbonyl carbon. PBOF not only meet the requirements of the high-performance packaging market but can quickly degrade after the end of their usage cycles, providing a new choice for green and environmental protection.

Theoretical Analysis of Coordination Geometries in Transition Metal-Histidine Complexes Using Quantum Chemical Calculations

A theoretical investigation utilizing density functional theory (DFT) calculations was conducted to explore the coordination complexes formed between histidine (His) ligands and various divalent transition metal ions (Mn2+, Fe2+, Co2+, Ni2+, Cu2+, and Zn2+). Conformational exploration of the His ligand was initially performed to assess its stability upon coordination. Both 1:1 and 1:2 of metal-to-ligand complexes were scrutinized to elucidate their structural features and the relative stability of the complexes. This study examined the ability of His to act as a bidentate or tridentate coordinating ligand, along with the differences in coordination geometry when solvent effects were incorporated. The reduced density gradient (RDG) analysis and local electron attachment energy (LEAE) analysis were employed to elucidate the interaction planes and the nucleophilic and electrophilic properties. The electronic properties were analyzed through electrostatic potential (ESP) maps and natural population analysis (NPA) of atomic charge distributions. This computational study provides valuable insights into the diverse coordination modes of His and its interactions with divalent transition metal ions, contributing to a better understanding of the role of this amino acid ligand in the formation of transition metal complexes. The findings can aid in the design and construction of self-assembled structures involving His-metal coordination.

Reaction rate constant: a theoretical description from local temperature

Application of various descriptors based on electron density and its associated quantities to quantify chemical reactivity within the conceptual density functional theory has recently come into spotlight. Among others and particularly relevant to our study, local temperature based on electron density as well as kinetic energy density, as a measure of the kinetic energy of an electron moving in the Kohn-Sham potential of systems, should be mentioned. In this work, we propose to use the local temperature for describing the reaction rate constant, where our main idea originates from the point that the smaller the local temperature at the reaction center, the easier the electron removal, leading to a larger rate constant. On the basis of theoretical considerations, it is proved that the rate constant variations caused by the substituent effects can well be proportional to the local temperature at the reaction center. In order to numerically validate our proposed approach, we have taken the phenol derivatives with the available experimental rate constants of their O-methylation reaction as working models. The reason for this choice is that one of the most versatile approaches for labeling biologically active compounds with the 11C nuclide for positron emission tomography (PET) is methylation by methyl iodide including 11C nuclide, [11C]MeI, where methylation of phenolic oxygen with [11C]MeI is utilized to label some important tracers for PET studies. Our results unveil that the local temperature changes at the reaction center of the aforementioned reaction are reasonably correlated with the rate constant variations. Hopefully, incorporating the proposed correlations between the local temperature and the kinetics data into a computer control algorithm not only provides a simple tool for predicting the rate constant of the O-methylation reaction for other substituted phenols, but also, as a part of the chemical artificial intelligence, the optimum [11C]MeI labeling conditions for a wide variety of phenol derivatives can be controlled.

Binding Strengths and Orientations in CO2 Adsorption on Cationic Scandium Oxides: Governing Factor Revealed by a Combined Infrared Spectroscopy and Theoretical Study

Carbon dioxide (CO2) adsorption is a critical step to curbing carbon emissions from fossil fuel combustion. Among various options, transition metal oxides have received extensive attention as promising CO2 adsorbents due to their affordability and sustainability for large-scale use. Here, the nature of binding interactions between CO2 molecules and cationic scandium oxides of different sizes, i.e., ScO+, Sc2O2+, and Sc3O4+, is investigated by mass-selective infrared photodissociation spectroscopy combined with quantum chemical calculations. The well-accepted electrostatic considerations failed to provide explanations for the trend in the binding strengths and variations in the binding orientations between CO2 and metal sites of cationic scandium oxides. The importance of orbital interactions in the driving forces for CO2 adsorption on cationic scandium oxides was revealed by energy decomposition analyses. A molecular surface property, known as the local electron attachment energy, is introduced to elucidate the binding affinity and orientation-specific reactivity of cationic scandium oxides upon the CO2 attachment. This study not only reveals the governing factor in the binding behaviors of CO2 adsorption on cationic scandium oxides but also serves as an archetype for predicting and rationalizing favorable binding sites and orientations in extended surface-adsorbate systems.

PB12+ and P2B12+/0/-: The Novel B12 Cage Doped by Nonmetallic P Atoms

A new kind of nonmetallic atom-doped boron cluster is described herein theoretically. When a phosphorus atom is added to the B12 motif and loses an electron, a novel B12 cage is obtained, composed of two B3 rings at both ends and one B6 ring in the middle, forming a triangular bifrustum. Interestingly, this B12 cage is formed by three B7 units joined together from three directions at an angle of 120°. When two P atoms are added to the B12 motif, this novel B12 cage is also obtained, and two P atoms are attached to the B3 rings at both ends of the triangular bifrustum, forming a triangular bipyramid (Johnson solid). Amazingly, the global minimums of neutral, monocationic, and monoanionic P2B12+/0/- have the same cage structure with a D3h symmetry; this is the smallest boron cage with the same structure. The P atom has five valence electrons, according to adaptive natural density partitioning bonding analyses of cage PB12+ and P2B12, in addition to one lone pair, the other three electrons of the P atom combine with an electron of each B atom on the B3 ring to form three 2c-2e σ bonds and form three electron sharing bonds with B atoms through covalent interactions, stabilizing the B12 cage. The calculated photoelectron spectra can be compared with future experimental values and provide a theoretical basis for the identification and confirmation of PnB12- (n = 1-2).

A broadly applicable quantitative relative reactivity model for nucleophilic aromatic substitution (SNAr) using simple descriptors

We report a multivariate linear regression model able to make accurate predictions for the relative rate and regioselectivity of nucleophilic aromatic substitution (S_NAr) reactions based on the electrophile structure. This model uses a diverse training/test set from experimentally-determined relative S_NAr rates between benzyl alcohol and 74 unique electrophiles, including heterocycles with multiple substitution patterns. There is a robust linear relationship between the experimental S_NAr free energies of activation and three molecular descriptors that can be obtained computationally: the electron affinity (EA) of the electrophile; the average molecular electrostatic potential (ESP) at the carbon undergoing substitution; and the sum of average ESP values for the ortho and para atoms relative to the reactive center. Despite using only simple descriptors calculated from ground state wavefunctions, this model demonstrates excellent correlation with previously measured S_NAr reaction rates, and is able to accurately predict site selectivity for multihalogenated substrates: 91% prediction accuracy across 82 individual examples. The excellent agreement between predicted and experimental outcomes makes this easy-to-implement reactivity model a potentially powerful tool for synthetic planning.

Machine learning meets mechanistic modelling for accurate prediction of experimental activation energies

Accurate prediction of chemical reactions in solution is challenging for current state-of-the-art approaches based on transition state modelling with density functional theory. Models based on machine learning have emerged as a promising alternative to address these problems, but these models currently lack the precision to give crucial information on the magnitude of barrier heights, influence of solvents and catalysts and extent of regio- and chemoselectivity. Here, we construct hybrid models which combine the traditional transition state modelling and machine learning to accurately predict reaction barriers. We train a Gaussian Process Regression model to reproduce high-quality experimental kinetic data for the nucleophilic aromatic substitution reaction and use it to predict barriers with a mean absolute error of 0.77 kcal mol-1 for an external test set. The model was further validated on regio- and chemoselectivity prediction on patent reaction data and achieved a competitive top-1 accuracy of 86%, despite not being trained explicitly for this task. Importantly, the model gives error bars for its predictions that can be used for risk assessment by the end user. Hybrid models emerge as the preferred alternative for accurate reaction prediction in the very common low-data situation where only 100-150 rate constants are available for a reaction class. With recent advances in deep learning for quickly predicting barriers and transition state geometries from density functional theory, we envision that hybrid models will soon become a standard alternative to complement current machine learning approaches based on ground-state physical organic descriptors or structural information such as molecular graphs or fingerprints.

Key activity descriptors of nickel-iron oxygen evolution electrocatalysts in the presence of alkali metal cations

Efficient oxygen evolution reaction (OER) electrocatalysts are pivotal for sustainable fuel production, where the Ni-Fe oxyhydroxide (OOH) is among the most active catalysts for alkaline OER. Electrolyte alkali metal cations have been shown to modify the activity and reaction intermediates, however, the exact mechanism is at question due to unexplained deviations from the cation size trend. Our X-ray absorption spectroelectrochemical results show that bigger cations shift the Ni^2+/(3+δ)+ redox peak and OER activity to lower potentials (however, with typical discrepancies), following the order CsOH > NaOH ≈ KOH > RbOH > LiOH. Here, we find that the OER activity follows the variations in electrolyte pH rather than a specific cation, which accounts for differences both in basicity of the alkali hydroxides and other contributing anomalies. Our density functional theory-derived reactivity descriptors confirm that cations impose negligible effect on the Lewis acidity of Ni, Fe, and O lattice sites, thus strengthening the conclusions of an indirect pH effect.

It is commonly accepted that electrolyte alkali metal cations modify the catalytic activity for oxygen evolution reaction. Here the authors challenge this assumption, showing that the activity is actually affected by a change in the electrolyte pH rather than a specific alkali cation.

tele-Substitution Reactions in the Synthesis of a Promising Class of 1,2,4-Triazolo[4,3-a]pyrazine-Based Antimalarials

We have discovered and studied a tele-substitution reaction in a biologically important heterocyclic ring system. Conditions that favor the tele-substitution pathway were identified: the use of increased equivalents of the nucleophile or decreased equivalents of base or the use of softer nucleophiles, less polar solvents, and larger halogens on the electrophile. Using results from X-ray crystallographic and isotope labeling experiments, a mechanism for this unusual transformation is proposed. We focused on this triazolopyrazine as it is the core structure of the in vivo active antiplasmodium compounds of Series 4 of the Open Source Malaria consortium.

F-Halogen Bond: Conditions for Its Existence

The question as to whether the F atom can engage in a halogen bond (XB) remains unsettled. This issue is addressed via density functional theory calculations which pair a wide range of organic and inorganic F-acids with various sorts of Lewis bases. From an energetic perspective, perfluorinated hydrocarbons with sp, sp², or sp³ C-hybridization are unable to form an XB with an N-base, but a very weak bond can be formed if electron-withdrawing C≡N substituents are added to the acid. There is little improvement for inorganic acids O₂NF, FOF, ClF, BrF, SiF₄, or GeF₄, but F₂ is capable of a stronger XB of up to 5 kcal/mol. These results are consistent with a geometric criterion, which compares the intermolecular equilibrium distance with the sum of atomic van der Waals radii. The intensity of the σ-hole on the F atom has predictive value in that a V_s,max of at least 10-15 kcal/mol is required for XB formation. Adding a positive charge to the Lewis acid enhances the strength of any XB and even more so if the base is anionic. The acid-base interaction induces a contraction of the r(AF) covalent bond in the acid in most cases and a deshielding of the NMR signal of the F nucleus.

Towards the critical understanding of selected vibrational features in biologically important dicyano aromatic conjugated molecules: Importance of electron donating/withdrawal groups and geometry associated with dicyano group

The Raman spectra of a series of synthesized DC molecules (benzylidene malononitrile derivatives) with different electron donating (EDG) and electron withdrawing (EWG) group have been presented and analyzed with DFT calculated spectra. In particular, different functional groups effect on cyano stretching (∼2200 cm^-1), phenyl ring breathing and alkenic double bond stretch which often appears mixed up (1475-1650 cm^-1) are studied systematically for several aromatic conjugated DC derivatives. Interestingly, symmetric stretching frequency of the DC compounds having two CN groups at geminal position appears at higher wavenumber (by 11-15 cm^-1) compared to their corresponding asymmetric stretch frequency. Angle (between dicyano group) dependent theoretical study indicates that the relative appearance of cyano symmetric/anti-symmetric stretching frequency depends on whether dicyano groups are at the geminal or vicinal position and the angle between them. Complete band assignments of observed Raman frequencies have been performed by potential energy distributions (PEDs) available in GAR2PED software. Our results will help to understand the vibrational feature of this important class of compounds in biological medium when used as probe.

Propagation of Holes and Electrons in Metal-Organic Frameworks

Charge transport in two zinc metal-organic frameworks (MOFs) has been investigated using periodic semiempirical molecular orbital calculations with the AM1* Hamiltonian. Restricted Hartree-Fock calculations underestimate the band gap using Koopmans theorem (ca. 2 eV compared to the experimental value of 2.8 eV). However, it almost doubles when the constraint on the wave function to remain spin-restricted is removed and the energies of the UHF Natural Orbitals are used. Charge-transport simulations using propagation of the electron- or hole-density in imaginary time allow charge-transport paths and mechanisms to be determined. The calculated relative mobilities in the directions of the three crystal axes agree with experimental expectations, but the absolute values are not reliable using the current technique. Hole-mobility along the crystal c-axis (along the metal stacks) is found to be 13 times higher in the zinc MOF with anthracene linker (Zn-ANMOF-74) than in the other directions, whereas the factor is far smaller (1.7) for electron mobility. Directional preferences are far less distinct in the equivalent structure with phenyl linkers (Zn-MOF-74). The imaginary-time simulation technique does not give quantitative mobilities. The simulations reveal a change in mechanism between the different directions: Coherent polaron migration is observed along the stacks but tunneling hops between them.

The local electron attachment energy and the electrostatic potential as descriptors of surface–adsorbate interactions

Local DFT-based properties are used for fast rationalization and accurate estimations of local surface reactivity of metal and oxide compounds.

Regium bonds between Mn clusters (M = Cu, Ag, Au and n = 2-6) and nucleophiles NH3 and HCN

The most stable geometries of the coinage metal (or regium) atom (Cu, Ag, Au) clusters Mn for n up to 6 are all planar, and adopt the lowest possible spin multiplicity. Clusters with even numbers of M atoms are thus singlets, while those with odd n are open-shell doublets. Examination of the molecular electrostatic potential (MEP) of each cluster provides strong indications of the most likely site of attack by an approaching nucleophile, generally one of two positions. A nucleophile (NH3 or HCN) most favorably approaches one particular M atom of each cluster, rather than a bond midpoint or face. In the closed-shell clusters, the interaction energies are highly dependent upon the intensity of the MEP, but this correlation fades for the open-shell systems studied in this work. The strength of the interaction is also closely related to the basicity of the nucleophile. Regium bond energies can be more than 30 kcal mol-1 and tend to follow the Au > Cu > Ag order. These interaction energies are in large part derived from Coulombic attraction, with a smaller orbital interaction contribution.

Electronegativity-a perspective

Electronegativity is a very useful concept but it is not a physical observable; it cannot be determined experimentally. Most practicing chemists view it as the electron-attracting power of an atom in a molecule. Various formulations of electronegativity have been proposed on this basis, and predictions made using different formulations generally agree reasonably well with each other and with chemical experience. A quite different approach, loosely linked to density functional theory, is based on a ground-state free atom or molecule, and equates electronegativity to the negative of an electronic chemical potential. A problem that is encountered with this approach is the differentiation of a noncontinuous function. We show that this approach leads to some results that are not chemically valid. A formulation of atomic electronegativity that does prove to be effective is to express it as the average local ionization energy on an outer contour of the atom's electronic density.

σ-Holes and σ-lumps direct the Lewis basic and acidic interactions of noble metal nanoparticles: introducing regium bonds

Using local DFT-based probes for electrostatic as well as charge transfer/polarization interactions, we are able to characterize Lewis basic and acidic sites on copper, silver and gold nanoparticles.