Phosphorene is a mono-elemental, two-dimensional (2D) substance with outstanding, highly directional properties and a bandgap that depends on the number of layers of the material1,2,3,4,5,6,7,8. Nanoribbons, meanwhile, combine the flexibility and unidirectional properties of one-dimensional nanomaterials, the high surface area of 2D nanomaterials and the electron-confinement and edge effects of both. The structures of nanoribbons can thus lead to exceptional control over electronic band structure, the emergence of novel phenomena and unique architectures for applications5,6,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24. Phosphorene’s intrinsically anisotropic structure has motivated numerous theoretical calculations of phosphorene nanoribbons (PNRs), predicting extraordinary properties5,6,12,13,14,15,16,17,18,19,20,21,22,23,24. So far, however, discrete PNRs have not been produced. Here we present a method for creating quantities of high-quality, individual PNRs by ionic scissoring of macroscopic black phosphorus crystals. This top–down process results in stable liquid dispersions of PNRs with typical widths of 4–50 nm, predominantly single-layer thickness, measured lengths of up to 75 μm and aspect ratios of up to 1,000. The nanoribbons are atomically flat single crystals, aligned exclusively in the zigzag crystallographic orientation. The ribbons have remarkably uniform widths along their entire lengths, and are extremely flexible. These properties—together with the ease of downstream manipulation via liquid-phase methods—should enable the search for predicted exotic states6,12,13,14,17,18,19,21, and an array of applications in which PNRs have been predicted to offer transformative advantages. These applications range from thermoelectric devices to high-capacity fast-charging batteries and integrated high-speed electronic circuits6,14,15,16,20,23,24.
Access optionsAccess options
Subscribe to Journal
Get full journal access for 1 year
only $3.90 per issue
All prices are NET prices.
VAT will be added later in the checkout.
Rent or Buy article
Get time limited or full article access on ReadCube.
All prices are NET prices.
The datasets generated and/or analysed during this study are available from the corresponding author on reasonable request.
Publisher’s note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Li, L. K. et al. Black phosphorus field-effect transistors. Nat. Nanotechnol. 9, 372–377 (2014).
Qiao, J. S., Kong, X. H., Hu, Z. X., Yang, F. & Ji, W. High-mobility transport anisotropy and linear dichroism in few-layer black phosphorus. Nat. Commun. 5, 4475 (2014).
Liu, H. et al. Phosphorene: an unexplored 2D semiconductor with a high hole mobility. ACS Nano 8, 4033–4041 (2014).
Xia, F. N., Wang, H. & Jia, Y. C. Rediscovering black phosphorus as an anisotropic layered material for optoelectronics and electronics. Nat. Commun. 5, 4458 (2014).
Lee, S. et al. Anisotropic in-plane thermal conductivity of black phosphorus nanoribbons at temperatures higher than 100 K. Nat. Commun. 6, 8573 (2015).
Sorkin, V., Cai, Y., Ong, Z., Zhang, G. & Zhang, Y. W. Recent advances in the study of phosphorene and its nanostructures. Crit. Rev. Solid State Mater. Sci. 42, 1–82 (2017).
Lewis, E. A., Brent, J. R., Derby, B., Haigh, S. J. & Lewis, D. J. Solution processing of two-dimensional black phosphorus. Chem. Commun. (Camb.) 53, 1445–1458 (2017).
Li, W. F., Yang, Y. M., Zhang, G. & Zhang, Y. W. Ultrafast and directional diffusion of lithium in phosphorene for high-performance lithium-ion battery. Nano Lett. 15, 1691–1697 (2015).
Cai, J. M. et al. Atomically precise bottom-up fabrication of graphene nanoribbons. Nature 466, 470–473 (2010).
Kosynkin, D. V. et al. Longitudinal unzipping of carbon nanotubes to form graphene nanoribbons. Nature 458, 872–876 (2009).
Yagmurcukardes, M., Peeters, F. M., Senger, R. T. & Sahin, H. Nanoribbons: from fundamentals to state-of-the-art applications. Appl. Phys. Rev. 3, 041302 (2016).
Guo, H., Lu, N., Dai, J., Wu, X. & Zeng, X. C. Phosphorene nanoribbons, phosphorus nanotubes, and van der Waals multilayers. J. Phys. Chem. C 118, 14051–14059 (2014).
Carvalho, A., Rodin, A. S. & Neto, A. H. C. Phosphorene nanoribbons. Europhys. Lett. 108, 47005 (2014).
Nourbakhsh, Z. & Asgari, R. Excitons and optical spectra of phosphorene nanoribbons. Phys. Rev. B 94, 035437 (2016).
Hu, W., Lin, L., Zhang, R. Q., Yang, C. & Yang, J. L. Highly efficient photocatalytic water splitting over edge-modified phosphorene nanoribbons. J. Am. Chem. Soc. 139, 15429–15436 (2017).
Das, P. M. et al. Controlled sculpture of black phosphorus nanoribbons. ACS Nano 10, 5687–5695 (2016).
Yang, Y. R., Zhang, Z. Q., Gu, L. & Fu, H. H. Spin-dependent Seebeck effect in zigzag black phosphorene nanoribbons. RSC Advances 6, 44019–44023 (2016).
Sisakht, E. T., Fazileh, F., Zare, M. H., Zarenia, M. & Peeters, F. M. Strain-induced topological phase transition in phosphorene and in phosphorene nanoribbons. Phys. Rev. B 94, 085417 (2016).
Wu, X. H., Zhang, X. L., Wang, X. L. & Zeng, Z. Spin density waves predicted in zigzag puckered phosphorene, arsenene and antimonene nanoribbons. AIP Adv. 6, 045318 (2016).
Poljak, M. & Suligoj, T. Immunity of electronic and transport properties of phosphorene nanoribbons to edge defects. Nano Res. 9, 1723–1734 (2016).
Yang, G., Xu, S. L., Zhang, W., Ma, T. X. & Wu, C. J. Room-temperature magnetism on the zigzag edges of phosphorene nanoribbons. Phys. Rev. B 94, 075106 (2016).
Wu, Q. Y. et al. Electronic and transport properties of phosphorene nanoribbons. Phys. Rev. B 92, 035436 (2015).
Zhang, J. et al. Phosphorene nanoribbon as a promising candidate for thermoelectric applications. Sci. Rep. 4, 6452 (2014); erratum 5, 11175 (2015).
Yao, Q. et al. Theoretical prediction of phosphorene and nanoribbons as fast-charging Li ion battery anode materials. J. Phys. Chem. C 119, 6923–6928 (2015).
Feng, X. et al. High mobility anisotropic black phosphorus nanoribbon field-effect transistor. Adv. Funct. Mater. 28, 1801524 (2018).
Abellán, G. et al. Exploring the formation of black phosphorus intercalation compounds with alkali metals. Angew. Chem. Int. Ed. 56, 15267–15273 (2017).
Cheng, Y. et al. Sodium-induced reordering of atomic stacks in black phosphorus. Chem. Mater. 29, 1350–1356 (2017).
Hembram, K. P. S. S. et al. Unraveling the atomistic sodiation mechanism of black phosphorus for sodium ion batteries by first-principles calculations. J. Phys. Chem. C 119, 15041–15046 (2015).
Cullen, P. L. et al. Ionic solutions of two-dimensional materials. Nat. Chem. 9, 244–249 (2017).
Ribeiro, H. B. et al. Edge phonons in black phosphorus. Nat. Commun. 7, 12191 (2016).
Wu, R. J. et al. Atomic and electronic structure of exfoliated black phosphorus. J. Vac. Sci. Technol. A 33, 060604 (2015).
Sun, J. et al. A phosphorene-graphene hybrid material as a high-capacity anode for sodium-ion batteries. Nat. Nanotechnol. 10, 980–985 (2015).
Peng, X. H. & Wei, Q. Chemical scissors cut phosphorene nanostructures. Mater. Res. Express 1, 045041 (2014).
Mayo, M., Griffith, K. J., Pickard, C. J. & Morris, A. J. Ab initio study of phosphorus anodes for lithium- and sodium-ion batteries. Chem. Mater. 28, 2011–2021 (2016).
Hembram, K. P. S. S. et al. A comparative first-principles study of the lithiation, sodiation, and magnesiation of black phosphorus for Li-, Na-, and Mg-ion batteries. Phys. Chem. Chem. Phys. 18, 21391–21397 (2016).
Yu, X. F., Ushiyama, H. & Yamashita, K. Comparative study of sodium and lithium intercalation and diffusion mechanism in black phosphorus from first-principles simulation. Chem. Lett. 43, 1940–1942 (2014).
Scheer, M., Balazs, G. & Seitz, A. P-4 activation by main group elements and compounds. Chem. Rev. 110, 4236–4256 (2010).
Evers, E. C. The alkali metal phosphides. 1. Reaction of alkali metals and white phosphorus in liquid ammonia. J. Am. Chem. Soc. 73, 2038–2040 (1951).
Zhang, R., Waters, J., Geim, A. K. & Grigorieva, I. V. Intercalant-independent transition temperature in superconducting black phosphorus. Nat. Commun. 8, 15036 (2017).
Mao, N., Zhang, S., Wu, J., Zhang, J. & Tong, L. Lattice vibration and Raman scattering in anisotropic black phosphorus crystals. Small Methods 2, 1700409 (2018).
Ribeiro, H. B., Pimenta, M. A. & de Matos, C. J. S. Raman spectroscopy in black phosphorus. J. Raman Spectrosc. 49, 76–90 (2017).
Phaneuf-L’Heureux, A.-L. et al. Polarization-resolved raman study of bulk-like and Davydov-induced vibrational modes of exfoliated black phosphorus. Nano Lett. 16, 7761–7767 (2016).
Abellán, G. et al. Fundamental insights into the degradation and stabilization of thin layer black phosphorus. J. Am. Chem. Soc. 139, 10432–10440 (2017).
Pénicaud, A., Poulin, P., Derre, A., Anglaret, E. & Petit, P. Spontaneous dissolution of a single-wall carbon nanotube salt. J. Am. Chem. Soc. 127, 8–9 (2005).
Pénicaud, A. & Drummond, C. Deconstructing graphite: graphenide solutions. Acc. Chem. Res. 46, 129–137 (2013).
Basma, N. S., Headen, T. F., Shaffer, M. S. P., Skipper, N. T. & Howard, C. A. Local structure and polar order in liquid n-methyl-2-pyrrolidone (NMP). J. Phys. Chem. B 122, 8963–8971 (2018).
Kneisl, P. & Zondlo, J. W. Vapor pressure, liquid density, and the latent heat of vaporization as functions of temperature for four dipolar aprotic solvents. J. Chem. Eng. Data 32, 11–13 (1987).
Clancy, A. J. et al. Charged carbon nanomaterials: redox chemistries of fullerenes, carbon nanotubes, and graphenes. Chem. Rev. 118, 7363–7408 (2018).
Reinhold, R. et al. Surface and electrochemical studies on silicon diphosphide as easy-to-handle anode material for lithium-based batteries—the phosphorus path. ACS Appl. Mater. Interfaces 10, 7096–7106 (2018).
Hanlon, D. et al. Liquid exfoliation of solvent-stabilized few-layer black phosphorus for applications beyond electronics. Nat. Commun. 6, 8563 (2015).
Brent, J. R. et al. On the stability of surfactant-stabilised few-layer black phosphorus in aqueous media. RSC Advances 6, 86955–86958 (2016).
Woomer, A. H. et al. Phosphorene: synthesis, scale-up, and quantitative optical spectroscopy. ACS Nano 9, 8869–8884 (2015).
Favron, A. et al. Photooxidation and quantum confinement effects in exfoliated black phosphorus. Nat. Mater. 14, 826–832 (2015).
Abate, Y. et al. Recent progress on stability and passivation of black phosphorus. Adv. Mater. 30, 1704749 (2018).
Zhang, T. et al. Degradation chemistry and stabilization of exfoliated few-layer black phosphorus in water. J. Am. Chem. Soc. 140, 7561–7567 (2018).
Zhou, Q., Chen, Q., Tong, Y. & Wang, J. Light-induced ambient degradation of few-layer black phosphorus: mechanism and protection. Angew. Chem. Int. Ed. 55, 11437–11441 (2016).
Edmonds, M. T. et al. Creating a stable oxide at the surface of black phosphorus. ACS Appl. Mater. Interfaces 7, 14557–14562 (2015).
Wood, J. D. et al. Effective passivation of exfoliated black phosphorus transistors against ambient degradation. Nano Lett. 14, 6964–6970 (2014).
Mikheikin, A. et al. DNA nanomapping using CRISPR-Cas9 as a programmable nanoparticle. Nat. Commun. 8, 1665 (2017).
Miller, T. S. et al. Single crystal, luminescent carbon nitride nanosheets formed by spontaneous dissolution. Nano Lett. 17, 5891–5896 (2017).
Picco, L. M. et al. Breaking the speed limit with atomic force microscopy. Nanotechnology 18, 044030 (2007).
Lui, C. H., Liu, L., Mak, K. F., Flynn, G. W. & Heinz, T. F. Ultraflat graphene. Nature 462, 339–341 (2009).
Shin, C. Variation-Aware Advanced CMOS Devices and SRAM Vol. 56 (Springer, Dordrecht, 2016).
Canet-Ferrer, P., Eugenio, C., Alicia, F.-A. & Elena, P.-C. Correction of the tip convolution effects in the imaging of nanostructures studied through scanning force microscopy. Nanotechnology 25, 395703 (2014).
Wang, X. et al. Highly anisotropic and robust excitons in monolayer black phosphorus. Nat. Nanotechnol. 10, 517–521 (2015).
This publication was funded, in part, by the Engineering and Physical Sciences Research Council (EPSRC). L.P. and O.D.P. thank the Royal Academy of Engineering for funding their time and for support, and also O. Ward for help with analysing the HS-AFM data. T.S.M. thanks the EPSRC for support via fellowship EP/P023851/1. P.L.C. thanks the EPSRC for support via fellowship EP/S001298/1. The HS-AFM data were collected using the University of Bristol’s Low Noise Labs. We thank M. Wood for early experiments, and A. Sella and A. Pakpour-Tabrizi for enlightening discussions.
Nature thanks David J. Lewis and the other anonymous reviewer(s) for their contribution to the peer review of this work.
Extended data figures and tables
Black phosphorus crystals were outgassed in a fused silica reaction tube, and lithium metal was added. The tube was cooled to −50 °C and high-purity ammonia (NH3) was condensed onto the black phosphorus and lithium, dissolving the alkali metal to form a Li:NH3 solution. The black phosphorus was left submerged in the Li:NH3 solution for about 24 hours, allowing intercalation, before the ammonia was slowly removed. Without air exposure, the intercalated black phosphorus crystal was transferred to a glovebox and placed into a vial, and NMP was added. This vial was then placed in an ultrasonic bath for 1 hour and centrifuged at low acceleration. This process resulted in solutions of phosphorene nanoribbons.
TEM micrographs showing PNRs (darker, ribbon-like features) deposited on ‘holey carbon’ TEM grids (lighter, webbed structure). a, Heterogeneous PNRs of different widths. b, Micrographs of high-aspect-ratio (greater than 100) PNRs with widths less than 5 nm. Dashed boxes demarcate composite images; solid boxes correspond to the extracted zoom insets. To enhance visibility, selected PNRs have been outlined by dashed orange lines; otherwise, red arrows indicate the locations of the ribbon. c, Images of PNRs that demonstrate their flexibility. Composite TEM images are demarcated by dashed black outlines.
a–f, SAED patterns (scale bar, 10 nm−1) taken at the locations corresponding to the boxes marked in the composite TEM image in g (scale bar, 1 μm). These images show the zigzag crystallographic alignment along the length of the PNR, as indicated by the red arrows in a–g. h, Indexed FFT transform of the HRTEM micrograph shown in Fig. 1h. Scale bar, 10 nm−1. i–k, STEM (i; scale bar, 2.5 µm) and TEM (j; scale bar, 2 µm) images of the same PNR, for which EELS spectra (l, m) were recorded. k, Left, closer view (scale bar, 1 µm) of the PNR shown in i; right; spectrum image chosen at a suspended-in-vacuum part of the ribbon, with the red square indicating the region from which the EEL spectra shown in l and m were averaged. l, Core-loss EELS at the phosphorus L2,3 edge (130 eV), showing features in regions I, II and III that correspond to features found in pristine few-layer black phosphorus31. m, Low-loss EELS spectrum, showing the expected plasmonic signature peak at roughly 19 eV and another feature at around 3.3 eV, the latter being similar to that found for few-layer black phosphorus and attributed to surface plasmons31. The broad feature at around 5 eV that is marked with an asterisk is from delocalized signal coming from the amorphous carbon support.
Extended Data Fig. 4 Heights and roughnesses of one- to five-layer PNRs from statistical analysis of HS-AFM images.
See Methods. Normalized histograms (black data points; error bars are proportional to the square root of the number of counts in each bin) and Gaussian fits (solid blue line) of independent, automated, pixel-height measurements from the HOPG substrate (far left) and for PNRs of one to five layers. Extracted (see Methods) average heights (z) and roughnesses (σ) are shown at the top. The numbers of individual pixel-height measurements (taken over multiple locations, PNRs and substrates) are: HOPG substrate, 886,908 measurements; one-layer PNRs, 61,181; two-layer PNRs, 50,637; three-layer PNRs, 18,740; four-layer PNRs, 78,514; and five-layer PNRs, 28,816. Scaled representations (red) of the crystal structures of the PNRs are overlaid.
a, Magnified image of the PNR shown in Fig. 2e. Scale bar, 5 µm. b, Magnified image of the bifurcation seen in a. Scale bar, 500 nm. c, d, Linecuts taken from b, before (c) and after (d) the bifurcation, demonstrating how the ribbon self-cleaves from its original height of roughly 3.5 nm (seven layers) to two branches with approximate heights of 1.5 nm (three layers) and 2 nm (four layers). e, f, HS-AFM images (scale bars, 200 nm and 500 nm, respectively) of areas over which the stability study in Extended Data Fig. 10a–f was undertaken. In e, the black box corresponds to a HOPG trench used to correct the stability data for tip effects (see Methods), and the blue box corresponds to the bilayer PNR in Extended Data Fig. 10d. In f, the blue box corresponds to the monolayer ribbon in Extended Data Fig. 10a.
Length versus width scatterplots of deposited PNRs, extracted from TEM data (Methods). The plots show the PNRs that have spontaneously dissolved from Li-BP crystals into NMP solvent (red circles) as well as those formed via bath sonication (blue circles). The magnified plot on the right corresponds to the dashed rectangle in the left-hand image and reveals that the spontaneously dissolved PNRs are on average wider than those formed by sonication.
a, XPS maps of the platinum-coated positive (left) and negative (right) terminal electrodes following PNR electrodeposition. The maps were constructed from the integrated area under those phosphorus 2p peaks assigned to the P–P environment. Scale bars, 500 µm. b, Schematic representations of the platinum pattern (grey) on the positive (left) and negative (right) electrodes. c, An example phosphorus 2p spectrum from a, fitted with mixed Gaussian/Lorentzian peaks assigned to P–P bonding and P–O surface contamination. d, Patterned electrodes, showing the grid over which the Raman maps in Fig. 3e, f were taken. Scale bar, 20 μm. e, Example Raman spectra. The blue and red dots in d show the specific locations where the example Raman were taken.
a, Example Raman spectrum (black circles) taken from PNRs deposited on a graphite substrate in an argon atmosphere (Methods). The spectrum is fitted with a sum (red curve) of three Lorentzian functions (blue curves) corresponding to the Ag1, B2g and Ag2 peaks of black phosphorus40,41. b–d, Histograms showing the positions (top) and widths (bottom) of the Raman modes. e, Integrated intensity (I) ratios of the Ag1/Ag2 peak (top) and the B2g/Ag2 peak (bottom), extracted by fitting 351 independent spectra taken from the same sample of solution-deposited PNRs. From the fits to the distributions in b–d, we determine the mean values (and standard deviations) of the Raman peak positions and widths to be: Ag1 position 362.5 (0.2) cm−1, width 3.3 (0.2) cm−1; B2g position 439.3 (0.2) cm−1, width 5.4 (0.3) cm−1; Ag2 position 467.5 (0.2) cm−1, width 3.7 (0.2) cm−1. f, Photoluminescence (PL) spectrum from PNR aggregates deposited on graphite. High-energy peaks are Raman features from the PNRs or the graphite substrate (the latter marked by asterisks). Mechanically exfoliated single-layer phosphorene emits at around 1.3 eV (ref. 65). g, UV/Vis absorption spectra from a PNR solution (shown in the inset). Red circles mark the positions of weak features below 500 nm.
a, b, Example Raman spectra for Li-BP, recorded with the incident laser polarization perpendicular (a) or parallel (b) to the vertical features that are visible in the bottom part of c. The top two spectra of a, b show the results obtained in argon; the bottom two spectra were obtained from the same region in the sample after 30 days in air. These spectra were fitted with a sum (red curve) of Lorentzian functions for black phosphorus (blue curves) and intercalation (green curves) Raman modes. See Methods. c, Top, optical micrograph of an Li-BP crystal (scale bar, 200 µm); bottom, a magnified image (scale bar, 18 µm). The bottom image also shows the area (red box) over which the Raman mapping data in f were taken from. d, Histograms showing the Li/P Ag2 peak position in argon (values taken from 396 individual spectra) and after 30 days in air (taken from 324 spectra). e, Histograms showing the Ag1/Ag2 peak intensity ratio for measurements performed with the incident laser polarization parallel (taken from 180 spectra) or perpendicular (taken from 396 spectra) to the vertical stripes in c. f, Map showing the position of the Ag2 peak for data taken over the red box in c. Scale bar, 5 µm. g, TEM image showing in-plane splitting of PNRs from a larger fragment. Scale bar, 500 nm. h, HS-AFM image (top; scale bar, 100 nm) and associated linescan (bottom) of the surface of an Li-BP crystal, showing a ribbon (line profile 2, black) liberated on the surface and the trench left behind (line profile 1, red).
a, HS-AFM images of the monolayer PNR shown in Extended Data Fig. 5f, taken immediately (top) and after six days of air exposure (bottom). Scale bars, 50 nm. b, Normalized histograms showing the width (left) and height (right) of the PNR in a, taken immediately (blue) and after six days in air (yellow). Each histogram has been fitted with a Gaussian function (red). c, Width (left) and height (right) of the PNR during six days in air, taken from fits to the histograms as shown in b. Data are means ± s.d. from the fit. For the histograms in b, the numbers of individual measurements are: 2,086 (day 0 width); 16,565 (day 6 width); 45,850 (day 0 height), and 281,980 (day 6 height). d, HS-AFM images of the bilayer PNR shown in in Extended Data Fig. 5e. Scale bars, 50 nm. e, f, As for b, c, but taken for the bilayer ribbon in d. For the histograms in e, numbers of individual measurements are: 8,343 (day 0 width); 4,169 (day 6 width); 65,614 (day 0 height), and 40,520 (day 6 height). g, Evolution of the Raman Ag1/Ag2 integrated intensity ratio, the Ag2/background intensity ratio, Ag2 position and Ag2 width (FWHM, full width at half-maximum), extracted by analysing maps of 90–100 independent spectra taken over the same location from a sample of PNRs drop-cast onto graphite (data points are means ± s.d. taken from fits; see Methods). h, Evolution of XPS spectra (black dots, data; red curves, fits to P peaks; blue curves, fits to P–O peaks) taken from a sample of PNRs drop-cast onto graphite as a function of increasing air exposure. i, Normalized absorbance at 465 nm (black dots) as a function of time for the same solution of PNRs in NMP over 50 days without exposure to air. A photograph of the solution is in the inset. The red curve is an exponential function fit to the data (see Methods).
A flyover visualisation of HS-AFM data taken of a phosphorene nanoribbon. Made by manually tracking a PNR over a 22 × 67 µm2 area using HS-AFM, this 3D flyover highlights the PNR’s compliance to the HOPG substrate and the relative heights changes of both the HOPG substrate and the PNR.