Aromaticity and Magnetic Behavior in Benzenoids: Unraveling Ring Current Combinations.

Nowadays, an active research topic is the connection between Clar's rule, aromaticity, and magnetic properties of polycyclic benzenoid hydrocarbons. In the present work, we employ a meticulous magnetically induced current density analysis to define the net current flowing through any cyclic circuit, connecting it to aromaticity based on the ring current concept. Our investigation reveals that the analyzed polycyclic systems display a prominent global ring current, contrasting with subdued semi-local and local ring currents. These patterns align with Clar's aromatic π-sextets only in cases where migrating π-sextet structures are invoked. The results of this study will enrich our comprehension of aromaticity and magnetic behavior in such systems, offering significant insights into coexisting ring current circuits in these systems.

Aromaticity in the Spectroscopic Spotlight of Hexaphyrins.

Spectroscopic properties are commonly used in the experimental evaluation of ground- and excited-state aromaticity in expanded porphyrins. Herein, we investigate if the defining photophysical properties still hold for a diverse set of hexaphyrins with varying redox states, topologies, peripheral substitutions, and core-modifications. By combining TD-DFT calculations with several aromaticity descriptors and chemical compound space maps, the intricate interplay between structural planarity, aromaticity, and absorption spectra is elucidated. Our results emphasize that the general assumption that antiaromatic porphyrinoids exhibit significantly attenuated absorption bands as compared to aromatic counterparts does not hold even for the unsubstituted macrocycles. To connect the spectroscopic properties to the hexaphyrins' aromaticity behaviour, we analyzed chemical compound space maps defined by the various aromaticity indices. The intensity of the Q-band is not well described by the macrocyclic aromaticity. Instead, the degeneracy of the frontier molecular orbitals appears to be a better indicator to identify hexaphyrins with enhanced light-absorbing abilities in the near-infrared region. Regions with highly planar (anti)aromatic hexaphyrin structuresare characterized by an intense B-band. Hence, we advise to use a combination of both global and local aromaticity descriptors rooted in different criteria to assess the aromaticity of expanded porphyrins instead of solely using the absorption spectrum.

New Perspectives on Delocalization Pathways in Aromatic Molecular Chameleons.

This study comprehensively analyzes the magnetically induced current density of polycyclic compounds labeled as "aromatic chameleons" since they can arrange their π-electrons to exhibit aromaticity in both the ground and the lowest triplet state. These compounds comprise benzenoid moieties fused to a central skeleton with 4n π-electrons and traditional magnetic descriptors are biased due to the superposition of local magnetic responses. In the S0 state, our analysis reveals that the molecular constituent fragments preserve their (anti)aromatic features in agreement with two types of resonant structures: one associated with aromatic benzenoids and the other with a central antiaromatic ring. Regarding the T1 state, a global and diatropic ring current is revealed. Our aromaticity study is complemented with advanced electronic and geometric descriptors to consider different aspects of aromaticity, particularly important in the evaluation of excited state aromaticity. Remarkably, these descriptors consistently align with the general features on the main delocalization pathways in polycyclic hydrocarbons consisting of fused 4n π-electron rings. Moreover, our study demonstrates an inverse correlation between the singlet-triplet energy difference and the antiaromatic character of the central ring in S0.

T-Shaped Palladium and Platinum {MNO}10 Nitrosyl Complexes.

The synthesis and characterization of a homologous series of T-shaped {MNO}10 nitrosyl complexes of the form [M(PR3)2(NO)]+ (M = Pd, Pt; R = tBu, Ad) are reported. These diamagnetic nitrosyls are obtained from monovalent or zerovalent precursors by treatment with NO and NO+, respectively, and are notable for distinctly bent M-NO angles of ∼123° in the solid state. Adoption of this coordination mode in solution is also supported by the analysis of isotopically enriched samples by 15N NMR spectroscopy. Effective oxidation states of M0/NO+ are calculated, and metal-nitrosyl bonding has been interrogated using DFT-based energy decomposition analysis techniques. While a linear nitrosyl coordination mode was found to be electronically preferred, the M-NO and P-M-P angles are inversely correlated to the extent that binding in this manner is prevented by steric repulsion between the bulky ancillary phosphine ligands.

Similarities and Differences in Benzene Reduction with Ca, Sr, Yb and Sm: Strong Evidence for Tetra-Anionic Benzene.

Inverse sandwich complexes of Yb and Sm stabilized by a bulky ß-diketiminate (BDI) ligand have been prepared: (BDI)Ln(η6,η6-C6H6)Ln(BDI); Ln=lanthanide. Coordinated benzene ligands can be neutral, di-anionic or, often controversially discussed, even tetra-anionic. The formal charge on benzene is correlated to assignment of the metal oxidation state which generally poses a problem. Herein, we take advantage of the structural similarities found when comparing CaII with YbII, and SrII with SmII complexes. In this work, we found an excellent overlap of the Ca/Yb inverse sandwich structures but a striking difference for the Sr/Sm pair. The much shorter Sm-N and Sm-C6H6 distances are strong evidence for a SmIII-benzene-4-SmIII assignment. This was further supported by NMR spectroscopy, magnetic susceptibility, reactivity and comprehensive computational investigation.

Quest for the Most Aromatic Pathway in Charged Expanded Porphyrins.

Despite the central role of aromaticity in the chemistry of expanded porphyrins, the evaluation of aromaticity remains difficult for these extended macrocycles. The presence of multiple conjugation pathways and different planar and nonplanar π-conjugation topologies makes the quantification of global and local aromaticity even more challenging. In neutral expanded porphyrins, the predominance of the aromatic conjugation pathway passing through the imine-type nitrogens and circumventing the amino NH groups is established. However, for charged macrocycles, the question about the main conjugation circuit remains open. Accordingly, different conjugation pathways in a set of neutral, anionic, and cationic expanded porphyrins were investigated by means of several aromaticity indices rooted in the structural, magnetic, and electronic criteria. Overall, our results reveal the predominance of the conjugation pathway that passes through all nitrogen atoms to describe the aromaticity of deprotonated expanded porphyrins, while the outer pathway through the perimeter carbon atoms becomes the most aromatic in protonated macrocycles. In nonplanar and charged macrocycles, a discrepancy between electronic and magnetic descriptors is observed. Nevertheless, our work demonstrates AVmin remains the best tool to determine the main conjugation pathway of expanded porphyrins.

Understanding the Reactivity of Supported Late Transition Metals on a Bare Anatase (101) Surface: A Periodic Conceptual DFT Investigation.

The rapidly growing interest for new heterogeneous catalytic systems providing high atomic efficiency along with high stability and reactivity triggered an impressive progress in the field of single-atom catalysis. Nevertheless, unravelling the factors governing the interaction strength between the support and the adsorbed metal atoms remains a major challenge. Based on periodic density functional theory (DFT) calculations, this paper provides insight into the adsorption of single late transition metals on a defect-free anatase surface. The obtained adsorption energies fluctuate, with the exception of Pd, between -3.11 and -3.80 eV and are indicative of a strong interaction. Depending on the considered transition metal, we could attribute the strength of this interaction with the support to i) an electron transfer towards anatase (Ru, Rh, Ni), ii) s-d orbital hybridisation effects (Pt), or iii) a synergistic effect between both factors (Fe, Co, Os, Ir). The driving forces behind the adsorption were also found to be strongly related to Klechkowsky's rule for orbital filling. In contrast, the deviating behaviour of Pd is most likely associated with the lower dissociation enthalpy of the Pd-O bond. Additionally, the reactivity of these systems was evaluated using the Fermi weighted density of states approach. The resulting softness values can be clearly related to the electron configuration of the catalytic systems as well as with the net charge on the transition metal. Finally, these indices were used to construct a model that predicts the adsorption strength of CO on these anatase-supported d-metal atoms. The values obtained from this regression model show, within a 95 % probability interval, a correlation of 84 % with the explicitly calculated CO adsorption energies.

Mechanochemical Felkin-Anh Model: Achieving Forbidden Reaction Outcomes with Mechanical Force.

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.

Wandering through quantum-mechanochemistry: from concepts to reactivity and switches.

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.

Designing hexaphyrins for high-potential NLO switches: the synergy of core-modifications and meso-substitutions.

Due to the enormous size of the chemical compound space, usually only small regions are traversed with traditional direct molecular design approaches making the discovery for novel functionalized molecules for nonlinear optical applications challenging. By applying inverse molecular design algorithms, we aim to efficiently explore larger regions of the compound space in search of promising hexaphyrin-based molecular switches as measured by their first-hyperpolarizability (ßHRS) contrast. We focus on the 28R → 30R switch with a functionalization pattern allowing for centrosymmetric OFF states yielding zero ßHRS response. This switch is particularly challenging as full meso-substitution with a single type of functional group or core-modifications result in almost no contrast enhancement. We carried out four inverse design procedures during which two sets of core-modifications and three sets of meso-substitutions sites were systematically optimized. All 4 optimal switches are characterized by a mix of meso-substitutions and core-modifications, of which the best performing switch yields a 10-fold improvement over the parent macrocycle. Throughout the inverse design procedures, we collected and analyzed a database biased towards high NLO contrasts that contains 277 different patterns for hexaphyrin-based switches. We derived three design rules to obtain highly functional 28R → 30R NLO switches: (I) a combination of 2 strong EWG and 1 EDG group is the ideal recipe for increasing the NLO contrast, though their position also plays an important role. (II) The type of core-modification is less important when only the diagonal positions are core-modified. Switches with 4 core-modifications show a clear preference for oxygen. (III) Keeping centrosymmetry in the OFF state remains highly beneficial given the investigated functionalization pattern. Finally, we have demonstrated that combining meso-substitutions with core-modifications can synergistically improve the NLO contrast.

Application of Inverse Design Approaches to the Discovery of Nonlinear Optical Switches.

Molecular switches, in which a stimulus induces a large and reversible change in molecular properties, are of significant interest in the domain of photonics. Due to their commutable redox states with distinct nonlinear optical (NLO) properties, hexaphyrins have emerged as a novel platform for multistate switches in nanoelectronics. In this study, we employ an inverse design algorithm to find functionalized 26R→28R redox switches with maximal ßHRS contrast. We focus on the role of core modifications, since a synergistic effect with meso-substitutions was recently found for the 30R-based switch. In contrast to these findings, the inverse design optima and subsequent database analysis of 26R-based switches confirm that core modifications are generally not favored when high NLO contrasts are targeted. Moreover, while push-pull combinations enhance the NLO contrast for both redox switches, they prefer a different arrangement in terms of electron-donating and electron-withdrawing functional groups. Finally, we aim at designing a three-state 26R→28R→ 30R switch with a similar NLO response for both ON states. Even though our best-performing three-state switch follows the design rules of the 30R-based component, our chemical compound space plots show that well-performing three-state switches can be found in regions shared by high-responsive 26R and 30R structures.

Modulation of Metal Carbonyl Stretching Frequencies in the Second Coordination Sphere through the Internal Stark Effect.

Spectroscopic and computational examination of a homologous series of rhodium(I) pybox carbonyl complexes has revealed a correlation between the conformation of the flanking aryl-substituted oxazoline donors and the carbonyl stretching frequency. This relationship is also observed experimentally for octahedral rhodium(III) and ruthenium(II) variants and cannot be explained through the classical, Dewar-Chatt-Duncanson, interpretation of metal-carbonyl bonding. Instead, these findings are reconciled by local changes in the magnitude of the electric field that is projected along the metal-carbonyl vector: the internal Stark effect.

Mechanochemically Triggered Topology Changes in Expanded Porphyrins.

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.

Reactivity of Single Transition Metal Atoms on a Hydroxylated Amorphous Silica Surface: A Periodic Conceptual DFT Investigation.

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

Designing Force Probes Based on Reversible 6π-Electrocyclizations in Polyenes Using Quantum Chemical Calculations.

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.

Ambident Nucleophilic Substitution: Understanding Non-HSAB Behavior through Activation Strain and Conceptual DFT Analyses.

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

Synthesis and Reactivity of Spirocarbocycles as Scaffolds for Nucleoside Analogues.

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

Stereoselective Reductions of 3-Substituted Cyclobutanones: A Comparison between Experiment and Theory.

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

Switching between Hückel and Möbius aromaticity: a density functional theory and information-theoretic approach study.

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.

Performance of Electronic Structure Methods for the Description of Hückel-Möbius Interconversions in Extended π-Systems.

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.