Diradical Organic One-Dimensional Polymers Synthesized on a Metallic Surface.

We report on the synthesis and characterization of atomically precise one-dimensional diradical peripentacene polymers on a Au(111) surface. By means of high-resolution scanning probe microscopy complemented by theoretical simulations, we provide evidence of their magnetic properties, which arise from the presence of two unpaired spins at their termini. Additionally, we probe a transition of their magnetic properties related to the length of the polymer. Peripentacene dimers exhibit an antiferromagnetic (S=0) singlet ground state. They are characterized by singlet-triplet spin-flip inelastic excitations with an effective exchange coupling (Jeff ) of 2.5 meV, whereas trimers and longer peripentacene polymers reveal a paramagnetic nature and feature Kondo fingerprints at each terminus due to the unpaired spin. Our work provides access to the precise fabrication of polymers featuring diradical character which are potentially useful in carbon-based optoelectronics and spintronics.

Introduction

Over the last decades, the precise structural control of organic polymers1 has become a primary subject of research in polymer science, since new structural topologies can endow polymers with unique properties and functionalities, paving the way to the development of fields such as plastic electronics and photonics.2, 3, 4 These emerging technologies are considered as a low‐cost alternative to conventional inorganic‐based optoelectronic devices, which have had an enormous impact in science and technology starting from the second half of the last century.5, 6, 7 Based on their structure, carbon‐based polymers can be classified as single‐ (one‐dimensional), double‐ (ladder and spiro) or multi‐strand (two‐dimensional) polymers, where the elementary constitutional units monomers that are connected to each other through two, three, four or more carbon atoms, respectively. Among the plethora of existing polymeric materials, polyacenes (trans‐polyene‐edge) and polyphenanthrene (cis‐polyene‐edge), graphene nanoribbons8 (GNRs) or graphene9 represent unique model systems to study the onset and evolution of diverse exotic physical and chemical properties, which arise from the difference in polymer size, edge structure and its termination.10, 11

In this context, the synthesis of radical organic materials by tailoring their topology12 has become a major scientific quest taking into account the envisioned scientific opportunity in fields such as magnetism,13 superconductivity,14 spintronics,15 and energy conversion.16, 17 Additionally, there is a renewed interest18 in such materials triggered by the recent success in the bottom‐up synthesis of polycyclic conjugated hydrocarbons (PCHs).19, 20 The topology of synthetized PCHs plays a crucial role in their electronic structure, which can exhibit closed‐shell or open‐shell behavior.21, 22, 23, 24 Polymers composed by PCHs as constitutional units may retain the open‐shell character of the former monomers, which can lead to the emergence of unpaired spins with potential for spintronic applications.12, 25, 26 However, the high reactivity of open‐shell PCHs and polymers makes the traditional solution‐based synthesis a great challenge and, consequently, successful experimental investigations on such compounds are limited.27, 28, 29

Inspired by these advances, here we report the on‐surface synthesis and characterization of the magnetic properties of one‐dimensional polymers based on cumulene‐bridged peripentacene units, which result from on‐surface reactions of the 13,13′‐bis (dibromomethylene)‐13H,13′H‐6,6′‐bipentacenylidene molecular precursor (1) on a bare Au(111) surface under ultra‐high‐vacuum (UHV) conditions. The characterization of the polymers via scanning tunneling microscopy (STM), non‐contact atomic force microscopy (nc‐AFM) and scanning tunneling spectroscopy (STS), complemented by multidisciplinary theoretical calculations, reveals that the polymers exhibit an experimental low band gap (E g≈0.8 eV) and pure diradical character, with single radical states at each terminus.

Results and Discussion

The solution synthesis of precursor 1 was achieved from pentacene‐6,13‐dione in three synthetic steps, following previously reported procedures that consist of: (i) reduction of pentacene‐6,13‐dione to pentacene‐6(13H)‐one;30 (ii) homocoupling of the resulting ketone to yield 13H,13′H‐[6,6′‐bipentacenylidene]‐13,13′‐dione; (iii) final Ramirez olefination that provides 1.31 Compound 1 was sublimed under UHV conditions onto an atomically clean Au(111) surface held at room temperature and subsequently annealed at 180 °C to induce debromination and homocoupling, followed by cyclodehydrogenation (see Figure 1 a for the on‐surface synthetic route toward the formation of polymer (3) and Figure S1 for the overview STM image of polymer (2) obtained after the debromination step). Overview STM images reveal the predominant presence of chain‐like nanostructures that follow the FCC regions of the Au(111) herringbone reconstruction, along with sporadically distributed individual molecules (Figure 1 b). Figure 1 c shows a STM image of a prototype polymer 3 acquired with a CO‐functionalized tip recorded in the Pauli repulsion regime,32 in this case constituted by five monomers. Interestingly, the longitudinal edges of the polymer show several pronounced lobes in the current channel, absent in the polymer backbone, when recorded at low bias voltages (3 mV), which indicates an enhancement of the local density of states (LDOS).

Figure 1

On‐surface synthesis of diradical peripentacene polymers on Au(111). a) Two‐step synthetic route toward the formation of polymer 3. b) Overview STM topography image after annealing precursor 1 on Au(111) at 180 °C for 180 minutes, revealing the predominant presence of chain‐like structures. V b=0.2 V, I t=10 pA, scale bar=10 nm. c) Constant‐height ultrahigh‐resolution STM image of 3 acquired with a CO‐functionalized tip showing intramolecular features attributed to peripentacene units together with enhancements of the LDOS observed at each longitudinal polymer edge. d) Constant‐height frequency‐shift nc‐AFM image of panel (c) acquired with a CO‐functionalized tip. Open feedback parameters for (c,d): V b=3 mV. Scale bars: 1 nm. e) Graph displaying the statistics of the polymer length as a function of the annealing time of the sample at 180 °C. The inset shows exemplary STM images of a dimer, a trimer and a polymer composed of 13 peripentacene units. Open feedback parameters: V b=3 mV. f) Open‐shell non‐Kekulé structure of polymer 3. Purple filled benzenoid rings denote Clar sextets.

State of the art structural details of polymer 3 are accessible by non‐contact atomic force microscopy (nc‐AFM) measurements using a CO‐functionalized tip.33, 34 Figure 1 d shows a constant‐ height frequency‐shift image where features assigned to the benzenoid molecular backbone of each peripentacene unit are clearly discerned. Remarkably, the linkage between peripentacene monomers appears as a sharp line with a homogeneous contrast, differently to the enhanced contrast at the central positions observed for triple bonds in nc‐AFM images.35, 36, 37 We attribute this linking bridge to three consecutive C−C double bonds, that is, cumulene‐like, with a C−C distance between cumulene bonds of 4.0 Å, in agreement with the experimentally measured length of 4.1±0.4 Å, as previously reported in recent studies.38, 39, 40 Hereby, the role of Kekulé and non‐Kekulé structures in the radical character of peripentacene units should be discussed. The presence of cumulene‐like linkages renders 3 a quinoidal character, maximizing the number of Clar sextets to five per peripentacene unit in the open shell non‐Kekulé structure (Figure 1 f) leading to a gain in the aromatic stabilization energy.41 On the contrary, the closed‐shell Kekulé structure presents only two migrating Clar sextets (see Figure S2 and explanation therein). The synthesis of substantially long polymer chains is critical in view of employing these materials as active components in electronic devices.42 We observe an increase in the polymer length with increased annealing times, as depicted in the graph in Figure 1 e. The difference in the average polymer length between 5 and 180 minutes of annealing time is ≈10 nm (i.e. the average length after 5 minutes is 4.5 nm and after 180 minutes is 14.5 nm, starting from fresh samples; see Figure S3 for overview STM images of samples prepared with different annealing times). Following this recipe, the termini of every polymer lose formally their :CBr2 moieties for all annealing times, being mostly passivated by a hydrogen atom (see Figure 1 d and Figure S4).

Next, we probed the electronic structure of 3 on Au(111) via STS measurements. The differential conductance dI/dV spectra show intense peaks in the density of states with onsets at around −0.4 V and +0.4 V, which we assign to the valence band maximum (VBM) and the conduction band minimum (CBM), respectively, corresponding to an STS‐band gap of ≈0.8 eV (Figure 2 a). The spatial distribution of the CB and the VB is revealed by mapping the dI/dV signal at their corresponding energies (Figure 2 b). The characteristic appearance of the VB in the middle of the polymer features lobes along the transversal edge of each peripentacene molecule, while the CB displays the same number of lobes along the transversal molecular edge together with bowtie shape features localized over the cumulene linkages. Simulated dI/dV maps using Density Functional Theory (DFT)‐calculated frontier molecular orbitals with a CO‐tip43 of a free‐standing polymer match very well the experimental results. Notably, the experimentally measured band gap is almost independent of the polymer length (see Figure S5). This is due to the non‐dispersive character of the flat localized bands revealed by the DFT band structure calculation for a free‐standing infinite polymer 3 (Figure 2 c).

Figure 2

Characterization of the electronic structure of polymer 3 on Au(111). a) dI/dV spectra acquired on a decamer (purple and green curves) at the positions marked with purple and green stars, in the inset STM image. The orange curve corresponds to the reference dI/dV spectrum acquired on Au(111). V b=0.2 V, I t=100 pA. b) Constant‐current differential conductance (dI/dV) maps (upper panels) and corresponding DFT calculated dI/dV maps (bottom panels) at the energetic positions corresponding to the onsets of the VB (left) and the CB (right). Tunneling parameters for the dI/dV maps: VB (V b=−0.4 V, I t=400 pA); CB (V b=0.7 V, I t=400 pA). c) Calculated band structure and PDOS of free‐standing infinite polymer 3.

To understand the origin of the intense features observed systematically at the longitudinal edges of 3 (Figure 1 c), we have recorded differential conductance dI/dV spectra on different oligomers at low bias voltages. Figure 3 a depicts such features observed for a trimer, which show pronounced zero‐bias peaks. Due to their characteristic temperature dependence, we attributed these peaks to a Kondo resonance,44 stemming from a localized magnetic moment at each edge screened by the conduction electrons of the Au(111) surface. The energy width of the Kondo resonances increases with temperature, following the characteristic trend of a Kondo screened state with a Kondo temperature of T K=30±1 K (see Figure S6 for the evolution of the peak width with increasing temperature for a hexamer).

Figure 3

Characterization of the magnetic properties of polymer 3 on Au(111). a) Zero‐bias features acquired at a trimer edges assigned to Kondo resonances (inset STM image). The orange curve corresponds to the reference dI/dV spectrum acquired on Au(111). V b=0.1 V, I t=100 pA. b) Split‐peak dI/dV features around zero bias acquired at a dimer edges assigned to singlet‐triplet inelastic excitations (inset STM image). The orange curve corresponds to the reference dI/dV spectrum measured on Au(111). V b=0.05 V, I t=100 pA. c) Schematic unrestricted spin diagram and corresponding frontier orbitals of a dimer plotted for the spin‐up and spin‐down configurations. d) Constant‐current dI/dV maps acquired at the approximate energies of the onsets of the VB and CB (compare to Figure 2) (top) with the corresponding DFT calculated maps (bottom). Tunneling parameters for the dI/dV maps: HOMO−1 (V b=−0.47 V, I t=300 pA); LUMO+1 (V b=0.9 V, I t=300 pA). e) Constant‐height dI/dV map acquired at the spin excitation bias voltage (top) and corresponding DFT‐calculated of the SOMO/SUMO (bottom). Tunneling parameters for the dI/dV map (V b=2 mV). Constant‐current dI/dV maps at the same bias voltage (Figure S9) do not reveal any feature, due to limitations intrinsic to the technique at biases very close to Fermi. Acquisition temperature=1.3 K.

In order to determine the coupling between the unpaired spins, we inspected the behavior of the termini as a function of the length of the polymer. For a dimer, the dI/dV spectra reveal a peak splitting around zero‐bias at ±2.5 meV (Figure 3 b). In order to assess its electronic and magnetic properties, we performed a multidisciplinary theoretical approach, which includes spin‐unrestricted calculations for a free‐standing dimer. Such methods, which comprise hybrid DFT,45 Mean Field Hubbard Tight Binding (MFH‐TB)12 and ab initio high‐level relativistic multireference configuration interaction (MRCI) calculations with CAS(8,8) active space,46 predict a singlet ground state for the system. The triplet state is higher in energy ranging from 3.3 to 53 meV depending on the employed method and parameters (see Figure S8). Namely, the quantum chemistry MRCI with CAS (8,8) active space method gives an energy difference between the singlet and triplet states for the dimer of about 3.5 and 6.7 meV at the MRCI and MRCI+Q(Pople) level, respectively, being reasonably close to the experimental values. Even though the results from the different theoretical approaches vary, they all point to a small antiferromagnetic coupling of the spins in the case of the dimer, which vanishes for the trimer and longer polymers. Thus, we attribute the split peak observed in the dI/dV spectra of the dimer to singlet‐triplet inelastic spin excitations that change the total spin of the two exchange coupled spins of the dimer induced by electrons tunneling from the metallic tip.

Usually, such singlet‐triplet inelastic excitations appear in dI/dV spectra as steps at the onset of spin excitations, when eV > J,47 from which, one can directly determine the strength of the exchange coupling J between the spins. Accordingly, for the case of the dimer, both spins are exchange‐coupled antiferromagnetically, with an experimental effective exchange parameter J eff=2.5 meV. Herein, it is important to highlight that magnetic‐coupling in a small nanographene (Clar's Goblet), triangulene dimers and in GNR junctions on metal surfaces have been recently reported, revealing how their antiferromagnetic ground state survives in contact with a metallic surface.48, 49, 50, 51 Figure 3 c shows a schematic energy diagram of frontier molecular orbitals obtained from spin unrestricted DFT calculations. The corresponding calculated dI/dV maps match well the experimental observation of the HOMO−1 and LUMO+1 (see Figure 3 d). In addition, DFT calculations predict the presence of two single occupied (SOMOs, each one hosting a single spin) and two single unoccupied molecular orbitals (SUMOs). Both SOMO and SUMO orbitals are energetically localized very close to HOMO−1 and LUMO+1, respectively. Thus, experimental dI/dV maps of a dimer acquired at energies of HOMO−1 (LUMO+1) contains contributions from the SOMO (SUMO). Such convolutions are also calculated and observed for longer oligomers (see Figure S9 for calculated and experimental dI/dV maps of a dimer and a trimer at distinct scanning tunneling conditions). Notably, the experimental dI/dV map recorded close to the spin excitation bias resembles the shape of the SOMO and SUMO (see Figure 3 e), as expected, revealing the sharp location at the edges of the oligomers of the unpaired spins (see Figure S9 for a trimer). Interestingly, the antiferromagnetic exchange coupling for the dimer appears despite the presence of the cumulene‐like bridging unit with sp‐bonding character. This indicates that the antiferromagnetic coupling is not limited to sp2‐hybridization. Comparing the above‐mentioned results, we conclude that there is a transition in magnetic properties going from antiferromagnetic dimers to paramagnetic longer oligomers and polymers. In addition to the majority of H‐terminated polymers, we observe, independently of polymer length, that the termini of the polymers can be occasionally pinned to the gold surface or terminated with H2, which induces the quenching of the spins (see Figure S4 for the characterization of a dimer of 3 with different edge termination).

Conclusion

In conclusion, we have demonstrated a synthetic pathway toward the atomically precise synthesis of one‐dimensional diradical cumulene‐bridged peripentacene polymers on a Au(111) surface. Based on precursor 1, functionalized with :CBr2 groups, polymer 3 has been synthetized as a result of on‐surface debromination and subsequent cyclodehydrogenation. The structure of polymer 3 is clearly elucidated by STM and nc‐AFM. STS measurements complemented with theoretical calculations reveal the low band gap character of the polymer, which is independent of the polymer size. Notably, the oligomers and polymers feature radical character at their terminal edges. Multidisciplinary theoretical calculations and high‐resolution STS display an antiferromagnetic ground state for peripentacene dimers, with a low exchange coupling of 2.5 meV, whereas longer oligomers and polymers behave as paramagnets, characterized by a Kondo resonance. Our work demonstrates the appearance of diradical character in π‐conjugated polymers linked by sp‐bonds stemming from the adaptation of the π‐topology to Clar's rules. We envision that our study will inspire future investigations of π‐conjugated systems with prospects in carbon‐based optoelectronics and spintronics.

Conflict of interest

The authors declare no conflict of interest.

As a service to our authors and readers, this journal provides supporting information supplied by the authors. Such materials are peer reviewed and may be re‐organized for online delivery, but are not copy‐edited or typeset. Technical support issues arising from supporting information (other than missing files) should be addressed to the authors.

Supplementary

Click here for additional data file.