Illuminating anti-hydrozirconation: controlled geometric isomerization of an organometallic species

A general strategy to enable the formal anti-hydrozirconation of arylacetylenes is reported that merges cis-hydrometallation using the Schwartz Reagent (Cp₂ZrHCl) with a subsequent light-mediated geometric isomerization at λ = 400 nm. Mechanistic delineation of the contra-thermodynamic isomerization step indicates that a minor reaction product functions as an efficient in situ generated photocatalyst. Coupling of the E-vinyl zirconium species with an alkyne unit generates a conjugated diene: this has been leveraged as a selective energy transfer catalyst to enable E → Z isomerization of an organometallic species. Through an Umpolung metal–halogen exchange process (Cl, Br, I), synthetically useful vinyl halides can be generated (up to Z : E = 90 : 10). This enabling platform provides a strategy to access nucleophilic and electrophilic alkene fragments in both geometric forms from simple arylacetylenes.

A general strategy to enable the formal anti-hydrozirconation of arylacetylenes is reported that merges cis-hydrometallation using the Schwartz Reagent (Cp₂ZrHCl) with a subsequent light-mediated geometric isomerization at λ = 400 nm.

The venerable Schwartz reagent (Cp₂ZrHCl) is totemic in the field of hydrometallation,¹ where reactivity is dominated by syn-selective M–H addition across the π-bond.^2,3 This mechanistic foundation can be leveraged to generate well-defined organometallic coupling partners that are amenable to stereospecific functionalization. Utilizing terminal alkynes as readily available precursors,⁴ hydrozirconation constitutes a powerful strategy to generate E-configured vinyl nucleophiles that, through metal–halogen exchange, can be converted to vinyl electrophiles in a formal Umpolung process.⁵ Whilst this provides a versatile platform to access the electronic antipodes of the E-isomer, the mechanistic course of addition renders access to the corresponding Z-isomer conspicuously challenging. To reconcile the synthetic importance of this transformation with the intrinsic challenges associated with anti-hydrometallation and metallometallation,⁶ it was envisaged that a platform to facilitate geometric isomerization⁷ would be of value. Moreover, coupling this to a metal–halogen exchange would provide a simple Umpolung matrix to access both stereo-isomers from a common alkyne precursor (Fig. 1).Open in a separate windowFig. 1The stereochemical course of alkyne hydrometallation using the Schwartz reagent and an Umpolung platform to generate both stereo-isomers from a common alkyne precursor.Confidence in this conceptual blueprint stemmed from a report by Erker and co-workers, in which irradiating the vinyl zirconium species derived from phenyl acetylene (0.5 M in benzene) with a mercury lamp (Philips HPK 125 and Pyrex filter) induced geometric isomerization.⁸ Whilst Hg lamps present challenges in terms of safety, temperature regulation, cost and wavelength specificity, advances in LED technology mitigate all of these points. Therefore, a process of reaction development was initiated to generalize the anti-hydrozirconation of arylacetylenes. Crucial to the success of this venture was identifying the light-based activation mode that facilitates alkene isomerization. Specifically, it was necessary to determine whether this process was enabled by direct irradiation of the vinyl zirconium species, or if the E → Z directionality results from a subsequent selective energy transfer process involving a facilitator. Several accounts of the incipient vinyl zirconium species reacting with a second alkyne unit to generate a conjugated diene have been disclosed.^9,10 It was therefore posited that the minor by-product diene may be a crucial determinant in driving this isomerization (Fig. 2).Open in a separate windowFig. 2A working hypothesis for the light-mediated anti-hydrozirconation via selective energy transfer catalysis.To advance this working hypothesis and generalize the formal anti-hydrozirconation process, the reaction of Cp₂ZrHCl with 1-bromo-4-ethynylbenzene (A-1) in CH₂Cl₂ was investigated (‡ for full details). This generates a versatile electrophile for downstream synthetic applications. Gratifyingly, after only 15 minutes, a Z : E-composition of 50 : 50 was reached (entry 1) and, following treatment with NBS, the desired vinyl bromide (Z)-1 was obtained in 76% yield (isomeric mixture) over the two steps. Further increasing the irradiation by 15 minute increments (entries 2–4) revealed that the optimum reaction time for the isomerization is 45 minutes (74%, Z : E = 73 : 27, entry 3). Extending the reaction time to 60 minutes (entry 4, 54%) did not lead to an improvement in selectivity and this was further confirmed by irradiating the reaction mixture for 90 minutes (entry 5). In both cases, a notable drop in yield was observed and therefore the remainder of the study was performed using the conditions described in entry 3. Next, the influence of the irradiation wavelength on the isomerization process was examined (entries 6–11). From a starting wavelength of λ = 369 nm, which gave a Z : E-ratio of 27 : 73 (entry 6), a steady improvement was observed by increasing the wavelength to λ = 374 nm (Z : E = 44 : 56, entry 7) and λ = 383 nm (Z : E = 53 : 47, entry 8). The selectivity reached a plateau at λ = 400 nm, with higher wavelengths proving to be detrimental (Z : E = 60 : 40 at λ = 414 nm, entry 9; Z : E = 26 : 74 at λ = 435 nm, entry 10). It is interesting to note that at λ = 520 nm, Z-1 was not detected by ¹H NMR (entry 11).Reaction optimization^a

Hydroxy-directed fluorination of remote unactivated C(sp3)–H bonds: a new age of diastereoselective radical fluorination

We report a photochemically induced, hydroxy-directed fluorination that addresses the prevailing challenge of high diastereoselectivity in this burgeoning field. Numerous simple and complex motifs showcase a spectrum of regio- and stereochemical outcomes based on the configuration of the hydroxy group. Notable examples include a long-sought switch in the selectivity of the refractory sclareolide core, an override of benzylic fluorination, and a rare case of 3,3′-difluorination. Furthermore, calculations illuminate a low barrier transition state for fluorination, supporting our notion that alcohols are engaged in coordinated reagent direction. A hydrogen bonding interaction between the innate hydroxy directing group and fluorine is also highlighted for several substrates with ¹⁹F–¹H HOESY experiments, calculations, and more.

We report a photochemical, hydroxy-directed fluorination that addresses the prevailing challenge of high diastereoselectivity. Numerous motifs showcase a range of regio- and stereochemical outcomes based on the configuration of the hydroxy group.

The hydroxy (OH) group is treasured and versatile in chemistry and biology.¹ Its ubiquity in nature and broad spectrum of chemical properties make it an attractive source as a potential directing group.² The exploitation of the mild Lewis basicity exhibited by alcohols has afforded several elegant pathways for selective functionalization (e.g., Sharpless epoxidation,³ homogeneous hydrogenation,⁴ cross-coupling reactions,⁵ among others⁶). Recently, we reported a photochemically promoted carbonyl-directed aliphatic fluorination, and most notably, established the key role that C–H⋯O hydrogen bonds play in the success of the reaction.⁷ Our detailed mechanistic investigations prompt us to postulate that other Lewis basic functional groups (such as –OH) can direct fluorination in highly complementary ways.⁸ In this communication, we report a hydroxy-directed aliphatic fluorination method that exhibits unique directing properties and greatly expands the domain of radical fluorination into the less established realm governing high diastereoselectivity.⁹Our first inclination that functional groups other than carbonyls may influence fluorination regiochemical outcomes was obtained while screening substrates for our published ketone-directed radical-based method (Scheme 1).^8a In this example, we surmised that oxidation of the tertiary hydroxy group on substrate 1 cannot occur and would demonstrate functional group tolerance (directing to C11, compound 2). Surprisingly, the two major regioisomers (products 3 and 4) are derivatized by Selectfluor (SF) on C12 and C16 – indicative of the freely rotating hydroxyl directing fluorination. Without an obvious explanation of how these groups could be involved in dictating regiochemistry, we continued the mechanistic study of carbonyl-directed fluorination (Scheme 2A). We established that the regioselective coordinated hydrogen atom abstraction occurs by hydrogen bonding between a strategically placed carbonyl and Selectfluor radical dication (SRD).⁷ However, we noted that the subsequent radical fluorination is not diastereoselective due to the locally planar nature of carbonyl groups. Thus, we posed the question: are there other directing groups that can provide both regio- and diastereoselectivity? Such a group would optimally be attached to a sp³ hybridized carbon; thus the “three dimensional” hydroxy carbon logically comes to mind as an attractive choice, and Scheme 1 illustrates the first positive hint.Open in a separate windowScheme 1Observed products for the fluorination of compound 1.Open in a separate windowScheme 2(A) Proposed mechanism, (B) β-caryophyllene alcohol hypochlorite derivative synthetic probe, (C) isodesmic relation of transition states showing the general importance of the hydroxy group to reactivity (ωB97xd/6-31+G*), and (D) ¹H NMR experiment with Selectfluor and various additives at different concentrations.We began our detailed study with a simple substrate that contains a tertiary hydroxyl group. Alcohol 5 was synthesized stereoselectively by the reaction of 3-methylcyclohexanone, FeCl₃, and 4-chlorophenylmagnesium bromide;¹⁰ the 4-chlorophenyl substituent allows for an uncomplicated product identification and isolation (aromatic chromophore). We sought to determine optimal reaction conditions by examination of numerous photosensitizers, bases, solvents, and light sources (7 Although we utilize cool blue LEDs (sharp cutoff ca. 400 nm), CFLs (small amount of UVB (280–315 nm) and UVA (315–400 nm)) are useable as well.¹¹ A mild base additive was also found to neutralize adventitious HF and improve yields in the substrates indicated (

Correction: Highly selective acid-catalyzed olefin isomerization of limonene to terpinolene by kinetic suppression of overreactions in a confined space of porous metal–macrocycle frameworks

Correction for ‘Highly selective acid-catalyzed olefin isomerization of limonene to terpinolene by kinetic suppression of overreactions in a confined space of porous metal–macrocycle frameworks’ by Wei He et al., Chem. Sci., 2022, 13, 8752–8758, https://doi.org/10.1039/d2sc01561g.

The authors regret that there were errors in Fig. 2, Fig. 5 and Fig. 6 in the original article and Fig. S18 of the ESI. The stereochemistry of the chemical structural formulas for (−)-α-pinene (6) and (−)-β-pinene (7) was incorrectly reversed. The correct versions of the figures are shown below, and in the updated version of the ESI.Open in a separate windowFig. 2Metal–macrocycle framework (MMF). (a) Self-assembly of asymmetrically twisted Pd^II-macrocycles into (b) a porous crystal MMF (sticks model) with five enantiomeric pairs of binding pockets (surface model). (c) Previously reported site-selective adsorption of (−)-α-pinene (6) (space-filling model) on the channel surface of the MMF.¹ Blue, yellow, red, or black dashed circles indicate the ceiling-, side-, bottom-, or tubular-pockets of the MMF, respectively. MMF: Pd, yellow; Cl, green; N, blue; C, grey. 6: C, pink; H, white. Hydrogen atoms attached to the MMF were omitted for clarity. Green or blue surface represents exposed Cl or N–H groups of the MMF, respectively.Open in a separate windowFig. 5Investigation of the inhibitory effects of additives on the isomerization reaction of 2 using 2-NBSA@MMF at 25 °C for 102 h.Open in a separate windowFig. 6Crystallographic study of MMFs soaked in (a) a CHCl₃ solution containing 1 (1.0 M), (b) a CHCl₃ solution containing 2 (1.0 M), and (c) a CH₃CN solution containing 7 (1.0 M). MMF: stick model or surface model; 1 and 7: space-filling model; water and CHCl₃: stick model. Red dashed circles indicate the bottom pocket of the MMF. MMF: Pd, yellow; Cl, green; N, blue; C, grey. 1: C, yellow; H, white. 7: C, pink; H, white. Water and CHCl₃: O, red; H, white; C, grey; Cl, green. Hydrogen atoms attached to the MMF were omitted for clarity. Green and blue surface represents exposed Cl and N–H groups of the MMF, respectively.The Royal Society of Chemistry apologises for these errors and any consequent inconvenience to authors and readers.     

Radical 1,2,3-tricarbofunctionalization of α-vinyl-β-ketoesters enabled by a carbon shift from an all-carbon quaternary center

Herein, we report an intermolecular, radical 1,2,3-tricarbofunctionalization of α-vinyl-β-ketoesters to achieve the goal of building molecular complexity via the one-pot multifunctionalization of alkenes. This reaction allows the expansion of the carbon ring by a carbon shift from an all-carbon quaternary center, and enables further C–C bond formation on the tertiary carbon intermediate with the aim of reconstructing a new all-carbon quaternary center. The good functional group compatibility ensures diverse synthetic transformations of this method. Experimental and theoretical studies reveal that the excellent diastereoselectivity should be attributed to the hydrogen bonding between the substrates and solvent.

Herein, we report an intermolecular, radical 1,2,3-tricarbofunctionalization of α-vinyl-β-ketoesters to achieve the goal of building molecular complexity via the one-pot multifunctionalization of alkenes.

A leading motive for the impressive achievements in the area of assembling molecular complexity is the transformation of simple feedstock chemicals into complex molecular skeletons with superior bioactive properties. In this respect, the direct functionalization of alkenes has been demonstrated as one of the most effective and simple strategies to meet this criterion at a high level. While the difunctionalization of alkenes in a one-pot process is the major theme of considerable interest in this field,¹ the multifunctionalization of alkenes,² for example, a 1,2,3-trifunctionalization of alkenes, has the power to simultaneously incorporate multifunctional groups. Therefore, this multifunctionalization reaction model can be regarded as an efficient and novel strategy to afford molecules with high structural diversity and complexity. However, such methods are elusive.During the last decades, radical alkene functionalizations have been revealed to be a powerful tool for building complex molecular frameworks by employing a radical initiator, a transition metal catalyst, or a photocatalyst.^1f–i However, only several successful methods for the radical multifunctionalization of alkenes have been achieved. For example, the Studer group reported an elegant 1,2-boryl shift-enabled radical 1,2,3-trifunctionalization of allylboronic esters using AIBN as the radical initiator (Fig. 1a).³ Shi et al. disclosed an excellent photocatalytic perfluoroalkylation of a vinyl-substituted all-carbon quaternary center through 1,2-aryl migration (Fig. 1b).⁴ Herein, we report a new one-pot protocol to realize an intermolecular, radical 1,2,3-tricarbofunctionalization of α-vinyl-β-ketoesters through a cascade process of deconstruction–reconstruction of the all-carbon quaternary center (Fig. 1c).⁵Open in a separate windowFig. 1Radical 1,2,3-trifunctionalization of alkenes. (a) Studer''s work; (b) Shi''s work; (c) This work.The direct incorporation of a fluorine atom or fluorinated moieties into organic compounds has been extensively investigated and proved to be an significant synthetic strategy in the field of discovering new pharmaceuticals.⁶ Recently, we are interested in the radical functionalization of alkenes with fluoroalkyl groups,⁷ and we envisioned that, different from the typical Dowd–Beckwith⁸ ring expansion reaction,⁹ the addition of a fluoroalkyl radical to the C C double bond would generate an adduct radical species I, which will transform into the radical intermediate II upon ring expansion (Fig. 1c). Finally, the cascade C–C coupling affords the product with a reconstructed all-carbon quaternary center. However, there are several challenging issues that need to be addressed: (1) the carbon shift from an all-carbon quaternary center to afford a tertiary carbon center which is bulkier than the tertiary carbon center formed in a typical Dowd–Beckwith ring expansion reaction; (2) the reconstruction of all-carbon quaternary center from tertiary carbon radical II will meet the associated conformational restriction and steric congestion; (3) side reactions, such as 1,2-radical addition to the alkenyl group, homolytic couplings of the carbon radical intermediates I and II, and direct H-atom abstraction;¹⁰ (4) how to control the diastereomeric ratio of the products. To meet these challenges, we developed a novel method for the 1,2,3-trifunctionalization of alkenes using alkynyl triflones as both the CF₃ (ref. 6) and alkynyl sources, providing the ring-expanded cyclic β-ketoesters with excellent diastereoselectivity and functional group diversity. In addition, good functional group compatibility of this method was observed, which ensures the diverse synthetic transformations. Moreover, hydrogen bonding between the substrates and 2,2,2-trifluoroethanol solvent was revealed to be the key factor for the excellent diastereoselectivity obtained in this reaction, and this result was confirmed by both experimental and theoretical studies.This study began by surveying radical initiators for 1,2,3-tricarbofunctionalizing α-vinyl-β-ketoester 1a with alkynyl triflone 2a¹¹ (12 (13 dramatically increased the diastereoselectivity and (±)-3a could be obtained in an identical yield with an even higher dr value (dr > 20 : 1) (14 Without the addition of a radical initiator, a reaction did not happen ( EntrySolventYield^b (%)1EA60 (dr = 13 : 1)^c2EA55 (dr = 11 : 1)^d3EA63 (dr = 12 : 1)4MTBE45 (dr = 10 : 1)5DCE63 (dr = 15 : 1)6TolueneTrace7DMFTrace8MeOHTrace9TFE63 (dr > 20 : 1)10^eTFE60 (dr > 20 : 1)11^fTFE56 (dr > 20 : 1)12^gTFE70 (dr > 20 : 1)13^hTFE76 (65)ⁱ (dr > 20 : 1)14^jTFE71 (dr > 20 : 1)15TFETraceOpen in a separate window^aReaction conditions: alkene 1a (0.2 mmol, 1 equiv.), 2a (0.6 mmol, 3.0 equiv.), and AIBN (0.3 equiv.) in 3 mL of solvent at 85 °C for 18 h in a sealed tube under a nitrogen atmosphere.^bCrude yield and crude diastereomeric ratio were determined by ¹⁹F NMR.^cLPO was used as the initiator.^dBPO was used as the initiator.^eThe reaction was performed at 100 °C.^fThe reaction was performed at 120 °C.^gAIBN (60 mol%) was used.^h2a (3.0 equiv.) and AIBN (60 mol%) were added as two equal portions with an interval of 9 h.ⁱIsolated yield in parentheses.^j2a (3.0 equiv.) and AIBN (60 mol%) were added as three equal portions with an interval of 6 h.Under optimal conditions, a diverse array of α-vinyl-β-ketoesters serve as substrates in this metal-free deconstruction–construction of all-carbon quaternary centers for the synthesis of carbon-ring expanded cyclic β-ketoesters (Fig. 2). In most of the cases, excellent diastereoselectivities (dr > 20 : 1) were observed by crude ¹⁹F NMR analysis. Substrates with the substituents at the 5- or 6-position of the α-vinyl-β-ketoesters generally produced the corresponding product (±)-3 in higher yields than those with the substituents at the 4-position. Apart from the carbonyl group and the ester group, functional groups such as chloride ((±)-3b and (±)-3f), fluoride ((±)-3c), a methoxyl group ((±)-3d and (±)-3h), a methyl group ((±)-3e and (±)-3g) and a phenyl group ((±)-3i) can be tolerated under the reaction conditions. Notably, the phenyl ring of the core structure with two substituents reacted smoothly to afford the corresponding products ((±)-3j and (±)-3k). When substrate 1l that lacks the fused benzene ring was used for this carbon-ring expansion reaction, a dramatical loss of diastereoselectivity was detected, presumably because of the feasible interconversion of the boat and chair conformations of the intermediate. Substrates with an ethyl ester or a benzyl ester group, as opposed to a methyl ester group, delivered the corresponding products ((±)-3m and (±)-3n) with moderate yields and excellent diastereoselectivity. When the CH₂ unit of the six membered-ring was replaced by a CMe₂ group, only a trace amount of the desired product (±)-3o was detected. A reaction with the purpose of realizing an extension from the six-membered ring was also carried out and (±)-3p was obtained, although with a low yield and low diastereoselectivity. Notably, the diastereochemistries of products (±)-3e and (±)-3h have been confirmed by X-ray crystallography.Open in a separate windowFig. 2Substrate scope of α-vinyl-β-ketoesters. ^aThe reaction was performed with 1p and 2b.The scope with respect to the alkynyl triflones was also investigated and the results are summarized in Fig. 3. Generally, substituents on the phenyl ring of the arylethynyl moiety have little impact on the yields of the corresponding products. The functional groups at the para-, meta-, or ortho-position of the phenyl ring produced the desired products ((±)-4a–(±)-4k) with excellent diastereoselectivities. Furthermore, the method is compatible with alkynyl triflones that have a thienyl group or a perfluorobutyl group and the reactions afforded the product ((±)-4l or (±)-4m) with an excellent dr value, respectively. However, when the arylethynyl moiety was replaced by an alkylethynyl or a silylethynyl part, the reaction failed to produce the targeted tricarbofunctionalization product ((±)-4n or (±)-4o).¹⁵ Moreover, when triflic azide or (Z)-TolCH CHSO₂CF₃ was used in place of the alkynyl triflone, the desired product was not obtained and most of the starting material was recovered. Notably, the diastereochemistry of product (±)-4a has been confirmed by X-ray crystallography.Open in a separate windowFig. 3Substrate scope of alkynyl triflones.This 1,2,3-trifunctionalization reaction not only allows the deconstruction and reconstruction of all-carbon quaternary centers, but features good functional group tolerance and excellent diastereoselectivity. Regarding the diverse reactivities of these functional groups, many valuable synthetic transformations have been successfully achieved (Fig. 4). For example, the C–C triple bond of (±)-4a can be completely reduced to a CH₂CH₂ unit ((±)-5) in the presence of hydrogen and a Pd/C catalyst,¹⁶ while the selective reduction of (±)-4a gives rise to a Z-alkene (±)-6 when quinoline is added as an additive for the Lindlar reduction.¹⁷ The diastereochemistry of (±)-6 has been confirmed by X-ray crystallography. The selective reducing methods afford formal approaches for radical 1,3-trifluoromethylalkylation and 1,3-trifluoromethylalkenylation of α-vinyl-β-ketoesters, respectively, to produce the corresponding products which are otherwise difficult to obtain. In addition, the C–C triple bond can be oxidized under oxidative conditions with RuCl₃/NaIO₄, and (±)-4a can be smoothly transformed into the trifluoromethylated triketone (±)-7 in 65% yield.¹⁸ With a large excess amount of reducing agent LiAlH₄, the carbonyl group and the ester group, together with the C–C triple bond, can be unexpectedly reduced simultaneously, affording the alkenyl diol (±)-8 in excellent regioselectivity. The hydrolysis process under basic conditions provided a reliable method for access to a free carboxylic acid (±)-9. Interestingly, when the reaction was performed under milder conditions compared to those for the synthesis of (±)-8, (±)-4a was successfully converted into an alkynyl diol (±)-10, which can be cyclized into a spiro compound (±)-11 (ref. 19) and an endocyclic compound (±)-12,²⁰ respectively. Notably, in the majority of these cases, the excellent diastereoselectivity was reserved. These synthetic applications can demonstrate the significant value of this method.Open in a separate windowFig. 4Synthetic transformations.In order to gain some mechanistic insights into this radical cascade reaction, subsequent efforts have been made (Fig. 5). First, the detection of trifluoromethylated toluene (with toluene as the solvent, Fig. 5a, see ESI^† for details). Second, we were curious about the excellent diastereoselectivity associated with the use of TFE as the solvent. As can be seen in Fig. 5b, ¹H NMR titration of 1a with increasing amounts of TFE showed a chemical shift of the resonance signal corresponding to protons. The 2D NOESY spectrum indicates the existence of an interaction between 1a and TFE (Fig. 5c). Moreover, Job plot studies by both ¹H NMR and ¹⁹F NMR imply a 1 : 1.5 stoichiometry of the complex adduct resulting from 1a and TFE (Fig. 5d). These mechanistic studies strongly suggest that the excellent diastereoselectivity of this reaction might be attributed to the hydrogen bonding between TFE and the α-vinyl-β-ketoester.Open in a separate windowFig. 5Mechanism studies. (a) Radical probe; (b) ¹H NMR titration; (c) 2D NOESY; (d) Job plot studies.On the other hand, density functional theory (DFT) calculations have also been performed at the B3LYP-D3(SMD)/Def2-TZVP//B3LYP-D3(SMD)/Def2-SVP level of theory in the TFE solvent model to further investigate the reaction pathways (Fig. 6). On the basis of the experimental results, herein, the radical pathway was considered. Initially, the CF₃ radical addition onto 1a was calculated, and a transition state, TS1, was located with a free energy barrier of 10.9 kcal mol⁻¹ to deliver the radical intermediate int1 with an exergonicity of 20.5 kcal mol⁻¹. Then, a bicyclic transition state, TS2,²¹ with a barrier of 11.0 kcal mol⁻¹ through a concerted 1,2-shift route was found to be the lower barrier TS for int2 formation than that of the addition to 2b for the byproduct (see Fig. S5 in ESI^†), which is consistent with the experimental results of the mainly hexacyclic products. Moreover, the intrinsic reaction coordinate (IRC) calculations and the root mean square (RMS) gradient of the potential energy surface from TS2 suggested that no transition state for the formation of the previously proposed strained alkoxyl radical was found. Next, the radical intermediate int2 attacking 2b was calculated. To understand the diastereoselectivity of this step, the transition states of the addition of 2b onto the Re and Si faces of C3 in int2 were located with barriers of 12.5 and 17.4 kcal mol⁻¹ (TS3 and TS3′), respectively. It is noteworthy that the torsion angle of C1–C2–C3–C4 in TS3′ is −62.3°, larger than that of −40.9° in int2 and −49.0° in TS3, indicating that the distortion factor in TS3′ is large due to the steric effect from the trifluoroethyl group in int2 and, therefore, increases the barrier. The transition states of 2b addition were also optimized in solvents DCE and EA, and the free energy barrier differences between TS3 and TS3′ [ΔG^‡ = G^‡(TS3′) − G^‡(TS3)] are 3.6 and 3.0 kcal mol⁻¹, respectively, in agreement with the experimental observations. Finally, dissociation of a SO₂ molecule with a CF₃ radical from int3 to deliver the product was conducted, and a transition state TS4 with a much lower barrier of only 7.1 kcal mol⁻¹ was located, which led to the major product (±)-4a with a relative free enthalpy of −51.6 kcal mol⁻¹.Open in a separate windowFig. 6Gibbs free energy profile for the synthesis of 4a in the TFE solvent model.     

Palladium/GF-Phos-catalyzed asymmetric carbenylative amination to access chiral pyrrolidines and piperidines

The cross-coupling of N-tosylhydrazones has emerged as a powerful method for the construction of structurally diverse molecules, but the development of catalytic enantioselective versions still poses considerable challenges and only very limited examples have been reported. We herein report an asymmetric palladium/GF-Phos-catalyzed carbenylative amination reaction of N-tosylhydrazones and (E)-vinyl iodides pendent with amine, which allows facile access to a range of chiral pyrrolidines and piperidines in good yields (45–93%) with up to 96.5 : 3.5 er. Moreover, mild conditions, general substrate scope, scaled-up preparation, as well as the efficient synthesis of natural product (−)-norruspoline are practical features of this method.

An efficient asymmetric palladium/GF-Phos-catalyzed carbenylative amination reaction to access structurally diverse chiral pyrrolidines and piperidines in good yields with high chemo-, regio- and enantioselectivities has been developed.

N-tosylhydrazones, readily prepared from aldehydes or ketones, served as a safe source of carbene precursors and have attracted much attention of chemists.¹N-tosylhydrazone-mediated applications have been continuously developed, such as cyclopropanation or cyclopropenation, X–H insertion, ylide formation, cycloaddition, aza-Wacker-type cyclization, asymmetric allylic substitution, etc.² Among them, transition-metal-catalyzed cross-coupling is one of the powerful protocols for C–X or C C bond formation in organic synthesis involving versatile intermediates, of which in situ generation of diazo compounds and carbene migratory insertion are considered key steps.^3–5 Over the past decades, considerable progress has been made in the asymmetric cross-coupling reactions of N-tosylhydrazones with various coupling partners, including cyclobutanols, terminal alkynes, silacyclobutanes and so on.⁴ Relatively, only a few examples focus on the cross-coupling reactions of aryl halides with N-tosylhydrazones involving benzyl metal intermediates [Scheme 1A, eqn. (a)].⁶ For example, Gu,^6a Wu,^6b Lassaletta^6c and coworkers have developed a palladium-catalyzed asymmetric synthesis of axial chiral compounds from aryl bromides and N-tosylhydrazones, ending with β-H elimination. Very recently, we realized palladium/GF-Phos catalyzed asymmetric three component cross-coupling reactions of aryl halides, N-tosylhydrazones, with terminal alkynes.^6f In contrast, much less progress has been made in N-tosylhydrazone-based carbenylative insertions from vinyl halides, which would generate a π-allylic metal intermediate followed by nucleophile attack, providing a unique approach for building C–X bonds, especially for N-heterocyclic compounds [Scheme 1A, eqn. (b)].⁷N-heterocycles are important structural motifs for the development of various types of valuable chemicals and materials.⁸ Importantly, optically active 2-substituted pyrrolidine and piperidine derivatives are privileged scaffolds in many natural products and pharmaceuticals with a wide range of biological activities,⁹ as well as the backbone of organocatalysts in asymmetric catalysis (Fig. 1).¹⁰Open in a separate windowScheme 1Asymmetric transition-metal-catalyzed carbenylative cross-coupling reactions.Open in a separate windowFig. 1Selected natural products and pharmaceuticals containing chiral 2-substituted pyrrolidine and piperidine units.Notably, Van Vranken and coworkers reported an elegant palladium-catalyzed carbenylative amination reaction of N-tosylhydrazones and (E)-vinyl iodides pendent with amine, providing facile access to pyrrolidine and piperidine ring systems that are common to alkaloid natural products (Scheme 1B).¹¹ Unfortunately, only up to 58.5 : 41.5 er was obtained after they made a lot of efforts to screen a series of chiral phosphine ligands, indicating that this asymmetric reaction indeed poses considerable challenges in addition to competitive side reactions such as the dimerization of vinyl iodides,¹² the formation of diene via the palladatropic rearrangement/β-H elimination or allene via β-H elimination from C_sp²,¹³ and the π-allylpalladium intermediate trapped by the byproduct sulfinic acid salt.¹⁴ Given the significance of chiral pyrrolidines and piperidines as core structures in alkaloid natural products, the development of an asymmetric version of this elegant carbenylative amination reaction is highly desirable. In recent years, our group has developed a series of chiral sulfinamide phosphine ligands (so-called Sadphos), which showed unique potential in asymmetric transition-metal catalysis,^6f,15 so we wondered whether Sadphos could address this challenging asymmetric carbenylative amination reaction (Scheme 1C).Initially, our study began with (E)-vinyl iodide 1a and N-tosylhydrazone 2a in the presence of Pd₂(dba)₃, t-BuOLi, Et₃N, and triethylbenzylammonium chloride (TEBAC) in THF at 30 °C. A series of commercially available chiral ligands were first screened (Fig. 2). Only (R, R)-DIOP (L1), (R)-DTBM-SegPhos (L3) and (R)-MOP (L4) provided the desired product 3aa with poor enantioselectivity and other ligands such as (R, R)-Ph-BPE (L2), (R, S)-Josiphos (L5) and (S, S)-ⁱPr-FOXAP (L6) showed low reactivity. We next turned to systematically investigate Sadphos, such as Wei-Phos,¹⁶ Xiao-Phos,^15d,17 Ming-Phos,^15a,18 Xu-Phos,^15b,19 Xiang-Phos²⁰ and PC-Phos^15c,21 (Fig. 2). To our delight, PC1 delivered 3aa in 32% yield and 85.5 : 14.5 er. Inspired by this result, we further screened PC2–PC5 which vary in the substituent of phenyl, but unfortunately none of them showed better results. Surprisingly, the reactivity of this reaction could be greatly improved with our recently developed GF-Phos GF1, delivering 71% yield. When steric hindered tert-butyl groups were introduced on the phenyl group (GF2), the product 3aa was obtained in 77% yield with 91.5 : 8.5 er. After screening different palladium catalysts and solvents (Open in a separate windowFig. 2Screened chiral ligands.Optimization of reaction conditions^a

0D–1D hybrid nanoarchitectonics: tailored design of FeCo@N–C yolk–shell nanoreactors with dual sites for excellent Fenton-like catalysis

Heterogeneous Fenton-like processes are very promising methods of treating organic pollutants through the generation of reactive oxygen containing radicals. Herein, we report novel 0D–1D hybrid nanoarchitectonics (necklace-like structures) consisting of FeCo@N–C yolk–shell nanoreactors as advanced catalysts for Fenton-like reactions. Each FeCo@N–C unit possesses a yolk–shell structure like a nanoreactor, which can accelerate the diffusion of reactive oxygen species and guard the active sites of FeCo. Furthermore, all the nanoreactors are threaded along carbon fibers, providing a highway for electron transport. FeCo@N–C nano-necklaces thereby exhibit excellent performance for pollutant removal via activation of peroxymonosulfate, achieving 100% bisphenol A (k = 0.8308 min⁻¹) degradation in 10 min with good cycling stability. The experiments and density-functional theory calculations reveal that FeCo dual sites are beneficial for activation of O–O, which is crucial for enhancing Fenton-like processes.

Novel 0D–1D hybrid nanoarchitectonics consisting of FeCo@N–C yolk–shell nanoreactors are developed for Fenton-like reaction. With the multilevel advantages of this design, FeCo@N–C nano-necklaces exhibit excellent performance for BPA removal.

Advanced oxidation processes (AOPs) are one of the most promising strategies to eliminate organic contaminants, sustainably generating reactive oxygen species (ROS) to ideally destroy all non-biodegradable, recalcitrant, toxic, or membrane-permeable organic impurities.^1–4 Among these AOPs, sulfate radical (SO₄˙⁻)-based Fenton-like processes have gained increasing attention as a water treatment strategy because of the strong oxidation potential of SO₄˙⁻ (3.1 V vs. normal hydrogen electrode) at wider pH ranges. SO₄˙⁻ is mainly produced by physical or chemical methods for activation of persulfate salts, such as peroxymonosulfate (PMS) and persulfate.^5–9 Over the past two decades, heterogeneous catalysis has emerged as the most effective approach to water treatment, with much effort dedicated to developing better catalysts, including transition metal-based and carbonaceous materials.^10,11 Unfortunately, most metal-based catalysts suffer from leaching of toxic metal ions, which can thwart their practical application,^12,13 and although carbonaceous catalysts produce no secondary pollution, their cycle performance is always depressed.¹⁴ There is therefore an urgent need to find robust catalysts with adequate activity and stability for Fenton-like processes.To achieve superior performance, an ideal Fenton-like catalyst should contain oxidants with favorably reactive centers for cleavage of peroxyl bonds (O–O), have structure optimized for target pollutant attraction, and have chainmail to protect the vulnerable active sites for long periods.^15–17 Recent studies have demonstrated Co–N–C active sites prefer to activate the O–O of PMS.¹⁸ Furthermore, introducing Fe-doping into the Co–N–C system not only suppresses Co²⁺ leaching, but also modulates the pyrrolic-N content, which is the adsorption site for capture of bisphenol A (BPA).¹⁹ We previously discovered that Co@C yolk–shell nanoreactors could enhance the catalytic activity because of the confinement effect in the nano-spaces between the core and shell, while the carbon shell acted like a chainmail protecting the Co active sites, keeping them highly reactive after five cycles.^20,21Combining different kinds of materials to generate novel hybrid material interfaces can enable the creation of new kinds of chemical and physical functionalities that do not currently exist. However, one cannot simply mix these materials in an uncontrolled manner, because the ensemble of interfaces created by random mixing tends to favour thermodynamically stable interfaces that are functionally less active. Therefore, to prepare new materials with high functionality, it is necessary to carefully control the hybridization of components in interfacial regions with nanometric or atomic precision. By further hybridization of different components e.g., zero to one dimension (0D–1D) hybrid structures, we can prepare the structure to increase not only the specific surface area but also the interfacial region between different materials.In this work, we report novel 0D–1D hybrid nanoarchitectonics (necklace-like structures) consisting of FeCo@N–C yolk–shell nanoreactors as a PMS activator for Fenton-like processes. This catalyst has multilevel advantages: (i) each FeCo@N–C unit is a well-formed yolk–shell nanoreactor, which can guarantee sufficient contact of reactants and active sites, as well as defend them for good durability; (ii) all single nanoreactors are threaded along the carbon fibers, providing a highway for electron transport; and (iii) all the carbon fibers constructed into a thin film with macroscopic structure, which overcomes the complex recyclability of powder catalysts. Benefiting from favorable composition and unique structure, the FeCo@N–C catalyst delivers excellent performance for BPA removal via activation of PMS accompanied with good stability.The synthesis processes of necklace-like nanoarchitecture containing FeCo@N–C yolk–shell nanoreactors are illustrated in Fig. 1a. First, uniform Fe–Co Prussian blue analogue (Fe–Co PBA) nanocubes with an average size of 800–900 nm (Fig. 1b) are encapsulated in polyacrylonitrile (PAN) nanofibers by electrospinning. The obtained necklace-like FeCo PBA–PAN fibers (Fig. 1c) are then pyrolyzed at 800 °C in N₂ atmosphere to produce FeCo@N–C nano-necklaces. The scanning electron microscopy (SEM) image (Fig. 1d) of the FeCo@N–C shows this necklace-like morphology with its large aspect ratio, with the FeCo@N–C particles strung along the PAN-derived carbon fibers. A broken particle (Fig. 1e) shows that the FeCo@N–C has a yolk–shell architecture, which is also identified by transmission electron microscopy (TEM). Fig. 1f and g show the well-defined space between the inner yolk and outer shell, which is attributed to the volume shrinkage of the original Fe–Co PBAs. During pyrolysis, Fe–Co PBA is reduced to FeCo (inner yolk) and PAN is carbonized (outer carbon shell), resulting in the unique necklace-like nanoarchitecture.^22–24 The high-resolution TEM in Fig. 1h shows a lattice fringe of 0.20 nm, which matches well with the (110) plane of FeCo alloy.²⁵ The scanning transmission electron microscopy (STEM) image (Fig. 1i) and corresponding elemental map (Fig. 1j) indicate that FeCo nanocrystals are well distributed in the inner core with some small FeCo nanocrystals located on external carbon shells. Furthermore, the control samples of Fe@N–C and Co@N–C nano-necklaces, prepared by only replacing the Fe–Co PBA nanocubes with Fe–Fe PB and Co–Co PBA (Fig. S1^†), also demonstrate the versatility of this synthetic strategy. The formation of hierarchical porous structure, beneficial to the PMS transportation on the surface of catalysts, could be determined by N₂ adsorption–desorption isotherms and corresponding pore volume analysis (Fig. S2 and Table S1^†).Open in a separate windowFig. 1(a) Preparation of FeCo@N–C necklace-like nanoarchitecture. SEM images of (b) Fe–Co PBA cubic particles and (c) the electrospun FeCo PBA–PAN fibers. (d and e) SEM, (f and g) TEM, and (h) high-resolution TEM images of FeCo@N–C nano-necklaces. (i) STEM and (j) the corresponding elemental mappings of C, N, Fe, and Co.The X-ray diffraction patterns of the as-prepared products are depicted in Fig. S3,^† with one prominent diffraction peak centered at 44.8° corresponding to the (110) lattice plane of FeCo alloy. All the products also have a characteristic signal at 26°, implying that graphite carbon is formed during pyrolysis. Raman spectroscopy further analyzed the crystal structures and defects of the FeCo@N–C nano-necklaces (Fig. S4^†), where peaks found at 1349 cm⁻¹ and 1585 cm⁻¹ index the disordered (D band) and graphitic carbon (G band), respectively.²⁶ X-ray photoelectron spectroscopy investigated the composition and valence band spectra of FeCo@N–C nano-necklaces. The survey spectrum (Fig. S5a^†) reveals the presence of Fe (1.4%), Co (1.2%), C (86.4%), N (4.5%), and O (6.5%) in the composite. The high-resolution N 1s spectrum (Fig. S5b^†) exhibits broad peaks at 398.1, 401.1, and 407.4 eV, corresponding to the pyridinic-N, graphitic-N, and σ* excitation of C–N, respectively.²⁷ The high-resolution Fe 2p spectrum (Fig. S5c^†) shows a broad peak at 707.4 eV, attributed to Fe⁰. Similarly, the 777.5 eV peak observed in the Co 2p spectrum (Fig. S5d^†) corresponds to Co⁰, implying that FeCo dual sites have formed.²⁸ The oxidation state of these sites was investigated by ⁵⁷Fe Mössbauer spectroscopy, which found a sextet in the Mössbauer spectrum of the FeCo@N–C nano-necklaces attributed to FeCo dual sites (Fig. 2a and Table S2^†).²⁹ The coordination environment of the FeCo dual sites was also verified by X-ray absorption fine structure (XAFS) spectroscopy. Fig. 2b shows that the X-ray absorption near-edge structure (XANES) spectra of the Fe K-edge, which demonstrates a similar near-edge structure to that of Fe foil, illustrating that the main valence state of Fe in FeCo@N–C nano-necklaces is Fe⁰. Furthermore, the extended-XAFS (EXAFS) spectra (Fig. 2c) displays a peak at 1.7 Å, which is ascribed to the Fe–N bond, and a remarkable peak at approximately 2.25 Å corresponding to the metal–metal band.^10,30 The Co K-edge and EXAFS spectra (Fig. S6^†) also confirm the presence of Co–N and the metal–metal band. These results provide a potential structure of the FeCo dual sites in the FeCo@N–C nano-necklaces, as illustrated in Fig. 2d.Open in a separate windowFig. 2(a) ⁵⁷Fe Mössbauer spectra of FeCo@N–C nano-necklaces at 298 K. (b) Fe K-edge XANES spectra of FeCo@N–C nano-necklaces and Fe foil. (c) Corresponding Fourier transformed k³-weighted of the EXAFS spectra for Fe K-edge. (d) Possible structure of the FeCo dual sites.This dual-metal center and necklace-like structure may be beneficial to enhance catalytic performance. Fig. 3a shows the Fenton-like performance for BPA degradation compared to Fe@N–C nano-necklaces, Co@N–C nano-necklaces, and FeCo@N–C particles (Fe–Co PBA directly carbonized without electrospinning). Here, the FeCo@N–C nano-necklaces display a higher catalytic performance, with BPA completely removed in 7 min. To clearly compare their catalytic behavior, the kinetics of BPA degradation was fitted by the first-order reaction. As shown in Fig. 3b, FeCo@N–C nano-necklaces exhibit the highest apparent rate constant (k = 0.83 min⁻¹), which is approximately 6.4, 2.6, and 1.2 times that of FeCo@N–C particles, Fe@N–C nano-necklaces, and Co@N–C nano-necklaces, respectively. The significantly enhanced performance of FeCo@N–C nano-necklaces suggests that the FeCo dual sites and necklace-like nanoarchitecture are crucial. Furthermore, the concentration of BPA and PMS in the solution is higher than that in yolk–shell nanoreactor, resulting a concentration gradient which helps to accelerate the diffusion rates of reactants (Fig. 3c).^31,32 For these nano-necklaces, the carbon shell acts like a chainmail protecting the FeCo active sites from attack by molecules and ions, and all the nanoreactors are threaded along the carbon fibers, providing a highway for electron transport, which is important for SO₄˙⁻ generation (SO₄˙⁻ production as eqn, HSO₅⁻ + e⁻ → SO₄˙⁻ + OH⁻). Electrochemical impedance spectroscopy further confirms the good conductivity of the FeCo@N–C nano-necklaces (Fig. 3d). In addition, the concentration of metal-ion leaching and cycling performance (Fig. 3e and f) reveal the high reusability of FeCo@N–C nano-necklaces, with 95% BPA removal in 20 min after five cycles, which is also proved by the SEM and TEM characterization (Fig. S7^†). The effect of other reaction parameters on the BPA degradation, such as pH, reaction temperature, PMS or catalysts dosage, and common anions, were investigated in detail (Fig. S8–S11^†). All the results demonstrate that FeCo@N–C nano-necklaces deliver a better performance for PMS catalysis. In addition, the turnover frequency (TOF) value of FeCo@N–C nano-necklaces is 5.5 min⁻¹ for BPA degradation, which is higher than many previously reported catalysts (detailed catalytic performance comparison as shown in Table S3^†).Open in a separate windowFig. 3(a) BPA degradation efficiency in different reaction systems and (b) the corresponding reaction rate constants. (c) Schematic illustration of PMS activation in FeCo@N–C nano-necklaces. (d) Nyquist plots of the catalysts. (e) The metal leaching in different reaction systems. (f) Cycling performance of FeCo@N–C nano-necklaces for BPA removal. Reaction conditions: [catalyst] = 0.15 g L⁻¹, [BPA] = 20 mg L⁻¹, [PMS] = 0.5 g L⁻¹, T = 298 K, and initial pH = 7.0.To examine the enhanced catalytic activity, radical quenching experiments were conducted. As shown in Fig. 4a, when NaN₃ is added to the reaction solution as a scavenger for ¹O₂, there is no significant reduction of BPA decomposition, implying that non-radicals are not the dominant reactive species. By comparison, when tert-butanol (TBA) (radical scavenger for ˙OH) is added, there is a slight (2.8%) decrease in BPA removal. However, if methanol (radical scavenger for SO₄˙⁻ and ˙OH) is added, the efficiency of BPA degradation declines by up to 59.2%, indicating that the major radicals generated from the PMS activation are SO₄˙⁻;³³ the presence of these radicals is also verified by electron paramagnetic resonance (EPR) (Fig. 4b). Furthermore, the significant inhibition ratio can be observed when KI (quencher for the surface) is added, demonstrating that BPA degradation is mainly attributed to reactions with SO₄˙⁻, which is produced by a surface catalytic process.³⁴Open in a separate windowFig. 4(a) Effects of the radical scavengers on BPA degradation. (b) EPR spectra of SO₄˙⁻ and ˙OH. (c) The energy profiles of PMS on FeCo@N–C nano-necklaces surface. (d) Optimized configurations of PMS adsorbed on FeCo@N–C nano-necklaces.Density-functional theory was applied to calculate the surface energy of PMS activation at FeCo dual sites (Fig. 4c, d and S12^†). The dissociation barrier of PMS into SO₄˙⁻ and OH⁻ is −2.25 eV, which is much lower than that on an Fe or Co single site, suggesting that cleavage of O–O bonds of PMS occurs more easily on FeCo dual sites. This is because FeCo dual sites provide two anchoring sites for the dissociated O atoms, leading to more efficient activation of O–O. The FeCo@N–C nano-necklaces can reduce the energy barrier of O–O bond breaking, which results in high activity for PMS activation and thus high productivity of SO₄˙⁻.     

Construction of vicinal 4°/3°-carbons via reductive Cope rearrangement

Herein reported is a strategy for constructing vicinal 4°/3° carbons via reductive Cope rearrangement. Substrates have been designed which exhibit Cope rearrangement kinetic barriers of ∼23 kcal mol⁻¹ with isoenergetic favorability (ΔG ∼ 0). These fluxional/shape-shifting molecules can be driven forward by chemoselective reduction to useful polyfunctionalized building blocks.

Herein reported is a strategy for constructing vicinal 4°/3° carbons via reductive Cope rearrangement.

Constructing sterically congested vicinal quaternary–tertiary carbons (4°/3° carbons) via Cope rearrangement is currently quite limited with only a handful of papers on the subject published over the past 40 years. This stands in stark contrast to the plethora of other methods for establishing sterically congested vicinal carbons.^1–5 Central to the challenge are kinetic and thermodynamic issues associated with the transformation. In the simplest sense, Cope rearrangements proceed in the direction that results in highest alkene substitution (Fig. 1).^6,7 To forge 4°/3° motifs by Cope rearrangement, additional driving forces must be introduced to reverse the [3,3] directionality and compensate for the energetic penalty associated with the steric and torsional strain of the targeted vicinal 4°/3° motif. With limited reports in all cases, oxy-Cope substrates (Scheme 1, eqn (1)),^8–14 divinylcyclopropanes (Scheme 1, eqn (2)),^15–20 and vinylidenecyclopropane-based 1,5-dienes²¹ (Scheme 1, eqn (3)) have demonstrated favourability for constructing vicinal 4°/3° carbons. Malachowski et al. put forth a series of studies on the construction of quaternary centers via Cope rearrangement driven forward by a conjugation event (Scheme 1, eqn (4)).^22–25 In their work, a single example related to the construction of vicinal 4°/3° centers was disclosed, though kinetic (180 °C) and thermodynamic (equilibrium mixtures) challenges are also observed.²³ And of particular relevance to this work, Wigfield et al. demonstrated that 3,3-dicyano-1,5-dienes with the potential to generate vicinal 4°/3° carbons instead react via an ionic mechanism yielding the less congested products (Scheme 1, eqn (5)).²⁶Open in a separate windowFig. 1Cope equilibrium of 1,1,6-trisubstituted 1,5-dienes.Open in a separate windowScheme 1(A) Cope rearrangements for constructing vicinal 4°/3°-centers (B) this report.Our group has been examining strategies to decrease kinetic barriers and increase the thermodynamic favourability of 3,3-dicyano-1,5-diene-based Cope substrates.^27–31 Beyond the simplest, unsubstituted variants, this class of 1,5-diene is not particularly reactive in both a kinetic and thermodynamic sense (e.g.Scheme 1, eqn (5)).^26,32 Reactivity issues aside, these substrates are attractive building blocks for two main reasons: (1) they have straightforward accessibility from alkylidenemalononitriles and allylic electrophiles by deconjugative allylic alkylation.³³ (2) The 1,5-diene termini are substantially different (malononitrile vs. simple alkene) thus allowing for orthogonal functional group interconversion facilitating target and analogue synthesis.³⁴ Herein we report that a combination of 1,5-diene structural engineering^28,31 and reductive conditions (the reductive Cope rearrangement^29,30) can result in the synthesis of building blocks containing vicinal gem-dimethyl 4°/3° carbons along with orthogonal malononitrile and styrene functional groups for interconversion (Scheme 1B). On this line, malononitrile can be directly converted to amides³⁴ yielding functionally dense β-gem-dimethylamides, important pharmaceutical scaffolds.³⁵This project began during the Covid-19 pandemic lockdown (ca. March–May 2020). As such, we were not permitted to use our laboratory out of an abundance of caution. We took this opportunity to first computationally investigate a Cope rearrangement that could result in vicinal 4°/3° carbons (Scheme 2). Then, when permitted to safely return to the lab, we would experimentally validate our findings (vide infra). From our previous work, it is known that by adding either a 4-aromatic group²⁸ or a 4-methyl group³¹ to a 3,3-dicyano-1,5-diene, low barrier (rt – 80 °C) diastereoselective Cope rearrangements can occur. Notably, the 4-substituent was found to destabilize the starting material (weaken the C3–C4 bond, conformationally bias the substrate for [3,3]), and stabilize the product side of the equilibrium via resonance (phenyl group) or hyperconjugation (methyl group). In this study, we modelled substrates 1, 3, and 5 that have variable 4-substitution and would result in vicinal gem-dimethyl- and phenyl-containing 4°/3° carbons upon Cope rearrangement to 2, 4, or 6, respectively. We chose to target this motif due to likely synthetic accessibility from simple starting materials but also because of the important and profound impact that gem-dimethyl groups impart on pharmaceuticals.³⁵ Substrate 1 lacking 4-substitution had an extremely unfavourable kinetic and thermodynamic profile (ΔG^‡ = 31.6; ΔG = +5.3 kcal mol⁻¹). When a 4-methyl group was added, the kinetic barrier (ΔG^‡) dropped appreciably to 28.2 kcal mol; however, the thermodynamics were still quite endergonic (ΔG = +4.4 kcal mol⁻¹). Most excitingly, it was uncovered that the 4-phenyl group dramatically impacted the kinetics and thermodynamics: the [3,3] has a barrier of 22.9 kcal mol⁻¹ (ΔG^‡) and is ∼isoenergetic (ΔG = +0.17 kcal mol⁻¹). Thus, the reaction appears to be fluxional/shape-shifting at room temperature.^36–40 For this substrate, we also modelled the dissociative pathway (Scheme 2D). It was found that bond breakage to two allylic radical intermediates is a higher energy process than the concerted transition state (Scheme 2Cvs.Scheme 2D). Specifically, the dissociative pathway was found to be kinetically less favourable (ΔG^‡ ∼ 27.6 kcal mol; ΔG = 26.2 kcal mol⁻¹) than the concerted process (ΔG^‡ = 22.9 kcal mol⁻¹). While the dissociative pathway is less favourable than the concerted transformation, we surmised that the two-step process becomes accessible at elevated temperature (vide infra). Finally, the ionic pathway was calculated to be significantly higher for this substrate (see the ESI^†).Open in a separate windowScheme 2Computational analysis of 3,3-dicyano-1,5-diene that in theory could result in vicinal 4°/3° carbons. (A) 4-Unsubstituted 3,3-dicyano-1,5-diene. (B) 4-Methyl 3,3-dicyano-1,5-diene. (C) 4-Phenyl 3,3-dicyano-1,5-diene. (D) The dissociative mechanism for substrate 5 is higher than the closed transition state. (E) visualization of the kinetic- and thermodynamic differences of transformations (A–D).The class of substrate uncovered from our computational investigation could be accessed from γ,γ-dimethyl-alkylidenemalononitrile (7a) and 1,3-diarylallyl electrophiles (such as 8a) by Pd-catalyzed deconjugative allylic alkylation (Scheme 3A).³³ As such, model 1,5-diene 5a was prepared to verify the computational results. It was found that upon synthesis of 5a, an inseparable 21 : 79 mixture of 1,5-diene 5a and the 1,5-diene 6a was observed. The predicted ratio of 5a to 6a was 57 : 43 (Scheme 2C). These two results are within the error of the calculations (predicted; slightly endergonic, observed; slightly exergonic). To determine whether the transformation was progressing through the predicted concerted pathway (Scheme 2C) over the dissociative pathway (Scheme 2D), substrate 5b was prepared by an analogous deconjugative allylic alkylation reaction. Similarly, two Cope equilibrium isomers 5b and 6b are observed at room temperature in a 12 : 88 ratio. Upon heating at 100 °C for 3 h, the 1,5-dienes “scramble” (e.g. iso-6b is observed; 0.2 : 1.0 : 1.5 ratio of 5b : 6b : iso-6b) indicating that the dissociative pathway is only accessible at elevated temperature. This is all in good agreement with the calculated kinetics and thermodynamics of this system (Scheme 2).Open in a separate windowScheme 3(A) Observation of fluxional [3,3] and confirmation of calculated predictions. (B) Optimization of a reductive Cope rearrangement protocol for constructing vicinal 4°/3° centers. (C) The Pd-catalyzed deconjugative allylic alkylation must be regioselective.With respect to the synthetic methodology, we aimed to increase the overall efficiency and applicability of the sequence (Scheme 3B). Specifically, we wanted to avoid [3,3] equilibrium mixtures and sensitive/unstable substates and intermediates. It was found that the direct coupling of 7a with diphenylallyl alcohol 9a could take place in the presence of DMAP, Ac₂O, and Pd(PPh₃)₄. When the coupling was complete, methanol and NaBH₄ were added to drive the Cope equilibrium forward, yielding the reduced Cope rearrangement product 10a in 76% isolated yield. In terms of practicality and efficiency, this method utilizes diphenylallyl alcohols, which are more stable and synthetically accessible than their respective acetates, and the [3,3] equilibrium mixture can be directly converted dynamically to a single reduced product.With an efficient protocol in hand for constructing malononitrile–styrene-tethered building blocks featuring central vicinal 4°/3° carbons, we next examined the scope of the transformation (Scheme 4). We chose diarylallyl alcohols with the propensity to react regioselectively via an electronic bias (Scheme 3C).^41,42 The combination of p-nitrophenyl and phenyl (10b) or p-methoxyphenyl (10c) yielded regioselective outcomes with the electron-deficient arene at the allylic position. This is consistent with the expected regiochemical outcome where the nucleophile reacts preferentially at the α-position and the electrophile reacts at the allylic position bearing the donor-arene (Scheme 3C).^41,42 Then, reductive Cope rearrangement occurs to position the electron-deficient arene adjacent to the gem-dimethyl quaternary center. This is an exciting outcome as many pharmaceutically relevant (hetero)arenes are electron deficient. Thus, fluorinated arenes were installed at the allylic position of products 10d–10k. While the phenyl group resulted in poor regioselectivity (1 : 1–3 : 1), the p-methoxyphenyl group enhanced the regiomeric ratios in all cases (3 : 1–15 : 1). The degree of selectivity is correlated with the number and position of fluorine atoms. N-Heterocycles could be incorporated with excellent regioselectivity, generally speaking (10l–10q). For example, 3-chloro-4-pyridyl (10l/10m) groups were installed at the allylic position with >20 : 1 rr. 4-Chloro-3-pyridyl was poorly regioselective (10n), but the combination of 4-trifluomethyl-3-pyridyl/p-methoxyphenyl (10o) gave good regioselectivity of 11 : 1. 2-Pyridyl/p-methoxyphenyl (10q) was also a regioselective combination. We also examined a few other heterocycles including quinoline (10s) and thiazole (10t and 10u) with excellent and modest regioselectivity observed, respectively. As a general trend, when the arenes on the allylic electrophile become less polarized, poor regioselectivity is observed in the Pd-catalyzed allylic alkylation. For example, the combination of p-chlorophenyl and p-methoxyphenyl (10v) or phenyl (10w) yields regioisomeric mixtures of products. This can be circumvented by utilizing symmetric electrophiles (to 10x).Open in a separate windowScheme 4Scope of the 4°/3°-center-generating reductive Cope rearrangement.The phenyl or the p-methoxyphenyl group is necessary to achieve the 4°/3° carbon-generating Cope rearrangement: it functions as an “activator” by lowering the kinetic barrier and increasing thermodynamic favourability. These activating groups can be removed through alkene C C cleavage reactions (e.g. metathesis (Scheme 5) and ozonolysis (Scheme 6B)). In this regard, highly substituted cycloheptenes 11 were prepared by allylation and metathesis (Scheme 4).^28,43 The yields were modest to excellent over this two-step sequence. In many cases, where 10 exists as a mixture of regioisomers, the major allylation/RCM products 11 could be chromatographically separated from their minor constituents. As shown in Scheme 6A, the malononitrile can be transformed via oxidative amidation³⁴ to products 12 containing a dense array of pharmaceutically relevant functionalities (amides, gem-dimethyl, fluoroaromatics, and heteroaromatics). Following this transformation, ozonolysis terminated with a NaBH₄ quench installs an alcohol moiety on small molecule 13a.Open in a separate windowScheme 5Removal of the “activating group” by ring-closing metathesis.Open in a separate windowScheme 6(A) oxidative amidation of malononitrile. (B) Removal of “activating group” by ozonolysis.These first computational and experimental studies utilizing 3,3-dicyano-1,5-dienes as substrates for constructing vicinal 4°/3° centers sets the stage for much further examination and application. For example, while we focused our efforts on gem-dimethyl-based quaternary carbons, it is likely that other functionality can be installed at this position. For example, while unoptimized, it appears the protocol is reasonably effective at incorporating a piperidine moiety in addition to heteroarenes from the allylic electrophile (7b + 9f → 14a; Scheme 7A). Similar functional group interconversion chemistry as described in Schemes 5 and ​and66 can thus yield functionally dense building blocks 15 and 16 in good yields.Open in a separate windowScheme 7(A) The construction of 4/3° centres on piperidines. (B) Promoting endergonic [3,3] rearrangements is possible, assuming the [3,3] kinetic barrier is sufficiently low.While the 4,6-diaryl-3,3-dicyano-1,5-dienes offered the most attractive energetic profile (low kinetic barrier, isoenergetic [3,3] equillibrium; Scheme 2C), the 4-methyl analogue is also intriguing to consider as a viable substrate class for reductive Cope rearrangement (Scheme 2B). The challenge here is that the kinetics and thermodynamics are quite unfavourable (not observable by NMR), but potentially not prohibitively so. It is extremely exciting to find that Cope equilibria that are significantly endergonic in the desired, forward direction (e.g.3a to 4a) can be promoted by a related reductive protocol (Scheme 7B). While unoptimized, we were able to isolate product 17 in xx% yield by heating at 90 °C in the presence of Hantzsch ester in DMF.     

Simplifying and expanding the scope of boron imidazolate framework (BIF) synthesis using mechanochemistry

Mechanochemistry enables rapid access to boron imidazolate frameworks (BIFs), including ultralight materials based on Li and Cu(i) nodes, as well as new, previously unexplored systems based on Ag(i) nodes. Compared to solution methods, mechanochemistry is faster, provides materials with improved porosity, and replaces harsh reactants (e.g. n-butylithium) with simpler and safer oxides, carbonates or hydroxides. Periodic density-functional theory (DFT) calculations on polymorphic pairs of BIFs based on Li⁺, Cu⁺ and Ag⁺ nodes reveals that heavy-atom nodes increase the stability of the open SOD-framework relative to the non-porous dia-polymorph.

Mechanochemistry enables rapid access to boron imidazolate frameworks (BIFs), including ultralight materials based on Li and Cu(i) nodes, as well as new, previously unexplored systems based on Ag(i) nodes.

Mechanochemistry^1–7 has emerged as a versatile methodology for the synthesis and discovery of advanced materials, including nanoparticle systems^8–10 and metal–organic frameworks (MOFs),^11–15 giving rise to materials that are challenging to obtain using conventional solution-based techniques.^16–18 Mechanochemical techniques such as ball milling, twin screw extrusion¹⁹ and acoustic mixing^20,21 have simplified and advanced the synthesis of a wide range of MOFs, permitting the use of simple starting materials such as metal oxides, hydroxides or carbonates,^22,23 at room temperature and without bulk solvents, yielding products of comparable stability and, after activation, higher surface areas than solution-generated counterparts.^24–29 The efficiency of mechanochemistry in MOF synthesis was recently highlighted by accessing zeolitic imidazolate frameworks (ZIFs)^30,31 that were theoretically predicted, but not accessible under conventional solution-based conditions.¹⁷The advantages of mechanochemistry in MOF chemistry led us to address the possibility of synthesizing boron imidazolate frameworks (BIFs),^32–34 an intriguing but poorly developed class of microporous materials analogous to ZIFs, comprising equimolar combinations of tetrahedrally coordinated boron(iii) and monovalent Li⁺ or Cu⁺ cations as nodes (Fig. 1A–C). Although BIFs offer an attractive opportunity to access microporous MOFs with lower molecular weights, particularly in the case of “ultralight” systems based on Li⁺ and B(iii) centers, this family of materials has remained largely unexplored – potentially due to the need for harsh synthetic conditions, including the use of n-butyllithium in a solvothermal environment.^32–34Open in a separate windowFig. 1Structures of previously reported BIFs with: (A) zni-, (B) dia-, or (C) SOD-topology (M = Li, Cu); (D) tetrakis(imidazolyl)boric acids used herein for mechanochemical BIF synthesis; and (E) schematic representation of the herein developed mechanosynthesis of dia- and SOD BIF polymorphs based on Li, Cu or Ag metal nodes.We now show how switching to the mechanochemical environment enables lithium- and copper(i)-based BIFs to be prepared rapidly (i.e., within 60–90 minutes), without elevated temperatures or bulk solvents, and from readily accessible solid reactants, such as hydroxides and oxides (Fig. 1D and E). While the mechanochemically-prepared BIFs exhibit significantly higher surface areas than the solvothermally-prepared counterparts, mechanochemistry allows for expanding this class of materials towards previously not reported Ag⁺ nodes. The introduction of BIFs isostructural with those based on Li⁺ or Cu⁺ but comprising of Ag⁺ ions, enables a periodic density-functional theory (DFT) evaluation of their stability. This reveals that switching to heavier elements as tetrahedral nodes improves the stability of sodalite topology (SOD) open BIFs with respect to close-packed diamondoid (dia) topology polymorphs.As a first attempt at mechanochemically synthesis of BIFs, we targeted the synthesis of previously reported zni-topology LiB(Im)₄ and CuB(Im)₄ frameworks (Li-BIF-1 and Cu-BIF-1, respectively, Fig. 1A) using a salt exchange reaction between LiCl or CuCl with commercially available sodium tetrakis(imidazolyl)borate (Na[B(Im)₄]) (Fig. 2A). Milling of LiCl and Na[B(Im)₄] in a 1 : 1 stoichiometric ratio for up to 60 minutes led to the appearance of Bragg reflections consistent with the target Li-BIF-1 (CSD MOXJEP) and the anticipated NaCl byproduct. The reaction was, however, incomplete, as seen by X-ray reflections of Na[B(Im)₄] starting material. In order to improve reactant conversion, we explored liquid-assisted grinding (LAG), i.e. milling in the presence of a small amount of a liquid phase (measured by the liquid-to-solid ratio η³⁵ in the range of ca. 0–2 μL mg⁻¹). Using LAG conditions with acetonitrile (MeCN, 120 μL, η = 0.5 μL mg⁻¹) led to the complete disappearance of reactant X-ray reflections, concomitant with the formation of Li-BIF-1 alongside NaCl within 60 minutes.Open in a separate windowFig. 2(A) Reaction scheme for the mechanochemical synthesis of Li-BIF-1 by a salt metathesis strategy. Selected PXRD patterns for: (B) Na[B(Im)₄] (C) LiCl, (D) simulated Li-BIF-1 (CSD MOXJPEP) and (E) synthesized BIF-1-Li by LAG for 60 minutes with MeCN (η = 0.5 μL mg⁻¹), (F) CuCl, (G) simulated Cu-BIF-1 (CSD MOXJIT), and (H) synthesized BIF-1-Cu by LAG for 60 minutes with MeOH (η = 0.50 μL mg⁻¹). Asterisks denote NaCl, a byproduct of the metathesis reaction. (Fig. 2B–E, also see ESI^†). The copper-based zni-CuB(Im)₄ (Cu-BIF-1) was readily obtained from CuCl within 60 minutes using similar LAG conditions. We also explored LAG with methanol (MeOH), revealing that the exchange reaction to form NaCl took place with both LiCl and CuCl starting materials. With LiCl, however, the PXRD pattern of the product could not be matched to known phases involving Li⁺ and B(Im)₄⁻ (see ESI^†). With CuCl as a reactant, LAG with MeOH (η = 0.5 μL mg⁻¹) cleanly produced Cu-BIF-1 alongside NaCl (see ESI^†).Next, we explored an alternative synthesis approach, analogous to that previously used to form ZIFs and other MOFs: an acid–base reaction between a metal oxide or hydroxide and the acid form of the linker: tetrakis(imidazolato)boric acid, HB(Im)₄ (Fig. 3A).^36–40 Neat milling LiOH with one equivalent of HB(Im)₄ in a stainless steel milling assembly led to the partial formation of Li-BIF-1, as evidenced by PXRD analysis (see ESI^†). Complete conversion of reactants into Li-BIF-1 was achieved in 60 minutes by LAG with MeCN (η = 0.25 μL mg⁻¹), as indicated by PXRD analysis (Fig. 3B–E), Fourier transform infrared attenuated total reflectance spectroscopy (FTIR-ATR), thermogravimetric analysis (TGA) in air, and analysis of metal content by inductively-coupled plasma mass spectrometry (ICP-MS) (see ESI^†).Open in a separate windowFig. 3(A) Reaction scheme for the mechanochemical synthesis of Li-BIF-1 using the acid–base strategy. Selected PXRD patterns for: (B) H[B(Im)₄] (C) LiOH, (D) simulated Li-BIF-1 (CSD MOXJPEP), (E) synthesized BIF-1-Li by LAG for 60 minutes with MeCN (η = 0.25 μL mg⁻¹), (F) Cu₂O, (G) simulated Cu-BIF-1 (CSD MOXJIT), and (H) synthesized Cu-BIF-1 by ILAG for 60 minutes with MeOH (η = 0.50 μL mg⁻¹) and NH₄NO₃ additive (5% by weight).Neat milling of HB(Im)₄ with Cu₂O under similar conditions gave a largely non-crystalline material, as evidenced by PXRD (see ESI^†). Switching to the ion- and liquid-assisted grinding (ILAG) methodology, in which the reactivity of a metal oxide is enhanced by a small amount of a weakly acidic ammonium salt, and which was introduced to prepare zinc and cadmium ZIFs from respective oxides,^37–40 enabled the synthesis of Cu-BIF-1 from Cu₂O. Specifically, PXRD analysis revealed complete disappearance of the oxide in samples obtained by ILAG with either MeOH or MeCN (η = 0.5 μL mg⁻¹) in the presence of NH₄NO₃ additive (5% by weight, see ESI^†). Notably, achieving complete disappearance of Cu₂O reactant signals also required switching from stainless steel to a zirconia-based milling assembly, presumably due to more efficient energy delivery.⁴¹ After washing with MeOH, the material was characterized by FTIR-ATR, TGA in air, and analysis of metal content by ICP-MS (see ESI^†).Whereas both the metathesis and acid–base approaches can be used to mechanochemically generate Li- and Cu-BIF-1, the latter approach has a clear advantage of circumventing the formation of the NaCl byproduct. Consequently, in order to further the development of mechanochemical routes to other BIFs, we focused on the acid–base strategy. As next targets, we turned to MOFs based on tetrakis(2-methylimidazole)boric acid H[B(Meim)₄],³⁶ previously reported³² to adopt either a non-porous diamondoid (dia) topology (BIF-2) or a microporous sodalite (SOD) topology (BIF-3) with either Li⁺ or Cu⁺ as nodes (Fig. 4). Attempts to selectively synthesize either Li-BIF-2 or Li-BIF-3 by neat milling or LAG (using MeOH or MeCN as liquid additives) with LiOH and a stoichiometric amount of HB(Meim)₄ were not successful. Exploration of different milling times and η-values produced only mixtures of residual reactants with Li-BIF-2, Li-BIF-3, and/or not yet identified phases (see ESI^†). Consequently, we explored milling in the presence of 2-aminobutanol (amb), which is a ubiquitous component of solvent systems used in the solvothermal syntheses of BIFs.^32,33 Gratifyingly, using a mixture of amb and MeCN in a 1 : 3 ratio by volume as the milling liquid led to an effective strategy for the selective synthesis of both the dia-topology Li-BIF-2 (CSD code MOXKUG), and the SOD-topology Li-BIF-3 (CSD code MUCLOM).^‡ The selective formation of phase-pure samples of Li-BIF-2 and Li-BIF-3 was confirmed by PXRD analysis, which revealed an excellent match to diffractograms simulated based on the previously reported structures (Fig. 4B–G). Systematic exploration of reaction conditions, including time (between 15 and 60 minutes) and η value (between 0.25 and 1 μL mg⁻¹) revealed that the open framework Li-BIF-3 is readily obtained at η either 0.75 or 1 μL mg⁻¹ after milling for 45 minutes or longer (Fig. 4B–G, also see ESI^†).^§ Lower η-values of 0.25 and 0.5 μL mg⁻¹ preferred the formation of the dia-topology Li-BIF-2, which was obtained as a phase-pure material upon 60 minutes milling at η = 0.5 μL mg⁻¹, following the initial appearance of a yet unidentified intermediate. The preferred formation of Li-BIF-2 at lower η-values is consistent with our previous observations that lower amounts of liquid promote mechanochemical formation of denser MOF polymorphs.³⁷Open in a separate windowFig. 4(A) Reaction scheme for the mechanochemical synthesis of Li-BIF-3. Comparison of selected PXRD patterns for the synthesis of Li-BIF-2 and Li-BIF-3: (B) H[B(Meim)₄] reactant; (C) LiOH reactant; (D) simulated for Li-BIF-3 (CSD MUCLOM); (E) simulated for Li-BIF-2 (CSD MOXKUG); (F) Li-BIF-3 mechanochemically synthesized by LAG for 60 minutes with a 1 : 3 by volume mixture of amb and MeCN (η = 1 μL mg⁻¹); and (G) Li-BIF-2 mechanochemically synthesized by LAG for 60 minutes with a 1 : 3 by volume mixture of amb and MeCN (η = 0.5 μL mg⁻¹). Comparison of selected PXRD patterns for the synthesis of Cu-BIF-2 and Li-BIF-3: (H) Cu₂O; (I) Cu-BIF-3 (CSD MOXJOZ); (J) Cu-BIF-2 (CSD MUCLIG); (K) Cu-BIF-3 mechanochemically synthesised by ILAG for 60 minutes using NH₄NO₃ ionic additive (5% by weight) and MeOH (η = 1 μL mg⁻¹); and (L) mechanochemically synthesised Cu-BIF-2 by ILAG for 90 minutes using NH₄NO₃ ionic additive (5% by weight) and MeOH (η = 0.5 μL mg⁻¹).Samples of both Li-BIF-2 and Li-BIF-3 after washing with MeCN were further characterized by FTIR-ATR, TGA in air, and analysis of metal content by ICP-MS (see ESI^†). Nitrogen sorption measurement on the mechanochemically obtained Li-BIF-3, after washing with MeCN and evacuation at 85 °C, revealed a highly microporous material with a Brunauer–Emmett–Teller (BET) surface area of 1010 m² g⁻¹ (Fig. 5A), which is close to the value expected from the crystal structure of the material (1200 m² g⁻¹, 32 For direct comparison with previous work,³² we also calculated the Langmuir surface area, revealing an almost 40% increase (1060 m² g⁻¹) compared to samples made solvothermally (762.5 m² g⁻¹) (Fig. 5A, inset).Experimental Brunauer–Emmett–Teller (BET) and Langmuir surface area (in m² g⁻¹) of mechanochemically synthesized SOD-topology BIFs, compared to previously measured and theoretically calculated values, along with average particle sizes (in nm) established by SEM and calculated energies (in eV) for all Li-, Cu-, and Ag-BIF polymorphs. The difference between calculated energies for SOD- and dia-polymorphs in each system is given as ΔE (in kJ mol⁻¹)

Do carbon nanotubes catalyse bromine/bromide redox chemistry?

The redox chemistries of both the bromide oxidation and bromine reduction reactions are studied at single multi-walled carbon nanotubes (MWCNTs) as a function of their electrical potential allowing inference of the electron transfer kinetics of the Br₂/Br⁻ redox couple, widely used in batteries. The nanotubes are shown to be mildly catalytic compared to a glassy carbon surface but much less as inferred from conventional voltammetry on porous ensembles of MWCNTs where the mixed transport regime masks the true catalytic response.

Schematic of a carbon nanotube impact in bromide solution.

The bromine–bromide redox couple plays an essential role in diverse energy storage devices including hydrogen–bromine, zinc–bromine, quinone–bromine, vanadium–bromide and bromide–polysulphide flow batteries.^1–5 The Br₂/Br⁻ redox couple is attractive as a cathode reaction due to its high standard potential, large solubility of both reagents, high power density and cost efficiency.⁶ The performance of such devices is generically limited by the thermodynamics and kinetics of the redox couple comprising the battery with fast (‘reversible’) electron transfer is essential. In many cases, including the Br₂/Br⁻ couple the electrode reaction involves more than one electron as given in the stoichiometric reaction:2Br⁻ − 2e⁻ ⇄ Br₂; E⁰ = 1.08 V vs. SHEwith, at high bromide concentrations, the possibility of the follow up chemical reaction⁷Br₂ + Br⁻ ⇄ Br₃⁻Since electrons are usually transferred sequentially this implies that the mechanism is multistep with any of the individual mechanistic steps in principle being rate limiting. For this reason catalysts are commonly required to enhance the electrode kinetics at otherwise favourable electrode materials. One type of catalyst which has seen wide usage, including for the Br₂/Br⁻ couple^8,9 are carbon nanotubes (CNTs) with suggested advantages which include high surface area and the inherent porosity of CNT composites.¹⁰ The deployment of CNTs as a porous composite presents a further level of complexity to the electrode reaction beyond its multistep character because of the ill-defined mass transport within the porous layer. In particular ascertaining the intrinsic electron transfer kinetics and hence the level of catalysis, if any, is essentially impossible since these are masked in the voltammetric response by diffusional mass transport effects.^11–14 Specifically the transport within the porous structure of CNT layers is dominated by thin-layer and other^15,16 effects which give the illusion of electrochemical reversibility. In order to unscramble possible electro-catalysis of the bromine/bromide couple a different approach is needed.In the following we study both the electro-oxidation of bromide (BOR) and the electro-reduction of bromine (BRR) at single MWCNTs via ‘nano-impact (aka ‘single entity’) electrochemistry’^17–20 in aqueous solution. In this approach a micro-wire electrode at a fixed potential is inserted in a suspension of CNTs in the solution of interest. From time to time a single CNT impacts the electrode, adopts the potential of the latter for the duration of the impact which in the case of CNTs can vary from 1–100 of seconds^21–23 and sustained catalytic currents flow if the oxidation/reduction of interest is faster at the nanotube in comparison with the micro-wire electrode. The catalytic currents are studied as a function of potential revealing the electron transfer kinetics. Fig. 1 shows the concept of the experiment.Open in a separate windowFig. 1Schematic representation of ‘nano-impact’ electrochemistry on a carbon micro wire electrode for the oxidation of aqueous bromide from which the kinetics of the BOR are inferred. Analogous experiments but showing negative impact currents allow the inference of the kinetics of the BRR.The BOR and BRR were studied first, however, voltammetrically at an unmodified glassy carbon (GC) electrode as shown in Fig. 2 (black line) using 5.0 mM solutions of either NaBr or Br₂ in 0.1 M HNO₃. The midpoint potential was 0.82 V versus the saturated calomel electrode (SCE) consistent with the literature values for the formal potential of the Br₂/Br⁻ couple.²⁴ The voltammograms were analysed to give transfer coefficients of 0.45 ± 0.01 and 0.33 ± 0.01 (ESI, Section 2^†) for the BOR and BRR respectively. Both processes were inferred to be diffusional and the diffusion coefficients D_Br⁻ and D_Br2 were calculated to be 2.05 (±0.04) × 10⁻⁵ cm² s⁻¹ and 1.50 (±0.04) × 10⁻⁵ cm² s⁻¹ (ESI, Section 3^†) using the Randles–Ševčík equation for an irreversible reaction the values are consistent with literature reports.²⁴ Then the electrodes were modified with 30 μg of MWCNTs consisting of ca. 125 monolayers (the calculation is given in the ESI, Section 9^†) of MWCNTs assuming that they are closely packed across the area of the GC electrode, and the resulting voltammograms are shown in Fig. 2 (red line). In comparison with the unmodified electrode, enhanced currents are seen for the Br₂/Br⁻ couple which partly reflects the enhanced capacitance of the interface reflecting in turn the large surface area of the deposited nanotubes (ca. 60–120 cm²). Larger signals are also seen indicating a thin layer contribution from the material occluded within the porous layer which also leads to the apparently quasi-reversible shape of the voltammograms obtained for both reactions. A log–log plot of peak current (I_p) vs. scan rate (ν) showed a gradient value of 0.68 (±0.01) and 0.66 (±0.03) for the BOR and BRR (ESI, Section 4^†) confirming a mixed mass transport regime^12,14 with a combination of semi-infinite diffusion and thin layer behaviour. The transition from the fully irreversible to the apparent quasi-reversible character is sometimes confused with electro-catalysis attributed to the CNTs rather than thin-layer diffusion. In order to ascertain the true catalytic response, single entity electrochemistry was measured to obtain the BOR and BRR responses at single CNTs.Open in a separate windowFig. 2Cyclic voltammograms at pristine GC (black line) and 30 μg MWCNTs dropcast on GC (red line) at a scan rate of 0.05 V s⁻¹ (a) for the bromide oxidation reaction (BOR) in 5.0 mM NaBr in 0.1 M HNO₃, (b) for the bromine reduction reaction (BRR) in 5.0 mM bromine in 0.1 M HNO₃.For single entity measurements, a clean carbon wire (CWE, length 1 mm and diameter 7 μm) working electrode was used. Chronoamperograms were recorded at a constant applied potential of 0.2 V vs. SCE and 1.3 V vs. SCE for the BOR and BRR respectively (5.0 mM solutions). These values were selected in the light of Fig. 2 to provide a large overpotential for each reaction. Clear oxidative and reductive current steps were observed (Fig. 3). These were ascribed to the arrival of a MWCNT at the electrode surface and the resulting catalytic electron transfer for the duration of the impact. No steps were observed in the absence of MWCNTs (ESI, Fig. S4^†). The average residence time of the MWCNT was 1.2 (±0.5) seconds and the frequency of the collisions was 0.3 (±0.1) impacts per second. The average impact current for the BOR at 1.3 V vs. SCE was 2.8 (±0.2) nA (65 impacts) and for the BRR at 0.2 V vs. SCE it was 3.8 (±0.1) nA (70 impacts). The impact currents were assumed to be entirely faradaic since control experiments in 0.1 M HNO₃ solution in the presence of 100 μg of MWCNTs (in the absence of Br⁻ and Br₂) showed no obvious impacts as shown in ESI Section 10.^†Open in a separate windowFig. 3Chronoamperograms showing the impact step current (a) for the BOR in 5.0 mM NaBr in 0.1 M HNO₃ at 1.3 V vs. SCE, (b) for the BRR in 5.0 mM bromine in 0.1 M HNO₃ at 0.2 V vs. SCE.Further, impacts for both the BOR and BRR were observed at various potentials (ESI, Section 11^†) and analysed to obtain the average faradaic current at each potential. The average impact step current was plotted against the applied potential (Fig. 4). Two sigmoidal curves were obtained reflecting the current–potential response for either the bromide oxidation (BOR) or the bromine reduction (BRR). The curves reflect the average voltammograms (current–potential characteristics) for the Br₂/Br⁻ redox reaction at single carbon nanotubes. The shape of the two sigmoidal curves reflects the onset of electrolysis followed by a diffusion controlled plateau at high over-potentials.²⁵ Mass transport corrected Tafel analysis (Fig. 4; inset) showed the transfer coefficients β to be ca. 0.42 and α to be ca. 0.20 from the impacts for the BOR and BRR respectively (ESI, Section 6^†). The length distribution of the MWCNTs was calculated (ESI, Section 6^†) from the currents recorded at potentials corresponding to the plateau in Fig. 4 assuming that the reactions are (Fickian) diffusion controlled at the potentials used and by modelling the CNTs as cylindrical electrodes²¹ assuming a nanotube radius of 15 (±5) nm and the diffusion coefficients reported above. Chronoamperometry was also conducted for the BOR and BRR in the absence of MWCNTs at 1.3 V and 0.2 V vs. SCE respectively to confirm that no impact currents were contributed by the redox species in the electrolyte (ESI, Section 5^†). Alongside, chronoamperograms in 0.1 M HNO₃ and 100 μg show that the impact current was contributed only by the Br⁻ and Br₂ redox reaction and the results are shown in the ESI, Section 10.^†Open in a separate windowFig. 4Average step currents observed as a function of applied potential (a) for the BOR in 5.0 mM NaBr in 0.1 M HNO₃ at, (b) for the BRR in 5.0 mM Bromine in 0.1 M HNO₃; insets in both the cases show mass transport corrected Tafel analyses.The lengths were found to be 5.4 (±3.4) μm (BOR) and 5.9 (±1.3) μm (BRR) and are given in Fig. 5 (see ESI, Section 7^† for calculations). These values were compared with previously reported dark-field optical microscopy data and good agreement was observed with the literature value of 5.3 (±2.1) μm.²⁶ The observed consistency provides strong support for the choice of modelling the single entity voltammetry by analogy with that of a cylindrical electrode.Open in a separate windowFig. 5The length of MWCNTs calculated from the impact currents for the BOR (at 1.3 V vs. SCE) and BRR (at 0.2 V vs. SCE).It is evident that the single entity measurements allow a clear analysis of the catalytic behaviour of the carbon nanotubes by providing a well-defined diffusional regime conducive to the extraction of the electrode kinetics of both the bromide oxidation and the bromine reduction process. In contrast, electrodes were formed by ensembles of carbon nanotubes in the form of a porous layer where the mixed transport regime is not amenable to ready modelling and the dissection of thin-layer effects from the measured voltammetry. The electron transfer kinetics for both the BOR and BRR at single MWCNTs was then obtained via full simulation of the two single entity ‘voltammograms’ using the above measured diffusion coefficients and again treating the impacted MWCNT as a cylindrical electrode with uniform diffusional access and further assuming Butler–Volmer kinetics. For the BOR, one electron transfer was considered as given below,For the BRR the two electron transfer was modelled as,Br₂ + 2e⁻ → 2Br⁻The set of parameters used for the analysis are given in the ESI, Section 8.^† By using the transfer coefficients deduced from Fig. 4, the only unknown is the standard electrochemical rate constant k which is determined by fitting the impact voltammogram measured relative to a formal potential for the Br₂/Br⁻ couple of 0.82 V vs. SCE obtained from the voltammogram at pristine GC. Fig. 6 shows the fitting for the BOR and the BRR with rate constants k_BOR of 1.0 (±0.1) × 10⁻³ cm s⁻¹ and k_BRR of 5.0 (±0.1) × 10⁻⁴ cm s⁻¹ respectively. The transfer coefficients and rate constants obtained from impacts were compared to the voltammograms obtained at pristine GC for the BOR and BRR and are given in Open in a separate windowFig. 6DIGISIM simulated curves (black line) for average impact currents obtained at different potentials (red circles) (a) for the BOR with a rate constant (k_BOR) of 1.0 (±0.1) × 10⁻³ cm s⁻¹; (b) BRR with a k_BRR of 5.0 (±0.1) × 10⁻⁴ cm s⁻¹.Transfer coefficients and rate constants for the BOR in 5.0 mM NaBr in 0.1 M HNO₃ and the BRR in 5.0 mM bromine in 0.1 M HNO₃ obtained at the glassy carbon macroelectrode GC, and single MWCNT impact current

Reversing electron transfer in a covalent triazine framework for efficient photocatalytic hydrogen evolution

Covalent triazine-based frameworks (CTFs) have emerged as some of the most important materials for photocatalytic water splitting. However, development of CTF-based photocatalytic systems with non-platinum cocatalysts for highly efficient hydrogen evolution still remains a challenge. Herein, we demonstrated, for the first time, a one-step phosphidation strategy for simultaneously achieving phosphorus atom bonding with the benzene rings of CTFs and the anchoring of well-defined dicobalt phosphide (Co₂P) nanocrystals (∼7 nm). The hydrogen evolution activities of CTFs were significantly enhanced under simulated solar-light (7.6 mmol h⁻¹ g⁻¹), more than 20 times higher than that of the CTF/Co₂P composite. Both comparative experiments and in situ X-ray photoelectron spectroscopy reveal that the strong interfacial P–C bonding and the anchoring of the Co₂P cocatalyst reverse the charge transfer direction from triazine to benzene rings, promote charge separation, and accelerate hydrogen evolution. Thus, the rational anchoring of transition-metal phosphides on conjugated polymers should be a promising approach for developing highly efficient photocatalysts for hydrogen evolution.

Reversing the electron transfer in a covalent triazine-based framework by Co₂P anchoring achieved highly efficient photocatalytic hydrogen evolution from water splitting.

Photocatalytic water splitting into hydrogen fuels has been considered as a promising technique for converting solar energy into chemical energy.^1–3 To achieve this target, it is necessary to design and construct photocatalysts with high solar-to-hydrogen (STH) conversion efficiency.^4,5 Recently, covalent triazine-based frameworks (CTFs),^6,7 as a new class of conjugated polymer materials, have attracted significant attention in the photocatalytic water splitting field owing to their visible-light response, organized architecture, adjustable pore-size, and controllable functionalization.^8–12 However, owing to the high charge recombination and fewer active sites for the hydrogen evolution reaction (HER), the photocatalytic activities of pristine CTFs are very low. To overcome this drawback, noble-metal platinum (Pt) is generally required in CTF-based photocatalytic systems as the HER cocatalyst for promoting charge separation and catalyzing hydrogen generation.^13–15 However, the high cost and scarcity of metallic Pt greatly limit the large-scale commercial applications. Thus, it is highly desirable to explore economical materials as noble-metal substitutes for achieving comparable or even superior hydrogen evolution activities.Recently, transition-metal phosphides (TMPs) have attracted intensive attention as HER cocatalysts for photocatalytic water splitting, owing to their unique structural and electronic properties.^16–20 Up to now, a variety of semiconductor materials, including metal oxides,^21,22 metal sulfides,^23,24 g-C₃N₄,^25,26 MOFs,^27,28etc., decorated with TMPs have been extensively investigated, and the hydrogen generation activity in some reports is even higher than that of Pt cocatalysts. However, related studies about the decoration of TMP cocatalysts on CTF photocatalysts for highly efficient hydrogen generation have not been reported so far. Taking into account the molecular structure of CTFs, the nitrogen sites in the triazine frameworks could easily coordinate with transition metal ions to form nitrogen–metal interactions.^29–31 In contrast, the coordination of TMPs with CTFs is relatively difficult owing to the high electronegativity of P sites which could draw electrons from metal atoms.¹⁷ Moreover, most reported TMPs were fabricated directly from precursor metal salts, oxides, etc.,^32–35 and the resultant large-dimensions usually lead to very limited contact-interface with CTFs, or even a physical mixture form. Accordingly, the charge transfer between TMPs and CTFs is significantly restrained. Thus, effective anchoring of TMPs on CTFs for achieving highly efficient photocatalytic hydrogen evolution still remains a great challenge.Herein, a one-step phosphidation strategy has been developed to achieve phosphorus bonding with the benzene rings of CTFs and the anchoring of well-defined dicobalt phosphide (Co₂P) nanocrystals (∼7 nm). The photocatalytic results clearly reveal that an excellent hydrogen evolution activity (7.6 mmol h⁻¹ g⁻¹) is achieved, which is much higher than that of the CTF/Co₂P composite (0.37 mmol h⁻¹ g⁻¹). More detailed studies confirm that the phosphidation strategy could effectively facilitate the interfacial bonding between CTFs and Co₂P nanocrystals. More importantly, the SI-XPS results clearly suggest that the charge transfer direction in CTFs is completely reversed, and the photo-generated electrons efficiently transferred from the triazine rings to P-bonded benzene rings, where the Co₂P cocatalyst attracted electrons through the interfacial P–C bonds for efficient hydrogen evolution. To our knowledge, this is the first report on incorporating TMPs on a CTF photocatalyst for enhancing the hydrogen evolution activity. Fig. 1A shows the basic procedures for the preparation of CTF polymer anchored Co₂P nanocrystals. Briefly, the CTFs with adsorbed cobalt ions were directly phosphatized by thermal decomposition of NaH₂PO₂ under an Ar atmosphere, which could simultaneously achieve P atom bonding with benzene rings and the anchoring of the Co₂P cocatalyst (marked as P-CTF-Co₂P). Fig. 1B shows the typical transmission electron microscopy (TEM) image of the obtained sample, clearly revealing that well-defined Co₂P nanocrystals with an average diameter of 5∼9 nm were uniformly dispersed on the CTF surface. Furthermore, the high-resolution TEM (HR-TEM) image (Fig. 1C) of the formed nanocrystals exhibited two lattice fringes with d-spacing values of 0.20 and 0.21 nm, respectively, which could be well indexed to the (211) and (121) planes of orthorhombic Co₂P crystals.^36,37Fig. 1D shows the X-ray diffraction (XRD) pattern of the obtained P-CTF-Co₂P sample. For comparison, the pristine CTF, and Co₂P nanocrystals have also been studied. It can be clearly seen that for the P-CTF-Co₂P sample, except for the two broad diffraction peaks at 7.3° and 26.1° attributed to the in-plane (100) facets and interlayer (001) stacking of CTFs,^38,39 the other XRD peaks well matched with those of the Co₂P crystals,^40,41 confirming the successful incorporation of the Co₂P cocatalyst in the CTFs. Furthermore, energy-dispersive X-ray (EDX) elemental mapping (Fig. 1E) clearly reveals the uniform distribution of C, N, Co and P elements in the whole detection region, further confirming the uniform dispersion of Co₂P on CTF polymers.Open in a separate windowFig. 1(A) Basic procedures and the ideal structure scheme for fabricating P-CTF-Co₂P catalysts; (B) TEM image and (C) HR-TEM image of the P-CTF-Co₂P catalyst; (D) XRD patterns of CTFs, Co₂P and P-CTF-Co₂P catalysts; (E) EDX elemental mapping images of the P-CTF-Co₂P catalyst.To further explore the P-bonding sites, high-resolution X-ray photoelectron spectroscopy (XPS) was performed on both P-CTF-Co₂P and pristine CTFs. As shown in Fig. 2A, for pristine CTF samples, the C 1s peak could be well fitted into two peaks located at 284.8 and 286.8 eV, which could be assigned to C–C C and N–C N bonds,^38,39,42 respectively. Notably, after the phosphidation treatment, a new peak at 285.3 eV attributed to P–C bonds was detected in the P-CTF-Co₂P sample,^43–45 and the C–C C to N–C N ratio was significantly decreased from 10.2 (pristine CTFs) to 7.1 (P-CTF-Co₂P). In contrast, compared with the N 1s spectrum of the pristine CTF sample, no evident change of peak shape and intensity could be detected in the P-CTF-Co₂P sample (Fig. 2B). These results clearly reveal that after the phosphidation treatment, P atoms should mainly bond with the carbon sites of the benzene rings in CTFs instead of the triazine rings.^46,47 To further confirm this inference, XPS studies on phosphatized CTFs without Co₂P nanocrystals (marked as P-CTFs) were also performed. It can be clearly seen from Fig. S9^† that similar changes for the C 1s and N 1s peaks of P-CTF-Co₂P (Fig. 2A and B) are observed in the P-CTF sample. Moreover, in both P-CTF-Co₂P and P-CTF samples, an evident P 2p peak corresponding to P–C bonds could be observed (Fig. 2C),⁴⁸ further indicating the bonding of P atoms with C atoms on the benzene rings in CTFs after the phosphidation treatment. However, note that the binding energy (BE) position of P–C peaks in P-CTF-Co₂P (133.6 eV) is slightly higher than that of P-CTF (133.2 eV), which should be attributed to the anchoring with Co₂P nanocrystals. In addition to the XPS investigations, solid-state cross-polarization magic angle spinning carbon-13 and phosphorus-31 nuclear magnetic resonance (¹³C-NMR and ³¹P-NMR) spectroscopy studies were further performed to explore the bonding sites of phosphorus atoms in the P-CTF (Fig. S10^†). The Fourier-transformed infrared (FTIR) spectrum was also employed (Fig. 2D). It can be clearly observed that the typical FTIR peaks of the triazine frameworks and the stretching vibrations of carbon–nitrogen (C–N) (1521 cm⁻¹ and 1354 cm⁻¹) in P-CTF-Co₂P are consistent with those of pristine CTFs.^38,39 However, the peaks at 1672 cm⁻¹ and 1014 cm⁻¹ corresponding to the stretching vibrations of C C and C–H in benzene rings were evidently decreased in P-CTF-Co₂P compared with those of pristine CTFs.^49,50 Combining the results of XPS, NMR, and FTIR, it can be concluded that in the obtained P-CTF-Co₂P sample, the carbon sites of the benzene rings in the CTFs were partially bonded with P-atoms, which should anchor Co₂P through P–Co bonding.Open in a separate windowFig. 2High-resolution XPS spectra of (A) C 1, (B) N 1s and (C) P 2p in CTF, P-CTF and P-CTF-Co₂P samples (the insets show the molecular structures); (D) Fourier transform infrared spectroscopy (FT-IR) of CTF and P-CTF-Co₂P samples.Furthermore, the hydrogen evolution activities of the P-CTF-Co₂P (2 wt% Co₂P) sample were evaluated under simulated solar light irradiation. For comparison, the CTF/Co₂P composite and pristine CTFs were also measured under the same conditions. As shown in Fig. 3A, the P-CTF-Co₂P sample exhibits a much higher H₂ evolution activity (15.2 mmol g⁻¹) than the CTF/Co₂P composite (0.7 mmol g⁻¹) at 2 h, while no evident hydrogen generation could be detected in the pristine CTF sample. To further explore the crucial roles of interfacial P-bonding in CTFs and the anchoring of the Co₂P cocatalyst, the hydrogen evolution rates of various samples were calculated and their comparison is shown in Fig. 3B. Obviously, compared with the excellent activity of P-CTF-Co₂P (7.6 mmol h⁻¹ g⁻¹), CTF/Co₂P and P-CTF/Co₂P only exhibit very low hydrogen evolution rates of 0.37 and 0.38 mmol h⁻¹ g⁻¹, respectively, while no evident activity could be detected for CTFs, Co₂P, and P-CTF samples. The above results clearly confirm that the significant enhancement of the hydrogen evolution activity of the P-CTF-Co₂P sample should be mainly attributed to the P-bonding and the anchoring of the Co₂P cocatalyst. Furthermore, the photocatalytic performances of CTFs decorated with Pt nanoparticles were also studied (the inset of Fig. 3B, Fig. S15^†), as Pt is generally recognized as the most active cocatalyst for hydrogen generation. Amazingly, the hydrogen production rates of the Pt cocatalyst at different amounts (1∼5 wt%) were all lower than that of the P-CTF-Co₂P sample, indicating that the rational bonding of the Co₂P cocatalyst on CTFs should be a promising strategy for achieving highly efficient photocatalytic hydrogen evolution.Open in a separate windowFig. 3(A) Photocatalytic H₂ evolution tests of the as-prepared photocatalysts; (B) comparative presentation of the hydrogen evolution rates; (C) ultraviolet-visible diffuse reflectance spectrum of the P-CTF-Co₂P photocatalyst (red solid line) and wavelength-dependent hydrogen production activities of the P-CTF-Co₂P photocatalyst within 2 h; light source: a 300 W Xe lamp equipped with various cut-off filters. (D) The cycling photocatalytic tests of P-CTF-Co₂P and the CTF/Co₂P composite; (E) I–t curves of various catalysts at −0.3 V (vs. SCE); (F) electrochemical impedance spectroscopy under light of the various catalysts. Measurements were conducted in 0.2 mol L⁻¹ Na₂SO₄ electrolyte solution under AM 1.5 G illumination.Furthermore, the relationship between wavelength and the hydrogen evolution activity of the P-CTF-Co₂P photocatalyst was studied and is shown in Fig. 3C. With increasing the wavelength from full-spectrum light irradiation (λ > 300) to 500 nm, the hydrogen evolution rate significantly decreased from 7.6 to 0.06 mmol h⁻¹ g⁻¹, which is generally consistent with their absorption spectrum. The highest value of the apparent quantum yields (AQYs) attained is 31.8% at 365 nm. Moreover, the photocatalytic stability and durability of P-CTF-Co₂P and CTF/Co₂P were examined by cycling experiments. As shown in Fig. 3D, the CTF/Co₂P composite demonstrated relatively poor stability, and there is no evident activity after only one cycling experiment, resulting from the incompact combination between CTFs and Co₂P (Fig. S5^†). In contrast, the P-CTF-Co₂P sample exhibits relatively high stability, and no evident inactivation was detected during the whole test. These results further confirm that in addition to enhancing activities, the interfacial P–C bonds and the anchoring of the Co₂P cocatalyst could also effectively promote the hydrogen evolution stability. To further confirm the efficient charge separation and electron transfer in P-CTF-Co₂P, photoelectrochemical (PEC) tests were performed. As shown in Fig. 3E, the amperometric I–t curves clearly reveal that P-CTF-Co₂P exhibits higher photocurrent density than pristine CTFs and CTF/Co₂P samples, confirming its more efficient charge separation capability. Moreover, electrochemical impedance spectroscopy (EIS) was also performed to explore the interfacial charge transfer process.^51,52 As shown in Fig. 3F, the P-CTF-Co₂P sample with the smallest arc radius revealed remarkably increased interfacial charge transport efficiency. These PEC results clearly demonstrate the highly efficient photogenerated charge separation and transfer in the P-CTF-Co₂P photocatalyst, which are highly consistent with the above photocatalytic hydrogen evolution results.The photo-induced charge separation and transfer in the excited state are crucial processes for determining the photocatalytic activity. Herein, we have demonstrated an in situ irradiation X-ray photoelectron spectroscopy (SI-XPS) technique for exploring the intrinsic charge separation and transfer mechanisms between Co₂P and CTFs. As shown in Fig. S20 and S21,^† no evident binding energy (BE) shift could be observed in both pristine CTF and Co₂P samples, confirming their relatively poor charge separation capability. Amazingly, for the P-CTF-Co₂P sample (Fig. 4), distinct variations for C 1s, N 1s, P 2p, and Co 2p peaks were detected in the excited state. More specifically, the C 1s and N 1s peaks in the CTFs shifted towards the high BE region by 0.3 eV accompanied by the broadening of peak shape (Fig. 4B and C), while the P 2p and Co 2p peaks of the Co₂P cocatalyst shifted towards the low BE direction by 0.3 and 0.4 eV (Fig. 4D and E), respectively.Interestingly, as shown in Fig. 4D, the P 2p peaks attributed to interfacial P–C bonds exhibited no evident BE change under dark and light irradiation. On the basis of the above SI-XPS results, it can be concluded that under light irradiation, the photo-excited electrons effectively transferred from CTFs to Co₂P through the P-bonding sites as well as the interfacial P–C bonds, leading to electron enrichment in Co₂P nanocrystals and hole enrichment in CTFs. To further confirm the crucial roles of interfacial P-bonding in CTFs and the anchoring of the Co₂P cocatalyst in promoting charge separation, SI-XPS studies for the CTF/Co₂P composite were also conducted. As shown in Fig. S22,^† no new peak or shape change of C 1s, N 1s, P 2p, and Co 2p peaks could be detected in CTF/Co₂P compared with pristine CTFs and Co₂P in the ground state, indicating no interfacial bonding in the CTF/Co₂P composite. Furthermore, in the excited state, no evident BE shifts and shape change could be observed, indicating the poor charge-separation capability of CTF/Co₂P, which is highly consistent with its fairly low hydrogen evolution activity. These SI-XPS results further confirm the crucial roles of interfacial P-bonding in CTFs and the anchoring of the Co₂P cocatalyst in efficiently promoting charge separation and enhancing the photocatalytic activity.Open in a separate windowFig. 4(A) Schematic illustration of the SI-XPS technique for direct observation of electron transfer in the excited state. The (B) C 1s, (C) N 1s, (D) P 2p and (E) Co 2p of SI-XPS spectra in the P-CTF-Co₂P sample tested under dark and light illumination.On the basis of the above results, it can be concluded that the P-bonding and Co₂P-cocatalyst anchoring should reverse the photo-induced electron transfer direction from the triazine to benzene rings in CTFs. To further confirm this speculation, the photocatalytic behavior of various CTF-based photocatalysts was studied and their electron transfer directions in excitation states have been proposed (Fig. 5A and B). First, when Co²⁺ ions were further bonded with the N sites of the triazine rings in the P-CTF-Co₂P sample (marked as P-CTF-Co₂P-Co), the hydrogen evolution activity remarkably increased from 7.6 up to 8.4 mmol h⁻¹ g⁻¹. This result clearly reveals that the Co–N bonding in P-CTF-Co₂P could further promote charge separation due to the effective hole trapping on Co sites. In contrast, when the Pt cocatalyst was decorated on the P-CTF-Co₂P sample (marked as P-CTF-Co₂P–Pt) to form Pt–N coordination, the photocatalytic activity significantly decreased to 4.6 mmol h⁻¹ g⁻¹, mainly resulting from the electron transfer competition between Pt and Co₂P. Furthermore, the decoration of the single Pt cocatalyst on CTFs promoted the photo-induced electron transfer from the benzene rings to triazine rings, and the photocatalytic activity achieved was up to 3.8 mmol h⁻¹ g⁻¹. Furthermore, in situ FTIR spectroscopy was performed for exploring the intermediate products in the photocatalytic process of the P-CTF-Co₂P sample with co-adsorption of H₂O vapor (Fig. S27^†), which also demonstrated the direction of electron transfer combined with the in situ XPS results. Accordingly, a possible mechanism has been proposed for clarifying the hydrogen evolution activities of the P-CTF-Co₂P photocatalyst (Fig. 5B). Owing to the P-bonding and the anchoring of Co₂P nanocrystals, the photo-generated electrons could effectively transfer from the triazine rings to benzene rings, where the electrons are trapped by the Co₂P cocatalyst through the interfacial P–C bonds for the hydrogen evolution reaction. Simultaneously, the holes left behind in the CTFs participated in oxidation reactions. More importantly, this mechanism may provide a new insight for understanding the crucial roles of interfacial P-bonding in CTFs and their bonding with Co₂P in promoting the charge separation for efficient hydrogen evolution.Open in a separate windowFig. 5(A) Photocatalytic hydrogen evolution activities of various samples and the scheme of electron-transfer direction in the excited state. (B) The ideal structural illustration of the interfacial bonding and charge transfer of the P-CTF-Co₂P photocatalyst for hydrogen evolution.     

Lithium achieves sequence selective ring-opening terpolymerisation (ROTERP) of ternary monomer mixtures

Heteroatom-containing degradable polymers have strong potential as sustainable replacements for petrochemically derived materials. However, to accelerate and broaden their uptake greater structural diversity and new synthetic methodologies are required. Here we report a sequence selective ring-opening terpolymerisation (ROTERP), in which three monomers (A, B, C) are selectively enchained into an (ABA′C)_n sequence by a simple lithium catalyst. Degradable poly(ester-alt-ester-alt-trithiocarbonate)s are obtained in a M_n range from 2.35 to 111.20 kDa which are not easily accessible via other polymerisation methodologies. The choice of alkali metal is key to achieve high activity and to control the terpolymer sequence. ROTERP is mechanistically compatible with ring-opening polymerisation (ROP) allowing switchable catalysis for blockpolymer synthesis. The ROTERP demonstrated in this study could be the first example of an entirely new family of sequence selective terpolymerisations.

A sequence selective ring-opening terpolymerisation of epoxides with CS₂ and phthalic thioanhydride yielding poly(ester-alt-ester-alt-trithiocarbonates) is reported.

Synthetic polymers are now more in demand than ever before and looking at their annually increasing production, a polymer free society is at best a vague memory rather than a vision for the future.¹ As most commodity polymers are based on chemically inert aliphatic –C–C– backbones, most polymer waste shows unappreciable degradation with respect to the polymer''s time in application.^2,3 In answer to the ever-increasing amount of plastic pollution, much effort focuses on the exploration of new heteroatom containing polymers which, because of their more polar bonds making up the backbone, are more susceptible to degradation via physical, chemical and biochemical pathways and even facilitate new chemical recycling methods.^4–8 Besides, there is also a constant demand for entirely new materials to enable technological innovation making new methodologies to synthesise heteroatom containing polymers necessary.Arguably one of the most popular methods to make such polymers is the ring-opening polymerisation (ROP) of a heterocycle A.^9–11 These polymerise under release of their associated ring strain energy to make polymers (A)_n such as poly(thio)ethers, poly(thio)ester and poly(thio)carbonates, only to name a few. Early on it has been realised that in some cases the ROP of three or four-membered heterocycles A can be coupled to the insertion of typically heteroallenes or cyclic anhydrides B to generate alternating copolymers (AB)_n.^12–14 The underlying requirements for (AB)_n polymerisations are a combination of kinetic factors, i.e. monomer A inserting orders of magnitude faster into the active catalyst growing chain-bond than monomer B, and chemoselectivities, i.e. insertion of monomer A resulting in type A active catalyst growing chain-bond that does not show appreciable reactivity with A but only with B (and vice versa). Prominent examples of this alternating (AB)_n ring-opening copolymerisation (ROCOP) include CO₂/epoxide ROCOP yielding polycarbonates and cyclic anhydride/epoxide ROCOP yielding polyesters.^15–17 Sulfur analogous are also accessible such as polythiocarbonate from CS₂/(epoxide or thiirane) and polythioesters from thioanhydride/(epoxide or thiirane) ROCOP.^18–31 Such sulfur rich polymers are attractive high-refractive index materials that can show improved crystallinity and degradability over their all-oxygen analogues and in some cases enable chemical polymer to monomer recycling.^32–39Most relevant to this study is a report by Werner and coworkers on lithium alkoxide catalysed CS₂/epoxide ROCOP yielding poly(monothio-alt-trithiocarbonate)s featuring R–O–C( S)–O–R and R–S–C( S)–S–R links (Fig. 1).^40,41 Such a polymer sequence is unexpected, as the formal product of alternating insertion would be a poly(dithiocarbonate) with R–O–C( S)–S–R links. Furthermore, the polymer shows an unusual “head-to-head-alt-tail-to-tail” selectivity meaning that the R–O–C( S)–O–R links sit adjacent to the tertiary CHR₃ positions of the ring opened epoxide (i.e. “head” position) and that the R–S–C( S)–S–R links sit adjacent to the secondary CH₂R₂ position of the ring opened epoxide (i.e. “tail” position). This sequence let the authors postulate a mechanism involving tail-selective epoxide ring-opening by a dithiocarbonate chain end which is formed by CS₂ insertion into an alkoxide intermediate. The resulting alkoxide intermediate was proposed to backbite into the adjacent dithiocarbonate link which after a rearrangement process resulted in an O/S exchange of the chain end. The rearrangement transforms the alkoxide into a thiolate chain end and the adjacent dithiocarbonate R–O–C( S)–S–R into a monothiocarbonate R–O–C( S)–O–R. CS₂ insertion of the thiolate generates a trithiocarbonate R–S–C( S)–S–R which after epoxide insertion regenerates the alkoxides. In contrast to alkoxides sitting adjacent to R–O–C( S)–S–R links, alkoxides next to R–S–C( S)–S–R links were not proposed to undergo backbiting and O/S exchange. Importantly the initial regioselectivity of the epoxide ring-opening was preserved throughout the rearrangement which explained the “head-to-head-alt-tail-to-tail” selectivity. Interestingly lithium appeared to be crucial as other alkali metal alkoxides failed to catalyse this ROCOP, while more sophisticated catalysts result in much more uncontrolled polymer sequences. Hence it appears that the Li controls which alkoxide intermediate precisely undergoes O/S exchange and to which degree this rearrangement occurs, but the reasons behind the special role of Li remains to be explored. Relatedly the ROCOP of thioanhydrides with epoxides has also been reported and similar exchange processes have been proposed as side reactions.²³ Although not directly proven, this mechanism seemed reasonable and let us hypothesise that lithium catalysts could grant a general access to control the O/S exchange process in ROCOP and even in the polymerisation of ternary monomer mixtures. Furthermore, we reasoned that the two distinct chain ends formed via O/S exchange, i.e. alkoxide and thiolate (type A vide supra), could enable discrimination between two different type B monomers and enable sequence selective terpolymerisations. It should be noted that reports exist in which mixtures of A, B and C (e.g. epoxide A, cyclic anhydride B and CO₂ C) either yield random terpolymers, (AB)_n-ran-(AC)_m poly(esters-ran-carbonate), or block-terpolymers, (AB)_n-b-(AC)_m, polyester-b-polycarbonate. In this case the monomer sequence depends on catalyst selection and reaction conditions, but sequence selective terpolymers, e.g. (ABC)_n or (ABAC)_n, are unknown.^42–48 The hypothesis of O/S exchange here led us to discover a new type of polymerisation, sequence selective ring-opening terpolymerisation (ROTERP) that we report in this contribution. ROTERP produces poly(ester-alt-ester-alt-trithiocarbonates), i.e. (ABA′C)_n sequences, from a mixture of the monosubstituted epoxides propylene oxide (PO) or butylene oxide (BO) A, phthalic thioanhydride (PTA) B and CS₂ C, employing simple lithium salts as the catalyst. Furthermore, model reactions proof the previously postulated O/S rearrangement that enable ROTERP and elucidate the role of the lithium catalyst. Finally, we employed ROTERP in so-called switchable catalysis, in which the onset of ROTERP stops the occurrence of ROP, for the construction of blockpolymers.Open in a separate windowFig. 1(Top) CS₂/epoxide ROCOP and postulated mechanism involving a central O/S exchange reaction, (middle) phthalic thioanhydride (PTA)/epoxide ROCOP and (bottom) PTA/CS₂/epoxide ROTERP reported in this study. R = Me, Et; [R_n] = polymer chain.ROTERP of mixtures of PTA, PO and CS₂ at loadings typically employed in ROCOP catalysis with lithium hexamethyldisilazide (LiHDMS) or lithium benzyloxide (LiOBn) as the catalyst at loadings in the range of 1 eq. LiX : (6.25–500 eq.) PTA : (31.25–2500 eq.) PO : (62.5–5000 eq.) CS₂ yield polymeric materials in 95% selectivity at 80 °C (see 40,47 Spectroscopic analysis of the isolated polymer reveals surprisingly clean NMR spectra given the potential for statistical terpolymerisation of these three monomers. The ¹H NMR spectrum (Fig. 2) shows two main aryl resonances for a symmetrically substituted terephthalate unit from ring opened PTA (δ = 7.72 and 7.56 ppm) as well as one main resonance for the CHMe (δ = 5.42 ppm; head position of the ring-opened epoxide) and CH₂ (δ = 3.90–3.40 ppm; tail position of the ring opened epoxide) groups respectively stemming from the ring opened PO. Correspondingly the ¹³C{¹H} NMR spectrum (ESI Fig. S1^†) reveals the almost exclusive presence of trithiocarbonate R–S–C( S)–S–R (δ = 222.8 ppm) and arylester R–C( O)–O–R (δ ∼166.6 ppm) resonances (94–98%) alongside minor thioester R–C( O)–S–R (δ = 192.7 ppm) resonances (2–6%). ¹H and 2D NMR spectra (ESI Fig. S3^†) show that trithiocarbonate units are positioned adjacent to CH₂ groups while arylesters are connected to the tertiary CHMe groups. The spectra remain unchanged after multiple precipitations from DCM : MeOH or THF : pentane and DOSY NMR shows that all ¹H NMR resonances diffuse at the same rate confirming that these are part of the same species. Furthermore, no other type of thiocarbonate R–(O/S)–C( O/S)–(O/S)–R are part of the polymer. Linkage identity was further substantiated by the ATR-IR spectrum (ESI Fig. S5^†) of these polymers showing an arylester C O stretch at = 1716 cm⁻¹ as well as a thiocarbonate C S stretch at = 1062 cm⁻¹ (ESI Fig. S5^†). Accordingly, the polymers are obtained as yellow solids due the presence of the C S chromophore (λ_abs = 435 nm, ESI Fig. S6^†). MALDI-TOF analysis unfortunately only led to decomposition of the materials and no signals could be identified as previously reported for sulfur-rich polymers.^26,38 Nevertheless ¹H NMR allows some conclusions regarding the topology as both for LiHMDS or LiOBn initiation, resonances for the HMDS and OBn groups can be identified to be part of the purified polymers. Insertion of alkalimetal alkoxides and amides into CS₂ yielding alkali dithiocarbonates and dithiocarbamates have been previously reported.^40,79 This makes initiation via CS₂ insertion likely which defines one end of the polymer and hence indicates that chains are linear rather than cyclic. As ROTERP is followed by CS₂/epoxide coupling once all PTA is consumed (vide infra) we infer that chains are terminated by CS₂ because heteroallene insertion products are generally established to be the resting states of heteroallene/heterocycle coupling reactions.¹²Data showing PTA/CS₂/epoxide ROTERP under different conditions

Ligand design for Rh(iii)-catalyzed C–H activation: an unsymmetrical cyclopentadienyl group enables a regioselective synthesis of dihydroisoquinolones

We report the regioselective synthesis of dihydroisoquinolones from aliphatic alkenes and O-pivaloyl benzhydroxamic acids mediated by a Rh(iii) precatalyst bearing sterically bulky substituents. While the prototypical Cp* ligand provides product with low selectivity, sterically bulky Cp^t affords product with excellent regioselectivity for a range of benzhydroxamic acids and alkenes. Crystallographic evidence offers insight as to the source of the increased regioselectivity.C–H activation mediated processes have provided a unique retrosynthetic approach to access a variety of substituted heterocycles.¹ One tactic that has received increased attention is the coupling of π-components with heteroatom containing molecules.² A variety of transition metals are capable of catalyzing this type of transformation, providing access to dozens of heterocyclic motifs.^1–3 A challenge for these methods is controlling the regioselectivity of migratory insertion across alkenes and alkynes after the metallacycle forming C–H activation (eqn 1).Steric and electronic effects are understood to control migratory insertion of unsymmetrical alkynes in Rh(iii) catalyzed isoquinolone syntheses (eqn 1). When the substituents are electronically similar, the larger group resides β- to Rh in the metallacycle to avoid unfavorable steric interactions (selectivity is generally >10 : 1).⁴ When the substituents are electronically different, the more electron-donating group prefers being α- to rhodium in the metallacycle, presumably to stabilize the electron poor metal.^5,6 The type of C–H bond being activated also plays an important role in the regioselectivity of migratory insertion; aromatic substrates typically provide synthetically useful regioselectivities when electronically different alkynes are used (>10 : 1) but alkenyl C–H activation leads to products with lower regioselectivities, presumably due to minimal steric interactions during migratory insertion.^7,8 We found that sterically bulky di-tert-butylcyclopentadienyl ligand (Cp^t) enhances the regioselectivity of the alkyne migratory insertion event in these cases, delivering regioselectivities (>10 : 1) modestly above those achievable by Cp* ligated Rh complexes (<6 : 1). However, when the alkyne migratory insertion was poorly selective with RhCp* (<3 : 1), RhCp^t complex was ineffective at providing synthetically useful levels of selectivity. Furthermore, the Cp^t ligand was only effective with aryl substituted alkynes, presumably because of strong steric interactions between the ligand and alkyne in the insertion event. Migratory insertion of alkenes to access heterocycles using C–H activation chemistry is still relatively rare, with seminal studies by Glorius and Fagnou reporting the synthesis of dihydroisoquinolones.^9–11 Similar to alkynes, alkenyl electron-donating groups favor the position adjacent to the metal in the metallacycle delivering high regioselectivity. In contrast to alkynes, aliphatic alkenes afford product with poor regioselectivity (2 : 1) (eqn 2).^5h,12 We hypothesized competing steric and electronic effects cause the low regioselectivity, with steric effects favoring the formation of a 4-substituted product and electronics favoring the formation of a 3-substituted product.¹³ As a temporary solution to this problem, our group and others have employed tethering strategies to increase the regioselectivity of the migratory insertion event (eqn 3).^14,15 Of course, regioselectivity controlled by the ligand on Rh would be the optimal solution to the selectivity problem (eqn 4).¹⁶ Consequently, we focused our attention toward developing an intermolecular variant of this reaction that would provide product with improved regioselectivity.As a model system, we explored the impact ligands have on the coupling of O-pivaloyl-benzhydroxamic acid 1a with 1-decene 2a to provide dihydroisoquinolones 3a and 3a′. When Cp* is used as a ligand, the desired products are isolated in excellent yield but poor selectivity (2.4 : 1 3a : 3a′) ( a

Modular assembly of MOF-derived carbon nanofibers into macroarchitectures for water treatment

The organized assembly of nanoparticles into complex macroarchitectures opens up a promising pathway to create functional materials. Here, we demonstrate a scalable strategy to fabricate macroarchitectures with high compressibility and elasticity from hollow particle-based carbon nanofibers. This strategy causes zeolitic imidazolate framework (ZIF-8)-polyacrylonitrile nanofibers to assemble into centimetre-sized aerogels (ZIF-8/NFAs) with expected shapes and tunable functions on a large scale. On further carbonization of ZIF-8/NFAs, ZIF-8 nanoparticles are transformed into a hollow structure to form the carbon nanofiber aerogels (CNFAs). The resulting CNFAs integrate the properties of zero-dimensional hollow structures, one-dimensional nanofibers, and three-dimensional carbon aerogels, and exhibit a low density of 7.32 mg cm⁻³, high mechanical strength (rapid recovery from 80% strain), outstanding adsorption capacity, and excellent photo-thermal conversion potential. These results provide a platform for the future development of macroarchitectured assemblies from nanometres to centimetres and facilitate the design of multifunctional materials.

A scalable strategy is established to generate macroarchitectures based on MOF-related nanofibers. The modular assembly of macroarchitectures with luffa-like structures exhibits high mechanical strength and low densities.

The assembly of simple nanoparticles (such as silica, polystyrene and metal–organic frameworks) into macroarchitectures has a unique attraction for engineering materials due to their variable sizes, shapes, and chemical and physical properties.^1–3 As a novel nanomaterial, the formed macroarchitecture with three-dimensional (3D) porous interconnected network structures has broad application prospects in various fields, including environment treatment, chemical sensing, energy storage, catalysis, and advanced electronic devices.^4–7 Moreover, the functions of macroarchitectures are mainly determined by the fundamental building blocks. On account of large surface area, high porosity and more exposure to active sites, the complex macroarchitectures, which are assembled by building blocks with hollow structures, possess greater advantages.^8,9In the past few decades, in order to seek high-performance hollow building blocks for macroarchitectures, much effort has been put into it. In particular, the emergence of one-dimensional (1D) carbon hollow nanostructures, including hollow porous carbon nanofibers (HPCNs)¹⁰ and carbon nanotubes (CNTs)¹¹ promotes the rapid development of this field. The HPCN- and CNT-based macroarchitectures realize the transformation from 1D nanomaterials to three-dimensional (3D) macroscopic materials with excellent properties (e.g., electrochemical energy storage and antimicrobial air filtration).^12,13 Furthermore, these macroarchitectures can not only retain the characteristics of the 1D material, but also generate many new kinds of features (e.g., high specific surface area, high mechanical strength, and low density) that the components do not possess. Unfortunately, though the HPCN- and CNT-based macroarchitectures exhibit improved conductivity and stability properties, the synthesis of their building blocks is usually expensive and complex.¹⁴ Besides, the assembly of building blocks into 3D macroarchitectures usually exhibits relatively poor mechanical properties and requires some adhesives or templates, which have to be eliminated by extra strategies.¹⁵ These complex synthetic procedures and less favourable structural stability largely hinder the scale-up production of carbon aerogels and their practical applications.To address these issues, we explore a novel and scalable method to synthesize functional macroarchitectures with robust mechanical properties fabricated from MOF-derived carbon nanofibers through manipulating nano-sized particles (MOFs) and micron-sized fibers. First, zeolitic imidazolate framework (ZIF-8) nanoparticles, which are spectacular for their large nitrogen content and surface area,^16–18 are embedded into polyacrylonitrile (PAN)/polyvinylpyrrolidone (PVP) nanofibers to form ZIF-8-PAN/PVP composite nanofibers. Subsequently, the ZIF-8-PAN/PVP nanofibers are assembled into centimetre-sized nanofiber aerogels (ZIF-8/NFAs) by a freeze-drying technique. After preoxidation and carbonization of ZIF-8/NFAs, the carbon nanofiber aerogels with hollow and porous interlayer structures are fabricated (named C-ZIF-8-CNFAs). The interlayer structure of C-ZIF-8-CNFAs is very similar to natural luffa consisting of a network of elastic frameworks. As with the interconnected nanofibers in the interlayer structure, cellulose skeletons are in the interior of luffa interconnected in a highly uniform manner to maximize strength, a porous structure. Therefore, we refer to our porous structure as luffa-like. The prepared C-ZIF-8-CNFAs exhibit a low density of 7.32 mg cm⁻³, high specific surface area (288 m² g⁻¹), large hierarchical pore volume (0.22 cm⁻³ g⁻¹), high mechanical strength (rapid recovery from 80% strain), outstanding adsorption capacity, and excellent photo-thermal conversion potential. Fig. 1a depicts the assembly strategy of macroarchitectures (C-ZIF-8-CNFAs) schematically. The preparation process begins with the fabrication of ZIF-8/nanofibers (Fig. S1b^†) using the method that ZIF-8 nanoparticles (Fig. S1a^†) are incorporated into nanofibers by electrospinning (specific preparation methods in the ESI^†). Upon homogenization in a mixed solution of ultrapure water and tert-butanol, the ZIF-8/nanofibers become wrapped around each other and dispersed uniformly. Subsequently, the homogenized nanofiber dispersion is frozen in a mold followed by freeze-drying into uncrosslinked ZIF-8/NFAs. To build further robust bonding among nanofibers, the obtained uncrosslinked ZIF-8/NFAs are preoxidized at 250 °C to form crosslinked ZIF-8/NFAs with a welding structure under the action of polyvinylpyrrolidone (PVP), providing elastic resilience to the resultant NFAs. Ultimately, the resulting preoxidized ZIF-8/NFAs are carbonized at 900 °C to form hollow C-ZIF-8-CNFAs under a N₂ atmosphere. As a result of the carbonization, the organic ligands of ZIF-8 and molecules of PAN and PVP are decomposed and transformed into N-doped carbon materials. In addition, the Zn²⁺ in ZIF-8 nanoparticles is reduced to metallic Zn and then evaporated at high temperature.^19–21 Nanoscale 0D MOFs are assembled into micron-scale 1D fibers, which are then assembled into centimeter-scale 3D carbon aerogels (Fig. 1b). This novel approach enables the super assembly on a multi-dimensional scale, which realizes the macroscopic application of nanoparticles and the functionalization of CAs. In order to verify that the performance of C-ZIF-8-CNFAs is improved after the introduction of ZIF-8, pure CNFAs (Fig. S2^†) without ZIF-8 (the synthetic process is shown in the ESI^†) are also prepared.Open in a separate windowFig. 1Preparation steps for C-ZIF-8-CNFAs. (a) More fabrication details for C-ZIF-8-CNFAs. (b) Schematic illustration of the fabrication of CNFAs.As illustrated in Fig. 2a, when the ZIF-8/nanofiber dispersion solution is frozen, the uniformly dispersed ZIF-8/nanofibers are extruded by the growth of ice crystals and assembled among ice crystals. After the sample is frozen completely, the nanofibers become lapped and locked into a 3D nanofibrous network. Subsequently, the ZIF-8/NFAs with a luffa-like structure are obtained after the sublimation of ice crystals through the freeze-drying process.²² Moreover, ZIF-8/NFAs can be made into diverse desired shapes such as cylinders, cubes, moon-like shapes, star-like shapes, heart-like shapes and intricate shapes of the letters (Fig. 2b). Fig. 2c shows the obvious reduction of the intensity of ZIF-8/NFA XRD patterns compared to the original ZIF-8, but the site hardly changes, which confirms that the introduced ZIF-8 nanoparticles are not destroyed during the electrospinning and the preparation process of aerogels. After the preoxidation and carbonization steps, the typical C-ZIF-8-CNFAs with an ultra-low density of 7.32 mg cm⁻³ can freely stand on the tip of a red maple leaf (Fig. 2b). The scanning electron microscopy (SEM) images in Fig. 2d–f show that C-ZIF-8-CNFAs have the hierarchical porous luffa-like structure with three kinds of pores (the picture of an actual luffa shown in the inset of Fig. 2d). The porous structure exhibits obvious rectangular pores of ∼25 μm, and the wall of these pores is made of interconnecting nanofibers (Fig. 2e). Meanwhile, the secondary pores of ∼1.5 μm are formed by the welded nanofibers that are interconnected with each other. The nanoscale pores of ∼200 nm also exist in these nanofibers, which come from the carbonization of ZIF-8 nanoparticles (Fig. 2f). As observed from the magnified SEM image (Fig. S3^†), the welded structure resulting from preoxdiation is still preserved through carbonization. Remarkably, ZIF-8 nanoparticles encapsulated by PAN and PVP are transformed into a hollow structure after carbonization. As can be seen from the TEM image (Fig. 2g), the ZIF-8-derived hollow structure is evenly dispersed in PAN/PVP-derived carbon nanofibers. Because Zn²⁺ ions coordinate with –C N groups existing on the surface of PAN/PVP nanofibers, the ZIF-8 particles become tightly encapsulated by the PAN/PVP layers. During carbonization, the PAN/PVP layers make ZIF-8 shrink from inside to outside, thereby leading to the generation of a hollow structure. Moreover, the confined carbonization process within the PAN/PVP matrix prevents the irreversible fusion and aggregation of carbonized ZIF-8 nanoparticles.¹⁹Open in a separate windowFig. 2(a) Schematical illustration of the formation principles for the hierarchical cellular structure. (b) Photographs of ZIF-8/NFAs with diverse shapes and the lightweight C-ZIF-8-CNFAs standing on the tip of a red maple leaf. (c) Wide-angle XRD patterns. (d–g) SEM (d–f) and TEM (g) images showing the microstructure of C-ZIF-8-CNFAs at various magnifications. (h and i) N₂ adsorption–desorption isotherm and pore-size distribution curve of CNFAs and C-ZIF-8-CNFAs.After ZIF-8 nanoparticles were introduced into CNFAs, the properties have been improved significantly. To comprehend the variation of the porous characteristics in C-ZIF-8-CNFAs, nitrogen (N₂) adsorption–desorption measurements were carried out. As observed from Fig. 2h and i, C-ZIF-8-CNFAs have a larger specific surface area of 288.3 m² g⁻¹ and pore volume of 0.22 cm⁻³ g⁻¹, while the specific surface area and pore volume of CNFAs are only 12.1 m² g⁻¹ and 0.01 cm⁻³ g⁻¹, respectively. Meanwhile, C-ZIF-8-CNFAs also have a hierarchical porous structure with micropores, mesopores, and macropores (Fig. 2i). Because of the existence of the hierarchical porous structure in C-ZIF-8-CNFAs, they have a lower density relative to CNFAs (20.73 mg cm⁻³) (Table S1^†). The chemical compositions and graphitic structure of C-ZIF-8-CNFAs are investigated by X-ray powder diffraction (XRD) and X-ray photoelectron spectroscopy (XPS). The XRD pattern (Fig. S4a^†) of C-ZIF-8-CNFAs only exhibits two broad peaks at about 25° and 44°, corresponding to the (002) and (101) diffraction facets of the graphitic structure, respectively.²³ The correlative XPS spectrum shows that C-ZIF-8-CNFAs consist of C, N and O (Fig. S4b^†). The high-resolution N 1s spectra (Fig. S4c^†) can be deconvoluted into four peaks: pyridinic N (398.6 eV), pyrrolic N (399.4 eV), graphitic N (400.9 eV), and oxidized N (403.7 eV), respectively. The percentage of nitrogen and its types are listed in Table S2.^† There is no Zn content thus indicating that it had evaporated during the carbonization process at high temperature (at 900 °C for 3 hours).²⁴ Furthermore the N content of C-ZIF-8-CNFAs increases with the introduction of ZIF-8 compared to CNFAs clearly (Table S3^†).In stark contrast to the hard and brittle characteristics of conventional carbon aerogels, C-ZIF-8-CNFAs show robust mechanical properties, sustaining large compressive strain without fracture (insets in Fig. 3a and movie S1^†). The compressive stress–strain (σ–ε) curves (Fig. 3a) show the compressive process of C-ZIF-8-CNFAs and two typical deformation regimes could be recognized: a Hookean or linear elastic regime of ε < 50% with a stable tangent modulus, and a densification regime of ε > 50% with σ and dσ/dε increasing sharply. When the maximum compressive strain increases from 30% to 80%, the maximum compressive stress steeply increases from 1.2 to 25.9 kPa, indicating that C-ZIF-8-CNFAs can bear over 25 000 times their own weight without cracking. Moreover, the cyclic compression test of the C-ZIF-8-CNFAs is performed to validate their durable cycling performance by applying 50 loading–unloading fatigue cycles at a large ε of 50% (Fig. 3b). C-ZIF-8-CNFAs hardly undergo plastic deformation at all after the 50th cycle, which is a huge improvement over traditional CAs with a crisp character. As shown in Fig. 3c, C-ZIF-8-CNFAs retain nearly 100% of the initial value of the Young''s modulus, maximum stress and energy loss coefficient, indicating that their strength or stiffness has no significant decline highlighting their exceptional structural robustness. Two reasons could account for the excellent mechanical elasticity of C-ZIF-8-CNFAs. On the one hand, the compressive strain is absorbed by the bending of the connected ZIF-8/nanofibers between layers; on the other hand, the strain that continues to grow (beyond 50%) is absorbed by the densification of minor pores, which is formed by the welding of PVP (Fig. 3d).²⁵ The two reasons can also be confirmed by the sharp increase of the compressive stress after ε > 50% (Fig. 3a).Open in a separate windowFig. 3(a) Compressive stress–strain curve of C-ZIF-8-CNFAs at different strains. The insets are photographs of C-ZIF-8-CNFAs under a compressing and releasing cycle (ε = 80%). (b) A 50-cycle compressive fatigue test with ε of 50%. (c) The Young^''s modulus, energy loss coefficient, and maximum stress versus compressive cycles. (d) Sketch of the changes in the hierarchical porous structure with compressive deformation.MOF-based derived CNFAs with hollow structures, ultra-low density, extraordinary mechanical elasticity, and large surface area are conducive to applications in various fields, including environment governance, energy absorption, and energy storage.^26–28 As a proof of concept, we evaluated the organic solvent absorption and photo-thermal conversion performance of C-ZIF-8-CNFAs. Excellent hydrophobicity is one of the important factors to ensure the absorption efficiency of organic solvents. As demonstrated in Fig. 4a, C-ZIF-8-CNFAs exhibit high hydrophobicity with a water contact angle of 142° and the water contact angle has no significant change after 120 s. The water droplet and the absorbed oil droplet are on the surface of C-ZIF-8-CNFAs, indicating the effective oil/water selectivity of C-ZIF-8-CNFAs (Fig. S5^†). As shown in Fig. 4b, C-ZIF-8-CNFAs demonstrate extraordinary absorption capacities for common oils and various organic solvents, 90–200 times their own weight, principally depending on the density and viscosity of the solvents. For example, the adsorption capacity for carbon tetrachloride (ρ = 1.595 g cm⁻³) is much higher than the adsorption capacity for n-hexane (ρ = 0.66 g cm⁻³). We also compared the maximum adsorption capacity of C-ZIF-8-CNFAs and CNFAs for several common organic solvents and oils. The adsorption performance of C-ZIF-8-CNFAs is also demonstrated to be better than that of CNFAs (Fig. S6^†). According to the cross-sectional diagram of the adsorption of organic solvents (Fig. 4c), CNFAs only absorb organic solvents through the capillary phenomenon of the channel, which is formed by nanofibers, while the absorbed organic solvent can also enter into the hollow structure formed by carbonization of ZIF-8 inside the nanofiber of C-ZIF-8-CNFAs. Therefore, the performance of C-ZIF-8-CNFAs is significantly improved after the addition of ZIF-8 nanoparticles. Moreover, the adsorption capacity of C-ZIF-8-CNFAs for organic solvents is greater than that of previously reported aerogels (Fig. S7^†). To verify the excellent adsorption performance of C-ZIF-8-CNFAs in practical application, the oil/water separation and cyclic adsorption experiments were carried out. As illustrated in Fig. S8,^† C-ZIF-8-CNFAs can quickly absorb heavy organic solvents such as carbon tetrachloride (dyed with oil red) sunk at the bottom of water, thus indicating their potential application for selectively removing oils from water. Recyclability and reusability are also crucial to evaluate the practical application possibility of adsorption materials. In view of the outstanding elasticity and structural robustness of C-ZIF-8-CNFAs, we chose ethanol as the absorption solvent for recycling tests. Through simple heating, ethanol absorbed by C-ZIF-8-CNFAs can be readily removed. As demonstrated in Fig. S9a and b,^† even after 10 cycles, C-ZIF-8-CNFAs still retained over 90% adsorption capacity for ethanol and their mass was reduced by less than 10%.Open in a separate windowFig. 4(a) Dynamic behaviors of a water droplet on the surface of C-ZIF-8-CNFAs. (b) Adsorption efficiency of C-ZIF-8-CNFAs towards commonly used organic solvents and oils. (c) Schematic illustration of the organic solvent adsorption process inside the CNFAs and C-ZIF-8-CNFAs. (d) Mass changes of evaporated water versus time under 1 sun illumination. (e) The temperature of C-ZIF-8-CNFAs and seawater under 1 sun illumination as a function of irradiation time. (f) The infrared images (IR) show the temperature distribution of C-ZIF-8-CNFAs and seawater under 1 sun illumination with an irradiation time of 0, 30, and 60 min.Given the abundant porous structure and blackbody characteristic, C-ZIF-8-CNFAs are also promising materials for interfacial solar steam generation (ISSG). To investigate the ISSG performance, the evaporation mass change of seawater and C-ZIF-8-CNFAs is measured under one sun illumination. As shown in Fig. 4d, the seawater in C-ZIF-8-CNFAs achieves the maximum evaporation rate of 3.74 kg m⁻² h⁻¹, which is 5.12 times the evaporation rate of bulk seawater (0.73 kg m⁻² h⁻¹) and 2.34 times faster than the evaporation rate of general 2D ISSG (∼1.6 kg m⁻² h⁻¹). An infrared camera is used to trace the surface temperatures of C-ZIF-8-CNFAs and bulk seawater under one-sun illumination to evaluate the photothermal behavior of ISSG (Fig. 4e). The surface temperature of C-ZIF-8-CNFAs presents a quick increase in 5 min and eventually reaches a stable state (∼48 °C) after 10 min, while the bulk seawater temperature stays unchanged (∼25 °C). The consecutive infrared images in Fig. 4f show the equilibrium temperature distribution and heat localization effect of C-ZIF-8-CNFAs and bulk seawater in 60 min under one-sun illumination.     

Tuning through-space interactions via the secondary coordination sphere of an artificial metalloenzyme leads to enhanced Rh(iii)-catalysis

We report computationally-guided protein engineering of monomeric streptavidin Rh(iii) artificial metalloenzyme to enhance catalysis of the enantioselective coupling of acrylamide hydroxamate esters and styrenes. Increased TON correlates with calculated distances between the Rh(iii) metal and surrounding residues, underscoring an artificial metalloenzyme''s propensity for additional control in metal-catalyzed transformations by through-space interactions.

We report computationally-guided protein engineering of monomeric streptavidin Rh(iii) artificial metalloenzyme to enhance catalysis of the enantioselective coupling of acrylamide hydroxamate esters and styrenes.

Artificial metalloenzymes (ArMs) can be made by anchoring a non-natural (metal) cofactor into a protein scaffold, with the goal of imbuing new-to-nature reactivity.¹ One of the most common ArM platforms is the biotin-tetrameric(strept)avidin (biotin-tSav) system pioneered by Whitesides and Ward.^2,3 These ArMs utilize high-affinity (up to K_D ∼10⁻¹⁴ M) interactions between tSav and biotin–metal conjugates. tSav-based ArMs have appeared in an increasing number of transition-metal catalyzed transformations.^4–6 In collaboration with the Ward group, we have previously described a tetrameric streptavidin (tSav) system containing a biotinylated Rh(iii) cofactor for the asymmetric synthesis of dihydroisoquinolones using benzhydroxamate esters and acrylate partners.⁷ Monomeric streptavidin (mSav), a streptavidin/rhizavidin hybrid designed to resist tetramerization, retains its high affinity for biotin (K_D ∼10⁻⁹ M).^8,9 We recently described the use of mSav as a new ArM,¹⁰ whose simpler topology encourages protein engineering via a site-directed mutagenesis approach.Traditional manipulation of a metal''s reactivity has been accomplished by modification of the electronic and steric properties of the bound ligands (Fig. 1a).^11,12 For example, we have documented and parsed the impact of Cp electronics and sterics on a number of Rh(iii) catalyzed transformations, by structural changes to the ligand in the primary coordination sphere of Rh.¹³ On the other hand, ArMs have traditionally been used as modifiers of a metal''s steric environment largely focusing on inducing asymmetry in the bond-forming events. Less broadly appreciated is the fact that any mutations in residues proximal to the active site may also impact the metal''s electronic properties via changes to the secondary coordination sphere (Fig. 1b), with the prospect of delivering more active catalysts for a given transformation.Open in a separate windowFig. 1Methods to modify the (a) primary and (b) secondary coordination sphere of a Rh(iii) catalyst.Previously, we described a mSav·Rh(iii) catalyst and demonstrated its use in the direct enantioselective coupling of acrylamide hydroxamate esters and styrenes.¹⁰ The reaction allows rapid access to piperidines – the most common N-heterocycle found in FDA-approved pharmaceuticals.¹⁴ One of the most interesting aspects of this reaction was our observation of a 7-fold increase in turnover number (TON) by embedding the cofactor into mSav''s active site.¹⁵ It has been a long-standing goal of ArMs to not only enable new-to-nature reactivity, but also for them to achieve the stellar kinetics of a native metalloenzyme. As these systems lack the evolutionary privilege of a natural metalloenzyme, extensive mutation of the protein scaffold may be required to find the optimal environment of the metal cofactor.Predicting the effects of specific mutations can prove very challenging, as any alterations to the protein conformation and charge distribution can impact reactivity regardless of the mutation''s distance from the active site.^16–19 In order to design a better mutant, we embarked on a collaborative experimental and computational study to define the role of the protein scaffold and how single point mutations affect reactivity. We identified two key residues that play a pivotal role in mSav·Rh(iii) ArM''s secondary coordination sphere, and have used this insight to design a more active mutant.For the purposes of this study, we focused on the mSav·Rh(iii) ArM-catalyzed coupling of methacrylamide with 4-methoxystyrene as our model reaction (Fig. 2a). Using a small model of the catalyst, the lowest energy pathway of this reaction''s proposed mechanism was generated (Fig. S9^†). The calculations were performed in Turbomole^20–32 with the M06 density functional.³³ Geometries were optimized with the def2-SVP basis set, and final electronic energies were calculated with the def2-TZVP basis set.³⁴ The conductor-like screening model (COSMO)³⁵ was used as implicit solvent with a dielectric of 80 to simulate water. These calculations predicted similar barriers for the N–H activation, the C–H activation, and the migratory insertion (differences less than 3 kcal mol⁻¹). Isotope-exchange experiments revealed that the C–H activation step is reversible, implicating the migratory insertion step as turnover-limiting.¹⁰Open in a separate windowFig. 2(a) Model transformation. (b) Snapshot of the transition state for alkene insertion illustrating key nearby residues Y112 (red), E124 (blue), and S119 (purple). (c) Computed barrier to alkene insertion in the presence and absence of phenol and acetate (shown in blue).The Cp* moiety of the Cp*^biotinRhX₂ cofactor is non-covalently localized in the active site likely due to a π–π stacking interaction with Y112 (Fig. 2b). This assignment is supported by the observation that mutant Y112A leads to lower yield and enantioselectivity.¹⁰ We hypothesized that we could further manipulate both the sterics and electronics of the Cp* moiety by either directly mutating Y112 or indirectly by mutating other residues that affect the Y112-Cp* interaction.To generate a model of mSav''s protein scaffold and active-site we used QM/DMD³⁶ – a hybrid quantum mechanics/molecular mechanics method that simulates proteins piecewise. Discrete molecular dynamics (DMD) equilibrates the entire system except for the metal and part of the substrate.³⁶ After a trajectory of ∼0.5 ns, quantum mechanics (QM) is used to optimize the metal region plus sidechains and residues immediately surrounding it. This process is repeated, providing efficient sampling of the entire protein scaffold while treating the metal environment quantum-mechanically. For this study, the migratory insertion transition state was modeled in WT by freezing the coordinates of the rhodium atom and the two carbon atoms forming a bond. For each system, five replicate simulations were run for ∼20 ns each.Residues E124 and S119 both hydrogen bond to Y112 and are in close proximity to the RhCp* catalytic site (Fig. 2b).³⁷ To estimate the electronic effects of these three residues on the reaction, an acetate ion, methanol molecule, and 4-methylphenol (p-cresol) molecule were added to a small catalyst model without constraints but initially positioned to mimic the sidechains of these residues (Fig. 2c). The migratory insertion energy barrier decreases by 2 kcal mol⁻¹ with incorporation of the three residues. However, this energy barrier decreases by an additional 3 kcal mol⁻¹ upon the deletion of the methanol molecule representing S119. Not only does this imply that these amino acid sidechains may be the primary reason for the increased activity of the protein-installed catalyst, but also suggest that a longer Y112–S119 distance is favorable, so long as no water can insert in this region and replace S119 in its H-bond with Y112. We hypothesize that the carboxylate group of E124 acts as a hydrogen bond acceptor, donating electron density to the Y112 phenol ring, which in turn donates electron density to the catalyst via π–π charge transfer. This could enhance the electron donation of the metal and decrease the energy barrier to the migratory insertion step. On the other hand, S119 acts as a hydrogen bond donor which would remove electron density from Y112 and subsequently the Rh(iii) moiety.Unfortunately, mutation of Y112 (Y112F and Y112W) results in negligible protein yields. We thus identified three flanking residues (T111, E113, H87) that may be expected to have a significant impact on Y112''s position, and one distal (T32) residue, chosen as distal mutations sometimes have significant impact (Fig. 3). Through this subset of mutants, we attempted to increase TON and establish a correlation between the Y112–Rh distance and Y112–S119 distance of the mutants and their reactivities.Open in a separate windowFig. 3Structure of mSav from two different views highlighting some of the mutated residues including their TON and enantioselectivity.We used QM/DMD to simulate a representative set of these mutants spanning a wide range of TONs measured in the experiment. The Y112–Rh and Y112–S119 distances were measured every ∼0.5 ps for every simulation. The results can be represented by a 3-dimensional plot with Y112–Rh distance on the X axis, Y112–S119 distance on the Y axis, and probability density on the Z axis (Fig. 4). We find the best correlation between TON and probability density in the conformational region where the Y112–Rh distance is the shortest and the Y112–S119 interaction is not energetically relevant.³⁸Open in a separate windowFig. 4Three-dimensional probability distributions from select mutants by simultaneous sampling of Rh–Y112 and S119–Y112 distances. Probabilities for the outlined regions are also shown.To clarify this correlation, we calculated the probability of having a Y112–S119 distance between 3.5–6 Å and a Y112–Rh distance less than 5.65 Å. This Y112–S119 distance corresponds to negligible hydrogen bonding.³⁹ Additionally, we constrained the small model catalyst shown in Fig. 5b (ref. 40) and calculated the corresponding energy barriers at different Y112–Rh distances (Fig. 5a). Since rate increases exponentially as the barrier decreases,³⁵ differences in probabilities in the region where the Y112–Rh distance is between 5.4–5.65 Å have the greatest impact on the relative TONs of our model methacrylamide styrene coupling. We conclude that mutants with increasing probability in this region provide increasing TON.Open in a separate windowFig. 5a) Theoretical dependence of migratory insertion barrier on Rh-phenol distance. (b) Small-model catalyst.Theoretically, a shorter Y112–Rh distance relative to WT would result in increased reactivity. Residue G49 is located under the Rh(iii) moiety (Fig. 6). We hypothesize that by mutating the glycine into an alanine, steric congestion would force the biotinylated Rh(iii) cofactor to shift upwards closer to the electron donating phenol side chain of residue Y112. Analyzing the critical portions of the Y112–Rh and Y112–S119 distances in tandem reveals that G49A has the highest probability density in this region (Fig. 4). Indeed, experimentally, this mutant gives 97 TON and 91% ee (Fig. 6). The combination of a short Y112–Rh distance and long Y112–S119 distance leads to an increase in reactivity. This is an approximate 3-fold improvement in the TON relative to WT. The G49A mutant serves as an experimental proof of concept that a computational analysis of an ArMs secondary coordination sphere can lead to the design of a more efficient ArM.Open in a separate windowFig. 6Snapshot of the transition state for alkene insertion highlighting the position of G49 (purple) relative to Rh. Y112 is shown in red and E124 is shown in blue.In summary, we have identified three key residues that contribute to accelerating the rate of a Rh(iii)-catalyzed reaction by electronic communication to the metal via the secondary coordination sphere. E124 hydrogen bonds to Y112 transferring electron density via π–π charge transfer, an effect that is attenuated by hydrogen bonding from S119. Optimal interaction of these residues can be described computationally by finding mutants that have multiple conformations bearing short Y112–Rh distances coupled with negligible bonding between Y112 and S119. This hypothesis was experimentally verified by a mutant that enforces a closer Y112–Rh distance leading to improved TON. This result demonstrates the use of a hypothesis-based site-directed mutagenesis of the secondary sphere residues, to optimize the metal''s electronic environment within the protein scaffold and enhance an ArM''s activity.     

HydroFlipper membrane tension probes: imaging membrane hydration and mechanical compression simultaneously in living cells

HydroFlippers are introduced as the first fluorescent membrane tension probes that report simultaneously on membrane compression and hydration. The probe design is centered around a sensing cycle that couples the mechanical planarization of twisted push–pull fluorophores with the dynamic covalent hydration of their exocyclic acceptor. In FLIM images of living cells, tension-induced deplanarization is reported as a decrease in fluorescence lifetime of the dehydrated mechanophore. Membrane hydration is reported as the ratio of the photon counts associated to the hydrated and dehydrated mechanophores in reconvoluted lifetime frequency histograms. Trends for tension-induced decompression and hydration of cellular membranes of interest (MOIs) covering plasma membrane, lysosomes, mitochondria, ER, and Golgi are found not to be the same. Tension-induced changes in mechanical compression are rather independent of the nature of the MOI, while the responsiveness to changes in hydration are highly dependent on the intrinsic order of the MOI. These results confirm the mechanical planarization of push–pull probes in the ground state as most robust mechanism to routinely image membrane tension in living cells, while the availability of simultaneous information on membrane hydration will open new perspectives in mechanobiology.

HydroFlippers respond to membrane compression and hydration in the same fluorescence lifetime imaging microscopy histogram: the responses do not correlate.

The detection and study of membrane mechanics in living cells is a topic of current concern.^1–14 To enable this research, appropriate chemistry tools, that is small-molecule fluorescent probes that allow imaging of membrane tension, are needed.¹⁵ With the direct imaging of physical forces being intrinsically impossible, design strategies toward such probes have to focus on the suprastructural changes caused by changes in membrane tension.¹⁵ These suprastructural changes are divers, often interconnected, and vary with the composition of the membrane.^15–25 Beyond the fundamental lipid compression and decompression, they include changes in membrane curvature, from rippling, buckling and budding to tubules extending from the membrane and excess lipid being ejected. Of similar importance are changes in membrane organization, particularly tension-induced phase separation and mixing, i.e. assembly and disassembly of microdomains. Consequences of these suprastructural changes include microdomain strengthening and softening and changes in membrane hydration and viscosity.^16–25The currently most developed fluorescent flipper probes have been introduced^26,27 to image membrane tension by responding to a combination of mechanical compression and microdomain assembly in equilibrium in the ground state.¹⁵ Extensive studies, including computational simulations,²⁸ have shown that flipper probes align non-invasively along the lipid tails of one leaflet and report changes in membrane order and tension as changes in fluorescent lifetimes and shifts of excitation maxima.¹⁵ Among other candidates, solvatochromic probes respond off-equilibrium in the excited state to changes in membrane hydration and have very recently been considered for the imaging of membrane tension in living cells.^29–36 So far not considered to image tension, ESIPT probes also report off equilibrium in the excited state on membrane hydration, but for different reasons.^37,38 Mechanosensitive molecular rotors respond off equilibrium in the excited state to changes in microviscosity.^{17,30,32,39–53} The same principle holds for the planarization of bent, papillon or flapping fluorophores.^54–57 The response of all possible probes to tension can further include less desired changes in positioning and partitioning between different domains, not to speak of more catastrophic probe aggregation, precipitation, disturbance of the surrounding membrane structure, and so on. Although the imaging of membrane tension is conceivable in principle with most of above approaches, the complex combination of parameters that has to be in place can thus far only be identified empirically, followed by much optimization.¹⁵The force-induced suprastructural changes are accompanied by the alteration in several unrelated physical properties of membranes. It is, for instance, well documented that membrane hydration increases with membrane disorder, from solid-ordered (S_o) to liquid-disordered (L_d) phases.^29,58 Increasing cholesterol content decreases membrane hydration in solid- and liquid-ordered membranes.⁵⁹ However, studies in model membranes also indicate that membrane hydration and membrane fluidity do not necessarily correlate.⁵⁹ The dissection of the individual parameters contributing to the response of fluorescent membrane tension probes would be important for probe design and understanding of their responses, but it remains a daunting challenge. In this study, we introduce fluorescent flipper probes that simultaneously report on mechanical membrane compression and membrane hydration at equilibrium in the ground state. Changes of both in response to changes in membrane tension and membrane composition are determined in various organelles in living cells.The dual hydration and membrane tension probes are referred to as HydroFlippers to highlight the newly added responsiveness to membrane hydration. The mechanosensing of lipid compression in bilayer membranes by flipper probes has been explored extensively.¹⁵ Fluorescent flippers²⁷ like 1 are designed as bioinspired⁶⁰ planarizable push–pull probes²⁶ (Fig. 1). They are constructed from two dithienothiophene fluorophores that are twisted out of co-planarity by repulsion of methyls and σ holes on sulfurs^61,62 next to the twistable bond. The push–pull system is constructed first from formal sulfide and sulfone redox bridges in the two twisted dithienothiophenes. These endocyclic donors and acceptors are supported by exocyclic ones, here a trifluoroketone acceptor and a triazole donor.⁶³ To assure stability, these endo- and exocyclic donors are turned off in the twisted ground state because of chalcogen bonding and repulsion, respectively.⁶²Open in a separate windowFig. 1The dual sensing cycle of HydroFlippers 1–5, made to target the indicated MOIs in living cells and responding to membrane compression by planarization and to membrane hydration by dynamic covalent ketone hydration. With indication of excitation maxima (ref. 63) and fluorescence lifetimes (this study).Mechanical planarization of the flipper probe establishes conjugation along the push–pull systems, electrons flow from endocyclic donors to acceptors, which turns on the exocyclic donors and acceptors to finalize the push–pull system.⁶² This elaborate, chalcogen-bonding cascade switch has been described elsewhere in detail, including high-level computational simulations.⁶² The planar high-energy conformer 1dp excels with red shifted excitation and increased quantum yield and lifetime compared to the twisted conformer 1dt because the less twisted Franck-Condon state favors emission through planar intramolecular charge transfer (PICT) over non-radiative decay through twisted ICT, or conical intersections.¹⁵Flipper probe 1 was considered for dual responsiveness to membrane tension and hydration because of the trifluoroketone acceptor.⁶³ Dynamic covalent hydration of 1dt yields hydrate 1ht.^64–76 Blue-shifted excitation and short lifetime of 1ht are not expected to improve much upon planarization because the hydrate is a poor acceptor and thus, the push–pull system in 1hp is weak. The dynamic covalent chemistry of the trifluoroketone acceptor has been characterized in detail in solution and in lipid bilayer membranes.⁶³To explore dual responsiveness to membrane tension in any membrane of interest (MOI) in living cells, HydroFlippers 2–5 were synthesized. While HydroFlipper 1 targets the plasma membrane (PM), HydroFlippers 2–4 were equipped with empirical targeting motifs.⁷⁷ HydroFlipper 5 terminates with a chloroalkane to react with the self-labeling HaloTag protein, which can be expressed in essentially any MOI.⁷⁸ Their substantial multistep synthesis was realized by adapting reported procedures (Schemes S1–S4^†).The MOIs labeling selectivity of HydroFlippers was determined in HeLa Kyoto (HK) cells by confocal laser scanning microscopy. Co-localization experiments of flippers 1–4 with the corresponding trackers gave Pearson correlation coefficients (PCCs) >0.80 for the targeting of mitochondria, lysosomes and the endoplasmic reticulum (ER, Fig. S4–S6^†). HydroFlipper 5 was first tested with stable HGM cells, which express both HaloTag and GFP on mitochondria (referred to as 5_M).^78,79 The well-established chloroalkane penetration assay demonstrated the efficient labeling of HaloTag protein by 5 as previously reported HaloFlippers (Fig. S3^†).⁷⁸ By transient transfection, HydroFlippers 5 were also directed to lysosomes (5_L), Golgi apparatus (GA, 5_G)⁸⁰ and peroxisomes (5_P) with HaloTag and GFP expressed on their surface.⁷⁸ PCCs >0.80 for co-localization of flipper and GFP emission confirmed that MOI labeling with genetically engineered cells was as efficient as with empirical trackers (Fig. S7–S11^†).Dual imaging of membrane compression and hydration was envisioned by analysis of fluorescence lifetime imaging microscopy (FLIM) images using a triexponential model (Fig. 2).⁸¹ FLIM images of ER HydroFlipper 4 in iso-osmotic HK cells were selected to illustrate the concept (Fig. 3a). Contrary to classical flipper probes, the fluorescence decay curve of the total FLIM image (Fig. 2a, grey) showed a poor fit to a biexponential model (Fig. 2a, cyan, b). Consistent with their expected dual sensing mode, a triexponential fit was excellent (Fig. 2a, dark blue, c). Lifetimes τ_1i = 4.3 ns (†) were obtained besides background. This three-component model was then applied to every pixel of FLIM images (Fig. 3c). The resulting reconvoluted FLIM histogram revealed three clearly separated populations for τ₁ (red), τ₂ (green), and background (τ₃, blue, Fig. 2d). Maxima of these three clear peaks were at the lifetimes estimated by triexponential fit of the global decay curve, thus demonstrating the validity of the methodology at necessarily small photon counts. Irreproducible fitting would give randomly scattered data without separated peaks.Open in a separate windowFig. 2(a) Fluorescence decay curve (grey, corresponding to the total image, not to a single pixel) with biexponential (cyan) and triexponential fit (dark blue). (b, c) Residual plots for bi- (b) and triexponential fit (c). (d) Histogram with the intensities associated with the τ₁ (red), τ₂ (green), and τ₃ (blue, background) components obtained by triexponential fit of the fluorescence decay curve of each pixel of the FLIM image, fit to Gaussian function (black solid curves).Open in a separate windowFig. 3FLIM images of HK cells labelled with ER flipper 4 before (a, c) and after (b, d) hyper-osmotic shock, showing average lifetimes τ_av (a, b) and τ₁ (c, d) from triexponential reconvolution; scale bars = 10 μm. (e) Distribution of the photon counts associated with the τ₁ component of 4 in HK cells after triexponential reconvolution of FLIM images before (c, τ_1i) and after (d, τ_1h) hyper-osmotic shock, showing decreasing lifetimes for τ₁ (4d). (f) The dehydration factor dh_i defined as total integrated photon counts for τ₁ (Στ₁) divided by Στ₂ (i.e., dh_i = area Στ_1i/area Στ_2i) for 4 in strongly hydrated ER (dh_i < 2, turquoise) and 1 in weakly hydrated plasma membrane (dh_i > 6, purple) of HK Kyoto cells under iso-osmotic conditions.Dual response of HydroFlippers to changes in membrane tension^a

Strong non-Arrhenius behavior at low temperatures in the OH + HCl → H2O + Cl reaction due to resonance induced quantum tunneling

The OH + HCl → H₂O + Cl reaction releases Cl atoms, which can catalyze the ozone destruction reaction in the stratosphere. The measured rate coefficients for the reaction deviate substantially from the Arrhenius limit at low temperatures and become essentially independent of temperature when T < 250 K, apparently due to quantum tunneling; however, the nature of the quantum tunneling is unknown. Here, we report a time-dependent wave packet study of the reactions on two newly constructed potential energy surfaces. It is found that the OH + HCl reaction possesses many Feshbach resonances trapped in a bending/torsion excited vibrational adiabatic potential well in the entrance channel due to hydrogen bond interaction. These resonance states greatly induce quantum tunneling of a hydrogen atom through the reaction barrier, causing the reaction rates to deviate substantially from Arrhenius behavior at low temperature, as observed experimentally.

The OH + HCl reaction possesses many Feshbach resonances trapped in the hydrogen bond well in the entrance channel, which substantially enhance the reaction rates at low temperatures.

In the classical picture, a chemical reaction with an energetic barrier can only occur at collision energies higher than the barrier, which leads to the well-known Arrhenius formula for chemical reaction rates. However, chemical reactions can happen at energies below the reaction barrier through quantum tunneling,^1–3 resulting in the deviation of the reaction rates from the Arrhenius behavior at low temperatures. The effect of quantum tunneling on the reaction rates increases with decreasing reaction temperature, and hence becomes especially important in low-temperature environments such as the interstellar medium and atmospheric processes.^4,5 Reaction resonances are quasi-trapped quantum states in the transition state region with some lifetime, and can substantially promote quantum tunneling through the reaction barrier. Over the past decades, great efforts have been devoted to detecting resonances in chemical reactions and to studying their structures and dynamics.^6–14 Theoretical studies on the O(³P) + HCl → OH + H reaction using the accurate ³A′′ potential energy surface revealed that the tunneling induced by the resonances trapped in the van der Waals well in the reactant channel can substantially enhance the thermal rate constants at low temperatures.^15–17 A combined experimental and theoretical investigation discovered that resonance-induced quantum tunneling dramatically enhances the reactivity of the F + p-H₂ → HF + H reaction,¹⁴ and is fully responsible for the unusually high chemical reactivity of the reaction in the low temperature interstellar medium. Recently, a quantum dynamics calculation showed that the presence of two resonance peaks strongly influence the rotational quenching of HF (j = 1, 2) with H, leading to an up to two-fold increase in the thermal rate coefficients at the low temperatures characteristic of the interstellar medium.¹⁸ Thus, understanding quantum tunneling, and in particular, resonance-induced quantum tunneling, in chemical reactions is of general interest and fundamental importance to low-temperature chemistry.The OH + HCl → H₂O + Cl reaction is of great importance in atmospheric chemistry because it releases Cl atoms from one of the principal chlorine-containing species in the stratosphere, HCl. The Cl atoms generated from the reaction can catalyze the ozone destruction reaction in the stratosphere, which was responsible for the formation of the ozone hole over Antarctica.¹⁹ Since the chlorine-catalyzed ozone destruction is proportional to the steady-state Cl atom concentration, which is directly controlled by the rate of the reaction, extensive studies have been carried out to measure the rate coefficient for the reaction with high accuracy over a wide range of temperatures.^20–29 The measured rate coefficients exhibit a small activation energy of a few hundred K, deviate substantially from the Arrhenius limit at low temperatures and become essentially independent of temperature when T < 250 K.^25,28,29 In addition, a large H/D kinetic isotope effect has also been found.^21,22,26,28 All these observations suggest the presence of an important quantum tunneling effect in the reaction.The dynamics of this reaction and its reverse have also attracted great attention in the past decades. In particular, the endothermic Cl + H₂O → HCl + OH reaction with a late barrier has been extensively investigated as a benchmark system for mode specificity and bond selectivity chemistry.^30–35 Recently, the construction of two high-quality potential energy surfaces (PESs) in its ground electronic state using the PIP-NN method have substantially advanced the theoretical study of the dynamics and kinetics of the system. The first PES was based on a large number of ab initio data points calculated at the multi-reference configuration interaction (MRCI) level of theory by Li, Dawes, and Guo (LDG),³³ and the second one was fitted to ab initio energy points obtained using an explicitly correlated unrestricted coupled-cluster method with single, double, and perturbative triple excitations (UCCSD(T)-F12b) and the augmented correlation-consistent polarized valence triple-zeta (aug-cc-pVTZ, or AVTZ) basis set by Zuo, Zhao, Guo, and Xie (ZZGX) with a fitting error of 6.9 meV.³⁶ The static barrier height is 2.86 and 2.23 kcal mol⁻¹, respectively, for these two PESs. On both PESs, there exists a well of about 3.5 kcal mol⁻¹ in the OH + HCl entrance channel due to the hydrogen bond (HB) interaction between OH and HCl. Extensive quantum dynamics studies on the PESs have revealed many interesting features of the reaction in both directions. In particular, time-dependent wave packet calculations for the title reaction with OH in the ground and vibrational excited states found one or two broad peaks in the total reaction probabilities, which are presumed to be the signature of the resonances supported by the reactant complex well.^36,37 However, the impact of these peaks on the reaction rates has not been investigated. Ring-polymer molecular dynamics (RPMD) calculations^38,39 were also carried out on both PESs to compute the thermal rate coefficients for the reaction.^40–42 It was found that the RPMD rates on the LDG PES underestimate the experimental data, while the RPMD rates on the latest ZZGX PES agree with the experimental results much better, and do not decrease further when the temperature falls below 300 K, apparently due to quantum tunneling. Unfortunately, RPMD calculations cannot provide any clue regarding the nature of the quantum tunneling.Therefore, despite the significant progress that has been made in theoretical studies of the system, some key issues still remain to be addressed: Are the broad peaks in the total reaction probabilities obtained from the time-dependent wave packet calculations indeed the signature of the resonances in the reaction? If not, do there exist resonances in the reaction? How do resonances affect the reaction rates at low temperature? Here, we report a quantum dynamics study of the reaction on two new and more accurate PESs. Good agreement is achieved between the rate coefficients calculated on these two PESs and the experimental data. Our calculations reveal that the HB well in the entrance channel of OH + HCl supports many low energy resonance states. These resonance states substantially enhance the quantum tunneling effect and have an important impact on the reaction rates at low temperatures.In order to improve the fitting accuracy of the ZZGX PES, we constructed two new PESs using the fundamental invariant neural network (FI-NN) method, fitted to ∼70 000 ab initio energy points calculated at the UCCSD(T)-F12a and UCCSD(T)-F12b levels of theory, respectively, both with the AVTZ basis set. The fitting RMSE is 3.07 and 3.12 meV, respectively, for the F12a and F12b PES, which is about half that for the ZZGX PES.³⁶ The spin–orbit coupling of the channels of both the reagent OH and the product Cl have been included using FI-NN fitting to about 38 000 points calculated at the MRCI/aug-cc-pVTZ level of theory with a fitting error of 0.19 meV. The CASSCF wave-function with an active space of (5e, 3o) was used as a reference for MRCI. Details of the new PESs are provided in the ESI^†. The static barrier height for PESa is 0.088 eV (0.095 eV with SO correction included), and that for PESb is 0.097 eV (0.104 eV). As can be seen from Tables S1 and S2,^† the geometries and energies of all the stationary points for the F12b PES without SO correction are in good agreement with the ZZGX PES. Table S2^† also shows the corresponding complete basis set (CBS) energies for these stationary points based on AVTZ, AVQZ and AV5Z calculations. Because the F12b energies are slightly closer to the CBS results, we will present the dynamical results obtained on the F12b PES in the main text and provide those for the F12a PES in the SI.On the new PESs, we carried out potential-averaged five-dimensional (PA5D) time-dependent wave packet^43,44 calculations to obtain the total reaction probabilities for the reaction by freezing the non-reacting OH bond in its ground vibrational state. Tests revealed that the PA5D treatment is capable of providing reaction probabilities for the ground rovibrational initial state that are essentially identical to those obtained using the full six-dimensional approach, as shown in Fig. S2.^†Fig. 1(A) shows the total reaction probabilities for the HCl + OH reaction as a function of collision energy calculated on the F12b PES with both reagents in the ground rovibrational state at propagation times of 60 000, 120 000, 360 000, and 2 400 000 a.u. At high collision energies, the reaction probabilities converge quickly with respect to the propagation time, and one can barely see any difference among the four reaction probability curves, which exhibit smooth increases with collision energy. However, in the low collision energy region, large differences appear for different propagation times. At t = 60 000 a.u., the reaction probability presents a smooth curve with some small and broad oscillations, as was observed in the wave packet calculations of Guo and coworkers.³⁶Open in a separate windowFig. 1(A) Total reaction probabilities for the ground initial state of the OH + HCl → Cl + H₂O reaction on the F12b PES at wave packet propagation times of T = 60 000, 120 000, 360 000, and 2 400 000 a.u. (B) Same as (A), except showing the collision energy between 0.0 and 0.04 eV. The crosses mark the points for which the wavefunctions are shown in Fig. 2. (C) Total reaction probabilities for some partial waves J = 0, 30, 60, and 90 as a function of the collision energy. (D) Same as (C) except showing the collision energy between 0.0 and 0.04 eV.When the propagation time is increased to t = 120 000, the reaction probabilities at collision energies below 0.05 eV increase substantially. With further increasing the propagation time to 360 000 a.u., many oscillatory structures emerge at collision energies below the barrier height of 0.104 eV, in particular in the very low collision energy region as shown in Fig. 1(B). These sharp structures become fully converged after around 2 400 000 a. u. of wave packet propagation (∼58 ps). The reaction probability even at a collision energy close to zero reaches 3–4%. The convergence of the total reaction probabilities on the F12a PES is very similar to that on the F12b PES (Fig. S3^†), except that the final converged reaction probabilities for these two PESs exhibit a small shift, apparently due to slightly different barrier heights. Therefore, it is clear that reaction resonances exist in the title reaction in the very low collision energy region, and that these resonance states substantially induce quantum tunneling and enhance the reactivity. The lifetimes for these resonance states are quite long, with many being longer than 6.5 ps and having corresponding widths smaller than 0.1 meV. Fig. 1(C) presents converged (t = 2 400 000 a.u. ≈ 58 ps) total reaction probabilities for the total angular momentum J = 0, 30, 60, and 90. With increasing J, the reaction probability curve shifts to higher energy. In the low collision energy region (<0.05 eV), the total reaction probabilities for J = 30 exhibit rich oscillatory structures as in the J = 0 case (Fig. 1(D)), which are expected to have a great influence on the rate constant at low temperature. With further increase of the total angular momentum, the influence of these resonances on the total reaction fades due to the centrifugal barrier, which prevents the low-energy wave function from entering the well. They only leave a small trace in the total reaction probabilities at low energies for J = 60, and do not have any effect for J = 90.To understand the nature of these resonances, we calculated scattering wave functions at two collision energies (1.08 and 4.26 meV) with the peak reaction probabilities indicated by x in Fig. 1(B). Fig. 2(A) shows the two dimensional (2D) contour at the collision energy of 1.08 meV in the Jacobi coordinates HCl bond length (r_H–Cl) and center of mass distance between OH and HCl (R_HCl–OH), with the other coordinates integrated. As can be seen, the wave function is localized in the HB well region in the entrance channel with a few nodes in the R coordinate and no node in the r_H–Cl coordinate. Inspection of the scattering wave function for the bending and torsion coordinates reveals nodes exist in these coordinates (Fig. S4 and S5^†). The 2D contour in the coordinates of R_HCl–OH and r_H–Cl at the collision energy of 4.26 meV shown in Fig. 2(B) looks similar to that shown in Fig. 2(A), except with more nodes in the R direction. Therefore, the observed resonance states in the reaction are Feshbach resonances trapped in a bending/torsion excited vibrationally adiabatic potential (VAP) well in the reactant complex region due to the HB interaction.Open in a separate windowFig. 2Reactive scattering wave functions for the OH + HCl → Cl + H₂O reaction on the F12b PES in the two Jacobi coordinates R(HCl–OH) and r(H–Cl) with other coordinates integrated at the collision energies of 1.08 (A) and 4.26 meV (B). The contour lines are the corresponding 2D PESs along the two reactive bonds R(HCl–OH) and r(H–Cl) with other coordinates optimized. The geometries for the saddle point and the HB minimum are displayed in (A). The coordinate units in the figures are a₀. Fig. 3(A) shows the accurate rate constants for the initial ground rovibrational state, k_g, based on the probabilities for J = 0, 30, 60, 90 using the uniform J-shifting approach with a temperature-dependent shifting constant.^45–47 A test shows that the J-shifting scheme based on these four individual J values only introduces a few percent error to the rate constants in the temperature region considered here (Fig. S6^†). As can be seen from the figure, with decreasing temperature, k_g first decreases rapidly from 1000 K to 700 K, then decreases slowly. It reaches a minimum at T ≈ 260 K, and increases slowly with further decrease of the temperature. At temperatures lower than 300 K, the rate constants for the ground rovibrational initial state are larger than the measured thermal rate coefficients, with k_g being larger than k_exp by ∼70% at T = 200 K.Open in a separate windowFig. 3(A) Accurate rate constants, k_g, for the initial ground rovibrational state of the HCl + OH → H₂O + Cl reaction calculated on the F12b PES, in comparison with k^JS_g (obtained using the J-shifting approximation) and k^NR_g (based on the background reaction probabilities for J = 0 with the resonance contribution removed shown below); (B) the background reaction probabilities up to E = 0.1 eV by connecting some valleys of the reaction probabilities marked by x.Also shown in Fig. 3(A) are the rate constants for the ground rovibrational initial state, k^JS_g, obtained from the J = 0 reaction probabilities using the J-shifting approximation (see ESI^† for details). As can be seen, the J-shifting approximation works very well at high temperatures around 1000 K, but begins to overestimate the rates with decreasing temperature. At T = 500 K, k^JS_g is about 10% higher than the true rate. At temperatures below 300 K, the J-shifting approximation underestimates the reaction rate, with k^JS_g being smaller than k_g by about 16% at T = 200 K. Overall, the J-shifting approximation works fairly well for the ground rovibrational initial state, although there are numerous resonance peaks in the reaction probabilities in the low collision energy region.Now we consider the issue of the effect of the resonance structures found in the reaction on the rate constant. For a reaction system with isolated resonances, it is rather straightforward to remove the resonance contributions from the reaction probabilities by fitting the resonance peaks to some Lorentzian functions and to obtain smooth background scattering probabilities, as demonstrated in the F + HD reaction and recently in the inelastic scattering of H + HF (ref. 6 and 18). However, the OH + HCl reaction possesses numerous highly overlapped resonances in the low energy region, as shown in Fig. 1(A), and it is impractical to fit these resonance peaks as accurately as the isolated resonances. Instead, we obtained an approximate background reaction probability curve, which is shown in Fig. 3(B), by smoothly connecting some resonance valleys as shown in the figure. The rate constants k^NR_g calculated using the background curve shown in Fig. 3(B) with the J-shifting approximation were compared with the original k_g in Fig. 3(A). Given the fact that the J-shifting approximation works fairly well for k^JS_g, as shown in Fig. 3(A), it is very reasonable to expect that will work even better for the k^NR_g values obtained from the background reaction probabilities with the resonance contribution removed. As can be seen, k^NR_g exhibits rate behavior typical for systems with a low barrier with some quantum tunneling effects. At T = 1000 K, k^NR_g is essentially identical to k_g, but decreases much faster than k_g as the temperature drops. At T = 200 K, k_g is larger than k^NR_g by a factor of about 5.6 (9.8 × 10⁻¹³vs. 1.74 × 10⁻¹³ cm³ s⁻¹), indicating that the reaction probabilities in the low collision energy region due to resonances substantially enhance the rate constants for the reaction at low temperatures. It is worthwhile to note that for systems with overlapped resonances like that shown in Fig. 3(B), the reaction probabilities at the valleys must be considerably higher than the true background reaction probabilities; therefore, the background curve shown in Fig. 3(B) is the upper limit of the background reaction probabilities and the rate based on the curve shown in Fig. 3(A) is also the upper limit of the rate for the ground rovibrational initial state without the resonance contributions. Therefore, the true enhancement due to the resonances must be larger than that shown in Fig. 3(A).For reliable comparison with the measured thermal rate coefficients, we must take into account the contributions from all the thermally populated initial states of the reagents. Due to the very large number of thermally populated rotational states for this reaction even at 200 K, we opted to calculate the cumulative reaction probabilities, NE(E), (the sum of the reaction probabilities for all the initial states with a fixed total energy) for the reaction from which the thermal rate constants can be reliably evaluated.^48–52 The transition state wave packet calculations were carried out to obtain NE(E) using the details given in the ESI^†. Due to huge computational efforts required to obtain the cumulative reaction probabilities for J > 0, we only calculated NE(E) for J = 0 and employed the J–K-shifting approximation^45,49,53 to obtain the thermal rate constant. NE(J = 0, E) as a function of total energy measured with respect to the ground rovibrational energy of OH and HCl is presented in Fig. S7.^†In Fig. 4, we present thermal rate constants for the reaction calculated on both the F12a and F12b PES, together with the rate constants for the ground rovibrational state (k_g) and the previous experimental measurements.^{21,22,24–29} As can be seen, the thermal rate coefficients are smaller than k_g over the entire temperature region, indicating that reagent rotation excitations diminish the reaction rates. Overall, the thermal rate coefficients calculated on both PESs agree with the experimental results rather well, with the F12a PES slightly overestimating and the other PES slightly underestimating compared to the experimental measurements. As observed in the experiments, the thermal rate coefficients decrease quickly with decreasing temperature in the high-temperature region, but decrease much more slowly at low temperatures, in particular at T < 300 K, substantially deviating from Arrhenius behavior.Open in a separate windowFig. 4Thermal rate constants of the HCl + OH → H₂O + Cl reaction calculated on both the F12a and F12b PES, compared with rate constants for the ground rovibrational state, the RPMD rates on ZZGX PES and the previous experimental measurements.^{21,22,24–29} Experimental data are taken from Husain et al. (in circles),^21,22 Molina et al. (in diamonds),²⁴ Ravishankara et al. (in upward-pointing triangles),²⁵ Smith et al. (in squares),²⁶ Sharkey et al. (in triangle left),²⁷ Battin-Leclerc et al. (in downward-pointing triangles),²⁸ and Bryukov et al. (in leftward-pointing triangles).²⁹Also shown in Fig. 4 are the RPMD rates calculated (see ESI^† for details) on the F12a, F12b and ZZGX PESs.⁴² As shown, the RPMD rates on the F12b PES agree with those on the ZZGX PES extremely well, although the barrier heights for the F12b PES are higher by 9 meV. The present RPMD rates agree with the J–K-shifting quantum rates rather well, except at 200 K, at which the RPMD rates are overestimated by about 40%. As a result, the RPMD rate is higher than the experimental value at T = 200 K, in particular for the F12a PES. Previous studies have shown that RPMD tends to underestimate reaction rates in the strong quantum tunneling region, even for systems with resonances such as the O(³P) + HCl → OH + H reaction;⁵⁴ however, it overestimates the rates for the F + H₂ reaction with pronounced post-barrier Feshbach resonances⁵⁵ and some insertion reactions.⁵⁶ The discrepancy between the RPMD and J–K-shifting quantum rate at T = 200 K could be caused by a possible underestimation of the rate by the J–K-shifting approximation, as in the J-shifting approximation of the ground rovibrational initial state shown in Fig. 3(A). On the other hand, the rapid increase in the rate with decreasing temperature (from 250 K to 200 K) could also be a feature of bimolecular reaction rates in the high-pressure limit for reactions with a pre-reactive minimum.⁵⁷ RPMD rate theory in its original form employs a free energy calculation,^38,39 which allows thermalization in the pre-reactive minimum, sampling tunneling pathways with energies not accessible in the low-pressure limit. Therefore, RPMD models a rate process that resembles rates in the high-pressure limit, which could lead to an overestimation of the rate at low temperature. A recent study also reported a similar issue with RPMD.⁵⁸ The authors attributed it to spurious resonances in RPMD favoring energy transfer in the pre-reactive minimum, which is similar to our argument. Certainly, more quantum dynamics calculations to provide rigorous quantum rates at low temperatures to assess the accuracy of the J–K-shifting approximation and the RPMD method on this reaction system with strong HB interaction in the entrance channel would be highly desirable.Therefore, the OH + HCl reaction possesses many long-lifetime Feshbach resonances trapped in a bending excited VAP well in the entrance channel due to HB interaction. These resonance states substantially induce quantum tunneling of the hydrogen atom through the reaction barrier and enhance the reactivity in the low-collision region. Consequently, the reaction rate for the reaction becomes essentially independent of temperature in the low-temperature region and deviates substantially from Arrhenius behavior, as observed experimentally. The resonance states in the reaction are very different in location from those in the F + H₂/H₂O reactions,^{8,11,12,14,59–63} which are trapped in the VAP well in the product channel. Furthermore, they are trapped in the bending/torsion excited VAP well with HCl in the ground vibrational state, unlike those in the F + H₂/H₂O reactions with HF in vibrationally excited states.^8,63,64 In nature, these resonance states in the reaction arise from HB interaction as in the F + H₂O (v = 0) reaction,^11,62,63 but are different from those in the F + H₂/HD/HOD (v = 1) reactions^{8,10,12,14,64–66} due to chemical bond softening. Because the resonance states are trapped in the bending/torsion excited VAP well, their lifetimes are much longer than those observed in the F + H₂/H₂O reactions. As a result, they should have an important impact on the differential cross sections. More efforts, and in particular more joint efforts involving theory and experiment, should be devoted to studying these Feshbach resonances in this reaction of great atmospheric importance in detail.     

Single-atom cobalt-hydroxyl modification of polymeric carbon nitride for highly enhanced photocatalytic water oxidation: ball milling increased single atom loading

Expediting the oxygen evolution reaction (OER) is the key to achieving efficient photocatalytic overall water splitting. Herein, single-atom Co–OH modified polymeric carbon nitride (Co-PCN) was synthesized with single-atom loading increased by ∼37 times with the assistance of ball milling that formed ultrathin nanosheets. The single-atom Co-N₄OH structure was confirmed experimentally and theoretically and was verified to enhance optical absorption and charge separation and work as the active site for the OER. Co-PCN exhibits the highest OER rate of 37.3 μmol h⁻¹ under visible light irradiation, ∼28-fold higher than that of common PCN/CoO_x, with the highest apparent quantum yields reaching 4.69, 2.06, and 0.46% at 400, 420, and 500 nm, respectively, and is among the best OER photocatalysts reported so far. This work provides an effective way to synthesize efficient OER photocatalysts.

Single-atom Co^II-OH modified polymeric carbon nitride synthesized with increased single-atom loading under the assistance of ball milling exhibits high photocatalytic water oxidation activity with Co-N₄OH as the highly active site.

Massive fuel energy consumption induced environmental and ecological problems, especially the greenhouse effect, and the resultant extreme climates and rise in sea level are threatening human life.¹ As a potential substitution for fuel energy, hydrogen energy conversion from solar energy via photocatalytic water splitting attracts great attention from scientists.^2–5 However, the photocatalytic hydrogen evolution efficiency from overall water splitting is still restricted by the sluggish oxygen evolution reaction (OER) that involves energy absorption, four-electron transfer, breakage of O–H bonds, and formation of O–O bonds,^6,7 and thus efficient OER photocatalysts become the key to achieving efficient overall water splitting. Though numerous hydrogen evolution photocatalysts have been reported, research on OER photocatalysts is mainly around a few semiconductors including BiVO₄, WO₃, Ag₃PO₄, α-Fe₂O₃, etc.8–11 and their activity is not high enough yet for practical applications. Therefore, exploring high-efficiency OER photocatalysts is still necessary.Polymeric carbon nitride (PCN) was first reported in 2009 (ref. 12) as a photocatalyst with a layered melon-type orthorhombic structure,¹³ and thereafter quickly became a “star” photocatalyst thanks to its advantages of being visible-light responsive and metal-free, non-toxic, and low cost, and its relatively high chemical stability.¹⁴ Because of several self-deficiencies including fast photogenerated charge recombination and a narrow optical absorption spectrum, PCN exhibits relatively low photocatalytic activity.¹⁵ Then, a series of strategies were put forward successively to enhance the photoactivity of PCN, such as enhancement of crystallinity,¹⁶ morphological control,¹⁷ structural modification¹⁸ (including extensively researched single atom modification in recent years^19,20), exfoliation,²¹ construction of hetero-(homo-)junctions,²² and loading of noble metals.²³ Though photocatalytic water splitting on PCN was extensively researched in the past, the research was mainly around the hydrogen evolution half-reaction used for exploring properties and the catalytic mechanism of photocatalysts, and little research was focused on the industrially useable overall water splitting process owing to the sluggish OER.¹⁵ Therefore, enhancing the photocatalytic OER activity of PCN becomes the key to practical applications.To increase OER rates of PCN, several kinds of methods were proposed, such as rational design of compound cocatalysts (e.g., CoO_x, IrO₂, CoP, CoPi, RhO_x, RuO_x, PtO_x, MnO_x, Co(OH)₂, Ni(OH)₂, and CoAl₂O₄ (ref. 24–30)), modification of carbon dots and carbon rings,^31,32 fabrication of special architectures of PCN (e.g., PCN quantum dot stacked nanowires³³), and single-atom (e.g., B, Co, and Mn^34–36) modification. For instance, Zhao and coauthors prepared B and N-vacancy comodified PCN that exhibits the highest OER rate of ∼28 μmol h⁻¹ (ref. 36) and recently their group further used these B doped PCN ultrathin nanosheets to fabricate a Z-scheme heterojunction for overall water splitting with a solar-to-hydrogen efficiency reaching ∼1.2%.³⁷ Comparatively, PCN loaded with compound cocatalysts can only enhance OER activity to a limited degree and there are finite methods for carbon modification and special architecture fabrication. Single-atom modification shows a bright prospect, on account of metal atoms capable of being inserted into the framework of PCN and effectively increasing the OER activity. However, reported single metal atom modification routes are all based on direct ion adsorption on PCN or calcination of mixtures of metal salts and PCN feedstocks.^34,35,38 New routes need be explored to increase effective loading of single atoms in PCN. Besides, the metal-OH structure is considered efficient for the OER,^30,39,40 and a single metal atom-OH structure has never been reported for modification of PCN, though Mn–OH was thought to play a key role in the OER process.³⁴Ball milling is an extensively used versatile and scalable way for preparation of heterogeneous catalysts and even single-atom catalysts,^41,42 but was rarely used in synthesis of PCN-based single-atom photocatalysts. In this work, we synthesized single-atom Co–OH modified PCN (Co-PCN) with the single-atom content in PCN highly increased with the assistance of ball milling. The simple synthetic route is shown in Fig. 1a. PCN was ball-milled to obtain BM-PCN that then adsorbed Co²⁺ till saturation to form BM-PCN/Co which was calcined to obtain BM-PCN/Co-c (Co-PCN). For comparison, PCN was directly used to adsorb Co²⁺ till saturation to form PCN/Co which was calcined to obtain PCN/Co-c. PCN mainly comprises large blocks with the size of several micrometers (Fig. S1^†), while BM-PCN contains massive irregular particles with the size reduced to several hundreds of nanometers (Fig. S2^†), indicative of high efficacy of ball milling. BM-PCN/Co-c exhibits a similar morphology as BM-PCN (Fig. 1b and S3^†) and PCN/Co-c exhibits a similar morphology to PCN (Fig. S4^†), but the surface area and mesopore volume of BM-PCN and BM-PCN/Co-c are not higher than those of PCN and PCN/Co-c (Fig. S5^†), manifesting that ball-milling and subsequent calcination did not form massive mesopores, which accords well with the particle size variation from several micrometers (before ball milling) to several hundreds of nanometers (after ball milling). However, the Co content in BM-PCN/Co-c, BM-PCN/Co, PCN/Co-c, and PCN/Co was measured to be 0.75, 0.50, 0.02, and ∼0.02 wt%, respectively, by inductively coupled plasma mass spectrometry (ICP-MS). The ∼37 times higher Co content in BM-PCN/Co-c than in PCN/Co-c suggests the ball-milling enhanced adsorption of Co²⁺ on surfaces of BM-PCN, which should arise mainly from the ball-milling induced increase of surface energy and adsorption sites.⁴³Open in a separate windowFig. 1(a) Schematic illustration for synthesis of single-atom Co^II-OH modified PCN (BM-PCN/Co-c); and (b) SEM, (inset in b) TEM, (c) AFM, (d) EDS elemental mapping, and (e) HAADF-STEM images of BM-PCN/Co-c.The TEM image shows the existence of small and ultrathin nanosheets in BM-PCN/Co-c (inset in Fig. 1b) which can also be observed in the atomic force microscopy (AFM) image with a thickness of ∼7–10 nm and lateral size of <70 nm (Fig. 1c), and formation of these ultrathin nanosheets results from the ball milling of PCN.⁴⁴ It should be noted that most formed ultrathin nanosheets with high surface energy may stack into compact particles upon ball milling, leading to no increase of the total surface area. Energy dispersive X-ray spectroscopy (EDS) elemental mapping images of BM-PCN/Co-c indicate homogeneous distribution of C, N, O, and Co elements in the sample (Fig. 1d). The high-angle annular dark-field scanning transmission electron microscopy (HAADF-STEM) image of BM-PCN/Co-c shows massive white spots (marked by circles) with a mean size of <1 Å dispersed in the sample (Fig. 1e and S6^†), which should correspond to single-atom Co.To further verify the single-atom Co structure in BM-PCN/Co-c, Co K-edge X-ray absorption near-edge structure spectroscopy (XANES) and extended X-ray absorption fine structure (EXAFS) analysis were performed. As shown in Fig. 2a, the absorption-edge position of BM-PCN/Co-c is quite close to that of CoO and their peak positions are similar and far from those of other reference samples, indicating that the valence of Co in BM-PCN/Co-c is about +2. The bonding structure around Co was determined by Fourier transformed (FT) k³-weighted EXAFS analysis. As shown in Fig. 2b, a distinct single Co-ligand peak at ∼1.6 Å for BM-PCN/Co-c is observed, which prominently differs from the Co–Co coordination peak at ∼2.2 Å for Co foil and the Co^II–O coordination peak at ∼1.7 Å for CoO. The wavelet transform (WT) contour plot of BM-PCN/Co-c shows only one intensity maximum (Fig. S7^†), and the Cl 2p core-level XPS spectrum of BM-PCN/Co-c reveals no residue of Cl (Fig. S8^†). These further indicate the single-atom dispersion of Co species. Apparently, the Co-ligand peak is almost consistent with the Co^II–N peak for Co porphyrin, suggesting that the single-atom Co in BM-PCN/Co-c mainly coordinates with N. Least-square EXAFS curve fitting was performed to confirm quantitative structural parameters of Co^II in BM-PCN/Co-c, as shown in Fig. 2c, S9, and S10 and Table S1.^† Simple Co–N single-shell fitting of BM-PCN/Co-c (Fig. S10^†) gave a coordination number of 5.6 ± 0.4 (Table S1^†), that is, Co^II coordinates with five atoms. Considering that the PCN monolayer provides four appropriate N coordination sites at most,⁴⁵ Co^II likely coordinates with four N atoms and one OH atom. Thus, we further performed Co–N₄/Co–O double-shell fitting (Fig. 2c) and the obtained R-factor (0.0011) remarkably reduces relative to that from Co–N single-shell fitting (0.0035), indicative of rationality of the proposed Co^II–N₄OH structure. Confirmed Co–N and Co–O bond lengths are 2.04 and 2.15 Å, respectively (Table S1^†).Open in a separate windowFig. 2(a) Co K-edge XANES and (b) EXAFS spectra of Co foil, Co porphyrin (Copr), CoO, Co₃O₄, Co₂O₃, and BM-PCN/Co-c; EXAFS (c) R space-fitting and (inset in c) K space-fitting curves of BM-PCN/Co-c; (d) optimized structure of PCN and Co-doped PCN with different doping configurations and calculated formation energies (e) of Co doped PCN; and (e) Co 2p and (f) O 1s core-level XPS spectra of samples.To further confirm rationality of the Co–N₄OH coordination structure, density functional theory (DFT) calculations were conducted. As shown in Fig. 2d, three possible Co^II coordination structures in the PCN monolayer were explored. The Co–N₄OH structure without removal of H from PCN exhibits a much lower formation energy (∼0.15 eV) than Co–N₄ and Co–N₃ structures with removal of two H atoms from PCN (∼2.51 and 3.55 eV), demonstrating a high probability of existence of the Co–N₄OH structure in BM-PCN/Co-c. This structure can also be evidenced by X-ray photoelectron spectroscopy (XPS). As shown in Fig. 2e, the Co 2p core-level XPS spectrum of BM-PCN/Co-c shows two distinct peaks at binding energies of 796.8 and 781.4 eV beside satellite peaks, corresponding to Co 2p_1/2 and 2p_3/2 of Co^II ions.⁴⁶ The spectrum of BM-PCN/Co also shows two Co 2p peaks but at a binding energy ∼1.1 eV higher, suggesting variation of the Co coordination structure from BM-PCN/Co to BM-PCN/Co-c. PCN/Co-c exhibits no peaks because of its low Co content. Fig. 2f shows O 1s core-level spectra of PCN, BM-PCN, BM-PCN/Co-c, and PCN/Co-c. All the samples exhibit one peak at a binding energy of ∼532.0 eV, ascribed to surface hydroxyl species,⁴⁷ but an additional peak could be obtained for BM-PCN or BM-PCN/Co-c after deconvolution. The peak at a binding energy of ∼531.3 eV for BM-PCN should be assigned to adsorbed H₂O at new active adsorption sites generated by ball milling. This peak can also be observed in the spectrum of BM-PCN/Co, but with a ∼0.1 eV shift to a higher binding energy (Fig. S11^†) owing to the influence of adsorbed Co^II ions. The peak at ∼531.2 eV for BM-PCN/Co-c should be assigned to Co–OH,⁴⁸ given that there is only one O 1s peak for BM-PCN-c (synthesized by direct calcination of BM-PCN) (Fig. S11^†). The calculated Co/O(–Co) molar ratio, based on the XPS data, is ∼1.07 (Table S2^†), close to 1, consistent with the Co–N₄OH coordination structure.In C 1s and N 1s core-level XPS spectra, BM-PCN, BM-PCN/Co-c, PCN/Co-c, and BM-PCN/Co exhibit similar peaks to PCN (Fig. S12a–d^†), indicative of their similar framework structure which can also be evidenced by their similar N/C molar ratios, 1.53 (Table S3^†), but the N–H peak of BM-PCN shifts ∼0.2 eV to a lower binding energy relative to that of PCN, likely arising from the ball-milling induced destruction of intralayer hydrogen bonds (Fig. S13^†). The Co content in BM-PCN/Co, BM-PCN/Co-c, and PCN/Co-c is too low to cause detectable variation of C 1s and N 1s peaks. Similar FT-IR absorption bands of the samples (Fig. S14a and b^†) also indicate their basic frame structure, but in enlarged spectra (Fig. S14c^†), ν(C–N) and ν(C N) absorption bands of BM-PCN shift 16 cm⁻¹ to a higher wavenumber and 19 cm⁻¹ to a lower wavenumber, respectively, relative to those of PCN at 1242 and 1640 cm⁻¹,⁴⁹ likely resulting from the ball-milling induced hydrogen bond destruction, and the shift of these two absorption bands turns smaller for BM-PCN/Co-c, suggesting calcination-induced reforming of the destroyed hydrogen bonds, which is consistent with the XPS results (Fig. S12c^†). Besides, BM-PCN exhibits a wider and relatively stronger ν(N–H)/ν(O–H) absorption band than PCN (Fig. S14a^†), probably owing to the hydrogen bond destruction and new adsorbed H₂O, while this absorption band for BM-PCN/Co-c becomes much weaker, suggesting hydrogen bond reforming and loss of new adsorbed H₂O (Fig. 2f). Zeta potentials of the samples dispersed in water reflect variation of surface adsorbed hydroxyl species. As shown in Fig. S15a,^† all the samples exhibit negative zeta potentials because of dissociation of surface hydroxyl species. The zeta potentials, following the order PCN (−24 mV) > BM-PCN (−41 mV) < BM-PCN/Co-c (−30 mV) ≈ PCN/Co-c (−28 mV), suggest the ball milling-induced increase of surface hydroxyls in BM-PCN and calcination-induced decrease in BM-PCN/Co-c, consistent with the FT-IR results.Solid-state ¹³C magic-angle-spinning nuclear magnetic resonance (NMR) spectra of PCN, BM-PCN, BM-PCN/Co-c, and PCN/Co-c show two similar peaks at chemical shifts of ∼164 and 156 ppm (Fig. S15b^†), ascribed to C−NH_x and N C–N, respectively,⁵⁰ indicating their similar molecular framework, but in enlarged spectra, BM-PCN exhibits ∼0.3° movement of the N C–N peak to a lower chemical shift compared with PCN, because of the ball-milling induced hydrogen bond destruction, and the C−NH_x peak of BM-PCN/Co-c moves ∼0.2° to a lower chemical shift, likely owing to formation of the C–N–Co structure whose peak lies close to the C−NH_x peak.⁵¹ The XRD patterns of the samples are shown in Fig. S15c.^† PCN and PCN/Co-c exhibit typical diffraction peaks of melon-type carbon nitride with a layered orthorhombic structure and peaks at 13.1° and 27.6° correspond to (210) and (002) facets, respectively,^13,52 but BM-PCN reveals remarkably decreased peak intensity and ∼0.2° shift of the (002) peak to a lower 2θ (indicative of the increased interlayer distance) relative to PCN, demonstrating the ball-milling induced hydrogen bond destruction and substantial decrease of crystallinity. The remarkable decrease of crystallinity and almost no change of the surface area of BM-PCN, compared with those of PCN, further suggest that ball milling may form massive thin nanosheets (Fig. 1c) most of which stack into compact particles (Fig. 1b) owing to their high surface energy. In comparison with BM-PCN, BM-PCN/Co-c exhibits a narrower (002) peak, suggesting enhanced crystallinity owing to the calcination-induced hydrogen bond reforming, consistent with the FT-IR results. On the whole, it is likely the ball-milling induced destruction of hydrogen bonds that contributes largely to the increase of surface energy and new active adsorption centers and thus Co²⁺ adsorption on BM-PCN.Optical absorption capability of samples was investigated by UV-vis diffuse reflectance spectroscopy (DRS). As shown in Fig. 3a, BM-PCN/Co-c, BM-PCN, and PCN/Co-c exhibit considerably higher, lower, and similar optical absorption than/to PCN, respectively. For BM-PCN/Co-c, the optical absorption enhancement at a wavelength of <400 nm may benefit from the electron-rich Co that enhances π–π* transitions in heptazine rings,⁵³ and the Urbach tail absorption should arise from the Co–OH doping.^54,55 Bandgaps (E_g) of PCN, BM-PCN, BM-PCN/Co-c, and PCN/Co-c were roughly confirmed as 2.70, 2.81, 2.56, and 2.73 eV, respectively, via the formula E_g/eV = 1240/(λ_ed/nm)⁵⁶ where λ_ed is the absorption edge determined by solid lines in the spectra. The wider E_g of BM-PCN probably results from the quantum size effect of massive ultrathin crystal nanosheets (Fig. 1c) formed by ball milling, and the narrower E_g of BM-PCN/Co-c arises from the Co–OH doping that was then verified by DFT calculations. As shown in Fig. S16,^† the calculated E_g of BM-PCN/Co-c, ∼1.90 eV, is much smaller than that of PCN (2.57 eV), in accordance with the experimental results. For PCN, the conduction band (CB) is contributed by C 2p and N 2p orbitals and the valence band (VB) mainly by N 2p orbitals, while for BM-PCN/Co-c, the CB is contributed by Co 3d, C 2p, and N 2p orbitals and the VB mainly by Co 3d and N 2p orbitals (Fig. S16c and d^†), effectively manifesting that the narrowing of E_g of BM-PCN/Co-c results from the Co–OH doping. In addition, there are prominent doping levels (E_d) in the bandgap of BM-PCN/Co-c, mainly contributed by Co 3d and O 2p orbitals (Fig. S16d^†), effectively proving the Co–OH doping effect in BM-PCN/Co-c. Similar calculation results have been reported for Pt–OH modified carbon nitride.⁵⁷ Given that the experimental Co content (0.75 wt%) is much lower than the theoretical (6.71 wt%), practical doping levels in the bandgap may approach more to the VB. CB edges of the samples (E_CB) could be roughly determined by using Mott-Schottky plots (Fig. S17^†) and their Fermi levels (E_f) were subsequently confirmed based on VB-XPS spectra (Fig. S18^†). Energy band levels of the samples are shown in Fig. 3b, and it seems that ball milling causes a slight downshift of the VB edge (E_VB) of BM-PCN, favorable for photocatalytic water splitting, but the Co–OH doping causes a slight downshift of E_CB and upshift of E_VB of BM-PCN/Co-c. It is noteworthy that the E_d close to the VB edge (E_VB) can capture photogenerated holes⁵⁸ and thus the single-atom Co–OH works as the active site for the OER (Fig. 3b).Open in a separate windowFig. 3(a) UV-vis diffuse reflectance spectra of PCN, BM-PCN, BM-PCN/Co-c, and PCN/Co-c; (b) energy band levels of the samples and schematic illustration for water oxidation on BM-PCN/Co-c; (c) photoluminescence spectra, (d) time-resolved fluorescence spectra, and (e) anodic photocurrent (J_a) response of the samples; and (f) EPR spectra of the samples in the dark and under visible light irradiation. Data in (d) are the results of fitting decay curves to a tri-exponential model. Dark J_a in (e) was set as zero for distinct comparison.Spectroscopy and photoelectrochemical tests were conducted to evaluate photogenerated charge separation and transfer performance. As shown in Fig. 3c, photoluminescence (PL) spectra of all the samples show one emission peak, basically corresponding to their bandgap emission. BM-PCN exhibits weaker PL intensity than PCN, revealing a decreased photogenerated charge recombination efficiency, which originates from faster charge transfer from the inside to the surface of ultrathin nanosheets (Fig. S19^†) and trapped by surface states.⁵⁹ BM-PCN/Co-c exhibits the lowest PL intensity and the PL intensity of PCN/Co-c is lower than that of PCN, which arises from the E_d capturing photogenerated holes to reduce their direct recombination with electrons beside the ultrathin nanosheet effect in BM-PCN/Co-c. Fig. 3d shows time-resolved fluorescence spectra of the samples. Decay curves were well fitted to a tri-exponential model (S3) and the obtained results are shown in Fig. 3d. Three lifetimes (τ₁–τ₃) and their mean lifetime (τ_m, 89.2 ns) of BM-PCN are all much longer than those of PCN (τ_m = 17.9 ns), further suggesting the faster charge transfer from the inside to the surface of ultrathin nanosheets in BM-PCN, decreasing the direct charge recombination efficiency, but with subsequent surface radiative recombination.⁶⁰ Interestingly, the τ₁–τ₃ and τ_m (10.8 ns) of BM-PCN/Co-c are much shorter than those of PCN, which should result from faster transfer of holes to E_d that effectively decreases the charge recombination efficiency, with subsequent nonradiative energy transformation.⁶¹ The Co–OH doping effect also makes PCN/Co-c exhibit shorter τ₁–τ₃ and τ_m (16.5 ns) than PCN. Fig. 3e shows the photocurrent response of the samples. Their anodic photocurrent density follows the order PCN < PCN/Co-c < BM-PCN < BM-PCN/Co-c, indicating gradually increased photogenerated charge separation efficiencies,⁶² basically consistent with the PL results. The relatively high photocurrent response of BM-PCN benefits from the applied bias that effectively inhibits surface recombination of photogenerated charge carriers.To assess charge mobility of the samples, their electrochemical impedance spectroscopy (EIS) spectra were tested with high-frequency data simply fitted to an equivalent circuit (Fig. S20^†). The obtained charge transfer resistance (R_ct) follows the order PCN (26 Ω) > BM-PCN (18 Ω) ≈ PCN/Co-c (19 Ω) > BM-PCN/Co-c (13 Ω). Apparently, BM-PCN/Co-c exhibits smaller R_ct than BM-PCN and PCN/Co-c, and PCN/Co-c exhibits smaller R_ct than PCN, indicating the highest charge transfer performance of BM-PCN/Co-c⁶³ which originates from the single-atom Co modification⁶⁴ that may increase the electron density to facilitate charge transport. The smaller R_ct of BM-PCN than that of PCN indicates the additional favorable effect of ultrathin nanosheets.⁶⁵Fig. 3f shows electron paramagnetic resonance (EPR) spectra of the samples. All reveal one single Lorentzian line centered at a g of 2.0039, attributed to unpaired electrons in heptazine rings.⁶⁶ In the dark, the EPR signal intensity follows the order PCN < BM-PCN < PCN/Co-c < BM-PCN/Co-c, and the stronger signal of BM-PCN than that of PCN results from formation of ultrathin nanosheets that enhances delocalization of unpaired electrons, while the stronger signal of BM-PCN/Co-c and PCN/Co-c mainly benefits from the Co doping that increases the delocalized electron density.⁶⁷ Under visible light irradiation, the samples exhibit remarkable signal enhancement, following the sequence PCN < BM-PCN < PCN/Co-c < BM-PCN/Co-c, similar to that of the signal intensity in the dark, suggesting that the increase in the delocalized electron density facilitates charge photoexcitation. The high delocalized electron density favors charge transport, consistent with the EIS results, and the high photoexcited charge density benefits enhancement of photocatalytic activity.Photocatalytic OER activity of various samples was well evaluated using Ag⁺ as the sacrificial agent (Fig. S21^†). The Co content in BM-PCN/Co-c was optimized according to the photocatalytic OER rates and BM-PCN-c exhibits no detectable OER activity (Fig. S22^†), indicating indispensability of the Co–OH structure for the OER. The influence of the calcination temperature (T_c °C) of BM-PCN/Co on OER rates of BM-PCN/Co-c (T_c = 460) and BM-PCN/Co-cT_c was investigated and BM-PCN/Co-c exhibits the highest photoactivity (Fig. 4a), manifesting that the optimal calcination temperature is 460 °C. Under both simulated solar light and visible light irradiation (λ ≥ 420 nm), BM-PCN/Co-c exhibits substantially higher OER activity than PCN/Co-c (Fig. 4b), further suggesting the significance of the single-atom Co loading amount, and remarkably higher activity than common PCN/CoO_x (with 0.75 wt% Co, obtained via photodeposition) and BM-PCN-c/Co(OH)₂ (with 0.75 wt% Co), demonstrating the high efficacy of the single-atom distribution of Co–OH in BM-PCN/Co-c. Besides, urea was used as the feedstock to synthesize carbon nitride (marked as PCN-urea) with a larger surface area (76 m² g⁻¹ (ref. 68)) than PCN, and PCN-urea was further used to synthesize PCN-urea/Co-c similar to the synthesis of BM-PCN/Co-c. The OER activity of BM-PCN/Co-c is prominently higher than that of PCN-urea/Co-c (with the optimized Co content and Co single atom distribution, Fig. S23^†), suggesting the significant role of ball milling in fabricating the single-atom Co–N₄OH structure. To quantitively compare photoactivity of the samples, their mean OER rates under visible light illumination for 2 h are shown in Fig. 4c. The OER rate of BM-PCN/Co-c can reach ∼37.3 μmol h⁻¹, about 13.8, 28.7, 2.6, and 2.0 times those of PCN/Co-c, PCN/CoO_x, BM-PCN-c/Co(OH)₂, and PCN-urea/Co-c, respectively. Comparatively, less N₂ was generated for BM-PCN/Co-c (Fig. S24^†), further demonstrating the significance of single-atom Co–OH modification.Open in a separate windowFig. 4(a) The influence of the calcination temperature (T_c °C) of BM-PCN/Co on photocatalytic OER activity of BM-PCN/Co-c (T_c = 460) and BM-PCN/Co-cT_c, under Xe-lamp illumination, with AgNo₃ as the sacrificial agent; (b) photocatalytic oxygen evolution on various samples under Xe-lamp illumination with or without using a 420-nm filter; (c) corresponding OER rates of the samples in 2 h; (d) photocatalytic OER rates of BM-PCN/Co-c under irradiation with various monochromatic light sources for 12 h; (e) apparent quantum yields (AQYs) of BM-PCN/Co-c at different wavelengths and reaction times and the highest AQY at every wavelength, along with the UV-DRS spectrum; and (f) proposed mechanism for photocatalytic water oxidation on the single-atom Co^II-OH structure.Photocatalytic oxygen evolution on BM-PCN/Co-c was also tested under monochromatic light irradiation (Fig. S25^†). Apparently, BM-PCN/Co-c can exhibit OER activity even at a wavelength of 500 nm. The mean OER rate in 12 h decreases from 1.85 to 0.54 μmol h⁻¹ with increasing wavelengths from 400 to 500 nm (Fig. 4d), independent of light intensity of the Xe lamp and is mainly dependent on optical absorption capability of BM-PCN/Co-c at various wavelengths (Fig. 3a). Fig. 4e shows apparent quantum yields (AQYs) of BM-PCN/Co-c at different reaction times and wavelengths. Basically, there are maxima of AQYs with increasing reaction time at every wavelength, suggesting the adverse effect of excessive photodeposited Ag on surfaces of samples. These maxima are shown in Fig. 4e and accord well with the UV-vis DRS spectrum with increasing wavelengths. The maxima of AQYs at 400, 420, 450, and 500 nm can reach 4.69, 2.06, 1.07, and 0.46%, respectively. Compared with the reported photocatalytic OER results for PCN (Table S4^†), BM-PCN/Co-c exhibits the top-class performance.To investigate chemical stability of BM-PCN/Co-c, the cyclic OER experiment was conducted. After five consecutive runs, OER rates of BM-PCN/Co-c decrease less (Fig. S26a^†), with the morphology similar to the original (Fig. S26b^†). Co single atoms in the sample could still be distinctly observed by HAADF-STEM (Fig. S26c and d^†). In addition, N 1s core-level XPS spectra of BM-PCN/Co-c are almost similar before and after the cyclic experiment (Fig. S26e^†). These indicate the high stability of the basic framework structure of the sample. However, Co 2p core-level spectra show remarkable differences before and after the experiment, not only the Co^II peak shift, probably owing to ion (e.g., IO₄⁻) adsorption, but also formation of a large amount of Co^III (Fig. S26f^†). Coexistence of Co^II/Co^III may suggest the photocatalytic OER mechanism.The proposed OER mechanism based on the Co–OH structure is shown in Fig. 4f, according to the reported results in Mn doped PCN.³⁴ Four holes are needed to complete four oxidation steps and obtain one O₂ molecule. The first step starting with one hole may involve formation of the Co^III O bond. The Co–N₄OH structure should facilitate the water oxidation more compared with that of Co–N₄ without OH coordination, by leaving out the initial adsorption process of H₂O molecules.³⁴ On the whole, the high photocatalytic OER activity of Co-PCN benefits from the Co–N₄OH structure that not only effectively enhances optical absorption, and charge separation and transport, but also works as the highly active site for the OER.     

One-electron bonds in copper–aluminum and copper–gallium complexes

Odd-electron bonds have unique electronic structures and are often encountered as transiently stable, homonuclear species. In this study, a pair of copper complexes supported by Group 13 metalloligands, M[N((o-C₆H₄)NCH₂PⁱPr₂)₃] (M = Al or Ga), featuring two-center/one-electron (2c/1e) σ-bonds were synthesized by one-electron reduction of the corresponding Cu(i) ⇢ M(III) counterparts. The copper bimetallic complexes were investigated by X-ray diffraction, cyclic voltammetry, electron paramagnetic spectroscopy, and density functional theory calculations. The combined experimental and theoretical data corroborate that the unpaired spin is delocalized across Cu, M, and ancillary atoms, and the singly occupied molecular orbital (SOMO) corresponds to a σ-(Cu–M) bond involving the Cu 4p_z and M ns/np_z atomic orbitals. Collectively, the data suggest the covalent nature of these interactions, which represent the first examples of odd-electron σ-bonds for the heavier Group 13 elements Al and Ga.

Hanging on by a thread. Formally zerovalent copper complexes with an Al(iii) or Ga(iii) support were investigated. The combined experimental and theoretical data corroborate the presence of an odd-electron σ-bond between Cu and the Group 13 center.

Odd-electron σ-bonds, where the electrons are delocalized between two atoms, can occur as two-center/one-electron (2c/1e) or two-center/three-electron (2c/3e) interactions. Proposed by Pauling in 1931,¹ odd-electron σ-bonds have garnered attention because of their fundamental importance to chemical bonding and their relationship to radical species generated during oxidative stress in biological systems.^2–14 Examples of compounds exhibiting odd-electron bonding are typically homonuclear (like H₂⁺, He₂⁺, and alkali metal dimers) and transiently stable, limiting them to spectroscopic characterization.^1,11,15–18The first solid-state structure of a formally one-electron σ-bond was a tetraphosphabenzene species (Fig. 1a) which was formed by the coupling of two diphosphirenyl radicals.¹⁹ Following this discovery, the formation of discrete 2c/1e σ-bonds, where the odd-electron is delocalized between two homonuclear main group centers, was reported for B·B and then extended to P·P.^8,17,20 Of note, the first solid-state structure of a B·B compound was reported in only 2014 (Fig. 1b).²¹ Examples of 2c/1e σ-bonds between the heavier Group 13 congeners are even more lacking because of the greater propensity for their unpaired spins to couple, forming larger more stable clusters.⁸ To our knowledge, there are only three structurally characterized examples of odd-electron bonds for the heavy Group 13 atoms,²² and these examples are all homonuclear π-radicals (Fig. 1c).^23–26Open in a separate windowFig. 1Select examples of structurally characterized molecules (a–d) featuring odd-electron bonds.Heteronuclear odd-electron σ-bonds are also rare. The Cu(TPB) complex, where TPB is a trisphosphinoborane, is the single structural example of a 2c/1e bond between heteroatoms (Fig. 1d).²⁷ The authors described the bonding as Cu·B, where the unpaired electron is heavily polarized toward B. A theoretical study predicted that such a bond would also exist between Cu and Al, but no heavier analogues of Cu(TPB) have been synthesized to date.²⁸ Furthermore, the heavier Group 13 elements by virtue of their lower electronegativity compared to B should facilitate greater covalent interactions with the Cu center.Hence, we sought to target formally zerovalent Cu complexes supported by Al(III) or Ga(III) as an extension of the previously reported isoelectronic nickelate species and Cu(TPB).²⁹ Herein, we describe the synthesis, structure, spectroscopic characterization, and DFT calculations of cationic [CuML]⁺ complexes (L = [N((o-C₆H₄)NCH₂PⁱPr₂)₃]³⁻; M = Al and Ga) as well as their one-electron reduced metalloradical counterparts that feature discrete 2c/1e bonds.     

Direct synthesis of pentasubstituted pyrroles and hexasubstituted pyrrolines from propargyl sulfonylamides and allenamides

Multisubstituted pyrroles are important fragments that appear in many bioactive small molecule scaffolds. Efficient synthesis of multisubstituted pyrroles with different substituents from easily accessible starting materials is challenging. Herein, we describe a metal-free method for the preparation of pentasubstituted pyrroles and hexasubstituted pyrrolines with different substituents and a free amino group by a base-promoted cascade addition–cyclization of propargylamides or allenamides with trimethylsilyl cyanide. This method would complement previous methods and support expansion of the toolbox for the synthesis of valuable, but previously inaccessible, highly substituted pyrroles and pyrrolines. Mechanistic studies to elucidate the reaction pathway have been conducted.

This method is a toolbox for the synthesis of valuable, but previously inaccessible, highly substituted pyrroles and pyrrolines.

Pyrroles are molecules of great interest in a variety of compounds including pharmaceuticals, natural products and other materials. Pyrrole fragments for example are key motifs in bioactive natural molecules, forming the subunit of heme, chlorophyll and bile pigments, and are also found in many clinical drugs, including those in Fig. 1a.¹ Although many classical methods of pyrrole synthesis, including the Paal–Knorr condensation,² the Knorr reaction,³ the Hantzsch reaction,⁴ transition metal-catalyzed reactions,⁵ and multicomponent coupling reactions,⁶ have been developed over many years, the efficient synthesis of multisubstituted pyrroles is still challenging. In condensation syntheses of pyrroles, the major problems lie in the extended syntheses of complex precursors and limited substitution patterns that are allowed. Multicomponent reactions are superior when building pyrrole core structures with more substituents. Among these, the [2+2+1] cycloaddition reaction of alkynes and primary amines is attractive because of the readily available alkyne and amine substrates and the ability to construct fully substituted pyrroles.⁷ However, with the exception of some rare examples,⁸ most [2+2+1] cycloaddition reactions afford pyrroles with two or more identical substituents. The synthesis of multisubstituted pyrroles with all different substituents from simple starting materials therefore remains a major challenge.Open in a separate windowFig. 1Previous reports and this work on propargylamides transformation.Easily accessible propargylamides are classical, privileged building blocks broadly utilized for the synthesis of a large variety of heterocyclic molecules such as pyrroles, pyridines, thiazoles, oxazoles and other relevant organic frameworks.⁹ For example, Looper¹⁰et al. reported the synthesis of 2-aminoimidazoles from propargyl cyanamides and Eycken¹¹ reported a method starting from propargyl guanidines which undergo a 5-exo-dig heterocyclization as shown in Fig. 1b. Subsequently, Wan¹²et al. revealed the cyclization of N-alkenyl propargyl sulfonamides into pyrroles via sulfonyl migration. Inspired by these transformations and multi-substituted pyrrole synthesis, we report herein a direct synthesis of pentasubstituted pyrroles and hexasubstituted pyrrolines with all different substituents from propargyl sulfonylamides and allenamides.Previously, Zhu,¹³ Ji¹⁴ and Qiu^13b,15 reported efficient syntheses of 2-aminopyrroles from isocyanides. Ye¹⁶ and Huang¹⁷ independently developed gold-catalyzed syntheses of 2-amino-pentasubstituted pyrroles with ynamides. Despite the many advantages of these methods, they all afford protected amines, rather than free amines. The deprotection of these amines may cause problems in further transformations of the products. Our method delivers pyrroles with an unprotected free amino group and are often complementary to the previously well-developed classical methods.Initially, the cyclization reaction of N-(1,3-diphenylprop-2-yn-1-yl)-N-ethylbenzenesulfonamide (1a) with trimethylsilyl cyanide (TMSCN) was carried out with Ni(PPh₃)₂Cl₂ as a catalyst, a base (Cs₂CO₃) and DMF as a solvent. Different metal catalysts, such as Ni(PPh₃)₂Cl₂, Pd(OAc)₂, Cu(OAc)₂, and Co(OAc)₂ provided the desired product with similar yields ( EntryCat.BaseSolventYield1Ni(PPh₃)₂Cl₂Cs₂CO₃DMF67%2Pd(OAc)₂Cs₂CO₃DMF65%3Cu(OAc)₂Cs₂CO₃DMF65%4Co(OAc)₂Cs₂CO₃DMF63%5Cs₂CO₃DMF65%6KFDMFTrace7K₃PO₄DMFTrace8K₂CO₃DMF48%9KOHDMF52%10KO^tBuDMF46%11Et₃NDMFTrace12Cs₂CO₃CH₃CN18%13Cs₂CO₃DME23%14Cs₂CO₃TolueneTrace15Cs₂CO₃DCETrace16Cs₂CO₃DioxaneTraceOpen in a separate window^aReaction conditions: 1a (0.1 mmol, 1 equiv.), TMSCN (0.3 mmol, 3 equiv.), cat. (0 or 10 mol%), base (0.3 mmol, 3 equiv.) and solvent (1 mL), at 80 °C for 10 h; isolated yield.With the optimal reaction conditions in hand, we investigated the scope of this reaction. As shown in Fig. 2, the transformation tolerates a broad variety of substituted propargylamides (1). The R¹ group could be an aryl group containing either electron-donating groups or electron-withdrawing groups, and the corresponding products (2b–2h) were obtained in yields of 62–80%. The substituent R¹ could also be an alkyl group such as 1-hexyl in which case the reaction provided the corresponding pyrrole (2i) in 53% yield. Exploration of the R² substituent was also conducted. Electron-rich and electron-deficient substituents in the aromatic ring of R² gave the desired products (2j–2o) with yields of 70–81%. The product bearing a furyl group (2p) can be produced in 61% yield. However, when R² group is an aliphatic group, the reaction failed to provide the desired product. Substituent groups R³, such as benzyl (2q) or 3,4-dimethoxyphenylethyl (2r) were also compatible in the reaction, providing the corresponding products in moderate yields. Significantly, this method has the potential to produce core structures (for example 2s) similar to that in Atorvastatin. Interestingly, when alkynyl substituted isoquinolines (1t–1v) were used as the substrates, the reactions smoothly afforded fused pyrrolo[2,1-α]isoquinoline derivatives (2t–2v), members of a class of compounds that are found widely in marine alkaloids and exhibit anticancer and antiviral activity.¹⁸Open in a separate windowFig. 2Substrate scope of propargylamides. Reaction conditions: 1 (0.20 mmol, 1 equiv.), TMSCN (0.60 mmol, 3 equiv.), Cs₂CO₃ (0.60 mmol, 3 equiv.) and DMF (2 mL), at 80 °C for 10 h; isolated yield.Allenes are key intermediates in the synthesis of many complex molecules.¹⁹ As a subtype of allenes, allenamines are also useful as reaction intermediates.²⁰ Although the transformation of allenamides to multisubstituted pyrroles has not been previously recorded, this reaction probably goes through the allenamide intermediates which can be derived from propargyl sulfonamides under basic conditions. To verify this hypothesis, the trisubstituted allenamide (3) was synthesized and subjected to the standard reaction conditions. A pyrrole (2a) was isolated in 82% yield from this reaction (Fig. 3). This result confirmed our assumption and raised a new question: is it possible to build hexasubstituted pyrrolines from tetrasubstituted allenamides? A range of tetrasubstituted allenamides²¹ was tested under the standard reaction conditions, and the hexasubstituted pyrrolines were obtained as is shown in Fig. 4. The R¹ group could be an aryl substituent or an alkyl chain, and the corresponding products (5a–5e) were obtained with good yields. Various aryl groups with either electron-donating groups or electron-withdrawing groups in the aromatic ring of R² provided the desired products (5f–5k) in 62–83% yields. In addition, the difluoromethyl group can also be replaced by a phenyl group, and the reaction provided the corresponding product 5l in 82% yield. It is worth noting that these pyrroline products are not easily accessible from other methods.Open in a separate windowFig. 3Synthesis of substituted pyrroles from allenes.Open in a separate windowFig. 4Substrate scope of tetrasubstituted allenamides. Reaction conditions: 4 (0.10 mmol, 1 equiv.), TMSCN (0.30 mmol, 3 equiv.), K₂CO₃ (0.30 mmol, 3 equiv.) and DMF (1 mL), at 80 °C for 10 h, isolated yield.Some synthetic applications of this method are shown in Fig. 5. The amide is a naturally occurring and ubiquitous functional group. When using benzoyl chloride to protect the free amino group of the fully-substituted pyrrole (2a), a bis-dibenzoyl amide (6) was obtained in the presence of a base, triethylamine while the monobenzoyl protected amide (7) was obtained in the presence of pyridine as the base (Fig. 5a). This method also provides a straightforward approach to pyrrole fused lactam structures (Fig. 5b). For examples, a five-membered lactam and a six-membered lactam were generated separately in a one pot reaction, directly from, (8 and 10), respectively. Taking advantage of this method, an analogue of the drug Atorvastatin was synthesized in 5 steps (Fig. 5c), demonstrating the synthetic value of the reaction.Open in a separate windowFig. 5Synthetic applications.Mechanistic experiments were performed (Fig. 6) to explore the mechanism of the reaction. When 3 equivalents of TEMPO were added, the reaction was not inhibited and the desired product (2a) was formed in 62% yield (Fig. 6a). This result suggested that the reaction might not involve a radical process. To probe the reaction further, a kinetic study was conducted (Fig. 6b). According to this study, the propargylamide (1a) was completely converted to an allenamide (3a) in 10 min under the standard conditions. The multi-substituted pyrrole (2a) was then gradually produced from the intermediate allenamide and no other reaction intermediates were observed or identified. On the other hand, DFT calculations of substrates 3b and 4a were carried out at the B3LYP-D3(SMD)/Def2-TZVP//B3LYP-D3/Def2-SVP level of theory to identify the natural bond orbital (NBO) charges on the carbons of the allene moieties. NBO charges on the internal carbon in both 3b and 4a are 0.11 and 0.18, respectively (Fig. 6c) indicating that the nucleophilic addition of cyanide anion onto the internal carbon should be reasonable as opposed to its addition onto the terminal carbon. Pathways of the cyano addition to 3b were also calculated (Fig. 6d). The transition state of cyano addition on the internal carbon (TS1), is indeed much lower than addition on the terminal carbon (TS2). The intermediate of internal carbon addition int1, is more stable than int2, implying that the internal carbon addition pathway is not only kinetically but also thermodynamically favoured.Open in a separate windowFig. 6Mechanistic studies and proposed mechanism.Based on the results of these mechanistic studies, a plausible reaction mechanism for the synthesis of pentasubstituted pyrroles and hexasubstituted pyrrolines is proposed and is shown in Fig. 6e. First, under basic conditions, the propargylamide isomerizes to an intermediate allenamide (A), which can be attacked nucleophilically by the cyanide anion to afford an intermediate imine (B) with release of the sulfonyl group. Then, the second cyanide anion attacks the imine to form an intermediate (C), which can undergo cyclization and protonation to afford the fully substituted pyrrole (2). Similarly, the hexasubstituted pyrroline product (5) can be obtained from double nucleophilic attack of the intermediate (A) by the cyanide ion.