Magnetism engineering of nanographene: An enrichment strategy by co-depositing diverse precursors on Au(111)

Hui Zhng,Jinhen Lu,Yong Zhng,Lei Go,Xin-Jing Zho,Yun-Zhi Tn,Jinming Ci,*

a Faculty of Materials Science and Engineering,Kunming University of Science and Technology,Kunming 650093,China

b Faculty of Science,Kunming University of Science and Technology,Kunming 650500,China

c State Key Laboratory for Physical Chemistry of Solid Surfaces,Collaborative Innovation Center of Chemistry for Energy Materials,and Department of Chemistry,College of Chemistry and Chemical Engineering,Xiamen University,Xiamen 361005,China

Keywords:Nanographene Magnetism enrichment Co-deposition Scanning tunneling microscope

ABSTRACT The magnetism of nanographene is dominated by the structure of its carbon skeleton.However,the magnetism engineering of nanographene is hindered due to finite precursors.Here,we demonstrate an ingenious synthetic strategy to engineer the magnetism of nanographene through hetero-coupling two precursors on Au(111) surface.Bond-resolved scanning tunneling microscopy and spectroscopy results show that two homo-coupled products host a closed-shell structure,while the products with five membered ring defects perform as an open-shell one with the total spin number of 1/2,confirmed by spin-polarized density functional theory calculations.While two hetero precursors on Au(111) substrate,the heterocoupled products both perform as the magnetic structure with total spin quantum numbers of 1/2 and 1,resulting from carbon skeleton transformations.Our work provides an effective way to engineer the magnetism of nanographene by enriching the magnetic products simultaneous,which could be extended into other controllable magnetic nanographene instruction.

The magnetism of carbon nanomaterials,such as nanographene(NG),has attracted extensive attention due to its high magnitudes of spin-wave stiffness [1],weak spin-orbit coupling [2],hyperfine couplings [3] and large spin coherence lifetimes [4],which hold promise for spintronics [5,6].The open-shell [7,8] structures of magnetic NGs with high reactivity results from the presence of unpair electrons (radicals),which make it challenging in solutionbased synthesis and magnetic ground state measurement [9,10].However,as an alternative method,on-surface synthesis under ultra-high vacuum has emerged as a powerful flexible synthetic toolkit in recent years [11].Current research focus of carbon magnetism inducing involve four aspects: NG or graphene nanoribbons (GNR) with spin-polarized zigzag edges states [12–15],sublattice imbalance as guided by Lieb’s theorem for bipartite lattices or defects induced systems [16–20],NGs with non-Kekulé structure formed by the topological frustrated [21] or size-induced Coulomb repulsion between valence electrons [22–26],and non-benzenoid rings imbedded NG [27,28].

The total spin quantum number (S) of NGs might be engineered in four ways.Firstly,artificially precursors transformed into open-shell structure by directly cyclodehydrogenation [16,17,19,21–24,29–33] to induce the magnetism.Secondly,the radical positions of NGs could be passivated by H atoms released from cyclodehydrogenation,which will modulateSfrom 0 to 1/2 [21] and 1 to 1/2[19].In this situation,theScould be recovered by applying bias to remove the redundant H atoms.Thirdly,quenching the spin by creating single C–metal bond,for example,bounding the radical positions to the elbow area of Au(111) surface [21] by applying positive voltage pulse [34].Lastly,the spin on/off could be switched due to the charge transfer by changing adsorption positions [35–37].Recently,an ingenious co-deposition method had been employed to induce magnetism in boron atoms doped graphene nanoribbon[38].

Fig.1.Synthetic routes toward NGs 1a,1b,2a and 2b.Route 1: precursor 1 undergoes homo-coupling into the defect-free NG 1a and defective 1b.Route 2: precursors 1 and 2 hetero-coupling into the defect-free NG 2a and defective 2b.

Here,we induce the co-deposition method into the engineering of nanographene and synthesize plenty products with different magnetic ground states by homo-coupling and hetero-coupling between precursors.We demonstrate this engineering strategy by co-depositing two precursors (precursor 1,9–bromo-10-(2,6-dimethylphenyl)anthracene,BDMPA [23]; precursor 2,10–bromo-9,9′-bianthracene,BBA [39]) on Au(111) substrate.As illustrated in route 1 of Fig.1,surface-catalyzed Ullmann homo-coupling and afterwards cyclodehydrogenation reactions (300°C,6°C/min) of precursor 1 lead to the formation of NG 1a (defect-free).Simultaneously,a NG 1b endowed with a pentagonal terminus resulting from losing a methyl group in the precursor 1 could also be formed during the annealing process.According to chemical structures and Lieb’s theorem [20],the total spin of 1a is 0.The pentagonal ring in NG 1b renders a non-Kekulé system with an unpaired electron,that means aS=1/2 grounds state of NG 1b.NG 1a can be described as either an open-shell non-Kekulé structure with 5 Clar sextets or a closed-shell Kekulé structure with 4 Clar sextets,respectively.This competition from the aromaticity of 5 Clar sextets is higher than 4 Clar sextets one,while the open-shell structure is more reactive than the closed-shell one.In addition,DFT calculation show that the real structure of NG 1a is still ambiguous(Table S1 in Supporting information) as the energies of open-shell and closed-shell structures are almost the same.When missing a methyl group in the precursor 1,a pentagonal topological defect in NG 1b raises a result of sublattice-imbalance,which manifests as a Kondo resonance.In another way,route 2,surface-catalyzed Ullmann hetero-coupling of the precursor 1 and the precursor 2 and afterwards cyclodehydrogenation cyclization reactions (300°C,6°C/min) lead to formations of NG 2a and NG 2b.NG 2a,a defectfree structure,possesses a triplet ground state (S=1).NG 2b,similar with NG 1b,formed by missing a methyl group,has a doublet ground state (S=1/2).

After depositing the precursor 1 on Au(111) surface and subsequent cyclodehydrogenation,the NGs 1a and 1b were obtained(Fig.2).Figs.2a and b show the STM images of NGs 1a and 1b,respectively.The NG 1a can be clearly recognized from the STM imageviathe two-lobe-like feature at both triangular termini (Fig.2a),while the NG 1b possesses such feature only at one of the triangular ends (Fig.2b).The structures are verified by our BR-STM images and corresponding chemical structure displayed in Figs.2c and d.The additional states are visible near non-defect area in Fig.2d.To elucidate the experimental electronic properties,we employed STS characterizations of the NG 1a (Fig.2e) and 1b (Fig.2f).For the NG 1a,the STS depicts several electronic resonances labeled as HOMO-1 (–1.35 V),LUMO+1 (+2.0 V) and a pair of electronic resonances marked as frontier molecular orbitals.The absence of any spin enhanced state near Fermi level from BR-STM as well as the inelastic spin excitation signal in a higher energy resolution dI/dVspectroscopy (Fig.S1 in Supporting information) alone with calculated energy (Table S1) which may illustrate a closedshell structure of NG 1a.The precise energy positions of LUMO and HOMO are identified to be–0.3 V and +0.21 V,respectively,which reveals a frontier gap of 0.51 eV for NG 1a.The closed-shell nature of NG 1a by the frontier maps (Fig.S1b in Supporting information) with different shapes and symmetries.The spatial maps of electronic state distribution of HOMO-1 (Fig.2g),and LUMO+1(Fig.2h) reveal similar results as reported before [18].For NG 1b,the STS mainly shows two obvious electronic resonances in Fig.S2b (Supporting information),and labeled as HOMO-1 (–1.1 V) and LUMO+1 (+1.6 V) in Fig.2f.Figs.2i and j show the spatial maps of electronic states distribution of HOMO-1 and LUMO+1 and the additional maps are displayed in Fig.S3 (Supporting information).The dI/dVspectra in Fig.2f depict the zero-energy states (ZESs)near the defect free terminus (the purple and green lines in Fig.2f) which correspond to the additional states in BR-STM (Fig.2d).To amplify on the ZESs of NG 1b,we performed low-energy range dI/dVspectra in Fig.S4 (Supporting information).The ZES appears when the tip position is close to defect free area and the resonance fits well with the Frota function [40].These results indicate that the ZES in NG 1b is attributed to the Kondo resonances[16–19,30,41].

Fig.2.STM images and electronic structures of NGs 1a and 1b.(a,b) STM images of NG 1a (c) and 1b (d) on Au(111) (V=-500 mV, I=200 pA,3×3 nm2).(c,d) BRSTM images (left panel) and Laplace filtered images (right panel) of NGs 1a and 1b(V=2 mV,2×2 nm2).dI/dV spectroscopy on NG 1a (e) and 1b (f) (V=-2 V, I=350 pA, Vrms=15 mV).(g,h) Experimental dI/dV maps of HOMO-1 and LUMO+1 for NG 1a (I=350 pA, Vrms=15 mV,3×3 nm2).(i,j) Experimental dI/dV maps of HOMO-1 and LUMO+1 for NG 1b (I=350 pA, Vrms=15 mV,3×3 nm2).

To engineer the NG magnetism,we co-deposited the precursor 1 and the precursor 2 onto the Au(111) substrate.After annealing process,two new products (NGs 2a and 2b) were observed,as shown in Fig.3a.NGs 2a and 2b are marked by solid white rectangles (Fig.3a).The difference between NGs 2a and 2b can be distinguished from zoom-in highly-resolved STM images (Figs.3b and c).Two-lobe like local density state distribution shown in Fig.3b corresponds to the triangular end of NG 2a,while the pentagonal defect terminus of NG 2b shows a faint feature (Fig.3c).To gain bond-resolved images of NGs 2a and 2b,we performed high resolution STM characterizations by attaching a CO molecule at the apex of tip.Figs.3d and e depict bond-resolved STM (BR-STM)images of NGs 2a and 2b,respectively,unambiguously confirming the defect-free feature of NG 2a and the pentagonal defect feature of NG 2b.Both BR-STM images show the feature with additional states,indicating the Kondo resonance enhanced state distributions[18,19].Moreover,the BR-STM image for NG 2a shows the feature with additional state around whole structure,while only the bottom part of BR-STM image for NG 2b shows the feature with additional state.The difference means that the Kondo resonance distributions of NG 2a and NG 2b [35,42–44].dI/dVspectra of NGs 2a and 2b shown in Figs.3f and g shed light on the ZESs existence as well as energy positions of HOMO-1 and LUMO+1.The distributions of ZESs for NGs 2a and 2b fit well with the electronic-state enhanced area of BR-STM images in Figs.3d and e.dI/dVmaps at–1.3 V,and +1.7 V for NG 2a (Figs.3h and i) and–1.0 V,and +1.4 V for NG 2b (Figs.3j and k) reveal spatial distributions of electronic states.

Fig.3.STM images and electronic structures of NGs 2a and 2b.(a) Large-scale STM image of NGs 2a and 2b on Au(111) (V=-1 V, I=30 pA.40×40 nm2).Highresolution STM images of NGs 2a (b) and 2b (c) (V=2 mV, I=200 pA.3×3 nm2).BR-STM images (left panel) with CO functionalized tip for NG 2a (d) and 2b (e)and corresponding Laplace-filtered images (V=2 mV,2.5×2.5 nm2).dI/dV spectra of NGs 2a (f) and 2b (g) (V=–2.0 V, I=350 pA).(h,i) Experimental dI/dV maps of HOMO-1 and LUMO+1 of NG 2a (I=350 pA, Vrms=15 mV).(j,k) Experimental dI/dV maps of HOMO-1 and LUMO+1 of NG 2b (I=350 pA,3×3 nm2).

To further investigate the electronic structures and magnetic properties of NGs 2a and 2b,we performed SP-DFT calculations.Compared with antiferromagnetic (AFM),the NG 2a is further verified to be ferromagnetic (FM) as preferable energy (Table S2 in Supporting information).Figs.4a and b present the SP-DFT calculated energy spectra of NGs 2a and 2b,which show single electron occupied frontier orbitals SOMOs and SUMOs of the nondegenerated energy levelψ2andψ3of NG 2a and the energy levelψ2of NG 2b respectively,triggering spin polarizations.According to the calculated energy level,magnetic ground states of NGs 2a and 2b areS=1 andS=1/2 respectively.The calculated spin-polarized wave functions of the SOMOs of NGs 2a and 2b are displayed in Figs.4c and d,respectively.ψ2↑andψ3↑of the NG 2a exhibit dominant localization at the top and the bottom terminus with slight overlap.ψ2↑of the NG 2b only exhibit dominant localization at the bottom area.The spin density distribution of the triplet ground state of NG 2a shown in Fig.4e is around the whole CBatoms,while the doublet ground state of NG 2b in Fig.4f is around at the bottom CBatoms of the NG 2b.Furthermore,the calculated spin density distributions fit well with the experimental BR-STM results (Figs.3d and e).Low-energy interval dI/dVspectra at the marked positions (Figs.4g and h) reveal the existence of Kondo resonances.The distributions of Kondo resonance fit well with the calculate spin density state distributed area.The Kondo peak positions are slight away from Fermi level and the shapes are slightly asymmetric,which can be attributed to the charge transfer and the hybridization with the substrate [35,36,45].As shown in Fig.4g,compared with the terminal ends,the Kondo peak acquired in middle positions of NG 2a are positively shifted (Table S3 in Supporting information),which may indicate a stronger charge transfer area at the middle part of NG 2a.

Fig.4.SP-DFT calculated magnetic properties of NGs 2a and 2b.(a,b) SP-DFT calculated energy spectra of NGs 2a and 2b.DFT calculated spin-polarized wave functions of the SOMO of NG 2a (c) and 2b (d).SP-DFT simulated spin density distributions of NG 2a (e) and 2b (f).gray and black balls represent carbon atoms from two different sublattice.The spin density marked by cyan and yellow represent two different spin directions.Experimental low-energy ranges dI/dV spectra of the Kondo resonances detected in 2a (g) and 2b (h),with fitting using the Frota function [40].

There are two obvious additional details,when taking a deeper analysis of the Kondo resonance by comparing NGs 2a and 2b.Firstly,the Kondo resonance acquired at both ends of NG 2a and around NG 2b shows obvious Kondo state difference,where NG 2a possesses smaller resonance amplitude and linewidth corresponding to the underscreen effect ofS=1 system [19,46] and NG 2b shows a complete wash out of the magnetic moment ofS=1/2.Secondly,the spatial variation of the Kondo peak when STS is measured at different lateral positions on the NG 2a which is different from a shorter NG 3 (Fig.S6 in Supporting information).Combining these Kondo resonance differences and calculated wave function results,we concluded two opening explanations for the NG 2a.As the longer distance of two spins compared with NG 3,they have negligible direct coupling leading theS=1/2 Kondo feature[19,21] in the middle area and indirect spin coupling or exchanging at both ends of NG 2a leading the under screened Kondo feature.Another possibility is that Kondo peak differences are result from relative position differences corresponding to the spin center and underneath gold substrate [47].

We demonstrate an efficient synthetic strategy to engineer and enrich the magnetic NG products by co-depositing diverse precursors on the Au(111) substrate.Two homo-coupled products (NG 1a and NG 1b) acquired by Ullmann coupling and cyclodehydrogenation reactions of precursor 1 and two hetero-coupled products (NG 2a and NG 2b) co-deposited by precursor 1 and precursor 2 have been successfully fabricated.Defective NG 1b hosts spins ofS=1/2 resulting from a sublattice imbalanced carbon skeleton.The magnetism for defect-free NG 2a and defective NG 2b host spins ofS=1/2 andS=1,respectively.Our work successfully realizes the magnetism engineering and enrichment by inducing another precursor and offers a valid opportunity via the co-deposition to extend the magnetism engineering of NG on metal substrate.

Declaration of competing interest

The authors declare no competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (No.61901200),the National Recruitment Program for Young Professionals (No.132310976002),the Yunnan Fundamental Research Projects (Nos.2019FD041,202101AV070008,202101AW070010 and 202101AU070043),the Strategic Priority Research Program of Chinese Academy of Sciences (No.XDB30000000),and the Dongguan Innovation Research Team Program.The numerical calculations in this paper have been done on Hefei advanced computing center.

Supplementary materials

Supplementary material associated with this article can be found,in the online version,at doi:10.1016/j.cclet.2022.04.048.