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.

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.

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.

Unravelling the Mechanism and Governing Factors in Lewis Acid and Non-Covalent Diels-Alder Catalysis: Different Perspectives.

In the current literature, many non-covalent interaction (NCI) donors have been proposed that can potentially catalyze Diels-Alder (DA) reactions. In this study, a detailed analysis of the governing factors in Lewis acid and non-covalent catalysis of three types of DA reactions was carried out, for which we selected a set of hydrogen-, halogen-, chalcogen-, and pnictogen-bond donors. We found that the more stable the NCI donor-dienophile complex, the larger the reduction in DA activation energy. We also showed that for active catalysts, a significant part of the stabilization was caused by orbital interactions, though electrostatic interactions dominated. Traditionally, DA catalysis was attributed to improved orbital interactions between the diene and dienophile. Recently, Vermeeren and co-workers applied the activation strain model (ASM) of reactivity, combined with the Ziegler-Rauk-type energy decomposition analysis (EDA), to catalyzed DA reactions in which energy contributions for the uncatalyzed and catalyzed reaction were compared at a consistent geometry. They concluded that reduced Pauli repulsion energy, and not enhanced orbital interaction energy, was responsible for the catalysis. However, when the degree of asynchronicity of the reaction is altered to a large extent, as is the case for our studied hetero-DA reactions, the ASM should be employed with caution. We therefore proposed an alternative and complementary approach, in which EDA values for the catalyzed transition-state geometry, with the catalyst present or deleted, can be compared one to one, directly measuring the effect of the catalyst on the physical factors governing the DA catalysis. We discovered that enhanced orbital interactions are often the main driver for catalysis and that Pauli repulsion plays a varying role.

The Halogen Bond in Weakly Bonded Complexes and the Consequences for Aromaticity and Spin-Orbit Coupling.

The halogen bond complexes CF3X⋯Y and C2F3X⋯Y, with Y = furan, thiophene, selenophene and X = Cl, Br, I, have been studied by using DFT and CCSD(T) in order to understand which factors govern the interaction between the halogen atom X and the aromatic ring. We found that PBE0-dDsC/QZ4P gives an adequate description of the interaction energies in these complexes, compared to CCSD(T) and experimental results. The interaction between the halogen atom X and the π-bonds in perpendicular orientation is stronger than the interaction with the in-plane lone pairs of the heteroatom of the aromatic cycle. The strength of the interaction follows the trend Cl < Br < I; the chalcogenide in the aromatic ring nor the hybridization of the C−X bond play a decisive role. The energy decomposition analysis shows that the interaction energy is dominated by all three contributions, viz., the electrostatic, orbital, and dispersion interactions: not one factor dominates the interaction energy. The aromaticity of the ring is undisturbed upon halogen bond formation: the π-ring current remains equally strong and diatropic in the complex as it is for the free aromatic ring. However, the spin-orbit coupling between the singlet and triplet π→π* states is increased upon halogen bond formation and a faster intersystem crossing between these states is therefore expected.

Hydrogen-Bond-Assisted Diels-Alder Kinetics or Self-Healing in Reversible Polymer Networks? A Combined Experimental and Theoretical Study.

Diels-Alder (DA) cycloadditions in reversible polymer networks are important for designing sustainable materials with self-healing properties. In this study, the DA kinetics of hydroxyl-substituted bis- and tetrafunctional furans with bis- and tris-functional maleimides, both containing ether-functionalized spacers, is investigated by modelling two equilibria representing the endo and exo cycloadduct formation. Concretely, the potential catalysis of the DA reaction through hydrogen bonding between hydroxyl of the furans and carbonyl of the maleimides or ether of the spacers is experimentally and theoretically scrutinized. Initial reaction rates and forward DA rate constants are determined by microcalorimetry at 20 °C for a model series of reversible networks, extended with (i) a hydroxyl-free network and hydroxyl-free linear or branched systems, and (ii) polypropylene glycol additives, increasing the hydroxyl concentration. A computational density-functional theory study is carried out on the endo and exo cycloadditions of furan and maleimide derivatives, representative for the experimental ones, in the absence and presence of ethylene glycol as additive. Additionally, an ester-substituted furan was investigated as a hydroxyl-free system for comparison. Experiment and theory indicate that the catalytic effect of H-bonding is absent or very limited. While increased concentration of H-bonding could in theory catalyze the DA reaction, the experimental results rule out this supposition.

Heteroaryl sulfonamide synthesis: scope and limitations.

Heteroaryl sulfonamides are important structural motifs in the medicinal and agrochemical industries. However, their synthesis often relies on the use of heteroaryl sulfonyl chlorides, which are unstable and toxic reagents. Herein, we report a protocol that allows direct oxidative coupling of heteroaryl thiols and primary amines, readily available and inexpensive commodity chemicals. The transformation proceeds under mild reaction conditions and yields the desired N-alkylated sulfonamides in good yields. N-alkyl heteroaryl sulfonamides can be further transformed using a microwave-promoted Fukuyama-Mitsunobu reaction to N,N-dialkyl heteroaryl sulfonamides. The developed protocols thus enable the preparation of previously difficult to prepare sulfonamides (toxic reagents, harsh conditions, and low yields) under mild conditions.

Can the Philicity of Radicals Be Influenced by Oriented External Electric Fields?

Herein, the effects of an electric field on radicals are investigated for a set of model radicals with varying complexity. An electric field impacts the intrinsic philicity of a radical, as quantified by the global electrophilicity index, ω. The extent of change in philicity depends on the directionality and strength of the applied electric field and the dipole moment and polarizability of the radical.

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.

A Combined Experimental/Quantum-Chemical Study of Tetrel, Pnictogen, and Chalcogen Bonds of Linear Triatomic Molecules.

Linear triatomic molecules (CO2, N2O, and OCS) are scrutinized for their propensity to form perpendicular tetrel (CO2 and OCS) or pnictogen (N2O) 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 sp1 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.

Preparation of 4'-Spirocyclobutyl Nucleoside Analogues as Novel and Versatile Adenosine Scaffolds.

Despite the large variety of modified nucleosides that have been reported, the preparation of constrained 4'-spirocyclic adenosine analogues has received very little attention. We discovered that the [2+2]-cycloaddition of dichloroketene on readily available 4'-exo-methylene furanose sugars efficiently results in the diastereoselective formation of novel 4'-spirocyclobutanones. The reaction mechanism was investigated via density functional theory (DFT) and found to proceed either via a non-synchronous or stepwise reaction sequence, controlled by the stereochemistry at the 3'-position of the sugar substrate. The obtained dichlorocyclobutanones were converted into nucleoside analogues, providing access to a novel class of chiral 4'-spirocyclobutyl adenosine mimetics in eight steps from commercially available sugars. Assessment of the biological activity of designed 4'-spirocyclic adenosine analogues identified potent inhibitors for protein methyltransferase target PRMT5.

Acceleration of Inverse Molecular Design by Using Predictive Techniques.

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

Metal-Free Cyclization of ortho-Nitroaryl Ynamides and Ynamines towards Spiropseudoindoxyls.

An efficient metal-free cascade reaction between 1-dibromovinyl-2-nitro-substituted arenes and secondary amines results in the formation of polycyclic pseudoindoxyls in a single step. The reaction mechanism leading to these fused ring systems was investigated, and is believed to involve the initial formation of nitroarylated ynamines/ynamides. These intermediates cycloisomerize towards N-alkenyl-tethered 2-aminoisatogens via a carbene intermediate as demonstrated by QTAIM (quantum theory of atoms in molecules) and ELF (electron localization function) analysis. A subsequent intramolecular dipolar cycloaddition afforded the title compounds.

Characterization of chalcogen bonding interactions via an in-depth conceptual quantum chemical analysis.

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.

Heterolytic Splitting of Molecular Hydrogen by Frustrated and Classical Lewis Pairs: A Unified Reactivity Concept.

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 H2 splitting reaction of various Lewis pairs are in good agreement with the experimentally observed feasibility of H2 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 H2 splitting characterizes activation processes corresponding to high-energy transition states and, accordingly, is not operative for H2-activating systems. One of the criteria of H2-activation, actually, is the availability of a low-energy transition state that represents a different H2 splitting mechanism, in which the electrostatic field generated in the cavity of Lewis pair plays a critical role: to induce a strong polarization of H2 that facilities an efficient end-on acid-H2 interaction and to stabilize the charge separated "H+-H-" moiety in the transition state.

Exploiting the σ-Hole Concept: An Infrared and Raman-Based Characterization of the S⋠⋠⋠O Chalcogen Bond between 2,2,4,4-Tetrafluoro-1,3-dithiethane and Dimethyl Ether.

In the last decade, halogen bonds, noncovalent interactions formed between positive regions in the electrostatic potential on halogen atoms, often referred to as σ-holes, and electron-rich sites, have gained a lot of interest. Recently, this interest has been expanded towards interactions with Group V and Group VI elements, giving rise to pnicogen and chalcogen bonds. Although chalcogen bonds have already shown some promising results for applications in crystallography and catalysis, experimental results characterising these noncovalent interactions remain scarce. In this combined experimental and theoretical study, original data allowing the characterization of S⋅⋅⋅O chalcogen bonds is obtained by studying the 1:1 molecular complexes between 2,2,4,4-tetrafluoro-1,3-dithiethane (C2 F4 S2 ) and dimethyl ether (DME). Ab initio calculations of the C2 F4 S2 ⋅DME dimer yield two stable chalcogen-bonded isomers, the difference being the presence or absence of secondary F⋅⋅⋅H interactions. Liquid-krypton solutions containing C2 F4 S2 and DME were studied using FTIR and Raman spectroscopy. Upon subtraction of rescaled monomer spectra, clear complex bands are observed. The observed complexation shifts agree favourably with the ab initio calculated shifts of the chalcogen-bonded complexes. The 1:1 stoichiometry of the complex is confirmed and a complexation enthalpy of -13.5(1) kJ mol-1 is found, which is in good agreement with the calculated values. A Ziegler-Rauk energy decomposition analysis revealed that electrostatic interactions prominently dominate over orbital interactions. Nevertheless, significant charge transfer occurs from the oxygen in DME to one of the sulfur atoms in C2 F4 S2 and the carbon along the extension of the chalcogen bond.

Tuning the HOMO-LUMO Energy Gap of Small Diamondoids Using Inverse Molecular Design.

Functionalized diamondoids show great potential as building blocks for various new optoelectronic applications. However, until now, only simple mono and double substitutions were investigated. In this work, we considered up to 10 and 6 sites for functionalization of the two smallest diamondoids, adamantane and diamantane, respectively, in search for diamondoid derivatives with a minimal and maximal HOMO-LUMO energy gap. To this end, the energy gap was optimized systematically using an inverse molecular design methodology based on the best-first search algorithm combined with a Monte Carlo component to escape local optima. Adamantane derivatives were found with HOMO-LUMO gaps ranging from 2.42 to 10.63 eV, with 9.45 eV being the energy gap of pure adamantane. For diamantane, similar values were obtained. The structures with the lowest HOMO-LUMO gaps showed apparent push-pull character. The push character is mainly formed by sulfur or nitrogen dopants and thiol groups, whereas the pull character is predominantly determined by the presence of electron-withdrawing nitro or carbonyl groups assisted by amino and hydroxyl groups via the formation of intramolecular hydrogen bonds. In contrast, maximal HOMO-LUMO gaps were obtained by introducing numerous electronegative groups.

Molecular Property Optimizations with Boundary Conditions through the Best First Search Scheme.

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.

Conceptual quantum chemical analysis of bonding and noncovalent interactions in the formation of frustrated Lewis pairs.

The contributions of covalent and noncovalent interactions to the formation of classical adducts of bulky Lewis acids and bases and frustrated Lewis pairs (FLPs) were scrutinized by using various conceptual quantum chemical techniques. Significantly negative complexation energies were calculated for fourteen investigated Lewis pairs containing bases and acids with substituents of various sizes. A Ziegler-Rauk-type energy decomposition analysis confirmed that two types of Lewis pairs can be distinguished on the basis of the nature of the primary interactions between reactants; dative-bond formation and concomitant charge transfer from the Lewis base to the acid is the dominant and most stabilizing factor in the formation of Lewis acid-base adducts, whereas weak interactions are the main thermodynamic driving force (>50 %) for FLPs. Moreover, the ease and extent of structural deformation of the monomers appears to be a key component in the formation of the former type of Lewis pairs. A Natural Orbital for Chemical Valence (NOCV) analysis, which was used to visualize and quantify the charge transfer between the base and the acid, clearly showed the importance and lack of this type of interaction for adducts and FLPs, respectively. The Noncovalent Interaction (NCI) method revealed several kinds of weak interactions between the acid and base components, such as dispersion, π-π stacking, C-H⋅⋅⋅π interaction, weak hydrogen bonding, halogen bonding, and weak acid-base interactions, whereas the Quantum Theory of Atoms in Molecules (QTAIM) provided further conceptual insight into strong acid-base interactions.

Pushing the boundaries of intrinsically stable radicals: inverse design using the thiadiazinyl radical as a template.

In this study, for the first time inverse design was applied to search for the intrinsically most stable radical system in a predefined chemical space of enormous size by scanning in a rational way that entire chemical space. The focus was predominantly on thermodynamic stabilization effects, such as stabilization through resonance. Two different properties were optimized: a newly introduced descriptor called the radical delocalization value and the intrinsic stability via a previously established bond dissociation enthalpy model. The thiadiazinyl radical was chosen as case study of this new approach of inverse design in stable radical chemistry. The resulting optimal structure is found to be highly stable, intrinsically more so than other well-known stable radicals, such as verdazyls and N,N-diphenyl-N'-picrylhydrazyl, and even rivaling the intrinsic stability of nitrogen monoxide.