Search

In this work we present the development of a novel, quantitative solubility descriptor based on the non-covalent interaction index. It is presented as a more insightful alternative to Hansen's solubility parameters and the COSMO model to assess and predict polymer solubility in different solvents. To this end, we studied the solvation behaviour as a function of the chain length of a single chain of arguably the most simple polymer, polyethylene, in anisole (solvent) and methanol (poor solvent) via molecular dynamics simulations. It was found that in anisole the solute maximized its interface with the solvent, whereas in methanol the macromolecule formed rod-like structures by folding on itself once the chain length surpassed a certain barrier. We assessed this behaviour - which can be related to solubility - quantitatively and qualitatively via well-known descriptors, namely the solvation free energy, and the solvent accessible surface area. In addition, we propose the non-covalent interaction (NCI) index as a versatile descriptor, providing information on the strength, as well as the nature, of the solute-solvent interactions, the solute's intramolecular interactions and on the solute's conformation, both qualitatively and quantitatively. Finally, as a quantitative measure for solubility, defined in this context as the solute's tendency to maximize its interactions with the solvent, we propose two new NCI-based descriptors: the relative integrated NCI density and the integrated NCI difference. The former represents the quantitative difference in solute-solvent interactions between a fully extended coil and the actual conformation during simulation and the latter the quantitative difference between the intermolecular (solute-solvent) and the intramolecular (in the solute) non-covalent interactions. The easy interpretation and calculation of these novel quantities open up the possibility of fast, reliable and insightful high-throughput screening of different (anti)solvent and solute combinations.

Linear triatomic molecules (CO 2 , N 2 O, and OCS) are scrutinized for their propensity to form perpendicular tetrel (CO 2 and OCS) or pnictogen (N 2 O) bonds with Lewis bases (dimethyl ether and trimethyl amine) as compared with their tendency to form end-on chalcogen bonds. Comparison of the IR spectra of the complexes with the corresponding monomers in cryogenic solutions in liquid argon enables to determine the stoichiometry and the nature of the complexes. In the present cases, perpendicular tetrel and pnictogen 1:1 complexes are identified mainly on the basis of the lifting of the degenerate ν 2 bending mode with the appearance of both a blue and a red shift. Van 't Hoff plots of equilibrium constants as a function of temperature lead to complexation enthalpies that, when converted to complexation energies, form the first series of experimental complexation energies on sp 1 tetrel bonds in the literature, directly comparable to quantum-chemically obtained values. Their order of magnitude corresponds with what can be expected on the basis of experimental work on halogen and chalcogen bonds and previous computational work on tetrel bonds. Both the order of magnitude and sequence are in fair agreement with both CCSD(T) and DFA calculations, certainly when taking into account the small differences in complexation energies of the different complexes (often not more than a few kJ mol -1 ) and the experimental error. It should, however, be noted that the OCS chalcogen complexes are not identified experimentally, most probably owing to entropic effects. For a given Lewis base, the stability sequence of the complexes is first successfully interpreted via a classical electrostatic quadrupole-dipole moment model, highlighting the importance of the magnitude and sign of the quadrupole moment of the Lewis acid. This approach is validated by a subsequent analysis of the molecular electrostatic potential, scrutinizing the σ and π holes, as well as the evolution in preference for chalcogen versus tetrel bonds when passing to "higher" chalcogens in agreement with the evolution of the quadrupole moment. The energy decomposition analysis gives further support to the importance/dominance of electrostatic effects, as it turns out to be the largest attractive term in all cases considered, followed by the orbital interaction and the dispersion term. The natural orbitals for chemical valence highlight the sequence of charge transfer in the orbital interaction term, which is dominated by an electron-donating effect of the N or O lone-pair(s) of the base to the central atom of the triatomics, with its value being lower than in the case of comparable halogen bonding situations. The effect is appreciably larger for TMA, in line with its much higher basicity than DME, explaining the comparable complexation energies for DME and TMA despite the much larger dipole moment for DME.

Supramolecular gels form a class of soft materials that has been heavily explored by the chemical community in the past 20 years. While a multitude of experimental techniques has demonstrated its usefulness when characterizing these materials, the potential value of computational techniques has received much less attention. This review aims to provide a complete overview of studies that employ computational tools to obtain a better fundamental understanding of the self-assembly behavior of supramolecular gels or to accelerate their development by means of prediction. As such, we hope to stimulate researchers to consider using computational tools when investigating these intriguing materials. In the concluding remarks, we address future challenges faced by the field and formulate our vision on how computational methods could help overcoming them.

The conjugated π-system in polyenes can be interrupted by electrocyclic ring-closure reactions. In this work, this 6π-electrocylization is shown by means of density functional calculations to be reversible by the application of an external mechanical pulling force at the terminal ends of the interrupted polyene chain. The test systems were constrained in a fused ring system, thus locking the orientation of three π-bonds and generally promoting 6π-electrocyclic ring-closure reactions. For several systems, the forward reaction is exergonic and the corresponding reaction barrier is comparable to those reported in the literature. The reverse reaction is triggered by an external pulling force of 2 nN (nano-Newton) or less and also becomes exergonic in all investigated polyenes under these force conditions. Moreover, it proceeds via a low reaction barrier when a pulling force of 2 nN is active, indicating that the mechanical force is an efficient stimulus for triggering ring-opening reactions. Analysis of the strain energy induced by this mechanical force confirms an optimal activation of the corresponding C-C σ-bond that breaks upon ring opening when the pulling positions are located on the polyene chain.

The halolactonization reaction is one of the most common electrophilic addition reactions to alkenes. The mechanism is generally viewed as a two-step pathway, which involves the formation of an ionic intermediate, in most cases a haliranium ion. Recently, an alternative concerted mechanism was proposed, in which the nucleophile of the reaction played a key role in the rate determining step by forming a pre-polarized complex with the alkene. This pathway was coined the nucleophile-assisted alkene activation (NAAA) mechanism. Metadynamics simulations on a series of model halolactonization reactions were used to obtain the full dynamic trajectory from reactant to product and investigate the explicit role of the halogen source and solvent molecules in the mechanism. The results in this work ratify the occasional preference of a concerted mechanism over the classic two-step transformation under specific reaction conditions. Nevertheless, as the stability of both the generated substrate cation and counter-anion increase, a transition towards the classic two-step mechanism was observed. NCI analyses on the transition states revealed that the activating role of the nucleophile is independent of the formation and stability of the intermediate. Additionally, the dynamic insights obtained from the metadynamics simulations and NCI analyses employed in this work, unveiled the presence of syn -directing noncovalent interactions, such as hydrogen bonding, between the alkenoic acid and the halogen source, which rationalized the experimentally observed diastereoselectivities. Explicit noncovalent interactions between the reactants and a protic solvent or basic additive are able to disrupt these syn -directing noncovalent interactions, affecting the diastereoselective outcome of the reaction., Competing Interests: There are no conflicts to declare., (This journal is © The Royal Society of Chemistry.)

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.)

The extension of the E = E[N, v] functional for exploring chemical reactivity in a conceptual DFT context to include external electric fields is discussed. Concentrating on the case of a homogeneous field the corresponding response functions are identified and integrated, together with the conventional response functions such as permanent dipole moment and polarizability, in an extended response function tree associated with the E = E[N, v, ε] functional. In a case study on the dihalogens F2, Cl2, Br2, I2 the sensitivity of condensed atomic charges (∂q/∂ε) is linked to the polarizability of the halogen atoms. The non-integrated (∂ρ(r)/∂ε) response function, directly related to the field induced density change, is at the basis of these features. It reveals symmetry breaking for a perpendicular field, not detectable in its atom condensed counterpart, and accounts for the induced dipole moment directly related to the molecular polarizability. The much higher sensitivity of the electronic chemical potential/electronegativity as compared to the chemical hardness is highlighted. The response of the condensed Fukui functions to a parallel electric field increases when going down in the periodic table and is interpreted in terms of the extension of the outer contours in the non-condensed Fukui function. In the case of a perpendicular field the (∂f(r)/∂ε) response function hints at stereoselectivity with a preferential side of attack which is not retrieved in its condensed form. In an application the nucleophilic attack on the carbonyl group in H2CO is discussed. Similar to the dihalogens, stereoselectivity is displayed in the Fukui function for nucleophilic attack (f+) in the case of a perpendicular electric field, and opposite to the one that would arise based on the induced density. Disentangling the expression for the evolution of the Fukui function in the presence of an electric field reveals that this difference can be traced back to local differences in the polarization or induced density between the anionic and the neutral system. This difference may be exploited, e.g. for an appropriately substituted H2CO, to generate enantioselectivity.

A novel class of substituted spiro[3.4]octanes can be accessed via a [2 + 2]-cycloaddition of dichloroketene on a readily prepared exo -methylene cyclopentane building block. This reaction sequence was found to be robust on a multigram scale and afforded a central spirocyclobutanone scaffold for carbocyclic nucleosides. The reactivity of this constrained building block was evaluated and compared to the corresponding 4'-spirocyclic furanose analogues. Density functional theory calculations were performed to support the observed selectivity in the carbonyl reduction of spirocyclobutanone building blocks. Starting from novel spirocyclic intermediates, we exemplified the preparation of an undescribed class of carbocyclic nucleoside analogues and provided a proof of concept for application as inhibitors for the protein methyltransferase target PRMT5.

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.)

Deep Eutectic Solvents (DESs) are hygroscopic liquids composed of a hydrogen bond donor (HBD) and acceptor (HBA). Their physical, chemical and electrochemical properties can be tailored to use them as solvents for different applications, i.e. electrodeposition, catalysis, extraction, etc. This can be done by changing the HBD, as well by adding water. However, the interrelated influence of H2O and HBD on the structure of the electrolyte, and on the behavior of the active species is not fully understood. In this work, we select nickel electrodeposition as a case study and we combine electrochemical techniques (cyclic voltammetry, chronoamperometry) with UV-vis spectroscopy and molecular dynamics to understand the influence of water and HBD on the electrochemical behaviour of DESs. The unique combination of these different experimental and modelling techniques provides new insights into the field. The addition of H2O changes, not only the interactions between the constituents of the liquid, but also the coordination of metal cations, which is reflected in the electrochemical performance of different DESs. More importantly, we show that, in the presence of very little (between 0.1 wt% and 2.8 wt%) and high (>4.5 wt%) water contents, DESs behave differently, and the changes in their electrochemical behavior are caused by both the complexation of metal cations and the electrolyte transport properties.

The stereoselective reduction of carbonyls is of key importance in the total synthesis of natural products and in medicinal chemistry. Nevertheless, models for rationalizing the stereoselectivity of the hydride reductions of cyclobutanones toward cyclobutanols are largely lacking, unlike cyclohexanone reductions. In order to elucidate the factors that control the stereoselectivity of these reductions, we have investigated the effect of the reaction temperature, solvent, substituent, and type of reducing agent using a synergistic experimental-computational approach. On the experimental side, the hydride reduction of 3-substituted cyclobutanones was proven to be highly selective for the formation of cis alcohol (>90%), irrespective of the size of the hydride reagent. The pronounced selectivity can be further enhanced by lowering the reaction temperature or decreasing the solvent polarity. On the computational side, density functional theory and noncovalent interaction analysis reveal that torsional strain plays a major role in the preference for the antifacial hydride approach, consistent with the Felkin-Anh model. In the presence of the benzyloxy substituent, the high selectivity for the cis isomer is also driven by repulsive electrostatic interactions in the case of a syn-facial hydride attack. The computed cis/trans ratios are in good agreement with the experimental ones and thus show the potential of computational chemistry for predicting and rationalizing the stereoselectivity of hydride reductions of cyclobutanones.

How do local reactivity descriptors, such as the Fukui function and the local spin density distribution, shape the potential energy surface (PES) associated with chemical reactions and thus govern reactivity trends and regioselective preferences? This is the question that is addressed here through a qualitative valence bond (VB) analysis. We demonstrate that common density functional theory (DFT)-based local reactivity descriptors can essentially be regarded-in one way or another-as indirect measures of delocalization, i.e., resonance stabilization, of the reactants within VB theory. The inherent connection between (spatial) delocalization and (energetic) resonance stabilization embedded in VB theory provides a natural and elegant framework for analyzing and comprehending the impact of individual local reactivity descriptors on the global PES. Our analysis provides new insights into the role played by local reactivity descriptors and illustrates under which conditions they can sometimes fail to predict reactivity trends and regioselective preferences, e.g., in the case of ambident reactivity. This treatment constitutes a first step toward a unification of VB theory and conceptual DFT .

Expanded porphyrins provide a versatile route to molecular switching devices due to their ability to shift between several π-conjugation topologies encoding distinct properties. DFT remains the workhorse for modeling such extended macrocycles, when taking into account their size and huge conformational flexibility. Nevertheless, the stability of Hückel and Möbius conformers depends on a complex interplay of different factors, such as hydrogen bonding, π···π stacking, steric effects, ring strain, and electron delocalization. As a consequence, the selection of an exchange-correlation functional for describing the energy profile of topological switches is very difficult. For these reasons, we have examined the performance of a variety of wave function methods and density functionals for describing the thermochemistry and kinetics of topology interconversions across a wide range of macrocycles. Especially for hexa- and heptaphyrins, the Möbius structures have a stronger degree of static correlation than the Hückel and twisted-Hückel structures, and as a result the relative energies of singly twisted structures are a challenging test for electronic structure methods. Comparison of limited orbital space full CI calculations with CCSD(T) calculations within the same active spaces shows that post-CCSD(T) correlation contributions to relative energies are very minor. At the same time, relative energies are weakly sensitive to further basis set expansion, as proven by the minor energy differences between the extrapolated MP2/CBS energies estimated from cc-pV{T,Q}Z, diffuse-augmented heavy-aug-cc-pV{T,Q}Z and explicitly correlated MP2-F12/cc-pVDZ-F12 calculations. Hence, our CCSD(T) reference values are reasonably well-converged in both 1-particle and n -particle spaces. While conventional MP2 and MP3 yield very poor results, SCS-MP2 and particularly SOS-MP2 and SCS-MP3 agree to better than 1 kcal mol -1 with the CCSD(T) relative energies. Regarding DFT methods, the range-separated double hybrids, such as ωB97M(2) and B2GP-PLYP, outperform other functionals with RMSDs of 0.6 and 0.8 kcal mol -1 , respectively. While the original DSD-PBEP86 double hybrid performs fairly poorly for these extended π-systems, the errors drop down to 1.9 kcal mol -1 for the revised revDOD-PBEP86-NL, which eliminates the same-spin correlation energy. Minnesota meta-GGA functionals with high fractions of exact exchange (M06-2X and M08-HX) also perform reasonably well, outperforming more robust and significantly less empirically parametrized functionals like SCAN0-D3.

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.)

Benziporphyrins are versatile macrocycles exhibiting aromaticity switching behaviors. The existence of both Hückel and Möbius (anti)aromaticity has been reported in these systems, whose validity is respectively governed by the [4n + 2] and [4n] π-electron rule on the macrocyclic pathway. Despite the experimental evidence on the floppiness of benziporphyrins, the switching mechanism between Hückel and Möbius structures is still not clear, as well as the factors influencing the stability of the different π-conjugation topologies. For these reasons, we performed a systematic study on A,D-di-p-benzihexaphyrins(1.1.1.1.1.1) with two redox states corresponding to [28] and [30] π-electron conjugation pathways. Whereas benzi[28]hexaphyrin obeys Möbius aromaticity, benzi[30]hexaphyrin follows Hückel aromaticity. The dynamic interconversion between Möbius and Hückel aromaticity is investigated through the rotation of a phenylene ring, which acts as the topology selector. Further analyses of the energy profiles using energy decomposition and information-theoretic approaches provide new insights into conformational stability, aromaticity and antiaromaticity for these species. Strong and opposite cross correlations between aromaticity indexes and information-theoretic quantities were found for the two macrocyclic systems with opposite global aromaticity and antiaromaticity behaviors. These results indicate that Hückel and Möbius aromaticity and antiaromaticity, though qualitatively different, are closely related and can be interchanged, and information-theoretic quantities provide a novel understanding about their relevance. Our present results should provide in-depth insights to appreciate the nature and origin about Möbius (anti)aromaticity and its close relationship with Hückel (anti)aromaticity.

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.)

In this publication, we propose a new set of reactivity/selectivity descriptors, derived within a Rayleigh-Schrödinger perturbation theory framework, for chemical systems undergoing an electrostatic (point-charge) perturbation. From the electron density polarization at first order, qualitative insight on reactivity is retrieved, while more quantitative information (noteworthy selectivity) can be obtained from either the second-order energy response or the number of shifted electrons under perturbation. Noteworthily, only a small number of excitations contribute significantly to the overall responses to perturbation, suggesting chemical reactivity could be foreseen by a careful scrutiny of the electron density reorganization upon excitation.

In our effort to implement the mechanical force used to activate single molecules in mechanochemistry in the context of conceptual density functional theory, we present a theoretical investigation of strained alkynes for rationalizing structural trends as well as the reactivity of cyclic alkynes that are of great importance in in vivo click reactions. The strain on the triple bond in cyclic alkynes is modeled by angular constraints in a 2-butyne fragment and the corresponding bending force is calculated by means of an extended COGEF (constrained geometries simulate external forces) model. In general, the force required to bend the triple bond is smaller with electron-withdrawing groups on the propargylic C-atom, which elegantly results in smaller angles around the triple bond in cyclic alkynes with such substitution pattern. By means of conceptual DFT descriptors, the electrophilic and nucleophilic character of bent triple bonds was investigated revealing moderate activation for small distortions from the linear geometry (0° to 15°) and a drastically more reactive π-space if the triple bond is bent further. This analysis of the intrinsic reactivity of the triple bond is in line with experimental observations, explaining the reactive nature of cyclooctynes and cycloheptynes, whereas larger cyclic systems do not drastically activate the triple bond., Competing Interests: There are no conflicts of interest to declare., (This journal is © The Royal Society of Chemistry.)

This review sets out to understand the reactivity of diradicals and how that may differ from monoradicals. In the first part of the review, we delineate the electronic structure of a diradical with its two degenerate or nearly degenerate molecular orbitals, occupied by two electrons. A classification of diradicals based on whether or not the two SOMOs can be located on different sites of the molecule is useful in determining the ground state spin. Important is a delocalized to localized orbital transformation that interchanges "closed-shell" to "open-shell" descriptions. The resulting duality is useful in understanding the dual reactivity of singlet diradicals. In the second part of the review, we examine, with a consistent level of theory, activation energies of prototypical radical reactions (dimerization, hydrogen abstraction, and addition to ethylene) for representative organic diradicals and diradicaloids in their two lowest spin states. Differences and similarities in reactivity of diradicals vs monoradicals, based on either a localized or delocalized view, whichever is suitable, are then discussed. The last part of this review begins with an extensive, comparative, and critical survey of available measures of diradical character and ends with an analysis of the consequences of diradical character for selected diradicaloids.

Alchemical derivatives have been used in recent years to obtain essentially qualitative information about transformations in which the number of electrons is unchanged. Within the context of Conceptual DFT, we present a systematic approach for combining changes in both the number of electrons and the nuclear charge so that for example one can navigate from one neutral atom to another. A general formalism is presented for transformations involving changes both in or , where Parr's parabolic approach for the dependence is considered as one particular case and the ensemble description in the 0 K limit as the second case. The B3LYP functional in its CAMB3LYP version combined with the aug-cc-pCVQZ basis has been chosen to perform Coupled Perturbed Kohn Sham calculations of the alchemical derivatives. The monotonic behaviour of the alchemical potential is scrutinised. The order of magnitude analysis of the derivatives preludes convergence at third order. These results are injected in two strategies for obtaining transmutation energies from neutral atoms to a neighbouring neutral atom: one road moving along the diagonal, the other one walking along a pure alchemical road after ionisation or electron attachment. Roads involving the anion of the reference atom are much less successful than those involving its cation. The transmutation energy for the cationic pathway displays chemical accuracy when the procedure is carried at third order in . The difficulties inherent to an accurate description of the anion and its response functions are responsible for the shortcomings along the anionic paths. As a direct application Ionization Energies (IE) and Electron Affinities (EA) are evaluated showing an almost perfect agreement with the direct evaluation and a difference with experimental values less than 0.5 eV for the IE. For the first EA reasonable agreement is obtained with direct and experimental values whereas the second EAs for atoms with stable mono-anions show a remarkable agreement with literature data. Besides proof of concept that with the information content of an atom one can get accurate energetics of its neighbours, the results indicate that alchemical derivatives are capable to yield quantitative information when navigating through Chemical Compound Space.

Proline is often found as a turn inducer in peptide or protein domains. Exploitation of its restricted conformational freedom led to the development of the d-Pro-l-Pro (corresponding to (R)-Pro-(S)-Pro) segment as a "templating" unit, frequently used in the design of β-hairpin peptidomimetics, in which conformational stability is, however, inherently linked to the cis-trans isomerization of the prolyl amide bonds. In this context, the stereoelectronic properties of the CF 3 group can aid in conformational control. Herein, the impact of α-trifluoromethylated proline analogues is examined for the design of enhanced β-turn inducers. A theoretical conformational study permitted the dipeptide (R)-Pro-(R)-TfmOxa (TfmOxa: 2-trifluoromethyloxazolidine-2-carboxylic acid) to be selected as a template with an increased trans-cis rotational energy barrier. NMR spectroscopic analysis of the Ac-(R)-Pro-(R)-TfmOxa-(S)-Val-OtBu β-turn model, obtained through an original synthetic pathway, validated the prevalence of a major trans-trans conformer and indicated the presence of an internal hydrogen bond. Altogether, it was shown that the (R)-Pro-(R)-TfmOxa template fulfilled all crucial β-turn-inducer criteria., (© 2019 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim.)

In this work, we report a systematic study on the global and local aromaticity of acenes using a series of model structures from 2-acene to 11-acene. A recently developed ansatz, an information-theoretic approach coached into density functional reactivity theory has been employed, which essentially provides different density functionals characterizing the molecular electron density distribution. Based on the correlation analysis of six conventional aromaticity indices with eight information-theoretic quantities, we examined the aromaticity of acenes from both global and local perspectives. From the global aromaticity viewpoint, our results suggest that different descriptors based on various physicochemical properties are intrinsically dependent. A novel laminated feature ruling local aromaticity of acenes has been unveiled, from which we found that the distance from the terminal rings plays the critical role. Based on the shape of the correlation plots between the conventional aromaticity indices and information-theoretic quantities, the latter could be separated into three subgroups. The seemingly contradictory results from global and local aromaticity perspectives not only present us the uniqueness of the acene systems but all demonstrate the effectiveness of the information-theoretic approach from density functional reactivity theory. Besides strengthening the validity of a series of new aromaticity descriptors, our results should lead to more clear insights into the chemical significance of the information-theoretic quantities.

Deep eutectic solvents (DESs) are a new class of solvents with wider potential window than that of water and high electrochemical stability, making them potential candidates for a wide range of electrochemical systems. However, due to the hygroscopic nature of DESs, the presence of latent water is unavoidable. Therefore, understanding the interfacial structure and the electrosorption and distribution of residual water at the electrified interface is of great importance for the use of these solvents in electrochemical systems. Using atomistic molecular dynamics, we explore the electrosorption and distribution of different amounts of water in 1 : 2 choline chloride-urea DES (Reline) at the electrified graphene interface. We found that both the water distribution and the interfacial structure are sensitive to the electrification of the graphene electrode. As a result, it is found that for moderately charged electrodes, water shows a preferential asymmetric adsorption in the vicinity of the positively charged electrode, partly due to strong intermolecular interactions with anions through hydrogen bonds. In contrast, for highly charged electrodes, water adsorbs at both electrodes due to a strongly enhanced external electrostatic interaction between the electrodes and the water dipoles., Competing Interests: There are no conflicts to declare., (This journal is © The Royal Society of Chemistry.)

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 performance of density functionals and wavefunction methods for describing the thermodynamics and kinetics of hydride reductions of 2-substituted cyclohexanones has been evaluated for the first time. A variety of exchange correlation functionals ranging from generalized gradient approximations to double hybrids have been tested and their performance to describe the facial selectivity of hydride reductions of cyclohexanones has been carefully assessed relative to the CCSD(T) method. Among the tested methods, an approach in which single-point energy calculations using the double hybrid B2PLYP-D3 functional on ωB97X-D optimized geometries provides the most accurate transition state energies for these kinetically-controlled reactions. Moreover, the role of torsional strain, temperature, solvation, noncovalent interactions on the stereoselectivity of these reductions was elucidated. Our results indicate a prominent role of the substituent on the cis / trans ratios driven by the delicate interplay between torsional strain and dispersion interactions., Competing Interests: The authors declare no conflict of interest.

This study addresses one of the most important drawbacks inherently related to molecular searches in chemical compound space by greedy algorithms such as Best First Search and Genetic Algorithm, i.e., the large computational cost required to optimize one or more quantum-chemical properties. Significant speed-ups are obtained by initial property screening via predictive techniques starting already from very small databases. It is shown that the attainable acceleration depends heavily on the molecular properties, the predictive model, the molecular descriptor, and the current size of the database. We discuss the implementation and performance of predictive techniques in molecular searches based on a fixed molecular framework with a selection of sites to be filled with groups from a chemical fragment library. It is shown that for some properties speed-ups of a factor of 5 to even 20 can be obtained, while inverse design procedures on more complex properties still reach speed-ups of a factor of 2 without losing performance.

This study investigates the rich mechanistic landscape of the iconic electrophilic aromatic substitution (EAS) reaction class, in the gas phase, in solvents, and under stimulation by oriented external electric fields. The study uses DFT calculations, complemented by a qualitative valence bond (VB) perspective. We construct a comprehensive and unifying framework that elucidates the many surprising mechanistic features, uncovered in recent years, of this class of reactions. For example, one of the puzzling issues which have attracted significant interest recently is the finding of a variety of concerted mechanisms that do not involve the formation of σ-complex intermediates, in apparent contradiction to the generally accepted textbook mechanism. Our VB modeling elucidates the existence of both the concerted and stepwise mechanisms and uncovers the root causes and necessary conditions for the appearance of these intermediates. Furthermore, our VB analysis offers insight into the potential applications of external electric fields as smart, green, and selective catalysts, which can control at will reaction rates, as well as mechanistic crossovers, for this class of reactions. Finally, we highlight how understanding of the electric fields effect on the EAS reaction could lead to the formulation of guiding principles for the design of improved heterogeneous catalysts. Overall, our analysis underscores the powerful synergy offered by combining molecular orbital and VB theory to tackle interesting and challenging mechanistic questions in chemistry.

In molecular mechanochemistry, the chemical properties of a reactant system are modified through the absorption of mechanical energy at the single molecular level, as opposed to thermochemistry (heat), electrochemistry (electrical current) and photochemistry (light). In this proof-of-concept paper, we present a theoretical approach to rationalize and predict the change in chemical properties and concepts when an external stretching force, Fext, is applied to a chemical bond. By considering perturbations to the total molecular energy due to Fext, a series of mechanochemical response indices are obtained within the framework of conceptual density functional theory (DFT). Trends in these indices monitoring the change in hardness and electrophilicity among others, were rationalized for a number of diatomic molecules in terms of the ground-state geometry and the frontier molecular orbitals of the molecule. Finally, we present a set of rules that can be easily adopted for polyatomic molecules to predict the changing chemical reactivity of single molecules when subjected to an external force.

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.)

Originated from the cyclic delocalization of electrons resulting in extra stability and instability, aromaticity and antiaromaticity are important chemical concepts whose appreciation and quantification are still much of recent interest in the literature. Employing information-theoretic quantities can provide us with more insights and better understanding about them, as we have previously demonstrated. In this work, we examine the triplet-state aromaticity and antiaromaticity, which are governed by Baird's 4 n rule, instead of Hückel's 4 n + 2 rule for the singlet state. To this end, we have made use of 4 different aromaticity indexes and 8 information-theoretic quantities, examined a total of 22 substituted fulvene derivatives, and compared the results both in singlet and triplet states. It is found that cross-correlations of these two categories of molecular property descriptors enable us to better understand the nature and propensity of aromaticity and antiaromaticity for the triplet state. Our results have not only demonstrated the existence and validity of Baird's rule but also shown that Hückel's rule and Baird's rule indeed share the same theoretical foundation because with these cross-correlation patterns we are able to distinguish them from each other simultaneously in both singlet and triplet states. Our results should provide new insights into the nature of aromaticity and antiaromaticity in the triplet state and pave the road toward new ways to quantify this pair of important chemical concepts., Competing Interests: The authors declare no competing financial interest.

In recent years, a considerable interest has grown in the design of molecular nanowires with an increasing conductance with length. The development of such nanowires is highly desirable because they could play an important role in future molecular-scale circuitry. Whereas the first experimental observation of this nonclassical behavior still has to be realized, a growing number of candidate wires have been proposed theoretically. In this Letter, we point out that all the wires with an anti-Ohmic increasing conductance with length proposed so far share a common characteristic: their diradical character increases with length. The conceptual connection between diradical character and conductance enables a systematic design of such anti-Ohmic wires and explains the difficulty in their syntheses. A strategy is proposed to balance the stability and conductance so that this nonclassical phenomenon can be observed.

Green, stable, and wide electrochemical window deep eutectic solvents (DESs) are ideal candidates for electrochemical systems. However, despite several studies of their bulk properties, their structure and properties under electrified confinement have barely been investigated, which has hindered widespread use of these solvents in electrochemical applications. In this Letter, we explore the electrical double layer structure of 1:2 choline chloride-urea (Reline), with a particular focus on the electrosorption of the hydrogen bond donor on a graphene electrode using atomistic molecular dynamics simulations. We discovered that the interface is composed of a mixed layer of urea and counterions followed by a mixed charged clustered structure of all of the Reline components. This interfacial structuring is strongly dependent on the balance between intermolecular interactions and surface polarization. These results provide new insights into the electrical double layer structure of a new generation of electrolytes whose interfacial structure can be tuned at the molecular level.

N,C,N-Chelated pnictinidenes ArE [where E = As, Sb or Bi; Ar = 2,6-(tBuN[double bond, length as m-dash]CH)2C6H3] were used as ligands for the coordination of various gold(i) complexes. Thus, the reaction of ArE with [AuCl(Me2S)] gave complexes [AuCl(ArE)] [where E = As (1) or Sb (2)] that exhibited only limited stability in solution. By contrast, the reaction of ArBi with [AuCl(Me2S)] led to the immediate deposition of gold metal and the oxidation of the bismuth atom giving ArBiCl2. The treatment of a tetrameric gold alkynyl complex [Au(C[triple bond, length as m-dash]CPh)]4 with ArAs and ArSb gave ionic compounds [Au(ArAs)2]+[Au2(C[triple bond, length as m-dash]CPh)3]- [denoted as 3+[Au2(C[triple bond, length as m-dash]CPh)3]-] and [Au(ArSb)2]+[Au(C[triple bond, length as m-dash]CPh)2]- [denoted as 4+[Au(C[triple bond, length as m-dash]CPh)2]-], respectively. Finally, the reaction of ArE with the carbene gold(i) complex [Au(IPr)(MeCN)]+[BF4]- [where IPr = 1,3-bis(2,6-diisopropylphenyl)imidazolin-2-ylidene, MeCN = acetonitrile] produced ionic complexes [Au(IPr)(ArE)]+[BF4]- [for cations: E = As (5+), Sb (6+) or Bi (7+)]. All complexes were characterized using 1H and 13C NMR, high mass accuracy electrospray ionization mass spectrometry (ESI-MS), IR and Raman spectroscopy and (except for 1) by single-crystal X-ray diffraction analysis. Furthermore, the structure and bonding of both neutral and ionic complexes with different coordination patterns have also been investigated in detail using a Density Functional Theory (DFT) computational approach.

With their versatile molecular topology and aromaticity, porphyrinoid systems combine remarkable chemistry with interesting photophysical properties and nonlinear optical properties. Hence, the field of application of porphyrinoids is very broad ranging from near-infrared dyes to opto-electronic materials. From previous experimental studies, aromaticity emerges as an important concept in determining the photophysical properties and two-photon absorption cross sections of porphyrinoids. Despite a considerable number of studies on porphyrinoids, few investigate the relationship between aromaticity, UV/vis absorption spectra and nonlinear properties. To assess such structure-property relationships, we performed a computational study focusing on a series of Hückel porphyrinoids to: (i) assess their (anti)aromatic character; (ii) determine the fingerprints of aromaticity on the UV/vis spectra; (iii) evaluate the role of aromaticity on the NLO properties. Using an extensive set of aromaticity descriptors based on energetic, magnetic, structural, reactivity and electronic criteria, the aromaticity of [4 n +2] π-electron porphyrinoids was evidenced as was the antiaromaticity for [4 n ] π-electron systems. In agreement with previous studies, the absorption spectra of aromatic systems display more intense B and Q bands in comparison to their antiaromatic homologues. The nature of these absorption bands was analyzed in detail in terms of polarization, intensity, splitting and composition. Finally, quantities such as the average polarizability and its anisotropy were found to be larger in aromatic systems, whereas first and second hyperpolarizability are influenced by the interplay between aromaticity, planarity and molecular symmetry. To conclude, aromaticity dictates the photophysical properties in porphyrinoids, whereas it is not the only factor determining the magnitude of NLO properties.

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.)

The chalcogen bond has been acknowledged as an influential noncovalent interaction (NCI) between an electron-deficient chalcogen (donor) and a Lewis base (acceptor). This work explores the main features of chalcogen bonding through a large-scale computational study on a series of donors and acceptors spanning a wide range in strength and character of this type of bond: (benzo)chalcogenadiazoles (with Ch = Te/Se/S) versus halides and neutral Lewis bases with O, N, and C as donor atoms. We start from Pearson's hard and soft acids and bases (HSAB) principle, where the hard nature of the chalcogen bond is quantified through the molecular electrostatic potential and the soft nature through the Fukui function. The σ-holes are more pronounced when going down in the periodic table and their directionality matches the structural orientation of donors and acceptors in the complexes. The Fukui functions point toward an n→σ*-type interaction. The initial conjectures are further scrutinized using quantum mechanical methods, mostly relating to the systems' electron density. A Ziegler-Rauk energy decomposition analysis shows that electrostatics plays a distinctly larger role for the soft halides than for the hard, uncharged acceptors, associated with the softness matching within the HSAB principle. The natural orbital for chemical valence analysis confirms the n→σ* electron donation mechanism. Finally, the electron density and local density energy at the bond critical point in the quantum theory of atoms in molecules study and the position of the spikes in the reduced density gradient versus density plot in the NCI theory situate the chalcogen bond in the same range as strong hydrogen bonds. © 2017 Wiley Periodicals, Inc., (© 2017 Wiley Periodicals, Inc.)

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.)

With the idea of using alchemical derivatives to explore in an efficient, computer- and cost-effective way Chemical Space was launched several years ago. In the context of Conceptual DFT response functions, these energies vs nuclear charge derivatives permit the estimatation of the energy of transmutants of a given starting or reference molecule showing different nuclear compositions. After an explorative study on small and planar molecules ( Balawender et al. J. Chem. Theory Comput. 2013 , 9 , 5327 ) by the present authors of this paper, the present study fully exploits the computational advantages of the alchemical derivatives in larger three-dimensional systems. Starting from a single reference calculation on C 60 , the complete BN substitution pattern, from single substituted C 58 BN via the belt (C 20 (BN) 20 and the ball C 12 (BN) 24 structures to the fully substituted (BN) 30 , is explored. Successive and simultaneous substitution strategies are followed and compared, indicating that both techniques yield identical results up to 13 substitutions but that for higher substitutions the simultaneous approach needs to be taken. Due to the cost-efficiency of the algorithm this path can indeed be followed as opposed to earlier work in the literature where for each step a full SCF calculation was at stake leading to prohibitively large computational demands for adopting the simultaneous approach. Previously formulated rules governing the substitution pattern by Kar and co-workers are scrutinized in this context and reformulated giving chemical insight in the gradual substitution process and the relative energies of the isomers. In its present form the method offers an interesting venue to study BN substitution patterns in higher fullerenes and graphene and in general paves the way for more efficient exploration of the Chemical Space.

Expanded porphyrins are flexible enough to switch between different π-conjugation topologies, namely Möbius, Hückel and twisted-Hückel, each with distinct electronic properties and aromaticity. Since these switches can be induced by different external stimuli, expanded porphyrins represent a promising platform to develop a novel type of molecular switch for molecular electronic devices. In this work, the feasibility of conductance switches based on topology and/or aromaticity changes in expanded porphyrins is assessed for the first time. In particular, the electron transport properties of penta-, hexa- and heptaphyrins with different π-conjugation topologies and aromaticity were carefully investigated using the nonequilibrium Green's function formalism in combination with density functional theory for various configurations of the gold contacts. Our results highlight the importance of the macrocyclic aromaticity and connectivity and, to a lesser extent, the molecular topology, in determining the transmission functions and local currents. When the electrodes are connected along the longitudinal axis of the macrocycle, we found that aromaticity of Hückel expanded porphyrins increases single-molecule junction conductance, contrary to the negative relationship between conductance and aromaticity found in single five-membered rings. For this particular connectivity, antiaromatic Hückel structures with [4n] π-electrons exhibit a sharp reduction in transmission near the Fermi level due to destructive quantum interference between the HOMO and LUMO. Belt-shaped Möbius aromatic structures exhibit a lower conductance as compared to the Hückel aromatic structures, and the current flow avoids the molecular twist. Importantly, we show that expanded porphyrins, upon redox and topology interconversions, could act as efficient three-level molecular switches with high ON/OFF ratio, up to 10 3 at low bias voltage.

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.)

Using a set of state-of-the-art quantum chemical techniques we scrutinized the characteristically different reactivity of frustrated and classical Lewis pairs towards molecular hydrogen. The mechanisms and reaction profiles computed for the H 2 splitting reaction of various Lewis pairs are in good agreement with the experimentally observed feasibility of H 2 activation. More importantly, the analysis of activation parameters unambiguously revealed the existence of two reaction pathways through a low-energy and a high-energy transition state. An exhaustive scrutiny of these transition states, including their stability, geometry and electronic structure, reflects that the electronic rearrangement in low-energy transition states is fundamentally different from that of high-energy transition states. Our findings reveal that the widespread consensus mechanism of H 2 splitting characterizes activation processes corresponding to high-energy transition states and, accordingly, is not operative for H 2 -activating systems. One of the criteria of H 2 -activation, actually, is the availability of a low-energy transition state that represents a different H 2 splitting mechanism, in which the electrostatic field generated in the cavity of Lewis pair plays a critical role: to induce a strong polarization of H 2 that facilities an efficient end-on acid-H 2 interaction and to stabilize the charge separated "H + -H - " moiety in the transition state.

In this work, electrical currents through nanographenes, a class of alternant hydrocarbons also known as polycyclic aromatic hydrocarbons, in molecular junctions under small bias are explored. We illustrate that when the π-current dominates, that is, when no quantum interference takes place, the current prefers the direction of the shortest bond (the bond with the highest double bond character) upon entering the molecule from the contacts. As such, the idea of electrons propagating through double bonds from contact to contact, originating from the curly arrow drawings used in a previously established selection rule for transmission, seems to be more deeply rooted in the actual physical process of electron transport than previously anticipated. Furthermore, this work confirms that the σ-current behaves completely differently than the π-current. When this type of current becomes important, that is, when quantum interference takes place, the current generally prefers the shortest path from contact to contact, irrespective of the length of the bonds constituting this path, in accordance with the strong distance dependency of the σ-current. Finally, it is demonstrated that keto groups (and cross-conjugating groups in general) can be used to seal off parts of the molecule for the current. No current flows through the sealed off part of the molecule under small bias and it does not influence the transmission spectrum of the considered system., (© 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.)

This corrects the article DOI: 10.1038/srep42494.

Even though the concept of aromaticity and antiaromaticity is extremely important and widely used, there still exist lots of controversies in the literature, which are believed to be originated from the fact that there are so many aromatic types discovered and at the same time there are many aromaticity indexes proposed. In this work, using seven series of substituted fulvene derivatives as an example and with the information-theoretic approach in density functional reactivity theory, we examine these concepts from a different perspective. We investigate the changing patterns of Shannon entropy, Fisher information, Ghosh-Berkowitz-Parr entropy, information gain, Onicescu information energy, and relative Renyi entropy on the ring carbon atoms of these systems. Meanwhile, we also consider variation trends of four representative kinds of aromaticity indexes such as FLU, HOMA, ASE and NICS. Statistical analyses among these quantities show that with the same ring structure of the derivatives, both information-theoretic quantities and aromaticity indexes obey the same changing pattern, which are valid across all seven systems studied. However, cross correlations between these two sets of quantities yield two completely opposite patterns. These ring-structure dependent correlations are in good agreement with Hückel's 4n + 2 rule of aromaticity and 4n rule of antiaromaticity. Our results should provide a novel and complementary viewpoint on how aromaticity and antiaromaticity should be appreciated and categorized. More studies are in progress to further our understanding about the matter.

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.)