Search

Conceptual Density Functional Theory (CDFT) has been extended beyond its traditional role in elucidating chemical reactivity to the development of density functional theory methods, e.g., the investigation of the delocalization error. This delocalization error causes the dependence of the energy on the number of electrons ( N ) to deviate from its exact piecewise linear behavior, an error which is the basis of many well-known limitations of commonly used density-functional approximations (DFAs). Following our previous work on the analytical hardness η ± for pure functionals, we extend its application to hybrid and range-separated functionals. A comparison is made between the analytical hardness and the slope of the delocalization function introduced by Hait and Head-Gordon. Our results show that there is a linear relationship between its slope and the analytical hardness. An approximate scheme is presented to construct the energy vs N curve without fractional occupation number calculations. The extension to densities is discussed.

Until quite recently, Conceptual DFT (CDFT) was mainly based on the energy functional, E [ N , v ], where the number of electrons N and the external potential v are state variables. One of the strengths of CDFT, however, is the ease with which additional and/or different state variables can be incorporated. Here, the incorporation of new variables-namely temperature and external fields-is discussed, outlining the motivation for these extensions, sketching their theoretical/computational context, and presenting some elucidative examples. Using the Grand Canonical Ensemble, finite temperature can be introduced, ameliorating the N -differentiability problem and leading to new well-behaved chemical reactivity concepts. Incorporating temperature as an additional fundamental variable enables us to formulate a novel suite of "thermodynamic" reactivity descriptors, including important concepts like the heat capacity of electronic systems. The mathematical structure underpinning the set of (well-behaved) finite-temperature reactivity indices can guide the formulation of plausible definitions for local analogs of global descriptors. This endeavor is especially significant in the case of local hardness, a concept that has remained elusive since the inception of CDFT. The ever-increasing portfolio of experimental reaction conditions to creatively synthesize new molecules, needs the introduction of various external "fields" like electric and magnetic fields ( ε and B ), mechanical forces ( F ) and pressure ( P ) to describe the state of the chemical system. The conventional energy functional can be expressed in a general form, E [ N , v , X ], where X denotes the "field". Response functions to changes in the field can then be defined in analogy to classical thermodynamics. The electric field results display a case of a field-induced selectivity in a reaction channel of the Fukui function. Remarkably, atomic electronegativity and hardness in magnetic fields display a piecewise behavior in magnetic fields, associated to configurational jumps upon increasing field strength. The overall compression of their ranges for stronger fields may be insightful when investigating chemistry in extremely high fields. The electronegativity and hardness of diatomics under mechanical force can be traced back to changes in equilibrium distances in the neutral, cationic and anionic state, parallel with the evolution of an intrinsic atomic volume under pressure., Competing Interests: There are no conflicts to declare., (This journal is © The Royal Society of Chemistry.)

While high-pressure chemistry has a well-established history, methods to simulate pressure at the single-molecule level have been somewhat lacking. The current work aims at comparing two static models (XP-PCM and GOSTSHYP) to apply isotropic pressure to single molecules, focusing on the equilibrium bond length and electric dipole moment of diatomic molecules. Numerical challenges arising in the potential energy surface using the XP-PCM method were examined, and a pragmatic approach was followed to mitigate these. The definition of the cavity was scrutinized, and two approaches to retrieve the isotropic character that could potentially be lost when using the standard methodology were suggested. Subsequently, equilibrium bond lengths under pressure were evaluated, showing reasonable agreement between GOSTSHYP and XP-PCM, but some discrepancies persist. A Taylor series analysis introduced elsewhere was then applied to rationalize the observed trends in terms of the bond surface. Finally, the dipole moment was shown to be highly sensitive to the cavity definition, and qualitative agreement necessitates the use of our adapted procedure.

The extension of conceptual density-functional theory (conceptual DFT) to include external electromagnetic fields in chemical systems is utilised to investigate the effects of strong magnetic fields on the electronic charge distribution and its consequences on the reactivity of π-systems. Formaldehyde, H 2 CO, is considered as a prototypical example and current-density-functional theory (current-DFT) calculations are used to evaluate the electric dipole moment together with two principal local conceptual DFT descriptors, the electron density and the Fukui functions, which provide insight into how H 2 CO behaves chemically in a magnetic field. In particular, the symmetry properties of these quantities are analysed on the basis of group, representation, and corepresentation theories using a recently developed automatic program for symbolic symmetry analysis, QSYM 2 . This allows us to leverage the simple symmetry constraints on the macroscopic electric dipole moment components to make profound predictions on the more nuanced symmetry transformation properties of the microscopic frontier molecular orbitals (MOs), electron densities, and Fukui functions. This is especially useful for complex-valued MOs in magnetic fields whose detailed symmetry analyses lead us to define the new concepts of modular and phasal symmetry breaking . Through these concepts, the deep connection between the vanishing constraints on the electric dipole moment components and the symmetry of electron densities and Fukui functions can be formalised, and the inability of the magnetic field in all three principal orientations considered to induce asymmetry with respect to the molecular plane of H 2 CO can be understood from a molecular perspective. Furthermore, the detailed forms of the Fukui functions reveal a remarkable reversal in the direction of the dipole moment along the CO bond in the presence of a parallel or perpendicular magnetic field, the origin of which can be attributed to the mixing between the frontier MOs due to their subduced symmetries in magnetic fields. The findings in this work are also discussed in the wider context of a long-standing debate on the possibility to create enantioselectivity by external fields.

Background & Aim: In hospitalized patients, daily protein intake remains far below WHO requirements for healthy adults (0.8 g·kg -1 ·d -1 ) as well as ESPEN guidelines for patients (1.2-1.5 g·kg -1 ·d -1 ). Providing access to a pre-sleep protein dense snack between dinner and going to bed may serve as a great opportunity to increase daily energy and protein intake in hospitalized patients. However, it remains to be assessed whether protein provision prior to sleep effectively increases protein intake, or may reduce food intake throughout the remainder of the day(s). The present study evaluated the impact of giving access to a pre-sleep snack on daily energy and protein intake in patients throughout their hospitalization., Methods: Patients admitted to the surgical wards of the Maastricht University Medical Centre+ were randomly allocated to usual care (n = 51) or given access to a pre-sleep snack (n = 50). The pre-sleep snack consisted of 103 g cheese cubes (30 g protein) provided between 7:30 and 9:30 PM, prior to sleep. All food provided and all food consumed was weighed and recorded throughout (2-7 days) hospitalization. Daily energy and protein intake and distribution were calculated. Data were analyzed by independent T-Tests with P < 0.05 considered as statistically significant., Results: Daily energy intake was higher in the pre-sleep group (1353 ± 424 kcal d -1 ) when compared to the usual care group (1190 ± 402 kcal·d -1 ; P = 0.049). Providing patients access to a pre-sleep snack resulted in a 17% (11 ± 9 g) higher daily protein intake (0.81 ± 0.29 g·kg -1 ·d -1 ) when compared to the usual care group (0.69 ± 0.28 g·kg -1 ·d -1 ; P = 0.045). Protein intake at breakfast, lunch, and dinner did not differ between the pre-sleep and usual care groups (all P > 0.05)., Conclusion: Providing access to a pre-sleep protein snack, in the form of protein dense food items such as cheese, represents an effective dietary strategy to increase daily energy and protein intake in hospitalized patients. Patients consuming pre-sleep protein snacks do not compensate by lowering energy or protein intake throughout the remainder of the days. Pre-sleep protein dense food provision should be implemented in hospital food logistics to improve the nutritional intake of patients., Trial Register No: NL8507 (https://trialsearch.who.int/)., (Copyright © 2024 The Authors. Published by Elsevier Ltd.. All rights reserved.)

A combined computational and interpretational DFT study is performed to investigate the regioselectivity of B 18 and B 30 leapfrog boron sheets upon reaction with XH 3 -type electrophiles and nucleophiles (X = N, P, As, B, Al). The M062X, B3LYP, and B3LYP-D3 functionals are used combined with the 6-31+G(d,p) basis. The molecular electrostatic potential (MEP), Fukui functions, and the dual descriptor are employed to predict the local reactivity of B 18 and B 30 . Our results reveal that both clusters are hard and prefer to react with hard bases and acids, such as NH 3 and BH 3 . Further, these leapfrog B 6n clusters can play the role of catalysts as they break B-H and Al-H bonds of BH 3 and AlH 3 in s -BH 3 -B 6n and s -AlH 3 -B 6 n complexes, respectively. Leapfrog B 6 n -XH 3 complexes (X = B and Al) can be considered as an interaction between two electron-deficient systems. Therefore, the chemical reactivity between these systems cannot be interpreted in terms of the Hard-Soft-Acid-Base principle.

In the context of the growing impact of conceptual density functional theory (DFT) as one of the most successful chemical reactivity theories, response functions up to second order have now been widely applied; in recent years, among others, particular attention has been focused on the linear response function and also extensions to higher order have been put forward. As the larger part of these studies have been carried using a finite difference approach to compute these concepts, we now embarked on (an extension of) an analytical approach to conceptual DFT. With the ultimate aim of providing a complete set of analytically computable second order properties, including the softness and hardness kernels, the hardness as the simplest second order response function is scrutinized again with numerical results highlighting the difference in nature between the analytical hardness (referred to as hardness condition) and the Parr-Pearson absolute chemical hardness. The hardness condition is investigated for its capability to gauge the (de)localization error of density functional approximations (DFAs). The analytical Fukui function, besides overcoming the difficulties in the finite difference approach in treating negatively charged systems, also showcases the errors of deviating from the straight-line behavior using fractional occupation number calculations. Subsequently, the softness kernel and its atom-condensed inverse, the hardness matrix, are accessed through the Berkowitz-Parr relation. Revisiting the softness kernel confirms and extends previous discussions on how Kohn's Nearsightedness of Electronic Matter principle can be retrieved and identified as the physicist's version of the chemist's "transferability of functional groups" concept. The accurate, analytical hardness matrix evaluation on the other hand provides further support for the basics of Nalewajski's charge sensitivity analysis. Based on Parr and Liu's functional expansion of the energy functional, a new energy decomposition is introduced with an order of magnitude analysis of the different terms for a series of simple molecules both at their equilibrium geometry and upon variation in bond length and dihedral angle. Finally, for the first time, the perturbation expansion of the energy functional is studied numerically up to second order now that all response functions and integration techniques are at hand. The perturbation expansion energies are in excellent agreement with those obtained directly from DFA calculations giving confidence in the convergence of the perturbation series and its use in judging the importance of the different terms in reactivity investigations.

Mechanochemistry has experienced a renaissance in recent years witnessing, at the molecular level, a remarkable interplay between theory and experiment. Molecular mechanochemistry has welcomed a broad spectrum of quantum-chemical methods to evaluate the influence of an external mechanical force on molecular properties. In this contribution, an overview is given on recent work on quantum mechanochemistry in the Brussels Quantum Chemistry group (ALGC). The effect of an external force was scrutinized both in fundamental topics, like reactivity descriptors in Conceptual DFT, and in applied topics, such as designing molecular force probes and tuning the stereoselectivity of certain types of reactions. In the conceptual part, a brief overview of the techniques introducing mechanical forces into a quantum-mechanical description of a molecule is followed by an introduction to conceptual DFT. The evolution of the electronic chemical potential (or electronegativity), chemical hardness and electrophilicity are investigated when a chemical bond in a series of diatomics is put under mechanical stress. Its counterpart, the influence of mechanical stress on bond angles, is analyzed by varying the strain present in alkyne triple bonds by applying a bending force, taking the strain promoted alkyne-azide coupling cycloaddition as an example. The increase of reactivity of the alkyne upon bending is probed by Fukui functions and the local softness. In the applied part, a new molecular force probe is presented based on an intramolecular 6π-electrocyclization in constrained polyenes operating under thermal conditions. A cyclic process is conceived where ring opening and closure are triggered by applying or removing an external pulling force. The efficiency of mechanical activation strongly depends on the magnitude of the applied force and the distance between the pulling points. The idea of pulling point distances as a tool to identify new mechanochemical processes is then tested in [28]hexaphyrins with an intricate equilibrium between Möbius aromatic and Hückel antiaromatic topologies. A mechanical force is shown to trigger the interconversion between the two topologies, using the distance matrix as a guide to select appropriate pulling points. In a final application, the Felkin-Anh model for the addition of nucleophiles to chiral carbonyls under the presence of an external mechanical force is scrutinized. By applying a force for restricting the conformational freedom of the chiral ketone, otherwise inaccessible reaction pathways are promoted on the force-modified potential energy surfaces resulting in a diastereoselectivity different from the force-free reaction.

The linear response kernel also referred to as linear response function (LRF) in the framework of conceptual density functional theory has gained tremendous success in time-dependent density functional theory. Comparatively less attention has been devoted to the LRF from a chemical reactivity perspective in its time- or frequency-independent context, although it has recently been used to qualitatively describe electron delocalization, (anti-)aromaticity, inductive and mesomeric effects, etc. Despite these successes, which were obtained by approximating the LRF using the independent particle approximation deriving from a coupled-perturbed Kohn-Sham computation, the robustness of this LRF approach needs to be assessed. In this work, we compute the LRF at four levels of approximation (independent particle approximation, random phase approximation, Hartree-Fock approximation, and the (exact) DFT (density functional theory) expression) using functionals from the first four rungs of Jacob's ladder of exchange-correlation energy functionals. To scrutinize the impact of these approximations, new visualization strategies are discussed and systematized. The overall conclusion is that the independent particle approximation yields qualitatively correct results (ergo previous conceptual applications of the LRF are trustworthy), but for quantitative results, LRF expressions including coulomb and exchange(-correlation) terms should be included. With respect to functionals, density-gradient contributions to the exchange-correlation kernel are less than 10% and may be omitted safely where that is preferable computationally.

Anti-Felkin-Anh diastereoselectivity can be achieved for nucleophilic additions to α-chiral ketones upon stretching the ketone with a mechanical pulling force. Herein, a mechanochemical Felkin-Anh model is proposed for predicting the outcome of a nucleophilic addition to an α-chiral ketone. Essentially, the fully stretched chiral ketone has one substituent shielding each side of the carbonyl, in contrast to the Felkin-Anh model, in which free rotation around a bond is required to achieve the two rotamers of the ketone. Depending on the pulling scenario, either Felkin-Anh or anti-Felkin-Anh diastereoselectivity is obtained. The model is entirely based on the distance between the pulling points, which is maximized in the anti-periplanar arrangement. The major diastereomer is associated with the approach with the least steric interactions. The intuitive model is validated by means of mechanochemical density functional theory calculations. Importantly, the ketone is fully stretched in the sub 1 nN force regime, thus minimizing the risk of undesired homolytic bond rupture. Moreover, the mechanical force is not used for lowering the reaction barriers associated with the nucleophilic addition; instead, it is solely applied for locking the conformation of a molecule and provoking otherwise inaccessible reaction pathways on the force-modified potential energy surface.

The necessity of the recent incorporation of new external variables in the context of conceptual DFT (CDFT) is discussed based on the ever-increasing portfolio of experimental reaction conditions in the endeavor of experimentalists to synthesize new molecules with unprecedented properties. Electric and magnetic fields (ε and B), mechanical forces (F), and confinement are proposed as valuable new variables, extending conventional CDFT and its associated response functions. A finite field approach is used to calculate the evolution of both global and local descriptors in a selected series of atomic and molecular applications, and from it derive new response function involving, with one exception, the first derivative to the field considered. The electric field results, displaying, for example, a case of a field-induced enantioselectivity in the Fukui function, may be instrumental in the recent upsurge of chemistry in oriented external electric fields. The study of atomic electronegativity and hardness in magnetic fields displays a piecewise behavior, associated to configurational jumps upon increasing field strength and reveals an overall compression of their ranges for stronger fields, which may be guiding upon investigating chemistry in extremely high fields like in white dwarfs. The evolution of the electronegativity and hardness of diatomics under mechanical force can elegantly be traced back to differences in their equilibrium distance in the neutral, cationic, and anionic state. The well-known reduction of the polarizability under confinement can be seen as a fore-runner of the increasing hardness of atoms under pressure, presently under investigation. Periodicity showing up in a spontaneous way in the variety of properties is a leitmotiv in this study, as well as the interconnections/analogies between the different response functions., (© 2022 Wiley Periodicals LLC.)

In this paper, the history, present status, and future of density-functional theory (DFT) is informally reviewed and discussed by 70 workers in the field, including molecular scientists, materials scientists, method developers and practitioners. The format of the paper is that of a roundtable discussion, in which the participants express and exchange views on DFT in the form of 302 individual contributions, formulated as responses to a preset list of 26 questions. Supported by a bibliography of 777 entries, the paper represents a broad snapshot of DFT, anno 2022.

An overview of mathematical properties of the non-local second order derivatives of the canonical, grand canonical, isomorphic, and grand isomorphic ensembles is given. The significance of their positive or negative semidefiniteness and the implications of these properties for atoms and molecules are discussed. Based on this property, many other interesting properties can be derived, such as the expansion in eigenfunctions, bounds on the diagonal and off-diagonal elements, and the eigenvalues of these kernels. We also prove Kato's theorem for the softness kernel and linear response and the dissociation limit of the linear responses as the sum of the linear responses of the individual fragments when dissociating a system into two non-interacting molecular fragments. Finally, strategies for the practical calculation of these kernels, their eigenfunctions, and their eigenvalues are discussed.

The position of conceptual density functional theory (CDFT) in the history of density functional theory (DFT) is sketched followed by a chronological report on the introduction of the various DFT descriptors such as the electronegativity, hardness, softness, Fukui function, local version of softness and hardness, dual descriptor, linear response function, and softness kernel. Through a perturbational approach they can all be characterized as response functions, reflecting the intrinsic reactivity of an atom or molecule upon perturbation by a different system, including recent extensions by external fields. Derived descriptors such as the electrophilicity or generalized philicity, derived from the nature of the energy vs. N behavior, complete this picture. These descriptors can be used as such or in the context of principles such as Sanderson's electronegativity equalization principle, Pearson's hard and soft acids and bases principle, the maximum hardness, and more recently, the minimum electrophilicity principle. CDFT has known an ever-growing use in various subdisciplines of chemistry: from organic to inorganic chemistry, from polymer to materials chemistry, and from catalysis to nanotechnology. The increasing size of the systems under study has been coped with thanks to methodological evolutions but also through the impressive evolution in software and hardware. In this flow, biosystems entered the application portfolio in the past twenty years with studies varying (among others) from enzymatic catalysis to biological activity and/or the toxicity of organic molecules and to computational peptidology. On the basis of this evolution, one can expect that "the best is yet to come".

High pressure chemistry offers the chemical community a range of possibilities to control chemical reactivity, develop new materials and fine-tune chemical properties. Despite the large changes that extreme pressure brings to the table, the field has mainly been restricted to the effects of volume changes and thermodynamics with less attention devoted to electronic effects at the molecular scale. This paper combines the conceptual DFT framework for analyzing chemical reactivity with the XP-PCM method for simulating pressures in the GPa range. Starting from the new derivatives of the energy with respect to external pressure, an electronic atomic volume and an atomic compressibility are found, comparable to their enthalpy analogues, respectively. The corresponding radii correlate well with major known sets of this quantity. The ionization potential and electron affinity are both found to decrease with pressure using two different methods. For the electronegativity and chemical hardness, a decreasing and increasing trend is obtained, respectively, and an electronic volume-based argument is proposed to rationalize the observed periodic trends. The cube of the softness is found to correlate well with the polarizability, both decreasing under pressure, while the interpretation of the electrophilicity becomes ambiguous at extreme pressures. Regarding the electron density, the radial distribution function shows a clear concentration of the electron density towards the inner region of the atom and periodic trends can be found in the density using the Carbó quantum similarity index and the Kullback-Leibler information deficiency. Overall, the extension of the CDFT framework with pressure yields clear periodic patterns., Competing Interests: There are no conflicts to declare., (This journal is © The Royal Society of Chemistry.)

An extension of conceptual DFT to include the influence of an external magnetic field is proposed in the context of a program set up to cope with the ever increasing variability of reaction conditions and concomitant reactivity. The two simplest global reactivity descriptors, the electronic chemical potential ( μ ) and the hardness ( η ), are considered for the main group atoms H-Kr using current density-functional theory. The magnetic field strength, | B |, is varied between 0.0 and 1.0 B 0 = ħe -1 a 0 -2 ≈ 2.3505 × 10 5 T, encompassing the Coulomb and intermediate regimes. The carbon atom is studied as an exemplar system to gain insight into the behaviour of the neutral, cationic and anionic species under these conditions. Their electronic configurations change with increasing | B |, leading to a piecewise behaviour of the ionization energy ( I ) and electron affinity ( A ) values as a function of | B |. This results in complex behaviour of properties such as the electronegativity χ = -1/2( I + A ) = - μ and hardness η = 1/2( I - A ). This raises an interesting question: to what extent are atomic properties periodic in the presence of a magnetic field? In the Coulomb regime, close to | B | = 0, we find the familiar periodicity of the atomic properties, and make the connections to response functions central to conceptual DFT. However, as the field increases in the intermediate regime configurational changes of the atomic species lead to discontinuous changes in their properties; fundamentally changing their behaviour, which is illustrated by constructing a periodic table of χ and η values at | B | = 0.5 B 0 . These values tend to increase for groups 1-2 and decrease for groups 16-18, leading to a narrower range overall and suggesting substantial changes in the chemistry of the main group elements. Changes within each group are also examined as a function of | B |. These are more complex to interpret due to the larger number of configurations accessible to heavier elements at high field. This is illustrated for group 17 where Cl and Br have qualitatively different configurations to their lighter cogener at | B | = 0.5 B 0 . The insight into periodic trends in strong magnetic fields may provide a crucial starting point for predicting chemical reactivity under these exotic conditions., Competing Interests: There are no conflicts to declare., (This journal is © The Royal Society of Chemistry.)

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.

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.

Aims: Previous research has shown the possibility to use the pre-operative period to improve a patient's tolerance for surgery. However, there is limited experience with prehabilitation in cardiac surgery. The aim of this study is to evaluate the effect of a comprehensive personalized teleprehabilitation programme on major adverse cardiac events (MACE) in patients scheduled for elective cardiac surgery. Secondary outcomes are post-operative complications, cardiovascular risk factors, quality of life, and cost-effectiveness., Methods and Results: In this single-centre randomized controlled trial, patients are eligible for inclusion when they are ≥18 years of age and cardiac surgery is scheduled at least 8 weeks from informed consent. Participants will be randomized to the teleprehabilitation group or the control group. After a digital baseline screening for perioperative risk factors, patients in the intervention arm can pre-operatively be referred to one or more of the prehabilitation modules (functional exercise training, inspiratory muscle training, psychological support, nutritional support, and/or smoking cessation). The programme is targeted at a duration of at least 6 weeks. It is executed by a multidisciplinary team using (video)calls and supported by a custom-made digital platform. During the pre-operative period, the platform is also used to inform patients about their upcoming surgery and for telemonitoring., Conclusion: Reducing perioperative risk factors might result in a reduction of MACE, post-operative complications, length of stay, and cardiovascular risk factors, as well as improved quality of life. Cost-effectiveness will be evaluated., (© The Author(s) 2021. Published by Oxford University Press on behalf of the European Society of Cardiology.)

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.

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

Background: Malnutrition is prevalent in hospitalized patients. To support muscle maintenance in older and chronically ill patients, a protein intake of 1.2-1.5 g/kg/d has been recommended during hospitalization. We assessed daily protein intake levels and distribution in older patients at risk for malnutrition during hospitalization., Methods: In this prospective, observational study, we measured actual food and food supplement consumption in patients (n = 102; age, 68 ± 14 years; hospital stay, 14 [8-28] days) at risk of malnutrition during hospitalization. Food provided by hospital meals, ONS, and snacks and the actual amount of food (not) consumed were weighed and recorded for all patients., Results: Hospital meals provided 1.03 [0.77-1.26] protein, whereas actual protein consumption was only 0.65 [0.37-0.93] g/kg/d. Protein intake at breakfast, lunch, and dinner was 10 [6-15], 9 [5-14], and 13 [9-18] g, respectively. The use of ONS (n = 62) resulted in greater energy (1.26 [0.40-1.79] MJ/d, 300 [100-430] kcal/d) and protein intake levels (11 [4-16] g/d), without changing the macronutrient composition of the diet., Conclusion: Despite protein provision of ∼1.0 g/kg/d, protein intake remains well below these values (∼0.65 g/kg/d), as 30%-40% of the provided food and supplements is not consumed. Provision of ONS may increase energy and protein intake but does not change the macronutrient composition of the diet. Current nutrition strategies to achieve the recommended daily protein intake in older patients during their hospitalization are not as effective as generally assumed., (© 2020 American Society for Parenteral and Enteral Nutrition.)

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 .

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.

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.

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.

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.

Objective: Hospitalization is generally accompanied by changes in food intake. Patients typically receive hospital meals upon personal preference within the framework of the food administration services of the hospital. In the present study, we assessed food provision and actual food and snack consumption in older patients admitted for elective hip or knee arthroplasty., Design: A prospective observational study., Setting: Orthopedic nursing ward of the Maastricht University Medical Centre+., Participants: In the present study, n=101 patients (age: 67±10 y; hospital stay: 6.1±1.8 d) were monitored during hospitalization following elective hip or knee arthroplasty., Measurements: Energy and protein provided by self-selected hospital meals and snacks, and actual energy and protein (amount, distribution, and source) consumed by patients was weighed and recorded throughout 1-6 days., Results: Self-selected meals provided 6.5±1.5 MJ•d-1, with 16, 48, and 34 En% provided as protein, carbohydrate, and fat, respectively. Self-selected hospital meals provided 0.75±0.16 and 0.79±0.21 g•kg-1•d-1 protein in males and females, respectively. Actual protein consumption averaged merely 0.59±0.18 and 0.50±0.21 g•kg-1•d-1, respectively. Protein consumption at breakfast, lunch, and dinner averaged 16±8, 18±9, and 20±6 g per meal, respectively., Conclusions: Though self-selected hospital meals provide patients with ~0.8 g•kg-1•d-1 protein during short-term hospitalization, actual protein consumption falls well below 0.6 g•kg-1•d-1 with a large proportion (~32%) of the provided food being discarded. Alternative strategies are required to ensure maintenance of habitual protein intake in older patients admitted for elective orthopedic surgery., Competing Interests: AAJV is an employee of Daily Fresh Food B.V., Geleen, the Netherlands. Daily Fresh Food B.V. is a manufacturer and supplier of fresh food products and convenient meal components for healthcare and consumer services. PE has received research grants of Zimmer-Biomet, Regentis, Active Implants, DSM Biomedical, DePuySynthes, and is a shareholder of Chondropeptix B.V. IWKK and LBV have received speaking honoraria from Nutricia Research. LJCvL has received research grants, consulting fees, and/or speaking honoraria from Friesland Campina and Nutricia Research. None of the other authors had any personal or financial conflict of 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.

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

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

We establish the physical origins of chemical transferability from the perspective of the nearsightedness of electronic matter. To do this, we explicitly evaluate the response of electron density to a change in the system, at constant chemical potential, by computing the softness kernel, [Formula: see text] The softness kernel is nearsighted, indicating that under constant-chemical-potential conditions like dilute solutions changing the composition of the molecule at [Formula: see text] has only local effects and does not have any significant impact on the reactivity at positions [Formula: see text] far away from point [Formula: see text] This locality principle elucidates the transferability of functional groups in chemistry., Competing Interests: The authors declare no conflict of interest., (Published under the PNAS license.)

The structure and morphology of three polymer/graphene nanocomposites have been studied using classical molecular dynamics (MD) simulations. The simulations use 10-monomer oligomeric chains of three polymers: polyethylene (PE), polystyrene (PS) and polyvinylidene fluoride (PVDF). The structure of the polymer chains at the graphene surface has been investigated and characterized by pair correlation functions (PCF), g(r), g(θ) and g(r,θ). In addition, the influence of the temperature on the graphene/polymer interactions has been analysed for each of the three polymer/graphene nanocomposite systems. The results indicate that graphene induces order in both the PE and PVDF systems by providing a nucleation site for crystallisation, steering the growth of oligomer crystals according to the orientation of the graphene sheet, whereas the PS system remains disordered in the presence of graphene. The overall results are in line with the findings in a recent quantumchemical study by some of the present authors.

Understanding the interaction between graphene and polymers is of essential interest when designing novel nanocomposites with reinforced mechanical and electrical properties. In this computational study, the interaction of pristine graphene (PG) and graphene oxide (GO) with a series of functional groups, representative of the functionalised buildings blocks occurring in different polymers, and attached to aliphatic and aromatic chains, is analyzed using dispersion-corrected semi-empirical methods (PM6-D3H4X) and density functional theory calculations with empirical dispersion corrections. Functional groups include alkyl, hydroxyl, aldehyde, carboxyl, amino and nitro groups, and the binding energies of these groups with graphene derivatives (PG and GO) are determined. Nitro- and carbonyl groups display stronger interactions in both aliphatic and aromatic chains. The importance of dispersion-type and non-covalent interactions (NCI) in general, which typically, double the interaction energies, is revealed. The results are interpreted in an extensive NCI analysis in order to characterize the different types of NCI, providing a better understanding of the nature of the interaction (π-π stacking, CH-π bonding, H-bonding and lone pair-π interaction) at stake. In order to highlight the influence of polymer structure/conformation on top of that of their functional groups, the binding of three polymers, polyethylene (PE), polystyrene (PS) and polyvinylidene fluoride (PVDF), on pristine graphene is also investigated. Our calculations indicate that, although all polymers exhibit evident attractive interactions with the graphene sheet, the overall interaction is strongly influenced by the specific polymer structure. Thus, three main conformations of PVDF (the so-called α, β and γ, ε conformations) are analyzed and we find that, although the α-conformer with a trans-gauche-trans-gauche (TGTG') conformation is the lowest energy conformer, the β-conformation of PVDF with the hydrogen atoms facing the graphene ("F-up") has the strongest interaction with the graphene surface among the polymers under consideration. Taken together, our computational approach sheds light on the character and importance of non-covalent graphene-polymer functional group interactions combined with the structural/conformational properties of the polymer, which are at stake in the design of novel nanocomposites with reinforced mechanical and electrical properties.

Electrophilic aromatic sulfonation of benzene with sulfur trioxide is studied with ab initio molecular dynamics simulations in gas phase, and in explicit noncomplexing (CCl 3 F) and complexing (CH 3 NO 2 ) solvent models. We investigate different possible reaction pathways, the number of SO 3 molecules participating in the reaction, and the influence of the solvent. Our simulations confirm the existence of a low-energy concerted pathway with formation of a cyclic transition state with two SO 3 molecules. Based on the simulation results, we propose a sequence of elementary reaction steps and a kinetic model compatible with experimental data. Furthermore, a new alternative reaction pathway is proposed in complexing solvent, involving two SO 3 and one CH 3 NO 2 .

There is an increasing interest in applying quantum chemistry to rationally design novel compounds with some desired characteristics. Furthermore, many applications require more than one property to be optimal. In this Concept, several inverse design strategies, based on the discrete best first search scheme, are introduced that allow for the simultaneous optimization of multiple properties or the optimization of the most vital target property with constraints for secondary properties. A detailed assessment of the different optimization techniques is carried out, and special attention is paid to improve the cost efficacy and performance by tuning the process parameters. Our suggested protocol allows for a more successful optimization routine when additional boundary conditions are desired., (© 2016 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Several expanded porphyrins switch between Hückel, Möbius and twisted-Hückel topologies, encoding different aromaticity and NLO properties. Since the topological switch can be induced by different external stimuli, expanded porphyrins represent a promising platform to develop molecular switches for molecular electronic devices. In order to determine the optimum conditions for efficient molecular switches from octaphyrins, we have carried out a comprehensive quantum chemical study focusing on the conformational preferences and aromaticity of [36]octaphyrins. Different external stimuli for triggering the topological switch have been considered in our work, such as protonation and redox reactions. Importantly, the structure-property relationships between the molecular conformation, the number of π-electrons and aromaticity in octaphyrins have been established by using energetic, magnetic, structural and reactivity descriptors. Remarkably, we found that the aromaticity of octaphyrins is highly dependent on the π-conjugation topology and the number of π-electrons and it can be modulated by protonation and redox reactions. A non-aromatic figure-eight conformation is strongly preferred by neutral [36]octaphyrins that switches to a Möbius aromatic conformation upon protonation. Such a change of topology involves an aromaticity switch in a single molecule and is accompanied by a drastic change in the NLO properties. In contrast, the twisted-Hückel topology remains the most stable one in the oxidized and reduced species, but the aromaticity is totally reversed upon redox reactions. Aromaticity is shown to be a key concept in expanded porphyrins, determining the electronic, magnetic and NLO properties of these macrocycles.

The role of the solvent and the influence of dynamics on the kinetics and mechanism of the SNAr reaction of several halonitrobenzenes in liquid ammonia, using both static calculations and dynamic ab initio molecular dynamics simulations, are investigated. A combination of metadynamics and committor analysis methods reveals how this reaction can change from a concerted, one-step mechanism in gas phase to a stepwise pathway, involving a metastable Meisenheimer complex, in liquid ammonia. This clearly establishes, among others, the important role of the solvent and highlights the fact that accurately treating solvation is of crucial importance to correctly unravel the reaction mechanism. It is indeed shown that H-bond formation of the reacting NH3 with the solvent drastically reduces the barrier of NH3 addition. The halide elimination step, however, is greatly facilitated by proton transfer from the reacting NH3 to the solvent. Furthermore, the free energy surface strongly depends on the halide substituent and the number of electron-withdrawing nitro substituents.

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