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.

Aromatic Character and Relative Stability of Pyrazoloporphyrin Tautomers and Related Protonated Species: Insights into How Pyrazole Changes the Properties of Carbaporphyrinoid Systems

Pyrazoloporphyrins (PzPs), which are porphyrin analogues incorporating a pyrazole subunit, are examples of carbaporphyrin-type structures with a carbon atom within the macrocyclic cavity. DFT calculations were used to assess a series of 17 PzP tautomers, nine monoprotonated species and four related diprotonated PzP dications. The geometries of the structures were optimized using M06-2X/6-311++G(d,p), and the relative stabilities computed with the cc-PVTZ functional. Nucleus independent chemical shifts, both NICS(0) and NICS(1)_zz, were calculated, and the anisotropy of the induced current density (AICD) plots were generated for all of the species under investigation. The results for free base PzPs show that fully aromatic PzP tautomers are not significantly more stable than weakly aromatic cross-conjugated species. In addition, strongly aromatic structures with internal CH_2's are much less stable, a feature that is also seen for protonated PzPs. The degree of planarity for the individual macrocycles does not significantly correlate with the stability of these structures. The results allow significant aromatic conjugation pathways to be identified in many cases, and provide insights into the aromatic properties of this poorly studied system. These investigations also complement experimental results for PzPs and emphasize the need for further studies in this area.

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 AV_min remains the best tool to determine the main conjugation pathway of expanded porphyrins.

Fine-Tuning of Nonlinear Optical Contrasts of Hexaphyrin-Based Molecular Switches Using Inverse Design

In the search for new nonlinear optical (NLO) switching devices, expanded porphyrins have emerged as ideal candidates thanks to their tunable chemical and photophysical properties. Introducing meso-substituents to these macrocycles is a successful strategy to enhance the NLO contrasts. Despite its potential, the influence of meso-substitution on their structural and geometrical properties has been scarcely investigated. In this work, we pursue to grasp the underlying pivotal concepts for the fine-tuning of the NLO contrasts of hexaphyrin-based molecular switches, with a particular focus on the first hyperpolarizability related to the hyper-Rayleigh scattering (β_HRS). Building further on these concepts, we also aim to develop a rational design protocol. Starting from the (un)substituted hexaphyrins with various π-conjugation topologies and redox states, structure-property relationships are established linking aromaticity, photophysical properties and β_HRS responses. Ultimately, inverse molecular design using the best-first search algorithm is applied on the most favorable switches with the aim to further explore the combinatorial chemical compound space of meso-substituted hexaphyrins in search of high-contrast NLO switches. Two definitions of the figure-of-merit of the switch performance were used as target objectives in the optimization problem. Several meso-substitution patterns and their underlying characteristics are identified, uncovering molecular symmetry and the electronic nature of the substituents as the key players for fine-tuning the β_HRS values and NLO contrasts of hexaphyrin-based switches.

Porphyrinoids, a unique platform for exploring excited-state aromaticity

Recently, Baird (anti)aromaticity has been referred to as a description of excited-state (anti)aromaticity. With the term of Baird's rule, recent studies have intensively verified that the Hückel aromatic [4n + 2]π (or antiaromatic [4n]π) molecules in the ground state are reversed to give Baird aromatic [4n]π (or Baird antiaromatic [4n + 2]π) molecules in the excited states. Since the Hückel (anti)aromaticity has great influence on the molecular properties and reaction mechanisms, the Baird (anti)aromaticity has been expected to act as a dominant factor in governing excited-state properties and processes, which has attracted intensive scientific investigations for the verification of the concept of reversed aromaticity in the excited states. In this scientific endeavor, porphyrinoids have recently played leading roles in the demonstration of the aromaticity reversal in the excited states and its conceptual development. The distinct structural and electronic nature of porphyhrinoids depending on their (anti)aromaticity allow the direct observation of excited-state aromaticity reversal, Baird's rule. The explicit experimental demonstration with porphyrinoids has contributed greatly to its conceptual development and application in novel functional organic materials. Based on the significant role of porphyrinoids in the field of excited-state aromaticity, this review provides an overview of the experimental verification of the reversal concept of excited-state aromaticity by porphyrinoids and the recent progress on its conceptual application in novel functional molecules.

Carbazole-Containing Carbadecaphyrins: Non-aromatic Expanded Porphyrins that Undergo Proton-Triggered Conformational Changes

A pair of meso-unsubstituted expanded carbaporphyrins containing two carbazole moieties were prepared in high isolated yields (82 and 76 %, respectively). The two macrocycles, namely 3 and 4, differ with respect to their substitution at the carbazole N-atoms i. e. by H and i-Bu, respectively. As prepared in their free-base forms, macrocycles 3 and 4 adopt figure-of-eight conformations and are best characterized as 40 π-electron, non-aromatic species possessing a decaphyrin(1.1.0.0.0.1.1.0.0.0) skeleton. Protonation of 3 with either trifluoroacetic acid (TFA) or perchloric acid (HClO₄ ) produces a parallelogram-shaped structure. A similar structure is produced when N-functionalized system 4 is treated with TFA. In contrast, protonation of 4 with HClO₄ leads it to adopt a twisted Möbius strip-like structure in the solid state, thus allowing access to three distinct conformational states as a function of the conditions.

Performance of Localized Coupled Cluster Methods in a Moderately Strong Correlation Regime: Hückel-Möbius Interconversions in Expanded Porphyrins

Localized orbital coupled cluster theory has recently emerged as a nonempirical alternative to DFT for large systems. Intuitively, one might expect such methods to perform less well for highly delocalized systems. In the present work, we apply both canonical CCSD(T) approximations and a variety of localized approximations to a set of flexible expanded porphyrins—macrocycles that can switch between Hückel, figure-eight, and Möbius topologies under external stimuli. Both minima and isomerization transition states are considered. We find that Möbius(-like) structures have much stronger static correlation character than the remaining structures, and that this causes significant errors in DLPNO-CCSD(T) and even DLPNO-CCSD(T₁) approaches, unless TightPNO cutoffs are employed. If sub-kcal mol^–1 accuracy with respect to canonical relative energies is required even for Möbius-type systems (or other systems plagued by strong static correlation), then Nagy and Kallay’s LNO-CCSD(T) method with “tight” settings is the suitable localized approach. We propose the present POLYPYR21 data set as a benchmark for localized orbital methods or, more broadly, for the ability of lower-level methods to handle energetics with strongly varying degrees of static correlation.

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.

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.

Synthesis of Azulitriphyrins(1.2.1) and Related Benzocarbatriphyrins

Bis(pyrrolylmethyl)azulene dialdehydes underwent intramolecular McMurry coupling, and following oxidation with DDQ and addition of trifluoroacetic acid, gave aromatic azulitriphyrin cations. The proton NMR spectra for these contracted carbaporphyrinoids showed the internal CH upfield at ca. 2 ppm, while the bridging methine units appeared downfield between 8 and 10 ppm, demonstrating that the macrocycles possess significant 14π electron diamagnetic ring currents. Although protonated azulitriphyrins(1.2.1) proved to be very stable, the free base forms were unstable and could not be isolated. Treatment of an azulitriphyrin with KOt-Bu-t-BuOOH gave rise to oxidative ring contractions that afforded the first examples of benzocarbatriphyrins. These contracted porphyrinoids also exhibit aromatic ring currents, and the internal CH resonances shifted upfield to 0.5-0.7 ppm. Protonation afforded cations with slightly enhanced diatropic properties. Nucleus independent chemical shift and anisotropy of induced current density calculations confirmed that benzocarbatriphyrins are 14π electron delocalized aromatic systems that bypass the fused arene units.

Azulichlorins and Benzocarbachlorins Derived Therefrom

Acid-catalyzed condensation of azulidipyrrane aldehydes with a dihydrodipyrrin carbaldehyde afforded the first examples of azulichlorins. These macrocyclic products were isolated in a monoprotonated form, and the free bases proved to be somewhat unstable. The monocations were strongly diatropic, and proton NMR spectroscopy showed the internal C-H at ca. -2 ppm. Addition of TFA gave the related dications, but these exhibited significantly reduced aromatic ring currents. Reaction of an azulichlorin with tert-butyl hydroperoxide and KOH in dichloromethane/methanol gave a benzocarbachlorin and two related aldehydes. The UV-vis spectrum for the benzocarbachlorin showed a split Soret band at 414 and 430 nm, together with a strong chlorin-like absorption at 684 nm. The proton NMR spectrum indicated that the carbachlorin is strongly aromatic and the internal C-H was observed at -4.64 ppm. Addition of TFA afforded a C-protonated dication with a significantly increased diatropic ring current. The proton NMR spectrum, NICS calculations, and AICD plots indicated that the system favors a 22π electron delocalization pathway that runs through the fused benzo unit. Addition of TFA to the benzocarbachlorin aldehydes primarily led to the formation of monocations, and the generation of C-protonated dications was no longer favored due to the presence of electron-withdrawing formyl moieties.

Aromaticity as a Guiding Concept for Spectroscopic Features and Nonlinear Optical Properties of Porphyrinoids

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 [4n+2] π-electron porphyrinoids was evidenced as was the antiaromaticity for [4n] π-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.

New electron delocalization tools to describe the aromaticity in porphyrinoids

There are several possible pathways in the macrocycle of large porphyrinoids and, among aromaticity indices, only AV_minis capable of recognizing the most aromatic one.

Computational Studies on Optoelectronic and Nonlinear Properties of Octaphyrin Derivatives

The electronic and nonlinear optical (NLO) properties of octaphyrin derivatives were studied by employing the DFT/TDFT at CAM-B3LYP/6-311++G (2d, 2p) level of the theory. Thiophene, phenyl, methyl and cyano moieties were substituted on the molecular framework of octaphyrin core, in order to observe the change in optoelectronic and nonlinear response of these systems. The frontier molecular orbital studies and values of electron affinity reveals that the studied compounds are stable against the oxygen and moisture present in air. The calculated ionization energies, adiabatic electron affinity and reorganization energy values indicate that octaphyrin derivatives can be employed as effective n-type material for Organic Light Emitting Diodes (OLEDs). This character shows an enhancement with the introduction of an electron withdrawing group in the octaphyrin framework. The polarizability and hyperpolarizability values of octaphyrin derivatives demonstrate that they are good candidates for NLO devices. The nonlinear response of these systems shows enhancement on the introduction of electron donating groups on octaphyrin moiety. However, these claims needs further experimental verification.

Understanding conductivity in molecular switches: a real space approach in octaphyrins

In recent years, expanded porphyrins have emerged as a promising class of π-conjugated switches whose conductance is studied from the electron density.

Understanding the molecular switching properties of octaphyrins

Triggering Hückel–Möbius topological and aromaticity switches in octaphyrins by protonation and redox reactions.

Tropone-fused carbaporphyrins

Previous attempts to prepare tropone-fused carbaporphyrins by reacting peroxides with azuliporphyrins under basic conditions afforded benzocarbaporphyrins instead. In this study, a methoxyazulitripyrrane condensed with a pyrrole dialdehyde in the presence of TFA, followed by oxidation with ferric chloride, to give a tropone-fused carbaporphyrin following a spontaneous demethylation. The porphyrinoid gave a modified UV-vis spectrum showing multiple bands in the Soret region, and the proton NMR spectrum showed that it has a reduced diamagnetic ring current in comparison to other carbaporphyrins. The tropone-fused derivative failed to react with tert-butyl hydroperoxide and potassium hydroxide, demonstrating that this type of structure is not an intermediate in the formation of benzocarbaporphyrins. However, the reaction with silver(I) acetate gave the corresponding silver(III) complex. Condensation of the methoxyazulitripyrrane with 2,5-thiophenedicarbaldehyde gave a related tropone-fused thiacarbaporphyrin together with a methoxythiaazuliporphyrin. Treatment of the carbaporphyrins with DBU resulted in the formation of anionic species, while addition of acid afforded dicationic structures. DFT studies were performed on a series of tautomers, protonated species, and anionic structures related to these tropone-fused carbaporphyrins, and NICS calculations were carried out. These results allowed favorable conjugation pathways to be identified. In addition, these studies predicted that protonation initially occurs on the carbonyl moiety rather than on the expected interior pyrrolenine nitrogen atom.

De Proft F. Exploring the structure–aromaticity relationship in Hückel and Möbius N-fused pentaphyrins using DFT

Fine-tuning of the molecular topology and macrocyclic aromaticity of N-fused pentaphyrins by changing the number of π electrons and meso-substituents.

Relative stability of benziporphyrin and naphthiporphyrin tautomers and the emergence of macrocyclic diatropicity

The conformations and relative stabilities of a series of benziporphyrin and naphthiporphyrin tautomers were calculated and the diatropic properties of each of these species were assessed. The results were in good agreement with experimental observations and allow favorable delocalization pathways to be identified.

Fused core-modified planar antiaromatic 32π heptaphyrins: unusual synthesis and structural diversity

The novel core-modified fused heptaphyrins were synthesized and their 4nπ Hückel antiaromatic electronic strucuture was exploited.

Stereoselectivity through a network of non-classical CH weak interactions: a prospective study of a bicyclic organocatalytic scaffold

The preferred stereoisomeric product of this catalytic Diels–Alder reaction is in part determined by noncovalent CH⋯π interactions.