Search

Appreciating reactivity in terms of physicochemical effects for chemical processes is one of the most important undertakings in chemistry. While transition state (TS) theory provides the framework enabling the reliable calculation of the barrier height for a given elementary step, analytical tools are necessary to gain insight into key factors governing the different processes during chemical reactions. In this contribution, we partition the potential energy surface of an elementary step along the intrinsic coordinate into three segments, the so-called Pre-TS, TS, and Post-TS regions, and then determine the most important factors dictating each segment. This analysis is based on the use of both reactivity descriptors from conceptual density functional theory and concepts from the information-theoretic approach in density functional theory. We found that in both Pre-TS and Post-TS regions, steric effects are the dominant factors, whereas in the TS region, it is the intrinsic electrophilic and nucleophilic propensity of the transition state structure that governs the reactivity. The wide applicability of our approach is shown by a validation for a total of 37 organic and inorganic reactions. Different partition approaches from density functional theory, energy decomposition analysis and wavefunction-based resonance theory are employed to support our main conclusions., (© 2024 Wiley-VCH GmbH.)

Rechargeable batteries are essential to the global shift towards renewable energy sources and their storage. At present, improvements in their safety and sustainability are of great importance as part of global sustainable development goals. A major contender in this shift are rechargeable solid-state sodium batteries, as a low-cost, safe, and sustainable alternative to conventional lithium-ion batteries. Recently, solid-state electrolytes with a high ionic conductivity and low flammability have been developed. However, these still face challenges with the highly reactive sodium metal electrode. The study of these electrolyte-electrode interfaces is challenging from a computational and experimental point of view, but recent advances in molecular dynamics neural-network potentials are finally enabling access to these environments compared to more computationally expensive conventional ab-initio techniques. In this study, heteroatom-substituted Na 3 PS 3 X 1 analogues, where X is sulfur, oxygen, selenium, tellurium, nitrogen, chlorine, and fluorine, are investigated using total-trajectory analysis and neural-network molecular dynamics. It was found that inductive electron-withdrawing and electron-donating effects, alongside differences in heteroatom atomic radius, electronegativity, and valency, influenced the electrolyte reactivity. The Na 3 PS 3 O 1 oxygen analogue was found to have superior chemical stability against the sodium metal electrode, paving the way towards high-performance, long lifetime and reliable rechargeable solid-state sodium batteries., (© 2023 The Authors. ChemSusChem published by Wiley-VCH GmbH.)

The drive to develop maximal atom-efficient catalysts coupled to the continuous striving for more sustainable reactions has led to an ever-increasing interest in single-atom catalysis. Based on a periodic conceptual density functional theory (cDFT) approach, fundamental insights into the reactivity and adsorption of single late transition metal atoms supported on a fully hydroxylated amorphous silica surface have been acquired. In particular, this investigation revealed that the influence of van der Waals dispersion forces is especially significant for a silver (98 %) or gold (78 %) atom, whereas the oxophilicity of the Group 8-10 transition metals plays a major role in the interaction strength of these atoms on the irreducible SiO 2 support. The adsorption energies for the less-electronegative row 4 elements (Fe, Co, Ni) ranged from -1.40 to -1.92 eV, whereas for the heavier row 5 and 6 metals, with the exception of Pd, these values are between -2.20 and -2.92 eV. The deviating behavior of Pd can be attributed to a fully filled d-shell and, hence, the absence of the hybridization effects. Through a systematic analysis of cDFT descriptors determined by using three different theoretical schemes, the Fermi weighted density of states approach was identified as the most suitable for describing the reactivity of the studied systems. The main advantage of this scheme is the fact that it is not influenced by fictitious Coulomb interactions between successive, charged reciprocal cells. Moreover, the contribution of the energy levels to the reactivity is simultaneously scaled based on their position relative to the Fermi level. Finally, the obtained Fermi weighted density of states reactivity trends show a good agreement with the chemical characteristics of the investigated metal atoms as well as the experimental data., (© 2020 Wiley-VCH GmbH.)

A hitherto unexplored class of molecules for molecular force probe applications are expanded porphyrins. This work proves that mechanical force is an effective stimulus to trigger the interconversion between Hückel and Möbius topologies in [28]hexaphyrin, making these expanded porphyrins suitable to act as conformational mechanophores operating at mild (sub-1 nN) force conditions. A straightforward approach based on distance matrices is proposed for the selection of pulling scenarios that promote either the planar Hückel topology or the three lowest lying Möbius topologies. This approach is supported by quantum mechanochemical calculations. Force distribution analyses reveal that [28]hexaphyrin selectively allocates the external mechanical energy to molecular regions that trigger Hückel-Möbius interconversions, explaining why certain pulling scenarios favor the Hückel two-sided topology and others favor Möbius single-sided topologies. The meso-substitution pattern on [28]hexaphyrin determines whether the energy difference between the different topologies can be overcome by mechanical activation., (© 2020 The Authors. Chemistry - A European Journal published by Wiley-VCH GmbH.)

Following the recent discovery of stable octa-coordinated alkaline earth metals with N 2 and CO, the role of group II metals in the catalytic reduction of these ligands by means of density functional theory (DFT) calculations and conceptual DFT-based reactivity indices is investigated. Cubic group IV and octahedral group VI transition metal complexes as well as the free ligands are computed for reference. The outer and most accessible atoms of N 2 and CO become much more nucleophilic and electrophilic in all complexes, relevant for N 2 fixation, as probed by the Fukui function and local softness. Within one row of the periodic table, the alkaline earth complexes often show the strongest activation. On the contrary, the electrostatic character is found to be virtually unaffected by complexation. Trends in the soft frontier orbital and hard electrostatic character are in agreement with calculated proton affinities and energy decomposition analyses of the protonated structures, demonstrating the dominance of the soft (HOMO-LUMO) orbital interactions., (© 2020 The Authors. Published by Wiley-VCH GmbH.)

The ability to understand and predict ambident reactivity is key to the rational design of organic syntheses. An approach to understand trends in ambident reactivity is the hard and soft acids and bases (HSAB) principle. The recent controversy over the general validity of this principle prompted us to investigate the competing gas-phase S N 2 reaction channels of archetypal ambident nucleophiles CN - , OCN - , and SCN - with CH 3 Cl (S N 2@C) and SiH 3 Cl (S N 2@Si), using DFT calculations. Our combined analyses highlight the inability of the HSAB principle to correctly predict the reactivity trends of these simple, model reactions. Instead, we have successfully traced reactivity trends to the canonical orbital-interaction mechanism and the resulting nucleophile-substrate interaction energy. The HOMO-LUMO orbital interactions set the trend in both S N 2@C and S N 2@Si reactions. We provide simple rules for predicting the ambident reactivity of nucleophiles based on our Kohn-Sham molecular orbital analysis., (© 2020 The Authors. Published by Wiley-VCH Verlag GmbH & Co. KGaA.)

Invited for this month's cover is the Quantum Chemistry group (ALGC) of the Free University of Brussels (VUB) in collaboration with the Molecular Design and Synthesis (MolDesignS) group of the KU Leuven. The image shows how a variety of computational techniques, such as density functional theory, molecular dynamics, and the noncovalent interaction index, can be used to virtually zoom in on the noncovalent interactions that are key for the supramolecular hydrogelation of bis-urea-based gelators. Read the full text of the article at 10.1002/cplu.201900551., (© 2020 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim.)

The current approach to designing low-molecular-weight gelators relies on a laborious trial-and-error process, mainly because of the lack of an accurate description of the noncovalent interactions crucial for supramolecular gelation. In this work, we report a multiscale bottom-up approach composed of several computational techniques to unravel the key interactions in a library of synthesized bis-urea-based gelators and rationalize their experimentally observed hydrogelation performance. In addition to density functional theory calculations and molecular dynamics, the noncovalent interaction index is applied as a tool to visualise and identify the different types of noncovalent interactions. Interestingly, as well as hydrogen bonds between urea moieties, hydrogen bonds between a urea moiety and a pyridine ring were shown to play a detrimental role in the early aggregation phase. These findings enabled us to explain the hydrogelation performance observed in a library of twelve bis-urea derivatives, which were synthesized with 58-95 % yields. From this library, three compounds were discovered to effectively gel water, with the most efficient hydrogelator only requiring a concentration of 0.2 w/v%., (© 2019 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim.)

The use of a new class of unsymmetrical cinchona-alkaloid-based, phthalazine-bridged organocatalysts enabled the highly enantioselective dichlorination of unfunctionalized alkenes. In combination with the electrophilic chlorinating agent 1,3-dichloro-5,5-dimethylhydantoin (DCDMH) and triethylsilyl chloride (TES-Cl) as the source of nucleophilic chloride, 1-aryl-2-alkyl alkenes were dichlorinated with enantioselectivities of up to 94:6 er. Initial mechanistic investigations suggest that no free chlorine is formed, and by replacement of the chloride by fluoride, enantioselective chlorofluorinations of alkenes are possible., (© 2019 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim.)

The importance of axially chiral biaryls has risen steeply in the recent decades. This structural motif proved to be successful in catalytic asymmetric synthesis and the configuration of the biaryl axis is decisive for the biological activity. A new approach for the atroposelective synthesis of biaryls would be through a cycloaddition between an enantiopure phenyl-substituted thiophene S-oxide and an alkyne. Importantly, the chiral center of the thiophene S-oxide needs to be stable enough to avoid pyramidal inversion during the cycloaddition. Considering that the racemization of thiophene monoxides has been scarcely investigated so far, we perform a thorough quantum chemical study on the inversion barriers of a large number of chiral thiophene S-oxide derivatives. Our main goal is to identify substitution patterns leading to stable atropisomers at room temperature. Appealingly, the role of stereoelectronic effects and the position of the substituents as well as the importance of aromaticity on the pyramidal inversion barrier are elucidated for the first time., (© 2019 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim.)

The solvent and autocatalytic effects of the electrophilic aromatic chlorination of benzene are studied using a combined approach of static calculations and ab initio metadynamics simulations. Different possible reaction pathways are investigated and the influence of the solvents (CCl 4 , acetonitrile and acetic acid) is thoroughly assessed. Our results show that the stability and lifetime of a charged σ-complex is increased by electrostatic stabilisation effects of the environment, which can originate from catalytic HCl, solvating effects of polar solvents (acetonitrile), or specific hydrogen bonding interactions with the solvent (acetic acid). Metadynamics simulations reveal a new chlorine addition mechanism explaining the autocatalytic effects of the reaction. The strength of combining static calculations and metadynamics simulations is highlighted, which provide complementary insight into chemical reactions in solvent., (© 2018 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim.)

In this study, an intuitive concept is derived, which explains the characteristic dependence of spin-state energetics on the exact exchange admixture of DFT functionals in the case of octahedral transition metal complexes. The change in electron density distributions upon varying the admixture, c 3 , in the B3LYP functional is analyzed for archetype ionic and covalent systems as well as for the Fe 2+ ion in an ideal octahedral field. An understanding of how the DFT description of the electronic structure of octahedral complexes changes as a function of c 3 is sought. A systematic spin-state energy analysis of 50 octahedral complexes of various metals and ligands with consistent experimental data is presented, allowing the derivation, in theory, of an optimal c 3 value for each system. The notion that the admixture dependence of DFT spin-state energetics stems from the treatment of nondynamic electrons arising from the mixing of (M-Lz2 ) 0 (dz2 ) 2 and (M-Lx2-y2 ) 0 (dx2-y2 ) 2 configurations into the dominant (M-Lx2-y2 ) 2 (dx2-y2 ) 0 and (M-Lx2-y2 ) 2 (dx2-y2 ) 0 ones in the low(er) spin states is put forward. That is, in the effort to mimic such electron-electron interactions, E x LDA overestimates, whereas exact exchange downplays the contribution of this type of electron correlation to the stability of low(er) spin states, leading to the widespread practical observation that the higher the exact exchange admixture, the more stable the high-spin-state configuration., (© 2018 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim.)

In the last decade, halogen bonds, noncovalent interactions formed between positive regions in the electrostatic potential on halogen atoms, often referred to as σ-holes, and electron-rich sites, have gained a lot of interest. Recently, this interest has been expanded towards interactions with Group V and Group VI elements, giving rise to pnicogen and chalcogen bonds. Although chalcogen bonds have already shown some promising results for applications in crystallography and catalysis, experimental results characterising these noncovalent interactions remain scarce. In this combined experimental and theoretical study, original data allowing the characterization of S⋅⋅⋅O chalcogen bonds is obtained by studying the 1:1 molecular complexes between 2,2,4,4-tetrafluoro-1,3-dithiethane (C 2 F 4 S 2 ) and dimethyl ether (DME). Ab initio calculations of the C 2 F 4 S 2 ⋅DME dimer yield two stable chalcogen-bonded isomers, the difference being the presence or absence of secondary F⋅⋅⋅H interactions. Liquid-krypton solutions containing C 2 F 4 S 2 and DME were studied using FTIR and Raman spectroscopy. Upon subtraction of rescaled monomer spectra, clear complex bands are observed. The observed complexation shifts agree favourably with the ab initio calculated shifts of the chalcogen-bonded complexes. The 1:1 stoichiometry of the complex is confirmed and a complexation enthalpy of -13.5(1) kJ mol -1 is found, which is in good agreement with the calculated values. A Ziegler-Rauk energy decomposition analysis revealed that electrostatic interactions prominently dominate over orbital interactions. Nevertheless, significant charge transfer occurs from the oxygen in DME to one of the sulfur atoms in C 2 F 4 S 2 and the carbon along the extension of the chalcogen bond., (© 2017 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim.)

The bonding in heavier Group 14 zero-valent complexes of a general formula L 2 E (E=Si-Pb; L=phosphine, N-heterocyclic and acyclic carbene, cyclic tetrylene and carbon monoxide) is probed by combining valence bond (VB) theory and maximum probability domain (MPD) approaches. All studied complexes are initially evaluated on the basis of the structural parameters and the shape of frontier orbitals revealing a bent structural motif and the presence of two lone pairs at the central E atom. For the VB calculations three resonance structures are suggested, representing the "ylidone", "ylidene" and "bent allene" structures, respectively. The influence of both ligands and central atoms on the bonding situation is clearly expressed in different weights of the resonance structures for the particular complexes. In general, the bonding in the studied E 0 compounds, the tetrylones, is best described as a resonating combination of "ylidone" and "ylidene" structures with a minor contribution of the "bent allene" structure. Moreover, the VB calculations allow for a straightforward assessment of the π-backbonding (E→L) stabilization energy. The validity of the suggested resonance model is further confirmed by the complementary MPD calculations focusing on the E lone pair region as well as the E-L bonding region. Likewise, the MPD method reveals a strong influence of the σ-donating and π-accepting properties of the ligand. In particular, either one single domain or two symmetrical domains are found in the lone pair region of the central atom, supporting the predominance of either the "ylidene" or "ylidone" structures having one or two lone pairs at the central atom, respectively. Furthermore, the calculated average populations in the lone pair MPDs correlate very well with the natural bond orbital (NBO) populations, and can be related to the average number of electrons that is backdonated to the ligands., (© 2017 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim.)

G protein-coupled receptors (GPCRs) play an important role in many cellular responses; as such, their mechanism of action is of utmost interest. To gain insight into the active conformation of GPCRs, the X-ray crystal structures of nanobody (Nb)-stabilized β 2 -adrenergic receptor (β 2 AR) have been reported. Nb80, in particular, is able to bind the intracellular G protein binding site of β 2 AR and stabilize the receptor in an active conformation. Within Nb80, the complementarity-determining region 3 (CDR3) is responsible for most of the binding interactions. Hence, we hypothesized that peptidomimetics of the CDR3 loop might be sufficient for binding to the receptor, inhibiting the interaction of β 2 AR with intracellular GPCR interacting proteins (e.g., G proteins). Based on previous crystallographic data, a set of peptidomimetics were synthesized that, similar to the Nb80 CDR3 loop, adopt a β-hairpin conformation. Syntheses, conformational analysis, binding and functional in vitro assays, as well as internalization experiments, were performed. We demonstrate that peptidomimetics can structurally mimic the CDR3 loop of a nanobody and its function by inhibiting G protein coupling as measured by partial inhibition of cAMP production., (© 2017 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Our attempts to synthesise N→M intramolecularly coordinated diorganometallic hydrides L2MH2 [M=Si (4), Ge (5), Sn (6)] containing the CH=N imine group (in which L is C,N-chelating ligand {2-[(2,6-iPr2C6H3)N=CH]C6 H4}(-)) yielded 1,1'-bis(2,6-diisopropylphenyl)-2,2'-spriobi[benzo[c][1,2]azasilole] (7), 1,1'-bis(2,6-diisopropylphenyl)-2,2'-spriobi[benzo[c][1,2]azagermole] (8) and C,N-chelated homoleptic stannylene L2Sn (10), respectively. Compounds 7 and 8 are an outcome of a spontaneous double hydrometallation of the two CH=N imine moieties induced by N→M intramolecular coordination (M=Si, Ge) in the absence of any catalyst. In contrast, the diorganotin hydride L2SnH2 (6) is redox-unstable and the reduction of the tin centre with the elimination of H2 provided the C,N-chelated homoleptic stannylene L2Sn (10). Compounds 7 and 8 were characterised by NMR spectroscopy and X-ray diffraction analysis. Because the proposed N→M intramolecularly coordinated diorganometallic hydrides L2MH2 [M=Si (4), Ge (5), Sn (6)] revealed two different types of reduction reactions, DFT calculations were performed to gain an insight into the structures and bonding of the non-isolable diorganometallic hydrides as well as the products of their subsequent reactions. Furthermore, the thermodynamic profiles of the different reaction pathways with respect to the central metal atom were also investigated., (© 2016 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

The heretofore unpredictable behavior of [26] and [28]hexaphyrins upon metalation has been elucidated through quantum chemical calculations. It is demonstrated that the molecular topology of Group 10 and Group 11 metal complexes of hexaphyrins depends on sensitive interplay between the intrinsic ligand strain and the metal-ligand interaction strength. As such, the aromaticity of the ligand and effective charge of the metal are revealed as key factors determining the binding mode and the preference for Möbius or Hückel structures. These findings offer a new perspective to rationalize experimental observations for metalated hexaphyrins. More importantly, the proposed guidelines could be useful for designing novel complexes of hexaphyrins, such as a hitherto unknown Möbius [26]hexaphyrin complex., (© 2015 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

The reduction of N,C,N-chelated bismuth chlorides [C6H3-2,6-(CH=NR)2]BiCl2 [where R = tBu (1), 2',6'-Me2C6H3 (2), or 4'-Me2NC6H4 (3)] or N,C-chelated analogues [C6H2-2-(CH=N-2',6'-iPr2C6H3)-4,6-(tBu)2]BiCl2 (4) and [C6H2-2-(CH2NEt2)-4,6-(tBu)2]BiCl2 (5) is reported. Reduction of compounds 1-3 gave monomeric N,C,N-chelated bismuthinidenes [C6H3-2,6-(CH=NR)2]Bi [where R = tBu (6), 2',6'-Me2C6H3 (7) or 4'-Me2NC6H4 (8)]. Similarly, the reduction of 4 led to the isolation of the compound [C6H2-2-(CH=N-2',6'-iPr2C6H3)-4,6-(tBu)2]Bi (9) as an unprecedented two-coordinated bismuthinidene that has been structurally characterized. In contrast, the dibismuthene {[C6H2-2-(CH2NEt2)-4,6-(tBu)2]Bi}2 (10) was obtained by the reduction of 5. Compounds 6-10 were characterized by using (1)H and (13)C NMR spectroscopy and their structures, except for 7, were determined with the help of single-crystal X-ray diffraction analysis. It is clear that the structure of the reduced products (bismuthinidene versus dibismuthene) is ligand-dependent and particularly influenced by the strength of the N→Bi intramolecular interaction(s). Therefore, a theoretical survey describing the bonding situation in the studied compounds and related bismuth(I) systems is included. Importantly, we found that the C3NBi chelating ring in the two-coordinated bismuthinidene 9 exhibits significant aromatic character by delocalization of the bismuth lone pair., (© 2015 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

The contributions of covalent and noncovalent interactions to the formation of classical adducts of bulky Lewis acids and bases and frustrated Lewis pairs (FLPs) were scrutinized by using various conceptual quantum chemical techniques. Significantly negative complexation energies were calculated for fourteen investigated Lewis pairs containing bases and acids with substituents of various sizes. A Ziegler-Rauk-type energy decomposition analysis confirmed that two types of Lewis pairs can be distinguished on the basis of the nature of the primary interactions between reactants; dative-bond formation and concomitant charge transfer from the Lewis base to the acid is the dominant and most stabilizing factor in the formation of Lewis acid-base adducts, whereas weak interactions are the main thermodynamic driving force (>50 %) for FLPs. Moreover, the ease and extent of structural deformation of the monomers appears to be a key component in the formation of the former type of Lewis pairs. A Natural Orbital for Chemical Valence (NOCV) analysis, which was used to visualize and quantify the charge transfer between the base and the acid, clearly showed the importance and lack of this type of interaction for adducts and FLPs, respectively. The Noncovalent Interaction (NCI) method revealed several kinds of weak interactions between the acid and base components, such as dispersion, π-π stacking, C-H⋅⋅⋅π interaction, weak hydrogen bonding, halogen bonding, and weak acid-base interactions, whereas the Quantum Theory of Atoms in Molecules (QTAIM) provided further conceptual insight into strong acid-base interactions., (© 2015 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Insight into the key factors driving the competition of halogen and hydrogen bonds is obtained by studying the affinity of the Lewis bases trimethylamine (TMA), dimethyl ether (DME), and methyl fluoride (MF) towards difluoroiodomethane (CHF(2) I). Analysis of the infrared and Raman spectra of solutions in liquid krypton containing mixtures of TMA and CHF(2) I and of DME and CHF(2) I reveals that for these Lewis bases hydrogen and halogen-bonded complexes appear simultaneously. In contrast, only a hydrogen-bonded complex is formed for the mixtures of CHF(2) I and MF. The complexation enthalpies for the C-H⋅⋅⋅Y hydrogen-bonded complexes with TMA, DME, and MF are determined to be -14.7(2), -10.5(5) and -5.1(6) kJ mol(-1), respectively. The values for the C-I⋅⋅⋅Y halogen-bonded isomers are -19.0(3) kJ mol(-1) for TMA and -9.9(8) kJ mol(-1) for DME. Generalization of the observed trends suggests that, at least for the bases studied here, softer Lewis bases such as TMA favor halogen bonding, whereas harder bases such as MF show a substantial preference for hydrogen bonding., (© 2014 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Novel organic-inorganic hybrids of various sizes were generated by reaction of 1,8-octanediphosphonic acid (ODP) and (NH4)6Mo7O24 in aqueous solution. The formation of rodlike hybrids with variable numbers of covalently bound ODP and polyoxomolybdate (POM) units can be tuned as a function of increasing (NH4)6Mo7O24 concentration at fixed ODP concentration. The chemical structure of the ODP/POM hybrids was characterized by (1)H, (31)P, and (95)Mo NMR spectroscopy. Heteronuclear (31)P DOSY (diffusion- ordered NMR spectroscopy) and molecular mechanics (MM) calculations were applied to determine the size and shape of the nanosized hybrids generated at various ODP/POM ratios. For this purpose, the structures of ODP/POM hybrids with variable numbers of ODP and POM units were optimized by MM and then approximated as cylinder-shaped objects by using a recently described mathematical algorithm. The thus-obtained cylinder length and diameter were further used to calculate the expected diffusion coefficients of the ODP/POM hybrids. Comparison of the calculated and experimentally determined diffusion coefficients led to the most probable ODP/POM hybrid length for each sample composition. The (31)P DOSY results show that the length of the hybrids increases with increasing POM concentration and reaches a maximum corresponding to an average of 8 ODP/7 POM units per chain at a sample composition of 20 mM ODP and 14 mM POM. With excess POM, above the latter concentration, the formation of shorter-chain hybrids terminated by Mo7 clusters at one or both ends was evidenced on further increasing the POM concentration. The results demonstrate that the combination of (31)P DOSY and MM, although virtually unexplored in POM chemistry, is a powerful innovative strategy for the detailed characterization of nanosized organic-inorganic POM-based hybrids in solution., (© 2014 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Noncovalent interactions involving aromatic rings, such as π-stacking and CH/π interactions, are central to many areas of modern chemistry. However, recent studies proved that aromaticity is not required for stacking interactions, since similar interaction energies were computed for several aromatic and aliphatic dimers. Herein, the nature and origin of π/π, σ/σ, and σ/π dispersion interactions has been investigated by using dispersion-corrected density functional theory, energy decomposition analysis, and the recently developed noncovalent interaction (NCI) method. Our analysis shows that π/π and σ/σ stacking interactions are equally important for the benzene and cyclohexane dimers, explaining why both compounds have similar boiling points. Also, similar dispersion forces are found in the benzene⋅⋅⋅methane and cyclohexane⋅⋅⋅methane complexes. However, for systems larger than naphthalene, there are enhanced stacking interactions in the aromatic dimers adopting a parallel-displaced configuration compared to the analogous saturated systems. Although dispersion plays a decisive role in stabilizing all the complexes, the origin of the π/π, σ/σ, and σ/π interactions is different. The NCI method reveals that the dispersion interactions between the hydrogen atoms are responsible for the surprisingly strong aliphatic interactions. Moreover, whereas σ/σ and σ/π interactions are local, the π/π stacking are inherently delocalized, which give rise to a non-additive effect. These new types of dispersion interactions between saturated groups can be exploited in the rational design of novel carbon materials., (© 2014 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Our attempts to synthesize the N→Si intramolecularly coordinated organosilanes Ph2 L(1) SiH (1 a), PhL(1) SiH2 (2 a), Ph2 L(2) SiH (3 a), and PhL(2) SiH2 (4 a) containing a CH=N imine group (in which L(1) is the C,N-chelating ligand {2-[CH=N(C6 H3 -2,6-iPr2)]C6 H4}(-) and L(2) is {2-[CH=N(tBu)]C6 H4}(-)) yielded 1-[2,6-bis(diisopropyl)phenyl]-2,2-diphenyl-1-aza-silole (1), 1-[2,6-bis(diisopropyl)phenyl]-2-phenyl-2-hydrido-1-aza-silole (2), 1-tert-butyl-2,2-diphenyl-1-aza-silole (3), and 1-tert-butyl-2-phenyl-2-hydrido-1-aza-silole (4), respectively. Isolated organosilicon amides 1-4 are an outcome of the spontaneous hydrosilylation of the CH=N imine moiety induced by N→Si intramolecular coordination. Compounds 1-4 were characterized by NMR spectroscopy and X-ray diffraction analysis. The geometries of organosilanes 1 a-4 a and their corresponding hydrosilylated products 1-4 were optimized and fully characterized at the B3LYP/6-31++G(d,p) level of theory. The molecular structure determination of 1-3 suggested the presence of a Si-N double bond. Natural bond orbital (NBO) analysis, however, shows a very strong donor-acceptor interaction between the lone pair of the nitrogen atom and the formal empty p orbital on the silicon and therefore, the calculations show that the Si-N bond is highly polarized pointing to a predominantly zwitterionic Si(+) N(-) bond in 1-4. Since compounds 1-4 are hydrosilylated products of 1 a-4 a, the free energies (ΔG298), enthalpies (ΔH298), and entropies (ΔH298) were computed for the hydrosilylation reaction of 1 a-4 a with both B3LYP and B3LYP-D methods. On the basis of the very negative ΔG298 values, the hydrosilylation reaction is highly exergonic and compounds 1 a-4 a are spontaneously transformed into 1-4 in the absence of a catalyst., (Copyright © 2014 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Homobimetallic metallophilic interactions between copper, silver, and gold-based [(NHC)MX]-type complexes (NHC=N-heterocyclic carbene, i.e, 1,3,4-trimethyl-4,5-dihydro-1H-1,2,4-triazol-5-ylidene; X=F, Cl, Br, I) were investigated by means of ab initio interaction energies, Ziegler-Rauk-type energy-decomposition analysis, the natural orbital for chemical valence (NOCV) framework, and the noncovalent interaction (NCI) index. It was found that the dimers of these complexes predominantly adopt a head-to-tail arrangement with typical M⋅⋅⋅M distance of 3.04-3.64 Å, in good agreement with the experimental X-ray structure determined for [{(NHC)AuCl}2 ], which has an Au⋅⋅⋅Au distance of 3.33 Å. The interaction energies between silver- and gold-based monomers are calculated to be about -25 kcal mol(-1) , whereas that for the Cu congener is significantly lower (-19.7 kcal mol(-1) ). With the inclusion of thermal and solvent contributions, both of which are destabilizing, by about 15 and 8 kcal mol(-1) , respectively, an equilibrium process is predicted for the formation of dimer complexes. Energy-decomposition analysis revealed a dominant electrostatic contribution to the interaction energy, besides significantly stabilizing dispersion and orbital interactions. This electrostatic contribution is rationalized by NHC(δ(+) )⋅⋅⋅halogen(δ(-) ) interactions between monomers, as demonstrated by electrostatic potentials and derived charges. The dominant NOCV orbital indicates weakening of the π backdonation in the monomers on dimer formation, whereas the second most dominant NOCV represents an electron-density deformation according to the formation of a very weak M⋅⋅⋅M bond. One of the characteristic signals found in the reduced density gradient versus electron density diagram corresponds to the noncovalent interactions between the metal centers of the monomers in the NCI plots, which is the manifestation of metallophilic interaction., (Copyright © 2014 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

The effect of non-polar and polar ligands and of monovalent cations on the one-electron reduction potential of the thiyl radical and the disulfide bond was evaluated. The reduction potentials E° for the CH3S(.)-nL/CH3S(-)-nL and CH3SSCH3-L/CH3SSCH3(.-)-L redox couples were calculated at the B3LYP, M06-2X and MP2 levels of theory, with n=1, 2 and L=CH4, C2H4, H2O, CH3OH, NH3, CH3COOH, CH3CONH2, NH4(+), Na(+), K(+) and Li(+). Non-polar ligands decrease the E° value of the thiyl radical and disulfide bond, while neutral polar ligands favour electron uptake. Charged polar ligands and cations favour electron capture by the thiyl radical while disfavouring electron uptake by the disulfide bond. Thus, the same type of ligand can have a different effect on E° depending on the redox couple. Therefore, properties of an isolated ligand cannot uniquely determine E°. The ligand effects on E° are discussed in terms of the vertical electron affinity and reorganization energy, as well as molecular orbital theory. For a given redox couple, the ligand type influences the nature of the anion formed upon electron capture and the corresponding reorganization process towards the reduced geometry., (Copyright © 2013 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

The switching of topology between "figure-eight", Möbius, and untwisted conformations in [32]heptaphyrins(1.1.1.1.1.1.1) has been investigated by using density functional theory calculations. Such a change is achieved by variation of one internal dihedral angle and, if properly controlled, can provide access to molecular switches with unique optical and magnetic properties. In this work, we have explored different conformational control methods, such as solvent, protonation and meso substituents. Despite its antiaromatic character, most of the [32]heptaphyrins (R=H, CH(3), CF(3), Ph, C(6)F(5)) adopt a figure-eight conformation in the neutral state, owing to their more-effective hydrogen-bonding interactions. The aromatic Möbius topology is only preferred with dichlorophenyl groups, which minimize the steric hindrance that arises from the bulky chlorine atoms. The conformational equilibrium is sensitive to the solvent, so polar solvents, such as DMSO, further stabilize the Möbius conformation. Protonation induces a conformational change into the Möbius topology, irrespective of the meso-aryl groups. In the triprotonated species, the conformational switch is blocked and a non-twisted conformer becomes much more stable than the figure-eight conformation. We have shown that the relative energies of the protonated [32]heptaphyrins are dominated by aromaticity. Importantly, this topology switching induces a dramatic change in the magnetic properties and reactivity of the macrocycles, as revealed by several energetic, magnetic, structural, and reactivity indices of aromaticity., (Copyright © 2013 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Halogen bonds between the trifluoromethyl halides CF(3)Cl, CF(3)Br and CF(3)I, and dimethyl ether, dimethyl sulfide, trimethylamine and trimethyl phosphine were investigated using Pearson's hard and soft acids and bases (HSAB) concept with conceptual DFT reactivity indices, the Ziegler-Rauk-type energy-decomposition analysis, the natural orbital for chemical valence (NOCV) framework and the non-covalent interaction (NCI) index. It is found that the relative importance of electrostatic and orbital (charge transfer) interactions varies as a function of both the donor and acceptor molecules. Hard and soft interactions were distinguished and characterised by atomic charges, electrophilicity and local softness indices. Dual-descriptor plots indicate an orbital σ hole on the halogen similar to the electrostatic σ hole manifested in the molecular electrostatic potential. The predicted high halogen-bond-acceptor affinity of N-heterocyclic carbenes was evidenced in the highest complexation energy for the hitherto unknown CF(3) I·NHC complex. The dominant NOCV orbital represents an electron-density deformation according to a n→σ*-type interaction. The characteristic signal found in the reduced density gradient versus electron-density diagram corresponds to the non-covalent interaction between contact atoms in the NCI plots, which is the manifestation of halogen bonding within the NCI theory. The unexpected C-X bond strengthening observed in several cases was rationalised within the molecular orbital framework., (Copyright © 2013 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Recently, hexaphyrins have emerged as a promising class of π-conjugated molecules that display a range of interesting electronic, optical, and conformational properties, including the formation of stable Möbius aromatic systems. Besides the Möbius topology, hexaphyrins can adopt a variety of conformations with Hückel and twisted Hückel topologies, which can be interconverted under certain conditions. To determine the optimum conditions for viable Möbius topologies, the conformational preferences of [26]- and [28]hexaphyrins and the dynamic interconversion between the Möbius and Hückel topologies were investigated by density functional calculations. In the absence of meso substituents, [26]hexaphyrin prefers a planar dumbbell conformation, strongly aromatic and relatively strain free. The Möbius topology is highly improbable: the most stable tautomer is 33 kcal mol(-1) higher in energy than the global minimum. On the other hand, the Möbius conformer of [28]hexaphyrin is only 6.5 kcal mol(-1) higher in energy than the most stable dumbbell conformation. This marked difference is due to aromatic stabilization in the Möbius 4n electron macrocycle as opposed to antiaromatic destabilization in the 4n+2 electron system, as revealed by several energetic, magnetic, structural, and reactivity indices of aromaticity. For [28]hexaphyrins, the computed activation barrier for interconversion between the Möbius aromatic and Hückel antiaromatic conformers ranges from 7.2 to 10.2 kcal mol(-1), in very good agreement with the available experimental data. The conformation of the hexaphyrin macrocycle is strongly dependent on oxidation state and solvent, and this feature creates a promising platform for the development of molecular switches., (Copyright © 2012 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

We analyse the electronic structure and aromaticity of graphene nanoribbons and carbon nanotubes through a series of delocalisation and geometry analysis methods. In particular, the six-centre index (SCI) is found to be in good agreement with the mean bond length (MBL) and ring bond dispersion (RBD) geometry descriptors. Based on DFT periodic calculations, three distinct classes of aromaticity patterns have been found for armchair graphene nanoribbons, appearing periodically as the width of the ribbon is increased. The periodicity in the band gap is found to be related to these aromaticity patterns. Also, the appearance of such distinct aromaticity distribution is explained within the framework of the Clar's sextet theory. Both delocalisation and geometry analysis methods are shown to be very fast and reliable tools for easily analysing the aromaticity in carbon nanosystems., (Copyright © 2012 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

In this work, we will show that the largest values of the local softness and hardness do not necessarily correspond to the softest and hardest regions of the molecule, respectively. Based on our results, we will argue that it is more useful to interpret the local softness and the local hardness as functions that measure the "local abundance" or "concentration" of the corresponding global properties. This new point of view helps reveal how and when these local reactivity indices are most useful.

The Woodward-Hoffmann rules for pericyclic reactions are explained entirely in terms of directly observable physical properties of molecules (specifically changes in electron density) without any recourse to model-dependent concepts, such as orbitals and aromaticity. This results in a fundamental explanation of how the physics of molecular interactions gives rise to the chemistry of pericyclic reactions. This construction removes one of the key outstanding problems in the qualitative density-functional theory of chemical reactivity (the so-called conceptual DFT). One innovation in this paper is that the link between molecular-orbital theory and conceptual DFT is treated very explicitly, revealing how molecular-orbital theory can be used to provide "back-of-the-envelope" approximations to the reactivity indicators of conceptual DFT.

Previous contributions concerning the computational approach to redox chemistry have made use of thermodynamic cycles and Car-Parrinello molecular dynamics simulations to obtain accurate redox potential values, whereas this article adopts a conceptual density functional theory (DFT) approach. Conceptual DFT descriptors have found widespread use in the study of thermodynamic and kinetic aspects of a variety of organic and inorganic reactions. However, redox reactions have not received much attention until now. In this contribution, we prove the usefulness of global and local electrophilicity descriptors for the prediction of the redox characteristics of first row transition metal ions (from Sc(3+) | Sc(2+) to Cu(3+) | Cu(2+)) and introduce a scaled definition of the electrophilicity based on the number of electrons an electrophile ideally accepts. This scaled electrophilicity concept acts as a good quantitative estimate of the redox potential. We also identify the first solvation sphere together with the metal ion as the primary active region during the electron uptake process, whereas the second solvation sphere functions as a non-reactive continuum region.

Current-density maps were calculated by the ipsocentric CTOCD-DZ/6-311G** (CTOCD-DZ=continuous transformation of origin of current density-diamagnetic zero) approach for three sets of inorganic monocycles: S(4) (2+), Se(4) (2+), S(2)N(2), P(5) (-) and As(5) (-) with 6 pi electrons; S(3)N(3) (-), S(4)N(3) (+) and S(4)N(4) (2+) with 10 pi electrons; and S(5)N(5) (+) with 14 pi electrons. Ipsocentric orbital analysis was used to partition the currents into contributions from small groups of active electrons and to interpret the contributions in terms of symmetry- and energy-based selection rules. All nine systems were found to support diatropic pi currents, reinforced by sigma circulations in P(5) (-), As(5) (-), S(3)N(3) (-), S(4)N(3) (+), S(4)N(4) (2+) and S(5)N(5) (+), but opposed by them in S(4) (2+), Se(4) (2+) and S(2)N(2). The opposition of pi and sigma effects in the four-membered rings is compatible with height profiles of calculated NICS (nucleus-independent chemical shifts).

The relative aromaticities of the three singlet benzyne isomers, 1,2-, 1,3-, and 1,4-didehydrobenzenes have been evaluated with a series of aromaticity indicators, including magnetic susceptibility anisotropies and exaltations, nucleus-independent chemical shifts (NICS), and aromatic stabilization energies (all evaluated at the DFT level), as well as valence-bond Pauling resonance energies. Most of the criteria point to the o-benzyne