Comment on “Natural Bond Orbitals and the Nature of the Hydrogen Bond”

Novel 2-mercaptobenzothiazole and 2-mercaptobenzimidazole-derived Ag16 and Ag18 nanoclusters: synthesis and optical properties

Silver and gold nanoclusters are promising nanomaterials for various applications such as sensing, catalysis, and bioimaging. However, their synthetic control and repeatability, and determination of their structures are highly complicated. Only a handful of crystal structures of silver nanoclusters (AgNCs) have been reported, while structures of a few others have been reported with the help of mass spectrometry. We synthesized two AgNCs, viz., Ag-MBTNC (Ag16 cluster) and Ag-MBINC (Ag18 cluster) respectively stabilized by 2-mercaptobenzothiazole (2-MBT) and 2-mercaptobenzimidazole (2-MBI) with excellent repeatability; determined their composition and plausible structures using XPS, TGA and MALDI-TOF mass spectrometry; and compared their optical properties. Interestingly, Ag-MBTNC is fluorescent while Ag-MBINC is not, although these are synthesized using stabilizing ligands that have difference in only one atom. The structural features of the clusters are found to be similar but they have contrasting optical behaviours due to the effect of one S atom (in 2-MBT) in place of one N atom (in 2-MBI).

MEDT analysis of mechanism and selectivities in non-catalyzed and lewis acid-catalyzed diels-alder reactions between R-carvone and isoprene

Within the context of Molecular Electronic Density Theory (MEDT), this study investigates the Diels-Alder reaction among isoprene (2) and R-carvone (1R) applying DFT simulations, with and without Lewis acid (LA) catalysis. The results show that carvone (1R) acts as an electrophile and isoprene (2) as a nucleophile in a polar process. LA catalysis increases the electrophilicity of carvone, thereby improving the reactivity and selectivity of the reaction by reducing the activation Gibbs free energy. Parr functions reveal that the C5=C6 double bond is more reactive than the C9=C10 double bond, indicating chemoselectivity. The examination of the Electron Localization Function (ELF) reveals high regio- and stereoselectivity, indicating an asynchronous mechanism for the LA-catalyzed DA reaction. Furthermore, it is suggested that cycloadduct 3 has great anti-HIV potential because it exhibits lower binding energies than azidothymidine (AZT) in the docking studies of cycloadducts 3 and 4 amongst a primary HIV-1protein (1A8O plus 5W4Q).

Halogen Bond via an Electrophilic π-Hole on Halogen in Molecules: Does It Exist?

This study reveals a new non-covalent interaction called a π-hole halogen bond, which is directional and potentially non-linear compared to its sister analog (σ-hole halogen bond). A π-hole is shown here to be observed on the surface of halogen in halogenated molecules, which can be tempered to display the aptness to form a π-hole halogen bond with a series of electron density-rich sites (Lewis bases) hosted individually by 32 other partner molecules. The [MP2/aug-cc-pVTZ] level characteristics of the π-hole halogen bonds in 33 binary complexes obtained from the charge density approaches (quantum theory of intramolecular atoms, molecular electrostatic surface potential, independent gradient model (IGM-δginter)), intermolecular geometries and energies, and second-order hyperconjugative charge transfer analyses are discussed, which are similar to other non-covalent interactions. That a π-hole can be observed on halogen in halogenated molecules is substantiated by experimentally reported crystals documented in the Cambridge Crystal Structure Database. The importance of the π-hole halogen bond in the design and growth of chemical systems in synthetic chemistry, crystallography, and crystal engineering is yet to be fully explicated.

Erythronium Bonds: Noncovalent Interactions Involving Group 5 Elements as Electron-Density Acceptors

Analyses of the Cambridge Structural Database and theoretical calculations (PBE0-D3/def2-TZVP level, atoms-in-molecules, natural bond orbital studies) prove the formation of net attractive noncovalent interactions between group 5 elements and electron-rich atoms (neutral or anionic). These kinds of bonding are markedly different from coordination bonds formed by the same elements and possess the distinctive features of σ-hole interactions. The term erythronium bond is proposed to denote these bonds. X-ray structures of vanadate-dependent bromoperoxidases show that these interactions are present also in biological systems.

Role of the solvent polarity on the optical and electronic characteristics of 1-iodoadamantane

The natural absorbance caused by the chromophore and chemical behavior of 1-iodoadamantane is highly influenced by the polarity of different solvent environments. This gives rise to the solvatochromatic shifts in the optical absorption and electronic structure and the experimentally measured UV-vis absorption spectra show significant solvatochromic shifts with respect to the solvent polarity. The absorption shift for both σ to σ*and n to σ* electronic transitions are more dominant in polar solvents than in nonpolar solvents. To obtain a better understanding of the impact of solvent polarity on the 1-iodoadamantane at the molecular level, computational calculations were carried out through implicit solvation. According to this, changes in the HOMO and LUMO energies and electron density distributions of various solvent continuums demonstrate the influence of solvent polarity on the HOMO and LUMO energy levels of the chemical system. This also shows an increment in the HOMO-LUMO gap with respect to the polarity of the solvent.

A molecular electron density theory study of asymmetric Diels-Alder [4 + 2] reaction's mechanism of furan with three substituted alkynes (5-R substituted-3-(3-(phenylsulfonyl)-propioloyl)-oxazolidin-2-one)

CONTEXT

The [4 +2 ] cycloaddition reactions between furan and three substituted alkynes (5-R-substituted-3-(3-(phenylsulfonyl)-propioloyl)-oxazolidin-2-one) have been investigated using the MEDT approach. Reactivity indices, reaction pathways, and activation energies are calculated. In an investigation of conceptual DFT indices, furan acts as a nucleophile, while the three substituted alkynes (5-R-substituted-3-(3-(phenylsulfonyl)-propioloyl)-oxazolidin-2-one) function as electrophiles in this reaction. The cycloaddition is regioiselective, as demonstrated by the activation and reaction energies, in clear agreement with the experiment's results. Hetero Diels-Alder [4 + 2] cycloadditions occur following a non-concerted two stages one-step molecular mechanism.

METHODS

For the purpose of this study, all calculations were performed using the Gaussian 09 software. Optimization was achieved through Berny's computational gradient optimization method, employing the B3LYP functional and the 6-31G(d) basis set. Analysis of both local and global reactivity indices provided insights into the reactivity tendencies of the reactants, distinguishing between electrophilic and nucleophilic characteristics via Parr functions. Frequency calculations were employed to identify and characterize stationary points, with transition states indicated by a single imaginary frequency and positive values of all frequencies for reactants and product. The electron localization function (ELF) was investigated using the Multiwfn software within the context of topological analyses.

Methylammonium Tetrel Halide Perovskite Ion Pairs and Their Dimers: The Interplay between the Hydrogen-, Pnictogen- and Tetrel-Bonding Interactions

The structural stability of the extensively studied organic-inorganic hybrid methylammonium tetrel halide perovskite semiconductors, MATtX₃ (MA = CH₃NH₃⁺; Tt = Ge, Sn, Pb; X = Cl, Br, I), arises as a result of non-covalent interactions between an organic cation (CH₃NH₃⁺) and an inorganic anion (TtX₃^-). However, the basic understanding of the underlying chemical bonding interactions in these systems that link the ionic moieties together in complex configurations is still limited. In this study, ion pair models constituting the organic and inorganic ions were regarded as the repeating units of periodic crystal systems and density functional theory simulations were performed to elucidate the nature of the non-covalent interactions between them. It is demonstrated that not only the charge-assisted N-H···X and C-H···X hydrogen bonds but also the C-N···X pnictogen bonds interact to stabilize the ion pairs and to define their geometries in the gas phase. Similar interactions are also responsible for the formation of crystalline MATtX₃ in the low-temperature phase, some of which have been delineated in previous studies. In contrast, the Tt···X tetrel bonding interactions, which are hidden as coordinate bonds in the crystals, play a vital role in holding the inorganic anionic moieties (TtX₃^-) together. We have demonstrated that each Tt in each [CH₃NH₃⁺•TtX₃^-] ion pair has the capacity to donate three tetrel (σ-hole) bonds to the halides of three nearest neighbor TtX₃^- units, thus causing the emergence of an infinite array of 3D TtX₆^4- octahedra in the crystalline phase. The TtX₄^4- octahedra are corner-shared to form cage-like inorganic frameworks that host the organic cation, leading to the formation of functional tetrel halide perovskite materials that have outstanding optoelectronic properties in the solid state. We harnessed the results using the quantum theory of atoms in molecules, natural bond orbital, molecular electrostatic surface potential and independent gradient models to validate these conclusions.

Importance of Noncovalent Interactions Involving Sulfur Atoms in Thiopeptide Antibiotics─Glycothiohexide α and Nocathiacin I

Noncovalent interactions involving sulfur atoms play essential roles in protein structure and function by significantly contributing to protein stability, folding, and biological activity. Sulfur is a highly polarizable atom that can participate in many types of noncovalent interactions including hydrogen bonding, sulfur-π interactions, and S-lone pair interactions, but the impact of these sulfur-based interactions on molecular recognition and drug design is still often underappreciated. Here, we examine, using quantum chemical calculations, the roles of sulfur-based noncovalent interactions in complex naturally occurring molecules representative of thiopeptide antibiotics: glycothiohexide α and its close structural analogue nocathiacin I. While donor-acceptor orbital interactions make only very small contributions, electrostatic and dispersion contributions are predicted to be significant in many cases. In pursuit of understanding the magnitudes and nature of these noncovalent interactions, we made potential structural modifications that could significantly expand the chemical space of effective thiopeptide antibiotics.

Stability and bonding of carbon(0)-iron-N2 complexes relevant to nitrogenase co-factor: EDA-NOCV analyses

The factors/structural features which are responsible for the binding, activation and reduction of N₂ to NH₃ by FeMoco of nitrogenase have not been completely understood well. Several relevant model complexes by Holland et al. and Peters et al. have been synthesized, characterized and studied by theoretical calculations. For a matter of fact, those complexes are much different than real active N₂ -binding Fe-sites of FeMoco, which possesses a central C(4-) ion having an eight valence electrons as an μ₆ -bridge. Here, a series of [(S₃ C(0))Fe(II/I/0)-N₂ ]^n- complexes in different charged/spin states containing a coordinated σ- and π-donor C(0)-atom which possesses eight outer shell electrons [carbone, (Ph₃ P)₂ C(0); Ph₃ P→C(0)←PPh₃ ] and three S-donor sites (i.e. ^- S-Ar), have been studied by DFT, QTAIM, and EDA-NOCV calculations. The effect of the weak field ligand on Fe-centres and the subsequent N₂ -binding has been studied by EDA-NOCV analysis. The role of the oxidation state of Fe and N₂ -binding in different charged and spin states of the complex have been investigated by EDA-NOCV analyses. The intrinsic interaction energies of the Fe-N₂ bond are in the range from -42/-35 to -67 kcal/mol in their corresponding ground states. The S₃ C(0) donor set is argued here to be closer to the actual coordination environment of one of the six Fe-centres of nitrogenase. In comparison, the captivating model complexes reported by Holland et al. and Peter et al. possess a stronger π-acceptor C-ring (S₂ C_ring donor, π-C donor) and stronger donor set like CP₃ (σ-C donor) ligands, respectively.

High-Density “Windowpane” Coordination Patterns of Water Clusters and Their NBO/NRT Characterization

Cluster mixture models for liquid water at higher pressures suggest the need for water clusters of higher coordination and density than those commonly based on tetrahedral H-bonding motifs. We show here how proton-ordered water clusters of increased coordination and density can assemble from a starting cyclic tetramer or twisted bicyclic (Möbius-like) heptamer to form extended Aufbau sequences of stable two-, three-, and four-coordinate “windowpane” motifs. Such windowpane clusters exhibit sharply reduced (~90°) bond angles that differ appreciably from the tetrahedral angles of idealized crystalline ice I_h. Computed free energy and natural resonance theory (NRT) bond orders provide quantitative descriptors for the relative stabilities of clusters and strengths of individual coordinative linkages. The unity and consistency of NRT description is demonstrated to extend from familiar supra-integer bonds of the molecular regime to the near-zero bond orders of the weakest linkages in the present H-bond clusters. Our results serve to confirm that H-bonding exemplifies resonance–covalent (fractional) bonding in the sub-integer range and to further discount the dichotomous conceptions of “electrostatics” for intermolecular bonding vs. “covalency” for intramolecular bonding that still pervade much of freshman-level pedagogy and force-field methodology.

Weakly Bound Complex Formation between HCN and CH3Cl: A Matrix-Isolation and Computational Study

The matrix-isolated infrared spectrum of a hydrogen cyanide-methyl chloride complex was investigated in a solid argon matrix. HCN and CH₃Cl were co-condensed onto a substrate held at 10 K with an excess of argon gas, and the infrared spectrum was measured using Fourier-transform infrared spectroscopy. Quantum chemical geometry optimization, harmonic frequency, and natural bonding orbital calculations indicate stabilized hydrogen- and halogen-bonded structures. The two resulting weakly bound complexes are both composed of one CH₃Cl molecule bound to a (HCN)₃ subunit, where the three HCN molecules are bound head-to-tail in a ring formation. Our study suggests that─in the presence of CH₃Cl─the formation of (HCN)₃ is promoted through complexation. Since HCN aggregates are an important precursor to prebiotic monomers (amino acids and nucleobases) and other life-bearing polymers, this study has astrophysical implications toward the search for life in space.

LModeA-nano: A PyMOL Plugin for Calculating Bond Strength in Solids, Surfaces, and Molecules via Local Vibrational Mode Analysis

The analysis of chemical bonding in crystal structures and surfaces is an important research topic in theoretical chemistry. In this work, we present a PyMOL plugin, named LModeA-nano, as implementation of the local vibrational mode theory for periodic systems (Tao et al. J. Chem. Theory Comput. 2019, 15, 1761) assessing bond strength in terms of local stretching force constants in extended systems of one, two, and three dimensions. LModeA-nano can also analyze chemical bonds in isolated molecular systems thus enabling a head-to-head comparison of bond strength across systems with different dimensions in periodicity (0-3D). The new code is interfaced to the output generated by various solid-state modeling packages including VASP, CP2K, Quantum ESPRESSO, CASTEP, and CRYSTAL. LModeA-nano is cross-platform, open-source and freely available on GitHub: https://github.com/smutao/LModeA-nano.

Ab Initio Study of Hydrostable Metal-Organic Frameworks for Postsynthetic Modification and Tuning toward Practical Applications

Metal-organic frameworks (MOFs), a subclass of nanoporous coordination polymers, have emerged as one of the most promising next-generation materials. The postsynthetic modification method, a strategy that provides tunability and control of these materials, plays an important role in enhancing its properties and functionalities. However, knowing adjustments which leads to a desired structure-function a priori remains a challenge. In this comprehensive study, the intermolecular interactions between 21 industrially important gases and a hydrostable STAM-17-OEt MOF were investigated using density functional theory. Substitutions on its 5-ethoxy isophthalate linker included two classes of chemical groups, electron-donating (-NH₂, -OH, and -CH₃) and electron-withdrawing (-CN, -COOH, and -F), as well as the effect of mono-, di-, and tri-substitutions. This resulted in 651 unique MOF-gas complexes. The adsorption energies at the ground state and room temperature, bond lengths, adsorption geometry, natural bond orbital analysis of the electric structure, HOMO-LUMO interactions, and the predicted zwitterionic properties are presented and discussed. This study provides a viable strategy for the functionalization, which leads to the strongest affinity for each gas, an insight into the role of different chemical groups in adsorbing various gas molecules, and identifies synthetic routes for moderating the gas adsorption capacity and reducing water adsorption. Recommendations for various applications are discussed. A custom Python script to assess and visualize the hypothetical separation of two equal gas mixtures of interest is provided. The methodology presented here provides new opportunities to expand the chemical space and physical properties of STAM-17-OEt and advances the development of other hydrostable MOFs.

Anti-Electrostatic Pi-Hole Bonding: How Covalency Conquers Coulombics

Intermolecular bonding attraction at π-bonded centers is often described as "electrostatically driven" and given quasi-classical rationalization in terms of a "pi hole" depletion region in the electrostatic potential. However, we demonstrate here that such bonding attraction also occurs between closed-shell ions of like charge, thereby yielding locally stable complexes that sharply violate classical electrostatic expectations. Standard DFT and MP2 computational methods are employed to investigate complexation of simple pi-bonded diatomic anions (BO-, CN-) with simple atomic anions (H-, F-) or with one another. Such "anti-electrostatic" anion-anion attractions are shown to lead to robust metastable binding wells (ranging up to 20-30 kcal/mol at DFT level, or still deeper at dynamically correlated MP2 level) that are shielded by broad predissociation barriers (ranging up to 1.5 Å width) from long-range ionic dissociation. Like-charge attraction at pi-centers thereby provides additional evidence for the dominance of 3-center/4-electron (3c/4e) nD-π*AX interactions that are fully analogous to the nD-σ*AH interactions of H-bonding. Using standard keyword options of natural bond orbital (NBO) analysis, we demonstrate that both n-σ* (sigma hole) and n-π* (pi hole) interactions represent simple variants of the essential resonance-type donor-acceptor (Bürgi-Dunitz-type) attraction that apparently underlies all intermolecular association phenomena of chemical interest. We further demonstrate that "deletion" of such π*-based donor-acceptor interaction obliterates the characteristic Bürgi-Dunitz signatures of pi-hole interactions, thereby establishing the unique cause/effect relationship to short-range covalency ("charge transfer") rather than envisioned Coulombic properties of unperturbed monomers.

Chameleonic Metal-bound Isocyanides: π-Donating CuI-center Imparts a Nucleophilicity to the Isocyanide Carbon toward Halogen Bonding

In the structures of the isostructural cocrystals [CuI3(CNXyl)3]·CHX3 (X = Br, I), two adjacent CuI-bound isocyanide groups, whose carbon lone pairs are blocked by the ligation, exhibit nucleophilic properties induced...

Electron Flow Characterization of Charge Transfer for Carbonic Acid to Strong Base Proton Transfer in Aqueous Solution

Protonation of the strong base methylamine CH₃NH₂ by carbonic acid H₂CO₃ in aqueous solution, HOCOOH···NH₂CH₃ → HOCOO^-···⁺HNH₂CH₃, has been previously studied ( J. Phys. Chem. B 2016, 109, 2271-2280; J. Phys. Chem. B 2016, 109, 2281-2290) via Car-Parinnello molecular dynamics. This proton transfer (PT) reaction within a hydrogen (H)-bonded complex was found to be barrierless and very rapid, with key reaction coordinates comprising the proton coordinate, the H-bond separation R_ON, and a solvent coordinate, reflecting the water solvent rearrangement involved in the neutral to ion pair conversion. In the present work, the reaction's charge flow aspects are analyzed in detail, especially a description via Mulliken charge transfer for PT (MCTPT). A natural bond orbital analysis and some extensions of them are employed for the complex's electronic structure during the reaction trajectories. Results demonstrate that consistent with the MCTPT picture, the charge transfer (CT) occurs from a methylamine base nonbonding orbital to a carbonic acid antibonding orbital. A complementary MCTPT reaction product perspective of CT from the antibonding orbital of the HN⁺ moiety to the nonbonding orbital of the oxygen in the H-bond complex is also presented. σ_OH and σ_HN⁺ bond order expressions show this CT to occur within the H-bond OHN triad, an aspect key for simultaneous bond-breaking and -forming in the PT reaction.

Bonding analysis of water clusters using quasi-atomic orbitals

The quasi-atomic orbital (QUAO) bonding analysis introduced by Ruedenberg and co-workers is used to develop an understanding of the hydrogen bonds in small water clusters, from the dimer through the hexamer (bag, boat, book, cyclic, prism and cage conformers). Using kinetic bond orders as a metric, it is demonstrated that as the number of waters in simple cyclic clusters increases, the hydrogen bonds strengthen, from the dimer through the cyclic hexamer. However, for the more complex hexamer isomers, the strength of the hydrogen bonds varies, depending on whether the cluster contains double acceptors and/or double donors. The QUAO analysis also reveals the three-center bonding nature of hydrogen bonds in water clusters.

Electrostatics, Charge Transfer, and the Nature of the Halide-Water Hydrogen Bond

Binary halide-water complexes X^-(H₂O) are examined by means of symmetry-adapted perturbation theory, using charge-constrained promolecular reference densities to extract a meaningful charge-transfer component from the induction energy. As is known, the X^-(H₂O) potential energy surface (for X = F, Cl, Br, or I) is characterized by symmetric left and right hydrogen bonds separated by a C_2v-symmetric saddle point, with a tunneling barrier height that is <2 kcal/mol except in the case of F^-(H₂O). Our analysis demonstrates that the charge-transfer energy is correspondingly small (<2 kcal/mol except for X = F), considerably smaller than the electrostatic interaction energy. Nevertheless, charge transfer plays a crucial role determining the conformational preferences of X^-(H₂O) and provides a driving force for the formation of quasi-linear X··· H-O hydrogen bonds. Charge-transfer energies correlate well with measured O-H vibrational redshifts for the halide-water complexes and also for OH^-(H₂O) and NO₂^-(H₂O), providing some indication of a general mechanism.

A reliable and efficient resonance theory based on analysis of DFT wave functions

Due to methodological difficulties and limitations of applicability, a quantitative bonding analysis based on the theory of resonance is presently not as convenient and popular as that based on the molecular orbital (MO) methods. Here, we propose an efficient quantitative resonance theory by expanding the DFT wave function in terms of a complete set of Lewis structures. By rigorously separating the resonance subsystem represented by a set of localized MOs, this approach is able to treat large molecules, nonplanar π-conjugate systems, and bonding systems mixing both σ and π electrons. Assessment in 2c-2e systems suggests a new projection-weighted symmetric orthogonalization method to evaluate the weights of resonance contributors, which overcomes the drawbacks of other weighting schemes. Applications to benzene, naphthalene and chlorobenzene show that the present method is insensitive to the basis set employed in the DFT calculations, and to the choices of the independent Lewis set determined by Rumer's rule. Advanced applications to diverse chemical problems provide unique and valuable insights into the understanding of hydrogen bonding, the π substituent effect on benzene, and the mechanism of Diels-Alder reactions.

Comparison of Chemical and Interpretative Methods: the Carbon-Boron π-Bond as a Test Case*

Quantum chemical calculations and NBO, ETS-NOCV, QTAIM and ELF interpretative approaches have been carried out on C-donor ligand-stabilized dihydrido borenium cations. Numerous descriptors of the C-B π-bond strength obtained from orbital localization, energy partitioning or topological methods as well as from structural and chemical parameters have been calculated for 39 C-donor ligands including N-heterocyclic carbenes and carbones. Comparison of the results allows the identification of relative and absolute descriptors of the π interaction. For both families of descriptors excellent correlations are obtained. This enables the establishment of a π-donation capability scale and shows that the interpretative methods, despite their conceptual differences, describe the same chemical properties. These results also reveal noticeable shortcomings in these popular methods, and some precautions that need to be taken to interpret their results adequately.

Cyclic Octamer of Hydroxyl-functionalized Cations with Net Charge Q=+8e Kinetically Stabilized by a 'Molecular Island' of Cooperative Hydrogen Bonds

Cyclic octamers are well-known structural motifs in chemistry, biology and physics. These include covalently bound cyclic octameric sulphur, cylic octa-alkanes, cyclo-octameric peptides as well as hydrogen-bonded ring clusters of alcohols. In this work, we show that even calculated cyclic octamers of hydroxy-functionalized pyridinium cations with a net charge Q=+8e are kinetically stable. Eight positively charged cations are kept together by hydrogen bonding despite the strong Coulomb repulsive forces. Sufficiently long hydroxy-octyl chains prevent "Coulomb explosion" by increasing the distance between the positive charges at the pyridinium rings, reducing the Coulomb repulsion and thus strengthen hydrogen bonds between the OH groups. The eightfold positively charged cyclic octamer shows spectroscopic properties similar to those obtained for hydrogen-bonded neutral cyclic octamers of methanol. Thus, the area of the hydrogen bonded OH ring represents a 'molecular island' within an overall cationic environment. Although not observable, the spectroscopic properties and the correlated NBO parameters of the calculated cationic octamer support the detection of smaller cationic clusters in ionic liquids, which we observed despite the competition with ion pairs wherein attractive Coulomb forces enhance hydrogen bonding between cation and anion.

Clusters of Hydroxyl-Functionalized Cations Stabilized by Cooperative Hydrogen Bonds: The Role of Polarizability and Alkyl Chain Length

We explore quantum chemical calculations for studying clusters of hydroxyl-functionalized cations kinetically stabilized by hydrogen bonding despite strongly repulsive electrostatic forces. In a comprehensive study, we calculate clusters of ammonium, piperidinium, pyrrolidinium, imidazolium, pyridinium, and imidazolium cations, which are prominent constituents of ionic liquids. All cations are decorated with hydroxy-alkyl chains allowing H-bond formation between ions of like charge. The cluster topologies comprise linear and cyclic clusters up to the size of hexamers. The ring structures exhibit cooperative hydrogen bonds opposing the repulsive Coulomb forces and leading to kinetic stability of the clusters. We discuss the importance of hydrogen bonding and dispersion forces for the stability of the differently sized clusters. We find the largest clusters when hydrogen bonding is maximized in cyclic topologies and dispersion interaction is properly taken into account. The kinetic stability of the clusters with short-chained cations is studied for the different types of cations ranging from hard to polarizable or exhibiting additional functional groups such as the acidic C(2)-H position in the imidazolium-based cation. Increasing the alkyl chain length, the cation effect diminishes and the kinetic stability is exclusively governed by the alkyl chain tether increasing the distance between the positively charged rings of the cations. With adding the counterion tetrafluoroborate (BF4-) to the cationic clusters, the binding energies immediately switch from strongly positive to strongly negative. In the neutral clusters, the OH functional groups of the cations can interact either with other cations or with the anions. The hexamer cluster with the cyclic H-bond motive and "released" anions is almost as stable as the hexamer built by H-bonded ion pairs exclusively, which is in accord with recent IR spectra of similar ionic liquids detecting both types of hydrogen bonding. For the cationic and neutral clusters, we discuss geometric and spectroscopic properties as sensitive probes of opposite- and like-charge interaction. Finally, we show that NMR proton chemical shifts and deuteron quadrupole coupling constants can be related to each other, allowing to predict properties which are not easily accessible by experiment.

Chemical Bonding: The Journey from Miniature Hooks to Density Functional Theory

Our modern understanding of chemistry is predicated upon bonding interactions between atoms and ions resulting in the assembly of all of the forms of matter that we encounter in our daily life. It was not always so. This review article traces the development of our understanding of bonding from prehistory, through the debates in the 19th century C.E. bearing on valence, to modern quantum chemical models and beyond.

Quantitative decomposition of resonance-assisted hydrogen bond energy in β-diketones into resonance and hydrogen bonding (π- and σ-) components using molecular tailoring and function-based approaches

Using the molecular tailoring and function-based approaches allows one to divide the energy of the O─H⋯O═C resonance-assisted hydrogen bond in a series of the β-diketones into resonance and hydrogen bonding components. The magnitude of the resonance component is assessed as about 6 kcal mol^-1 . This value increases by ca. 1 kcal mol^-1 on going from the weak to strong resonance-assisted hydrogen bonding. The magnitude of the hydrogen bonding component varies in the wide range from 2 to 20 kcal mol^-1 depending on the structure of the β-diketone in question.

Innovative Organic Electroluminescent Materials with a Doublet Ground State: A Theoretical Investigation

The displaced and distorted harmonic oscillator model, which has been proven to be appropriate in calculating vibronic spectra, is employed to treat the emission spectrum of title molecules in combination with a thermal vibration correlation function. The calculated results indicate that the main peak of the emission spectrum is visibly impacted by the normal modes with lower frequencies and that the shoulder peak is originated from the middle-frequency modes. On the level of time-dependent density functional theory (TDDFT), the calculated fluorescence lifetimes of TTM-3NCz and TTM-3PCz are 22.1 and 26.0 ns, respectively, which happen to coincide with the observed values of TTM-3NCz (17.2 ns) and TTM-3PCz (21.2 ns). The above data indicate that both the calculated radiative decay rates are reasonable at room temperature. Furthermore, we investigate the influence of the Duschinsky effect on the fluorescence quantum efficiency (FQE). When it is considered, the predicted FQE of the TTM-3NCz molecule is only 0.11%, and the observed value (49% in toluene) deviates significantly. If we ignore the Duschinsky effect, the FQE of TTM-3NCz increases dramatically to 41.8%. For the TTM-3PCz molecule (the FQE is 46% in toluene), the calculated FQE is 0.042% with the Duschinsky effect and increases to 45.2% without the Duschinsky effect. This phenomenon might be related to external factors and the nature of the TDDFT only considering a single configuration. In addition, the fluorescent properties of the fluorinated TTM-3NCz molecules are studied predictably. The obtained results show that the perfluorinated TTM-3NCz shows better luminous performance due to larger oscillator strength. Finally, the dimers, which are composed of both single title molecules, are explored theoretically to determine how they impact the fluorescent property; however, the effect can be nearly eliminated because of the small binding energies.

Nine questions on energy decomposition analysis

The paper collects the answers of the authors to the following questions: Is the lack of precision in the definition of many chemical concepts one of the reasons for the coexistence of many partition schemes? Does the adoption of a given partition scheme imply a set of more precise definitions of the underlying chemical concepts? How can one use the results of a partition scheme to improve the clarity of definitions of concepts? Are partition schemes subject to scientific Darwinism? If so, what is the influence of a community's sociological pressure in the "natural selection" process? To what extent does/can/should investigated systems influence the choice of a particular partition scheme? Do we need more focused chemical validation of Energy Decomposition Analysis (EDA) methodology and descriptors/terms in general? Is there any interest in developing common benchmarks and test sets for cross-validation of methods? Is it possible to contemplate a unified partition scheme (let us call it the "standard model" of partitioning), that is proper for all applications in chemistry, in the foreseeable future or even in principle? In the end, science is about experiments and the real world. Can one, therefore, use any experiment or experimental data be used to favor one partition scheme over another? © 2019 Wiley Periodicals, Inc.

Electrostatics and Dispersion in X-H···Y (X = C, N, O; Y = N, O) Hydrogen Bonds and Their Role in X-H Vibrational Frequency Shifts

The frequency shifts of donor stretching vibration in X-H···Y (X = C, N, O; Y = N, O) hydrogen-bonded complexes of phenylacetylene, indole, and phenol are linearly correlated with the electrostatic component of the interaction energy. This linear correlation suggests that the electrostatic component, which is the first-order perturbative correction to the stabilization energy, is essentially localized on the X-H group. The linear correlation suggests that the electrostatic tuning rate, which is a measure of the X-H oscillator to undergo shifts upon hydrogen bonding per unit increase in the electrostatic component of the stabilization energy, was found to be in the order of O-H > N-H > C-H. Interestingly, for each of the donor groups, viz., C-H, N-H, and O-H, the vibrational frequency shifts were inversely correlated to the dipole moment of the acceptor separately, which is counterintuitive vis-à-vis the electrostatic component. This implies that extrapolation to zero dipole moment of the acceptor will yield very large shifts in the hydrogen-bonded X-H stretching frequencies. The trends in the variation of the dispersion and exchange-repulsion components and the total interaction energy vis-à-vis frequency shifts of donor stretching vibration are similar for hydrogen-bonded complexes of phenylacetylene, indole, and phenol. Furthermore, it was observed that the vibrational frequency shifts of all of the complexes are linearly correlated with the charge transfer from the filled orbital of the hydrogen acceptor to the vacant antibonding (σ*) orbital of the X-H donor group on the basis of natural bonding orbital calculations.

Continuum of covalent to intermolecular bonding in the halogen-bonded complexes of 1,4-diazabicyclo[2.2.2]octane with bromine-containing electrophiles

The Br⋯N bonds in a series of halogen-bonded complexes change gradually from the typical intermolecular to the traditional covalent bond.