Synthesis of Functionalized Tetrasubstituted Allenes by the Addition of Bis(trimethylsilyl)ketene Acetals to Ynones Catalyzed by Gold(I).

We have developed a straightforward and rapid methodology for the synthesis of tetrasubstituted allenes bearing carboxylic acids in the 1,3-position through the gold(I)-catalyzed nucleophilic addition of bis(trimethylsilyl)ketene acetals to ynones. The reaction was evaluated with several substrates, and 21 allenes were obtained in moderate to good yields. Using DFT calculations, we studied the mechanism of the reaction, which suggested a nucleophilic 1,4-addition pathway. The potential of allenes to act as a source of highly functionalized lactones was also explored.

On the structure and electronic properties of Ptn clusters: new most stable structures for n = 16-17.

The lowest energy structures and electronic properties of Ptn clusters up to n = 17 are investigated by using a genetic algorithm in combination with density functional theory calculations. There are several putative global minimum structures for platinum clusters which have been reported by using different approaches, but a comprehensive study for n = 15-17 has not been carried out so far. Herein, we perform a consensus using GGA (PBE), meta-GGA (TPSS) and hybrid (B3PW91, PBE0, PBEh-3c, M06-L) functionals in conjunction with the Def2-TZVP basis set. New most stable structures are found for Pt16 and Pt17, which are slightly lower in energy than the previously reported global minima. Molecular dynamics simulations show that the clusters are rigid at room temperature. We analyze the structural, electronic, energy and vibrational data of the investigated clusters in detail.

Wave function analyses of scandium-doped aluminium clusters, AlnSc (n = 1-24), and their CO2 fixation abilities.

Nanoclusters represent a connection between (i) solid state systems and (ii) species in the atomic and molecular domains. Additionally, nanoclusters can also have very interesting electronic, optical and magnetic properties. For example, some aluminium clusters behave as superatoms and the doping of these clusters might strengthen their adsorption capabilities. Thus, we address herein the structural, energetic and electronic characterisation of scandium-doped aluminium clusters (AlnSc (n = 1-24)) by means of density functional theory calculations and quantum chemical topology wave function analyses. We studied the effect of Sc-doping on the structure and charge distribution by considering pure Al clusters as well. The quantum theory of atoms in molecules (QTAIM) reveals that interior Al atoms have large negative atomic charges (≈2a.u.) and hence the atoms surrounding them are considerably electron deficient. The Interacting Quantum Atoms (IQA) energy partition allowed us to establish the nature of the interaction between the Al13 superatom and the Al12Sc cluster with Al to form the complexes Al14 and Al13Sc, respectively. We also used the IQA approach to examine (i) the influence of Sc on the geometry of the AlnSc complexes along with (ii) the cooperative effects in the binding of AlnSc and Aln+1 clusters. We also exploited the QTAIM and IQA approaches to study the interaction of the electrophilic surface of the examined systems with CO2. Overall, we observe that the investigated Sc-doped Al complexes with a marked stability towards disproportionation reactions exhibit strong adsorption energies with CO2. Concomitantly, the carbon dioxide molecule is considerably distorted and destabilised, conditions which might prepare it for further chemical reactions. Altogether, this paper gives valuable insights on the tuning of the properties of metallic clusters for their design and exploitation in custom-made materials.

Atoms in molecules in real space: a fertile field for chemical bonding.

In this perspective, we review some recent advances in the concept of atoms-in-molecules from a real space perspective. We first introduce the general formalism of atomic weight factors that allows unifying the treatment of fuzzy and non-fuzzy decompositions under a common algebraic umbrella. We then show how the use of reduced density matrices and their cumulants allows partitioning any quantum mechanical observable into atomic or group contributions. This circumstance provides access to electron counting as well as energy partitioning, on the same footing. We focus on how the fluctuations of atomic populations, as measured by the statistical cumulants of the electron distribution functions, are related to general multi-center bonding descriptors. Then we turn our attention to the interacting quantum atom energy partitioning, which is briefly reviewed since several general accounts on it have already appeared in the literature. More attention is paid to recent applications to large systems. Finally, we consider how a common formalism to extract electron counts and energies can be used to establish an algebraic justification for the extensively used bond order-bond energy relationships. We also briefly review a path to recover one-electron functions from real space partitions. Although most of the applications considered will be restricted to real space atoms taken from the quantum theory of atoms in molecules, arguably the most successful of all the atomic partitions devised so far, all the take-home messages from this perspective are generalizable to any real space decompositions.

A QCT View of the Interplay between Hydrogen Bonds and Aromaticity in Small CHON Derivatives.

The somewhat elusive concept of aromaticity plays an undeniable role in the chemical narrative, often being considered the principal cause of the unusual properties and stability exhibited by certain π skeletons. More recently, the concept of aromaticity has also been utilised to explain the modulation of the strength of non-covalent interactions (NCIs), such as hydrogen bonding (HB), paving the way towards the in silico prediction and design of tailor-made interacting systems. In this work, we try to shed light on this area by exploiting real space techniques, such as the Quantum Theory of Atoms in Molecules (QTAIM), the Interacting Quantum Atoms (IQA) approaches along with the electron delocalisation indicators Aromatic Fluctuation (FLU) and Multicenter (MCI) indices. The QTAIM and IQA methods have been proven capable of providing an unbiased and rigorous picture of NCIs in a wide variety of scenarios, whereas the FLU and MCI descriptors have been successfully exploited in the study of diverse aromatic and antiaromatic systems. We used a collection of simple archetypal examples of aromatic, non-aromatic and antiaromatic moieties within organic molecules to examine the changes in π delocalisation and aromaticity induced by the Aromaticity and Antiaromaticity Modulated Hydrogen Bonds (AMHB). We observed fundamental differences in the behaviour of systems containing the HB acceptor within and outside the ring, e.g., a destabilisation of the rings in the former as opposed to a stabilisation of the latter upon the formation of the corresponding molecular clusters. The results of this work provide a physically sound basis to rationalise the strengthening and weakening of AMHBs with respect to suitable non-cyclic non-aromatic references. We also found significant differences in the chemical bonding scenarios of aromatic and antiaromatic systems in the formation of AMHB. Altogether, our investigation provide novel, valuable insights about the complex mutual influence between hydrogen bonds and π systems.

Binding Energy Partition of Promising IRAK-4 Inhibitor (Zimlovisertib) for the Treatment of COVID-19 Pneumonia.

The technique of Fragment-Based Drug Design (FBDD) considers the interactions of different moieties of molecules with biological targets for the rational construction of potential drugs. One basic assumption of FBDD is that the different functional groups of a ligand interact with a biological target in an approximately additive, that is, independent manner. We investigated the interactions of different fragments of ligands and Interleukin-1 Receptor-Associated Kinase 4 (IRAK-4) throughout the FBDD design of Zimlovisertib, a promising anti-inflammatory, currently in trials to be used for the treatment of COVID-19 pneumonia. We utilised state-of-the-art methods of wave function analyses mainly the Interacting Quantum Atoms (IQA) energy partition for this purpose. By means of IQA, we assessed the suitability of every change to the ligand in the five stages of FBDD which led to Zimlovisertib on a quantitative basis. We determined the energetics of the interaction of different functional groups in the ligands with the IRAK-4 protein target and thereby demonstrated the adequacy (or lack thereof) of the changes made across the design of this drug. This analysis permits to verify whether a given alteration of a prospective drug leads to the intended tuning of non-covalent interactions with its protein objective. Overall, we expect that the methods exploited in this paper will prove valuable in the understanding and control of chemical modifications across FBDD processes.

Partition of the electronic energy of the PM7 method via the interacting quantum atoms approach.

Partitions of the electronic energy such as that provided by the Interacting Quantum Atoms (IQA) approach have given valuable insights for numerous chemical systems and processes. Unfortunately, this kind of analysis may involve the integration of scalar fields over very irregular volumes, a condition which leads to a large and often prohibitive computational effort. These circumstances have limited the use of these energy partitions to systems comprising a few tens of atoms at most. On the other hand, semiempirical methods have proved useful in the study of systems of several thousands of atoms. Therefore, the goal of this work is to carry out partitions of the semiempirical method PM7 in compliance with the IQA approach. For this purpose, we computed one- and two-atomic energetic contributions whose sum equals the PM7 electronic energy. We illustrate how one might exploit the partition of electronic energies computed via the PM7 method by considering small organic and inorganic molecules and the energetics of individual hydrogen bond interactions within several water clusters which include (H2O)30, (H2O)50 and (H2O)100. We also considered the solvation of the amphiphilic caprylate anion to exemplify how to exploit the energy partition proposed in this paper. Overall, this investigation shows how the approach put forward herein might give further insights of the interactions occurring within complex systems in physical and biological chemistry.

Towards an atomistic understanding of polymorphism in molecular solids.

Understanding and controlling polymorphism in molecular solids is a major unsolved problem in crystal engineering. While the ability to calculate accurate lattice energies with atomistic modelling provides valuable insight into the associated energy scales, existing methods cannot connect energy differences to the delicate balances of intra- and intermolecular forces that ultimately determine polymorph stability ordering. We report herein a protocol for applying Quantum Chemical Topology (QCT) to study the key intra- and intermolecular interactions in molecular solids, which we use to compare the three known polymorphs of succinic acid including the recently-discovered γ form. QCT provides a rigorous partitioning of the total energy into contributions associated with topological atoms, and a quantitative and chemically intuitive description of the intra- and intermolecular interactions. The newly-proposed Relative Energy Gradient (REG) method ranks atomistic energy terms (steric, electrostatic and exchange) by their importance in constructing the total energy profile for a chemical process. We find that the conformation of the succinic acid molecule is governed by a balance of large and opposing electrostatic interactions, while the H-bond dimerisation is governed by a combination of electrostatics and sterics. In the solids, an atomistic energy balance emerges that governs the contraction, towards the equilibrium geometry, of a molecular cluster representing the bulk crystal. The protocol we put forward is as general as the capabilities of the underlying quantum-mechanical model and it can provide novel perspectives on polymorphism in a wide range of chemical systems.

Stronger-together: the cooperativity of aurophilic interactions.

Crystallographic distances and the electron density of bi- and tri-nuclear gold(I) compounds reveal that the existence of multiple Au⋯Au interactions increases their individual strength in the order of 0.9-2.9 kcal mol-1. We observed this behaviour both experimentally and theoretically in multinuclear systems, confirming a novel important cooperative character in aurophilic contacts.

Water clusters as bifunctional catalysts in organic chemistry: the hydrolysis of oxirane and its methyl derivatives.

This contribution explores the bifunctional catalytic activity of water clusters ((H2O)n with n = 1-5) in organic chemistry similar to that observed in the formation of H2SO4 in acid rain (Chem. Commun., 53, 3516, (2017)). We considered for this purpose the Hydrolysis of Epoxides (HE), in particular, that of oxirane and its methyl derivatives. Surrounding water molecules with H-bond cooperative effects decrease the activation energy of the rate-limiting step of HE in condensed phase, especially when they lead to an anti-periplanar attack on the alkoxide leaving group. Furthermore, the water molecules have a bifunctional catalytic role in HE by (i) increasing the nucleophilic and electrophilic character of the attacking oxygen atom and the leaving group of the reaction, respectively, and (ii) placing the reactants in a suitable disposition for the substitution to occur. Overall, this investigation provides relevant insights into the collective action of water molecules on organic reactions in neutral, basic and acid media.

On the Relationship between Hydrogen Bond Strength and the Formation Energy in Resonance-Assisted Hydrogen Bonds.

Resonance-assisted hydrogen bonds (RAHB) are intramolecular contacts that are characterised by being particularly energetic. This fact is often attributed to the delocalisation of π electrons in the system. In the present article, we assess this thesis via the examination of the effect of electron-withdrawing and electron-donating groups, namely -F, -Cl, -Br, -CF3, -N(CH3)2, -OCH3, -NHCOCH3 on the strength of the RAHB in malondialdehyde by using the Quantum Theory of Atoms in Molecules (QTAIM) and the Interacting Quantum Atoms (IQA) analyses. We show that the influence of the investigated substituents on the strength of the investigated RAHBs depends largely on its position within the π skeleton. We also examine the relationship between the formation energy of the RAHB and the hydrogen bond interaction energy as defined by the IQA method of wave function analysis. We demonstrate that these substituents can have different effects on the formation and interaction energies, casting doubts regarding the use of different parameters as indicators of the RAHB formation energies. Finally, we also demonstrate how the energy density can offer an estimation of the IQA interaction energy, and therefore of the HB strength, at a reduced computational cost for these important interactions. We expected that the results reported herein will provide a valuable understanding in the assessment of the energetics of RAHB and other intramolecular interactions.

The nature of the intermolecular interaction in (H2X)2 (X = O, S, Se).

Hydrogen bonds (HBs) are crucial non-covalent interactions in chemistry. Recently, the occurrence of an HB in (H2S)2 has been reported (Arunan et al., Angew. Chem., Int. Ed., 2018, 57, 15199), challenging the textbook view of H2S dimers as mere van der Waals clusters. We herein try to shed light on the nature of the intermolecular interactions in the H2O, H2S, and H2Se dimers via correlated electronic structure calculations, Symmetry Adapted Perturbation Theory (SAPT) and Quantum Chemical Topology (QCT). Although (H2S)2 and (H2Se)2 meet some of the criteria for the occurrence of an HB, potential energy curves as well as SAPT and QCT analyses indicate that the nature of the interaction in (H2O)2 is substantially different (e.g. more anisotropic) from that in (H2S)2 and (H2Se)2. QCT reveals that the HB in (H2O)2 includes substantial covalent, dispersion and electrostatic contributions, while the last-mentioned component plays only a minor role in (H2S)2 and (H2Se)2. The major contributions to the interactions of the dimers of H2S and H2Se are covalency and dispersion as revealed by the exchange-correlation components of QCT energy partitions. The picture yielded by SAPT is somewhat different but compatible with that offered by QCT. Overall, our results indicate that neither (H2S)2 nor (H2Se)2 are hydrogen-bonded systems, showing how the nature of intermolecular contacts involving hydrogen atoms evolves in a group down the periodic table.

Cooperativity and Anticooperativity in Ion-Water Interactions: Implications for the Aqueous Solvation of Ions.

Non-additive effects in hydrogen bonds (HB) take place as a consequence of electronic charge transfers. Therefore, it is natural to expect cooperativity and anticooperativity in ion-water interactions. Nevertheless, investigations on this matter are scarce. This paper addresses the interactions of (i) the cations Li+ , Na+ , K+ , Be2+ , Mg2+ , and Ca2+ together with (ii) the anions F- , Cl- , Br- , NO3- and SO42- with water clusters (H2 O)n , n=1-8, and the effects of these ions on the HBs within the complete molecular adducts. We used quantum chemical topology tools, specifically the quantum theory of atoms in molecules and the interacting quantum atoms energy partition to investigate non-additive effects among the interactions studied herein. Our results show a decrease on the interaction energy between ions and the first neighbouring water molecules with an increment of the coordination number. We also found strong cooperative effects in the interplay between HBs and ion-dipole interactions within the studied systems. Such cooperativity affects considerably the interactions among ions with their first and second solvation shells in aqueous environments. Overall, we believe this article provides valuable information about how ion-dipole contacts interact with each other and how they relate to other interactions, such as HBs, in the framework of non-additive effects in aqueous media.

Photochemistry in Real Space: Batho- and Hypsochromism in the Water Dimer.

Invited for the cover of this issue is Alberto Fernández-Alarcón and co-workers at The Institute of Chemistry of the National Autonomous University of Mexico and The School of Chemistry of the University of Oviedo. The image depicts the real space analysis of the excitation energies in the double blue and red shift of the water dimer. Read the full text of the article at 10.1002/chem.202002854.

Interacting Quantum Atoms-A Review.

The aim of this review is threefold. On the one hand, we intend it to serve as a gentle introduction to the Interacting Quantum Atoms (IQA) methodology for those unfamiliar with it. Second, we expect it to act as an up-to-date reference of recent developments related to IQA. Finally, we want it to highlight a non-exhaustive, yet representative set of showcase examples about how to use IQA to shed light in different chemical problems. To accomplish this, we start by providing a brief context to justify the development of IQA as a real space alternative to other existent energy partition schemes of the non-relativistic energy of molecules. We then introduce a self-contained algebraic derivation of the methodological IQA ecosystem as well as an overview of how these formulations vary with the level of theory employed to obtain the molecular wavefunction upon which the IQA procedure relies. Finally, we review the several applications of IQA as examined by different research groups worldwide to investigate a wide variety of chemical problems.

On the strength of hydrogen bonding within water clusters on the coordination limit.

Hydrogen bonds (HB) are arguably the most important noncovalent interactions in chemistry. We study herein how differences in connectivity alter the strength of HBs within water clusters of different sizes. We used for this purpose the interacting quantum atoms energy partition, which allows for the quantification of HB formation energies within a molecular cluster. We could expand our previously reported hierarchy of HB strength in these systems (Phys. Chem. Chem. Phys., 2016, 18, 19557) to include tetracoordinated monomers. Surprisingly, the HBs between tetracoordinated water molecules are not the strongest HBs despite the widespread occurrence of these motifs (e.g., in ice Ih ). The strongest HBs within H2 O clusters involve tricoordinated monomers. Nonetheless, HB tetracoordination is preferred in large water clusters because (a) it reduces HB anticooperativity associated with double HB donors and acceptors and (b) it results in a larger number of favorable interactions in the system. Finally, we also discuss (a) the importance of exchange-correlation to discriminate among the different examined types of HBs within H2 O clusters, (b) the use of the above-mentioned scale to quickly assess the relative stability of different isomers of a given water cluster, and (c) how the findings of this research can be exploited to indagate about the formation of polymorphs in crystallography. Overall, we expect that this investigation will provide valuable insights into the subtle interplay of tri- and tetracoordination in HB donors and acceptors as well as the ensuing interaction energies within H2 O clusters.

Photochemistry in Real Space: Batho- and Hypsochromism in the Water Dimer.

The development of chemical intuition in photochemistry faces several difficulties that result from the inadequacy of the one-particle picture, the Born-Oppenheimer approximation, and other basic ideas used to build models. It is shown herein how real-space approaches can be efficiently used to gain valuable insights in photochemistry through a simple example of red and blue shift effects: the double hypso- and bathochromic shifts in the low-lying valence excited states of (H2 O)2 . It is demonstrated that 1) the use of these techniques allows the perturbative language used in the theory of intermolecular interactions, even in the strongly interacting short-range regime, to be maintained; 2) one and only one molecule is photoexcited in each of the addressed excited states and 3) the electrostatic interaction between the in-the-cluster molecular dipoles provides a fairly intuitive rationalisation of the observed batho- and hypsochromism. The methods exploited and illustrated herein are able to maintain the individuality and properties of the interacting entities in a molecular aggregate, and thereby they allow chemical intuition in general states, at any geometry and using a broad variety of electronic structure methods to be kept and built.

Directing the Crystal Packing in Triphenylphosphine Gold(I) Thiolates by Ligand Fluorination.

We explore herein the supramolecular interactions that control the crystalline packing in a series of fluorothiolate triphenylphosphine gold(I) compounds with the general formula [Au(SRF)(Ph3P)] in which Ph3P = triphenylphosphine and SRF = SC6F5, SC6HF4-4, SC6F4(CF3)-4, SC6H3F2-2,4, SC6H3F2-3,4, SC6H3F2-3,5, SC6H4(CF3)-2, SC6H4F-2, SC6H4F-3, SC6H4F-4, SCF3, and SCH2CF3. We use for this purpose (i) DFT electronic structure calculations and (ii) the quantum theory of atoms in molecules and the non-covalent interactions index methods of wave function analyses. Our combined experimental and computational approach yields a general understanding of the effects of ligand fluorination in the crystalline self-assembly of the examined systems, in particular, about the relative force of aurophilic contacts compared with other supramolecular interactions. We expect this information to be useful in the design of materials based on gold coordination compounds.

Efficient implementation of the interacting quantum atoms energy partition of the second-order Møller-Plesset energy.

We describe an efficient implementation of the partition of the second-order Møller-Plesset (MP2) correlation energy within the interacting quantum atoms (IQA) energy decomposition. We simplify the IQA integration bottleneck by considering only the occupied to virtual elements of the second order reduced density matrix, a procedure that reduces substantially the size of the two-electron matrix, which has to be addressed. The algorithmic improvements described herein allow to perform the decomposition of the MP2 correlation energy for medium size molecular systems using moderate computational resources. We expect that the methods developed in this investigation will prove useful to understand electron correlation effects through a real space perspective.

The effect of chiral N-substituents with methyl or trifluoromethyl groups on the catalytic performance of mono- and bifunctional thioureas.

We evaluated thiourea organocatalysts that incorporate a chiral group which includes a trifluoromethyl moiety and contrasted their performance with non-fluorinated analogs. The comparison between such systems allows the direct study of the NH acidity of a thiourea bonded to an aliphatic substituent. In principle, -CF3 systems feature an enhanced hydrogen bond (HB) donor capacity that is undoubtedly beneficial for HB-catalysis applied to the Baylis-Hillman reaction. We found that the thiourea substituted on both nitrogens with this group accelerates this reaction like Schreiner's thiourea. On the other hand, we observed a different behavior in reactions promoted by bifunctional catalysts (thiourea-primary amine). In the Michael addition of isobutyraldehyde to methyl benzylidenepyruvate, the -CF3 containing catalysts were better than the -CH3 systems, whereas the conjugate addition to N-phenylmaleimide showed the opposite behavior. Theoretical calculations of the transition states indicated that the phenylethyl group in fluorinated and non-fluorinated compounds have different kinds of interactions with the electrophile. These interactions are responsible for a different arrangement of the electrophile and thereby the selectivity of the catalyst. Therefore, it cannot be generalized that in all cases NH acidity correlates with the performance of the catalyst, particularly, with aliphatic substituents that unlike the aromatic ones possess groups that are outside the plane of the thiourea.