Search

Our quantum chemical analyses elucidated how the replacement of O in the amide bonds of benzene-1,3-5-tricarboxamides (OBTAs) with the larger chalcogens S and Se enhances the intermolecular interactions and thereby the stability of the obtained hydrogen-bonded supramolecular polymers due to two unexpected reasons: i) the SBTA and SeBTA monomers have a better geometry for self-assembly and ii) induce stronger covalent (hydrogen-bond) interactions besides enhanced dispersion interactions. In addition, it is shown that the cooperativity in benzene-1,3,5-triamide (BTA) self-assembly is caused by charge separation in the σ-electronic system following the covalency in the hydrogen bonds.

The redox-conversion reaction of cobalt(II)-disulfide to cobalt(III)-thiolate complexes triggered by addition of the bidentate ligand 2,2'-bipyridine has been investigated. Reaction of the cobalt(II)-disulfide complex [Co2(L1SSL1)(X)4] (L1SSL1 = di-2-(bis(2-pyridylmethyl)amino)-ethyldisulfide; X = Cl or Br) [1X] with 2,2'-bipyridine (bpy) resulted in the formation of two different products, namely the cobalt(III)-thiolate complex [Co(L1S)(bpy)]X2 and the unexpected side product [Co2(L1SSL1)(bpy)2(X)2]X2. Crystals of [Co2(L1SSL1)(bpy)2(Cl)2](BPh4)2 [2Cl](BPh4)2 obtained after anion exchange showed the cobalt(II) ions to be in octahedral geometries with the nitrogen donors of the disulfide ligand arranged in a facial conformation and the chloride ion trans to the tertiary amine nitrogen. Remarkably, this side product cannot be converted to the cobalt(III)-thiolate compound [Co(L1S)(bpy)](SbF6)2 [3](SbF6)2 by removal of the chloride ion with use of a silver salt, as this causes scrambling of the ligands, resulting in the formation of [Co(bpy)3]n+. [Co(L1S)(bpy)](SbF6)2 was obtained in a pure form by addition of bpy to a solution in acetonitrile of the compound [Co(L1S)(MeCN)2]2+ [4]2+. Addition of NEt4Cl to [3](SbF6)2 regenerates the cobalt(II)-disulfide complex [1Cl] as confirmed spectroscopically. DFT studies revealed that the conversion from [1Cl] to [3]2+ most likely occurs via the hypothetical intermediate species [2Cl]2+mer, in which the nitrogen atoms of the disulfide ligand are arranged in a meridional conformation. Interestingly, the estimated d-orbital splitting energy of [3]2+ is lower than that of [4]2+, indicating that the ligand-field strength of bpy is lower than anticipated, which hampers clean conversion in the redox-conversion reaction. This study shows that the redox-conversion reaction between cobalt(II)-disulfide and cobalt(III)-thiolate complexes is intricate rather than straightforward.

We discuss the fundamental aspects of the EDA-NOCV method and address some critical comments that have been made recently. The EDA-NOCV method unlike most other methods focuses on the process of bond formation between the interacting species and not just only on the analysis of the finally formed bond. This is demonstrated using LiF as an example. There is a difference between the interactions between the initial species which form the bond and are also the final product of bond cleavage, and the interactions between the fragments in the eventually formed molecule. The flexibility of the method allows the choice of the interacting fragments which helps to identify the charge and electron configuration of the fragments which describe the bond. This is very helpful in cases where the bond may be described with several Lewis structures. We reject the idea that it would be a disadvantage to have “bond path functions” as the energy components in the EDA, which actually indicate the variability of the method. The bonding analysis in a different sequence of the bond formation gives important results for the various questions that can be asked. This is demonstrated by using CH2, CO2 and the formation of a guanine quartet as examples. The fact that a bond is always defined by the bound molecule, the fragments, and their states is universal and deeply physical, as we show here again for various examples. The results of the EDA-NOCV method are in full accordance with the physical mechanism of the chemical bond as revealed by Ruedenberg.

Hydrogen-bonded supramolecular systems are known to obtain extra stabilization from the complexation with ions, like guanine quadruplex (GQ). They experience strong hydrogen bonds due to cooperative effects. To gain deeper understanding of the interplay between ions and hydrogen-bonding cooperativity, relativistic dispersion-corrected density functional theory (DFT-D) computations were performed on triple-layer hydrogen-bonded rosettes of ammeline interacting with alkali metal cations and halides. Our results show that when ions are placed between the stacks, the hydrogen bonds are weakened but, at the same time, the cooperativity is strengthened. This phenomenon can be traced back to the shrinkage of the cavity as the ions pull the monomers closer together and therefore the distance between the monomers becomes smaller. On one hand this results in a larger steric repulsion, but on the other hand, the donor-acceptor interactions are enhanced due to the larger overlap between the donating and accepting orbitals leading to more charge donation and therefore an enhanced electrostatic attraction.

The synthesis of the selenium-based ligand L1SeSeL1 (2,2’-diselanediylbis(N,N-bis(pyridin-2-ylmethyl)ethan-1-amine) is described along with its reactivity with cobalt(II) salts. The cobalt(II)-diselenide complex [Co2(L1SeSeL1)Cl4] was obtained in good yield, and its spectroscopic properties closely resemble that of its sulfur analog. Reaction of L1SeSeL1 with Co(II) thiocyanate results in the formation of the cobalt(III) compound [Co(L1Se)(NCS)2], similar to reaction of L1SSL1. The redox-conversion reactions from the Co(II)-diselenide compound [Co(L1SeSeL1)Cl4] using external triggers such as removal of the halide ions or the addition of the strong-field ligand 8-quinolinolate resulted in good yields of the Co(III)-selenolate complexes [Co(L1Se)(MeCN)2](SbF6)2 and [Co(L1Se)(quin)]Cl. Our computational studies show that the ligand-field splitting energy of the selenium compounds is smaller than their sulfur analogs, indicating that redox-conversion of cobalt(II)-diselenide to cobalt(III)-selenolate complexes may be more arduous than that for the related sulfur compounds.

The gauche conformer in 1-X,2-Y-disubstituted ethanes, that is, the staggered orientation in which X and Y are in closer contact, is only favored for relatively small substituents that do not give rise to large X···Y steric repulsion. For more diffuse substituents, weakly attractive orbital interactions between antiperiplanar bonds (i.e., hyperconjugation) cannot overrule the repulsive forces between X and Y. Our quantum chemical analyses of the rotational isomerism of XCH2CH2Y (X = F, OH; Y = I) at ZORA-BP86-D3(BJ)/QZ4P reveal that indeed the anti conformer is generally favored due to a less destabilizing I···F and I···O-H steric repulsion. The only case when the gauche conformer is preferred is when the hydroxyl hydrogen is oriented toward the iodine atom in the 2-iodoethanol. This is because of the significantly stabilizing covalent component of the I···H-O intramolecular hydrogen bond. Therefore, we show that strong intramolecular interactions can overcome the steric repulsion between bulky substituents in 1,2-disubstituted ethanes and cause the gauche effect. Our quantum chemical computations have guided nuclear magnetic resonance experiments that confirm the increase in the gauche population as X goes from F to OH.

The Tc-99m nucleus is the most used nuclide in radiopharmaceuticals designed for imaging diagnosis. The metal can exist in nine distinct oxidation states and forms distinct coordination complexes with a variety of chelating agents and geometries. These complexes are usually characterized through Tc-99 NMR that is very sensitive to the Tc coordination sphere. Therefore, predicting Tc-99 NMR might be useful to assist experimentalists in structural characterization. In the present study, we propose three computational protocols for predicting Tc-99 NMR chemical shifts based on density functional theory calculations using relativistic and nonrelativistic Hamiltonians: the relativistic Model 1, the nonrelativistic Model 2, and the empirical nonrelativistic Model 3. In Models 2 and 3, the NMR-DKH basis set was used for all atoms, including the Tc, for which it was developed here. All models were applied for a set of 41 Tc-complexes with metal oxidation states 0, I, and V, for which the Tc-99 chemical shift was available experimentally. The mean absolute deviation and the mean relative deviation were 67 ppm and 4.8% (Model 1), 92 ppm and 6.2% (Model 2), and 65 ppm and 4.9% (Model 3), respectively. Last, the effect of the explicit solvent was evaluated for the [TcO2(en)2]+─Tc(V) complex. The calculated results for the Tc-99 NMR chemical shift at SO-ZORA-SSB-D/TZ2P-ZORA/COSMO//TPSS/def2-SVP/IEF-PCM(UFF) show that the inclusion of 14 water molecules (first solvation shell) together with the implicit solvation model leads to an absolute deviation of only 7 ppm (0.3%) from the experimental value, indicating that the solvent effects play a key role in predicting Tc-99 NMR.

We have investigated 44 crystal structures, found in the Cambridge Structural Database, containing the X3 synthon (where X = Cl, Br, I) in order to verify whether three type II halogen-halogen contacts forming the synthon exhibit cooperativity. A hypothesis that this triangular halogen-bonded motif is stabilized by cooperative effects is postulated on the basis of structural data. However, theoretical investigations of simplified model systems in which the X3 motif is present demonstrate that weak synergy occurs only in the case of the I3 motif. In the present paper we computationally investigate crystal structures in which the X3 synthon is present, including halomesitylene structures, that are usually described as being additionally stabilized by a synergic interaction. Our computations find no cooperativity for halomesitylene trimers containing the X3 motif. Only in the case of two other structures containing the I3 synthon a very weak or weak synergy, i.e. the cooperative effect being stronger than -0.40 kcal mol-1, is found. The crystal structure of iodoform has the most pronounced cooperativity of all investigated systems, amounting to about 10% of the total interaction energy.

The amino groups of thio- and selenoamides can act as stronger hydrogen-bond donors than of carboxamides, despite the lower electronegativity of S and Se. This phenomenon has been experimentally explored, particularly in organocatalysis, but a sound electronic explanation is lacking. Our quantum chemical investigations show that the NH2 groups in thio- and selenoamides are more positively charged than in carboxamides. This originates from the larger electronic density flow from the nitrogen lone pair of the NH2 group towards the lower-lying π*C=S and π*C=Se orbitals than to the high-lying π*C=O orbital. The relative energies of the π* orbitals result from the overlap between the chalcogen np and carbon 2p atomic orbitals, which is set by the carbon-chalcogen equilibrium distance, a consequence of the Pauli repulsion between the two bonded atoms. Thus, neither the electronegativity nor the often-suggested polarizability but the steric size of the chalcogen atom determines the amide's hydrogen-bond donor capability.

The design of molecular-level communication tools engages the field of nanobiotechnology. In this issue of Chem, Clayden and co-workers report a novel mechanism of spatial molecular communication through a single strand of hydrogen-bonded oligomers. These “minimalistic” polarity switches distinguish themselves from conventional strategies by not employing chiral inversion.

We have quantum chemically investigated the bonding between archetypical Lewis acids and bases. Our state‐of‐the‐art computations on the X3B−NY3 Lewis pairs have revealed the origin behind the systematic increase in B−N bond strength as X and Y are varied from F to Cl, Br, I, H. For H3B−NY3, the bonding trend is driven by the commonly accepted mechanism of donor−acceptor [HOMO(base)−LUMO(acid)] interaction. Interestingly, for X3B−NH3, the bonding mechanism is determined by the energy required to deform the BX3 to the pyramidal geometry it adopts in the adduct. Thus, Lewis acids that can more easily pyramidalize form stronger bonds with Lewis bases. The decrease in the strain energy of pyramidalization on going from BF3 to BI3 is directly caused by the weakening of the B−X bond strength, which stems primarily from the bonding in the plane of the molecule (σ‐like) and not in the π system, at variance with the currently accepted mechanism., That's strainful! The stability of Lewis pairs involving boron trihalides and ammonia is driven by the energy required to deform the Lewis acids upon complexation, while the stability of Lewis pairs between borane and nitrogen trihalides is determined by the strength of orbital interactions.

The hydroaminomethylation (HAM) reaction converts alkenes into N-alkylated amines and has been well studied for rhodium- and ruthenium-based catalytic systems. Cobalt-based catalytic systems are able to perform the essential hydroformylation reaction, but are also known to form very active hydrogenation catalysts, therefore we examined such a system for its potential use in the HAM reaction. Thus, we have quantum-chemically explored the hydrogenation activity of [HCo(CO)3] in model reactions with ethene, methyleneamine, formaldehyde, and vinylamine using dispersion-corrected relativistic density functional theory at ZORA-BLYP-D3(BJ)/TZ2P. Our computations reveal essentially identical overall barriers for the catalytic hydrogenation of ethene, formaldehyde, and vinylamine. This strongly suggests that a cobalt-based catalytic system will lack hydrogenation selectivity in experimental HAM reactions. Our HAM experiments with a cobalt-based catalytic system (consisting of Co2(CO)8 as cobalt source and P(n-Bu)3 as ligand) resulted in the formation of the desired N-alkylated amine. However, significant amounts of hydrogenated starting material as well as alcohol (hydrogenated aldehyde) were always formed. The use of cobalt-based catalysts in the HAM reaction to selectively form N-alkylated amines seems therefore not feasible. This confirms our computational prediction and highlights the usefulness of state-of-the-art DFT computations for guiding future experiments.

For many years non-covalently bonded complexes of nucleobases have attracted considerable interest. However, there is a lack of information about the nature of hydrogen bonding between nucleobases when the bonding is affected by metal coordination to one of the nucleobases, and how the individual hydrogen bonds and aromaticity of nucleobases respond to the presence of the metal cation. Here we report a DFT computational study of nucleobase pairs interacting with alkali metal cations. The metal cations contribute to the stabilization of the base pairs to varying degrees depending on their position. The energy decomposition analysis revealed that the nature of bonding between nucleobases does not change much upon metal coordination. The effect of the cations on individual hydrogen bonds were described by changes in VDD charges on frontier atoms, H bond length, bond energy from NBO analysis, and delocalization index from QTAIM calculations. The aromaticity changes were determined by HOMA index M. Solà and O.A.S. are grateful to the Ministerio de Economía y Competitividad (MINECO) of Spain (project CTQ2017-85341-P and Juan de la Cierva formación contract FJCI‐2017‐ 32757 to O.A.S.) and the Generalitat de Catalunya (project 2017SGR39). HS and TMK thank the National Science Centre of Poland for supporting this work under the grant no. UMO2016/23/B/ST4/00082. M. Swart acknowledges MICINN/MINECO (projects CTQ2014 59212 P, CTQ2015-70851-ERC, CTQ2017-87392-P), GenCat (2014SGR1202, 2017SGR1434 and XRQTC network) and European Fund for Regional Development (FEDER, UNGI104E801). C.F.G. thanks the Netherlands Organization for Scientific Research (NWO) for financial support

We have evaluated the performance of various density functionals, covering generalized gradient approximation (GGA), global hybrid (GH) and range‐separated hybrid (RSH), using time dependent density functional theory (TDDFT) for computing vertical excitation energies against experimental absorption maximum (λmax) for a set of 10 different core‐substituted naphthalene diimides (cNDI) recorded in dichloromethane. The computed excitation in case of GH PBE0 is most accurate while the trend is most systematic with RSH LCY‐BLYP compared to λmax. We highlight the importance of including solvent effects for optimal agreement with the λmax. Increasing the basis set size from TZ2P to QZ4P has a negligible influence on the computed excitation energies. Notably, RSH CAMY‐B3LYP gave the least error for charge‐transfer excitation. The poorest agreement with λmax is obtained with semi‐local GGA functionals. Use of the optimally‐tuned RSH LCY‐BLYP* is not recommended because of the high computational cost and marginal improvement in results., We have evaluated the performance of various exchange‐correlation functionals for reproducing the experimental absorption maximum in core‐substituted naphthalene diimides through quantum‐chemical TDDFT and statistical calculations. Results are analyzed and an appropriate TDDFT protocol has been developed.

Herein, we have investigated the effect of an endocyclic group (forming the N–C–C–F fragment) on the conformational preferences of 2-fluorocyclohexanone analogs. A combined approach of nuclear magnetic resonance and density functional theory calculations was employed to assess the conformational equilibrium in several media. In turn, natural bond orbital analysis and the conformational behavior of other 2-halocyclohexanone analogs were used to get more insights about the intramolecular interactions governing the conformer stabilities. Our results reveal that any stabilization from interactions featured in the gauche effect is overcome by a short-range interaction of the fluorine substituent with the carbonyl group. Consequently, the gauche effect in heterocyclic compounds is not as stabilizing as in their acyclic counterparts. Only the electrostatic gauche effect takes place even in polar solvents owing to an attraction between the axial fluorine and an endocyclic quaternary ammonium group.

Ion recognition is still an emerging topic in supramolecular chemistry and has aroused great attention in the last few years. In this work, we have examined the assemblies of selected hexameric rosettes of melamine and ammeline and their capacities to host halide and alkali ions in the gas phase and in water. Using relativistic dispersion-corrected density functional theory (DFT-D), we first studied the stability and the effect of introducing monovalent anions (Cl-, Br-, and I-) and cations (Na+, K+, and Rb+) in the center of the rosette´s cavity. Finally, we explored the interactions in two stacked rosettes with an interlayer ion. Our computations reveal that amine-substituted triazines are promising candidates for anion and cation recognition either in self-assembled monolayers or pillar array structures. The anion recognition process is governed by both the electrostatic and charge-transfer (donor-acceptor) interactions, while the cation recognition is governed by electrostatic and polarization. In addition, melamine and ammeline could constitute a potent mixture for dual-ion recognition strategies. Fil: Petelski, Andre Nicolai. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Nordeste; Argentina. Vrije Universiteit Amsterdam; Países Bajos. Universidad Tecnológica Nacional. Facultad Regional Resistencia. Departamento de Ingeniería Química. Laboratorio de Química Teórica y Experimental; Argentina Fil: Fonseca Guerra, Célia. Vrije Universiteit Amsterdam; Países Bajos

When carbonyl ligands coordinate to transition metals, their bond distance either increases (classical) or decreases (nonclassical) with respect to the bond length in the isolated CO molecule. C−O expansion can easily be understood by π‐back‐donation, which results in a population of the CO's π*‐antibonding orbital and hence a weakening of its bond. Nonclassical carbonyl ligands are less straightforward to explain, and their nature is still subject of an ongoing debate. In this work, we studied five isoelectronic octahedral complexes, namely Fe(CO)6 2+, Mn(CO)6 +, Cr(CO)6, V(CO)6 − and Ti(CO)6 2−, at the ZORA‐BLYP/TZ2P level of theory to explain this nonclassical behavior in the framework of Kohn–Sham molecular orbital theory. We show that there are two competing forces that affect the C−O bond length, namely electrostatic interactions (favoring C−O contraction) and π‐back‐donation (favoring C−O expansion). It is a balance between those two terms that determines whether the carbonyl is classical or nonclassical. By further decomposing the electrostatic interaction ΔV elstat into four fundamental terms, we are able to rationalize why ΔV elstat gives rise to the nonclassical behavior, leading to new insights into the driving forces behind C−O contraction., Classical or nonclassical? We show that there are two competing forces that affect the CO bond length in transition metal complexes, namely electrostatic interactions (favoring CO contraction) and π‐back‐donation (favoring CO expansion). The origin is elucidated by quantitative MO‐theory, leading to new insights into the driving forces behind nonclassical behavior.

It has been experimentally observed that water-ice-embedded polycyclic aromatic hydrocarbons (PAHs) form radical cations when exposed to vacuum UV irradiation, whereas ammonia-embedded PAHs lead to the formation of radical anions. In this study, we explain this phenomenon by investigating the fundamental electronic differences between water and ammonia, the implications of these differences on the PAH-water and PAH-ammonia interaction, and the possible ionization pathways in these complexes using density functional theory (DFT) computations. In the framework of the Kohn-Sham molecular orbital (MO) theory, we show that the ionic state of the PAH photoproducts results from the degree of occupied-occupied MO mixing between the PAHs and the matrix molecules. When interacting with the PAH, the lone pair-type highest occupied molecular orbital (HOMO) of water has poor orbital overlap and is too low in energy to mix with the filled π-orbitals of the PAH. As the lone-pair HOMO of ammonia is significantly higher in energy and has better overlap with filled π-orbitals of the PAH, the subsequent Pauli repulsion leads to mixed MOs with both PAH and ammonia character. By time-dependent DFT calculations, we demonstrate that the formation of mixed PAH-ammonia MOs opens alternative charge-transfer excitation pathways as now electronic density from ammonia can be transferred to unoccupied PAH levels, yielding anionic PAHs. As this pathway is much less available for water-embedded PAHs, charge transfer mainly occurs from localized PAH MOs to mixed PAH-water virtual levels, leading to cationic PAHs.

The formation of Co(III)-thiolate complexes from Co(II)-disulfide complexes using the anionic ligand 8-quinolinolate (quin-) has been studied experimentally and quantum chemically. Two Co(II)-disulfide complexes [Co2(LxSSLx)(Cl)4] (x = 1 or 2; L1SSL1 = 2,2'-disulfanediylbis(N,N-bis(pyridin-2-ylmethyl)ethan-1-amine; L2SSL2 = 2,2'-disulfanedylbis (N-((6-methylpyridin-2-yl)methyl)-N-(pyridin-2-ylmethyl) ethan-1-amine) have been successfully converted with high yield to their corresponding Co(III)-thiolate complexes upon addition of the ligand 8-quinolinolate. Using density functional theory (DFT) computations the d-orbital splitting energies of the cobalt-thiolate compounds [Co(L1S)(quin)]+ and [Co(L2S)(quin)]+ were estimated to be 3.10 eV and 3.07 eV, indicating a slightly smaller ligand-field strength of ligand L2SSL2 than of L1SSL1. Furthermore, the orientation of the quin- ligand in the thiolate compounds determines the stability of the thiolate complex. DFT computations show that the thiolate structure benefits from more electrostatic attraction when the oxygen atom of the quin- ligand is positioned trans to the sulfur atom of the [Co(L1S)]2+ fragment. Quin- is the first auxiliary ligand with which it appeared possible to induce the redox-conversion reaction in cobalt(II) compounds of the relatively weak-field ligand L2SSL2.

We have quantum chemically analyzed the cooperative effects and structural deformations of hydrogen-bonded urea, deltamide, and squaramide linear chains using dispersion-corrected density functional theory at BLYP-D3(BJ)/TZ2P level of theory. Our purpose is twofold: (i) reveal the bonding mechanism of the studied systems that lead to their self-assembly in linear chains; and (ii) rationalize the C-C bond equalization in the ring moieties of deltamide and squaramide upon polymerization. Our energy decomposition and Kohn-Sham molecular orbital analyses reveal cooperativity in all studied systems, stemming from the charge separation within the sigma-electronic system by charge transfer from the carbonyl oxygen lone pair donor orbital of one monomer towards the sigma* N-H antibonding acceptor orbital of the neighboring monomer. This key orbital interaction causes the C=O bonds to elongate, which, in turn, results in the contraction of the adjacent C-C single bonds that, ultimately, makes the ring moieties of deltamide and squaramide to become more regular. Notably, the pi-electron delocalization plays a much smaller role in the total interaction between the monomers in the chain.

Some nitrile-boron halide adducts exhibit a double-well potential energy surface with two distinct minima: a “long bond” geometry (LB, a van der Waals interaction mostly based on electrostatics, but including a residual charge transfer component) and a “short bond” structure (SB, a covalent dative bond). This behavior can be considered as a “weak” form of bond stretch isomerism. Our computations reveal that complexes RCN−BX3 (R=CH3, FCH2, BrCH2, and X=Cl, Br) exhibit a fast interconversion from LB to SB geometries even close to the absolute zero thanks to a boron atom tunneling mechanism. The computed half-lives of the meta-stable LB compounds vary between minutes to nanoseconds at cryogenic conditions. Accordingly, we predict that the long bond structures are practically impossible to isolate or characterize, which agrees with previous matrix-isolation experiments.

The formation of guanine quadruplexes (GQ) in DNA is crucial in telomere homeostasis and regulation of gene expression. Pollution metals can interfere with these DNA superstructures upon coordination. In this work, we study the affinity of the internal GQ channel site towards alkaline earth metal (Mg2+, Ca2+, Sr2+, and Ba2+), and (post-)transition metal (Zn2+, Cd2+, Hg2+, and Pb2+) cations using density functional theory computations. We find that divalent cations generally bind to the GQ cavity with a higher affinity than conventional monovalent cations (e. g. K+). Importantly, we establish the nature of the cation-GQ interaction and highlight the relationship between ionic and nuclear charge, and the electrostatic and covalent interactions. The covalent interaction strength plays an important role in the cation affinity and can be traced back to the relative stabilization of cations’ unoccupied atomic orbitals. Overall, our findings contribute to a deeper understanding of how pollution metals could induce genomic instability.

We have carried out an extensive quantum chemical exploration of gas-phase alkali metal cation affinities (AMCAs) of archetypal neutral bases across the periodic system using relativistic density functional theory. One objective of this work is to provide an intrinsically consistent set of values of the 298 K AMCAs of all neutral maingroup-element hydrides XHn of groups 15-18 along the periods 1-6. Our main purpose is to understand these trends in terms of the underlying bonding mechanism using Kohn-Sham molecular orbital theory together with a canonical energy decomposition analysis (EDA). We compare the trends in XHn AMCAs with the trends in XHn proton affinities (PAs). We also examine the differences between the trends in AMCAs of the neutral XHn bases with those in the corresponding anionic XHn-1 - bases. Furthermore, we analyze how the cation affinity of our neutral Lewis bases changes along the group-1 cations H+, Li+, Na+, K+, Rb+, and Cs+.

Direct synthetic routes to amidines are desired, as they are widely present in many biologically active compounds and organometallic complexes. N-Acyl amidines in particular can be used as a starting material for the synthesis of heterocycles and have several other applications. Here, we describe a fast and practical copper-catalyzed three-component reaction of aryl acetylenes, amines, and easily accessible 1,4,2-dioxazol-5-ones to N-acyl amidines, generating CO2 as the only byproduct. Transformation of the dioxazolones on the Cu catalyst generates acyl nitrenes that rapidly insert into the copper acetylide Cu-C bond rather than undergoing an undesired Curtius rearrangement. For nonaromatic dioxazolones, [Cu(OAc)(Xantphos)] is a superior catalyst for this transformation, leading to full substrate conversion within 10 min. For the direct synthesis of N-benzoyl amidine derivatives from aromatic dioxazolones, [Cu(OAc)(Xantphos)] proved to be inactive, but moderate to good yields were obtained when using simple copper(I) iodide (CuI) as the catalyst. Mechanistic studies revealed the aerobic instability of one of the intermediates at low catalyst loadings, but the reaction could still be performed in air for most substrates when using catalyst loadings of 5 mol The herein reported procedure not only provides a new, practical, and direct route to N-acyl amidines but also represents a new type of C-N bond formation.

The secondary electrostatic interaction (SEI) model is often used to predict and explain relative hydrogen bond strengths of self-assembled systems. The SEI model oversimplifies the hydrogen-bonding mechanisms by viewing them as interacting point charges, but nevertheless experimental binding strengths are often in line with the model's predictions. To understand how this rudimentary model can be predictive, we computationally studied two tautomeric quadruple hydrogen-bonded systems, DDAA-AADD and DADA-ADAD. Our results reveal that when the proton donors D (which are electron-donating) and the proton acceptors A (which are electron-withdrawing) are grouped together as in DDAA, there is a larger accumulation of charge around the frontier atoms than when the proton donor and acceptor groups are alternating as in DADA. This accumulation of charge makes the proton donors more positive and the proton acceptors more negative, which enhances both the electrostatic and covalent interactions in the DDAA dimer. The SEI model is thus predictive because it provides a measure for the charge accumulation in hydrogen-bonded monomers. Our findings can be understood from simple physical organic chemistry principles and provide supramolecular chemists with meaningful understanding for tuning hydrogen bond strengths and thus for controlling the properties of self-assembled systems.

We discuss the concept of cation affinities (CA) and provide an overview of topical trends throughout the periodic table, inferred from state-of-the-art relativistic quantum-chemical computations. The CA of a base B (−) for a cation Y + is the energy or enthalpy required to dissociate the complex BY (+) into molecular fragments B (−) and Y + . Probably the best-known CA is the proton affinity (PA) of Lewis bases. We extend this concept here to include methyl cation (MCA) and other alkyl cation affinities (ACA) as well as alkali metal cation affinities (AMCA). The Lewis bases covered herein are the anionic and neutral element hydrides B − = XH n−1 − and B = XH n , respectively, as well as methyl-substituted variants thereof. The element “X” in our model Lewis bases covers the maingroup elements of groups 14–18 in rows 1–6 of the periodic table. Emerging trends are analyzed and explained in terms of quantitative molecular orbital (MO) theory as contained in Kohn–Sham density functional theory (KS-DFT). Making the often implicitly used idea of a CA explicit serves a more rational design of compounds with a particular affinity for, or reactivity toward, other species throughout the molecular sciences, from inorganic via organic to biological chemistry.

For a series of organonitrile [RCN (R = Me, CF3, Ph, CH3Ph, CF3Ph)] ligands, the nature of the N–Pt bond in the related cis-/trans-(RCN)2PtCl2 complexes has been computationally investigated by Density Functional Theory. A fragment based bond analysis has been performed in the canonical Kohn–Sham molecular orbitals framework, and it has been ultimately assessed that this bond is characterized both by N→Pt σ and by N←Pt π contributions. Voronoi Deformation Density charges further confirms the occurrence of N←Pt π interactions. Moreover, the Energy Decomposition Analysis-Natural Orbital for Chemical Valence (EDA-NOCV) method shows that the strength of the N←Pt π interaction is not negligible by contributing to about 30–40% of the total orbital interaction. Finally, the well-known νCN blue-shift occurring upon coordination to PtII, has been thoroughly investigated by exploiting the EDA-NOCV and by evaluating νCN and force constants. The origin of the νCN blue-shift in these systems has been discussed on the basis of the CN bond polarization. N←Pt π backbonding causes only a systematic decrease of the observed νCN blue-shift when compared to the one calculated for RCN–X (X = H+, alkaline, Lewis acids) herein reported (X = purely σ acceptors).