Search

Non-covalent interactions involving aromatic rings play a central role in many areas of modern chemistry. In medicinal and bioorganic chemistry, the intermolecular interactions between the aromatic side chains of amino acids, such as phenylalanine and tyrosine, are of great interest. To enhance the affinity between such aromatic side chains, halogenation is a promising modification strategy. In the current work, the nature and strength of halogenated p-p stacked phenylalanine (Phe) dimers have been investigated using density functional theory, energy decomposition analyses and the non-covalent interaction (NCI) method. Our analysis shows that increasing the degree of halogenation enhances the strength of the stacking interactions and, moreover, the heavier halides (Cl, Br and I) lead to stronger interactions compared to the lighter F. This effect was traced back to local secondary interactions of the halide with the aliphatic C-H bonds of the phenylalanine side chain. Based on the computational findings, a set of peptide hydrogelators was synthesized, and the resulting hydrogel properties were further investigated via dynamic rheometry. Experimental observations can be correlated to the trends found in the theoretical analysis, suggesting that local interactions indeed play a noticeable role in enhancing peptide-based hydrogel strength. This journal is

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.

This chapter illustrates the good performance of the recently introduced noncovalent interactions (NCI) method in understanding molecular reactivity. This method is not only helpful in identifying the nature of the NCIs but can be also used to gain a deeper insight into the influence of such interactions on the outcome of different fundamental transformations in chemistry, including catalysed processes. To this end, representative catalysed transformations were selected where the NCI method was key to rationalizing different aspects such as reactivity trends and selectivity. The catalysed reactions chosen range from relatively simple transformations such as Diels–Alder cycloadditions to more intricate transition metal- and organo-catalysed processes.

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.

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 nature and origin of ion-π and ion-σ interactions has been systematically investigated using dispersion-corrected density functional theory and the recently developed noncovalent interaction (NCI) method. A detailed analysis of these interactions is performed with the aim to identify the requirements that have to be fulfilled by the molecular system for strong ion-ligand interactions. Interestingly, our results indicate that aliphatic systems, such as cyclohexane, can interact as strong as aromatic ones with both cations and anions, despite of having a negligible quadrupole moment. In fact, cyclohexane binds anions stronger than benzene itself but slightly weaker that hexafluorobenzene. The NCI method reveals that the interaction between the ions and three C–H bonds of the saturated fragment are responsible for the surprisingly strong ion-σ interaction. A weakening of the ion-σ interactions is observed in the order: Li + > F − > Na + > Cl − > Br − ≈ K + . In addition, a complete Ziegler–Rauk type energy decomposition analysis has been carried out in order to reveal the origins of the thermodynamic driving force for complex formations. The electron density deformation upon complex formation has been scrutinized with a complementary NOCV analysis allowing the identification of molecular orbital interaction contributions to the stabilization. Based on these analysis, it is shown that the formally anion-π interaction is rather an anion-σ ∗ interaction.

International audience; The role of aromaticity in porphyrinoids is a current subject of debate due to the intricate structure of these macrocycles, which in some cases adopt Hückel, Möbius and even figure-eight conformers. One of the main challenges in these large π-conjugated structures is identifying the most conjugated pathway because, among aromaticity descriptors, there are very few that can be applied coherently to this variety of conformers. In this paper, we have investigated the conjugated pathways in nine porphyrinoid compounds using several aromaticity descriptors, including BLA, BOA, FLU and HOMA, as well as the recently introduced AV1245 and AVmin indices. All the indices agree on the general features of these compounds, such as the fulfillment of Hückel's rule or which compounds should be more or less aromatic from the series. However, our results evince the difficulty in finding the most aromatic pathway in the macrocycle for large porphyrinoids. In fact, only AVmin is capable of recognizing the annulene pathway as the most aromatic one in the nine studied structures. Finally, we study the effect of the exchange in DFT functionals on the description of the aromaticity of the porphyrinoids. The amount of exact exchange quantitatively changes the picture for most aromaticity descriptors, AVmin being the only exception that shows the same qualitative results in all cases.

International audience; Describing non-covalent interactions (NCIs) has shown to be of paramount importance in many areas of theoretical chemistry and related disciplines, such as biochemistry and material science. However, non-covalent interactions are subtle effects, very difficult to reproduce from most common computational approaches. Electron density studies have shown to provide a good semiquantitative visual approach to such interactions, which are much less prone to method dependency. But to which extent? This is the question addressed in this contribution. The NCI approach based on the reduced density gradient is given the third degree so as to provide the user with a benchmark on how it is affected by the computational method and the basis set of choice. We have assessed the dependence of the NCI results on the geometry. This last question is addressed in detail to dissect how, why and when the NCI method can be used to understand dispersion interactions. Along various examples, we will show that the NCI index is very little dependent on the method and basis set used in the calculation of the electron density as long as the geometry is kept fixed. Indeed, the biggest variations in NCI come from changes in the geometry. Thus, methods which provide descriptions of a given interaction type of different accuracies will yield different electron density organizations. This gives no qualitative variations in the NCI 3D picture. But it is reflected in quantitative NCI measures even in very subtle cases. Moreover, in the case of a failure of the calculation method, NCI can also reveal the sources of its error. NCI volumes are able to locate the energetic ordering in various conformational situations, but always in a relative manner. Absolute values should not be used in comparisons, nor between compounds that do not belong to the same family.

International audience; In recent years, expanded porphyrins have emerged as a promising class of π-conjugated molecules that display unique electronic, optical and conformational properties. Several expanded porphyrins can switch between planar and twisted conformations, which have different photophysical properties. Such a change of topology involves a Hückel–Möbius aromaticity switch in a single molecule and it can be induced by solvent, pH and metallation. These features make expanded porphyrins suitable for the development of a novel type of molecular switches for molecular electronic devices. Octaphyrins consisting of eight pyrrole rings, exhibit twisted-Hückel, Möbius and Hückel π-conjugation topologies depending on the oxidation and protonation state, with distinct electronic structures and aromaticity. Our working hypothesis is that a significant change in the conductance of expanded porphyrins will be observed after the topology switching. Despite the potential of Hückel–Möbius systems as conductance switches, the relationship between the conductance and the molecular topology is not yet understood. We have explored the performance of local descriptors of conductivity in simple molecules, as well as the relationship with conductance. Since these indexes provide a qualitative measure of delocalization and conductance in the probe molecules, we have carried out a local analysis of electrical conductance changes as a function of the π-conjugation in two examples. In one of them, the locality of the electronic changes ensures the ability of these indexes to describe the conductance as local. Moreover, it enables to identify which conformational switch would be more efficient from an electronic device perspective. However, we also show that it is not always possible to reduce conductance changes to one bond, and in those molecules where a deep rearrangement occurs far from the structural perturbation, local measures show a limited efficiency. This is a first step for the description of the connection between the molecular structure and conductance in molecular switches.

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

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. ispartof: Chemistryopen vol:8 issue:6 pages:788-806 ispartof: location:Germany status: Published online