## Enhanced ambient stability of exfoliated black phosphorus by passivation with nickel nanoparticles

Maria Caporali, a* Manuel Serrano-Ruiz,  Francesca Telesio,  Stefan Heun,  Alberto Verdini,  Albano a b b c Cossaro,  Matteo Dal Miglio,  Andrea Goldoni,  Maurizio Peruzzini c d d a*

a CNR ICCOM, Via Madonna del Piano 10, 50019 Sesto Fiorentino, Italy.

b NEST, Istituto Nanoscienze-CNR and Scuola Normale Superiore, Piazza S. Silvestro 12, 56127 Pisa, Italy.

c TASC   CNR IOM, Area Science Park, Basovizza, 34149 Trieste, Italy.

d Elettra-Sincrotrone Trieste, Area Science Park, S.S. 14 km 163.5, I-34149 Trieste, Italy.

E-mail : maria.caporali@iccom.cnr.it, mperuzzini@iccom.cnr.it

Keywords : phosphorene, nickel nanoparticles, passivation, Raman, high resolution XPS

## Abstract

Since its discovery, the environmental instability of exfoliated black phosphorus (2D bP) has emerged as a challenge that hampers its wide application in chemistry, physics, and materials science. Many studies have been carried out to overcome this drawback. Here we show a relevant enhancement of ambient stability in few-layer bP decorated with nickel nanoparticles as compared to pristine bP. In detail, the behavior of the Ni-functionalized material exposed to ambient conditions in the dark is accurately  studied  by  TEM  (Transmission  Electron  Microscopy),  Raman  Spectroscopy,  and  high resolution X-ray Photoemission and Absorption Spectroscopy. These  techniques provide a morphological and quantitative insight of the oxidation process taking place at the surface of the bP flakes. In the presence of Ni NPs, the decay time of 2D bP to phosphorus oxides is more than three time slower compared to pristine bP, demonstrating an improved structural stability within twenty months of observation.

## Introduction

Within  the  family  of  two-dimensional  (2D)  materials,  exfoliated  black  phosphorus  (2D  bP)  has emerged as a hot and highly intriguing research topic over the last few  years [1]. This novel 2D material is endowed with a unique structure and outstanding physical properties, thus it has a broad prospect in the construction of electronic devices [2]. In this respect, a prominent feature of 2D bP is that it is a semiconductor with a thickness-dependent direct band gap, ranging from 0.3 eV (bulk material) to approximately 2.0 eV (monolayer, called phosphorene) [3, 4]. Besides, the high mobility

of charge carriers makes it a solid candidate for nanoelectronics [5]. Given the presence of a lone electron pair on each P-atom available to chemically interact with electron-poor atoms or moelcules, 2D bP suffers from a higher chemical reactivity in comparison to other 2D materials, which implies a lower stability in ambient conditions, being reactive with water and oxygen, and its reactivity is strongly  accelerated by light [6-8]. Thus, exfoliated black phosphorus needs to be protected by a capping layer [9-14] in order to avoid the formation of various phosphorus oxides species, which inevitably  cause  surface  degradation.  On  the  other  hand,  the  reactivity  of  few-layer  bP  can  be exploited  for  various  functionalizations,  with  molecules  [15,  16],  metal  adatoms  [17]  or  metal nanoparticles [18-20]. In particular, functionalization of bP with metal nanoparticles provided novel hybrids active in catalysis [18-20], electrocatalysis ( i.e. hydrogen evolution reaction) [21], Li storage [21] and in biomedicine [22].

## Materials and Methods

## Functionalization of bP with nickel nanoparticles

Bulk black phosphorus (bP) [23] and colloidal nickel nanoparticles [24] were prepared as previously described in literature. Liquid phase exfoliation of the bulk material was carried out under the action of  ultrasound  by  using  dimethylsulfoxide  as  solvent  in  inert  atmosphere  following  a  procedure optimized  in  our  laboratories  [25]. The  synthesis  of  the  nanohybrid  Ni/2D  bP  was  performed according  to  our  published  procedure  [26]. To  a  freshly  prepared  suspension  of  few-layer  bP  in tetrahydrofuran (1.2 mL, 1.2 mg, 0.0387 mmol P) a black colloidal suspension of nickel nanoparticles dispersed in toluene (205 µL, 0.018 M, 0.00368 mmol, P/Ni = 10.5) was added drop-wise under nitrogen at room temperature. After stirring for 25 minutes, degassed acetone (5.0 mL) was added, and the mixture was centrifuged at 8000 RPM for 15 minutes. The black residue was washed once more with acetone (5.0 mL) and dried under a stream of nitrogen. The Ni content of the isolated nanohybrid was measured by ICP-MS (Inductively Coupled Plasma-Mass Spectrometry) and resulted equal to 15.2 wt.%. The average size of Ni NPs is (11.9±0.8) nm, see Figure S1.

The NiNPs functionalized bP solution and a solution of pristine bP were drop-casted on different substrates, suitable for the different characterization techniques. The samples of pristine 2D bP and Ni/2D bP were aged storing them under ambient conditions. In detail, TEM samples were kept at T = (23 ± 3)°C, in the dark with a relative humidity h = (40 ± 10)%, see Figure S2. Raman measurements of the samples were made in another laboratory and the Raman samples were kept at T = (23.4 ± 1.5)°C and at relative humidity h = (30 ± 8)%, see Figure S3.

## Characterization

A transmission electron microscope (Philips) was used to obtain images of the materials. The incident energy of the electron beam was set to 80 kV. Dispersed fresh suspensions of 2D bP and Ni/2D bP in tetrahydrofuran (0.1 mg bP/mL in both cases) were drop-casted on carbon coated TEM grids and dried under a stream of nitrogen.

Raman spectroscopy was carried out on the same 2D bP and Ni/2D bP samples used above for TEM. For this purpose, the suspensions in tetrahydrofuran were drop-casted on a SiO2/Si wafer, kept for few minutes, then rinsed with ethanol and dried under a stream of nitrogen. A Renishaw InVia microRaman was used, equipped with an optical microscope and a motorized stage. A laser power of 0.546 mW was used, with a spot size of 2 µm diameter, to avoid any photon-induced damage to the sample. High  resolution  X-ray  Photoemission  and  Absorption  Spectroscopy  was  performed  at  the  Aloisa beamline [27] at the Elettra Synchrotron Radiation facility in Trieste (Italy), using suspensions of 2D bP  and  Ni/2D  bP  in  tetrahydrofuran  (3  mg  bP/mL  in  both  cases)  that  were  drop-casted  on Au(111)/mica, kept for 10 minutes, rinsed with ethanol, and dried under a stream of nitrogen.

## Results and discussion

## TEM study

Both on the pristine 2D bP sample and on the Ni/2D bP sample stored under ambient conditions in the dark, several flakes were chosen for TEM inspection and measured as prepared and after regular time intervals, to look for signatures of degradation. Figure 1 shows a TEM series for a pristine bP flake. Most of the pristine bP flakes were partially oxidized and hydrolyzed after one week, and within 15 days they lost completely their initial morphology and appeared as drops, see Figure 1 and Figure S4. The proposed degradation mechanism [28] goes via the leading role of water that being highly polar, interacts with oxygen adsorbed on the surface of bP and facilitates the electron transfer from bP to oxygen thus starting the oxidation process. Once P-O and P=O bonds are formed, being the Poxides  highly  hygroscopic,  they  react  with  ambient  water  yielding  molecular  compounds, i.e. phosphorus oxyacids, H3PO2, H3PO3 and H3PO4. This explains why the flake morphology is lost, and only drops are present.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00237v1-figure-3750b1a5-a5c2-4c27-b927-92da0011b13a.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00237v1-figure-c777aac0-d34e-4ef4-837f-d28911cb5b6a.jpg)

(c)

Figure 1 . TEM images of a pristine BP flake: a) freshly prepared, b) after one week, c) after two weeks in ambient conditions. All images show the same flake.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00237v1-figure-e0a0f7af-2dc2-40a6-a63c-2fbded960819.jpg)

Due to the nature of bP, Ni/2D bP flakes are  also expected to react with  oxygen and water, but surprisingly, after one week of ambient exposure, the flakes were completely preserved. Figure 2 shows the topography of the evolution of a selected flake. With prolonged time of observation, the flake appears covered by bubbles (Figures 2c-2f) due to the formation of P-oxides and P-oxyacids. Nevertheless, the morphology is preserved intact during seven months of aging, as shown in Figure 2, panels (c)-(f), suggesting that the oxidation process involves only the surface layers without going in depth.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00237v1-figure-332e3cff-763a-46f4-8c3b-47faf00278b1.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00237v1-figure-d27895f4-57af-4f56-8c96-a06b0031c6ea.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00237v1-figure-8e501b7a-2718-40aa-9ed9-ebe6c3b99066.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00237v1-figure-2848f5a2-92af-4ecc-8926-448c2f5c3fdb.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00237v1-figure-346ce129-8ecb-4eb6-b04c-4900c27742fe.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00237v1-figure-900e96ef-1bb6-4129-9699-5ccf3ca44e73.jpg)

Figure 2 . TEM images of a Ni/2D bP flake: a) freshly prepared; b) after one week; c) after 18 days; d) after one month; e) after three months, f) after seven months. The aging is performed keeping the samples under ambient conditions in the dark. All images show the same flake.

The lateral dimensions, length and width, of the observed flake (named afterwards flake 1) were obtained from TEM images and are plotted in Fig. 3. It clearly appears that the flake keeps roughly its  size  for  two  months, which indicates that the presence of the nickel nanoparticles stabilized it during  this  period.  Afterwards  there  is  a  dramatic  increase  in  size,  which  is  then  approximately constant  during  the  subsequent  three  months  of  observation.  In  other  words,  the  degraded  area increases fast between the second and third month of aging, and afterwards it reaches a saturation, following an S-shaped growth curve. A similar behavior has been observed for bP flakes completely coated with 5 nm thick Al2O3 and exposed to ambient conditions [29]. It is remarkable that in our case the surface coverage with Ni NPs ranges from 8% to 15% (as calculated from TEM images) far away from being a thick layer covering 100% of the surface of the bP flake.

Figure 3 .  Variation of the lateral dimensions of the Ni/2D bP flake, length (blue) and width (red) during six months of ambient exposure.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00237v1-figure-9c5d3c32-cf88-4387-84a7-c4f6899166a3.jpg)

Other Ni/2D bP flakes were studied by TEM, as well, within the six months of ambient exposure. All of them preserved their morphology. Their lateral dimensions were found to be enlarged up to 60 % in comparison to the fresh sample, following the same trend discussed above for flake 1, see Table 1 and Figures S5-S7.

Table 1 . Variation of the lateral dimensions (in nm) of Ni/2D bP flakes kept in ambient conditions.

| Entry   | dimension   |   day 0 |   2 months |    % 2 months |   3 months |    % 3 months |   6 months |    % 6 months |
|---------|-------------|---------|------------|----------------|------------|----------------|------------|----------------|
| Flake 1 | lenght      |     515 |        536 |            4.1 |        667 |           29.5 |        705 |           36.9 |
|         | width       |     503 |        520 |            3.4 |        628 |           24.9 |        659 |           31   |
| Flake 2 | lenght      |    1107 |       1040 |           -6.1 |       1290 |           16.5 |       1307 |           18.1 |
|         | width       |     768 |        758 |           -1.3 |       1000 |           30.2 |        972 |           26.6 |
| Flake 3 | lenght      |     704 |        742 |            5.4 |        910 |           29.3 |        970 |           37.8 |
|         | width       |     165 |        166 |            0.6 |        260 |           57.6 |        265 |           60.6 |
| Flake 4 | lenght      |     737 |        760 |            3.1 |       1000 |           35.7 |       1012 |           37.3 |
|         | width       |     381 |        394 |            3.4 |        539 |           41.5 |        550 |           44.4 |

A comparison between the crystallographic structure of 2D bP [30] and P2O5 (which is also a layered material featuring an in-plane anisotropic structure) [31] shows a remarkable expansion of the unit cell going from pristine bP to the oxidized species: the size of the unit cell in the zig-zag direction increases from 3.316(1) Å to 4.890(4) Å (  = 47.5%), and in the armchair direction from 4.379(1) Å to  7.162(7)  Å  (  =  63.5%),  as  shown  in  Figure  S8.  Thus,  the  observed  increase  of  the  lateral dimensions of the flakes is an indication that O2 reacted with bP forming P-O bonds. As will be discussed  later,  the  presence  of  phosphorus  oxides  is  further  confirmed  by  x-ray  photoelectron spectroscopy (XPS) measurements.

## Raman study

The Raman signature of bP is very clear and consists of three main peaks [6] at 361 cm -1 , 438 cm -1 , and 466 cm -1 . They are labeled as A 1 g, B2g, and A 2 g, and they are related to the out-of-plane, in-plane along  zigzag,  and  in-plane  along  armchair  vibrational  modes,  respectively  [32].  Since  these vibrational modes are strongly suppressed by oxidation [6], Raman spectroscopy is another powerful tool to study bP aging. These measurements were performed on dropcasted samples. Because of dropcasting, the small 2D bP flakes (both pristine and functionalized) tend to aggregate during solvent evaporation: a first morphological information that we can extract from the freshly prepared samples is that aggregation is much stronger in the case of pristine bP, with respect to the functionalized bP. In Figure 4 three typical areas for each sample are shown.

Figure 4 . Typical morphology of aggregates of a freshly prepared 2D bP sample [(a) - (c)] and of a Ni NPs functionalized 2D bP sample [(d) - (f)]. Aggregation is more pronounced for the pristine sample. Scale bar = 10 µm.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00237v1-figure-6093c83e-8019-42e6-abbe-b582cc458d14.jpg)

Since aggregates of various morphology and size are present on the substrate, ranging from less than one µm to several 10 µm in lateral size, a statistical analysis of the Raman data has been performed, measuring several flakes and aggregates of different morphology - up to 45 at a time - to capture as much as possible the aging trend of the overall  sample.  Over  all  the  experiment,  more  than  200 aggregates have been measured for each of the two samples.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00237v1-figure-f67e7646-a7be-4404-9ba4-35c5ecf6cdaf.jpg)

Figure  5 . Typical  morphology  of  the  pristine  sample  -  panels  (a)-(c)  -  and  of  the  Ni  NPs functionalized sample - panels (d)-(f) - after 15 months of aging. A more pronounced hydrolysis is

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00237v1-figure-0dddd2e9-e27e-486a-98fc-1691c8681ec4.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00237v1-figure-8a092661-bdbc-4bfb-b44b-8acacb10076a.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00237v1-figure-ef2502b3-c05a-4713-aace-09f3f9122480.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00237v1-figure-0aa06df1-fedd-4891-9339-7ef37c313932.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00237v1-figure-5bfc8bc7-250f-4fbc-b33a-73eb5a9f8d75.jpg)

present on the pristine 2D bP sample, while even small aggregates display a quite clear structure for the Ni NPs functionalized sample. Scale bar = 10 µm.

The morphology at the end of month 15 is shown in Figure 5. While the pristine sample [images (a)(c)] shows severe marks of degradation, the functionalized sample [images (d)-(f)] shows still several small aggregates which are just slightly degraded.

The averaged Raman spectrum from all aggregates measured at the beginning of the experiment, normalized to the Si peak intensity of the substrate, is presented in Figure 6(a) for the pristine bP sample and in Fig. (b) for the functionalized one. The Si peak at 520 cm -1 is more pronounced for the functionalized sample because the aggregates have a much smaller lateral size, often smaller than the diameter of the laser spot. Averaged normalized spectra  after 15 months of exposure to  ambient conditions are presented in Figure 6(c) for the not-functionalized sample and in Figure 6(d) for the functionalized one. The relative decrease of the bP peaks intensity is much stronger for the nonfunctionalized sample, for which many of the measured aggregates were not Raman active anymore, while the Ni NPs functionalized sample displays still all Raman peaks with higher intensity.

Figure 6 . Averaged Raman spectra for (a) the pristine bP sample and (b) the Ni/2D bP sample at the beginning of the experiment. (c) and (d), the same after 15 months of exposure to ambient conditions.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00237v1-figure-528a25b7-4e8e-4c83-94b8-7cdce4dc0422.jpg)

All spectra were acquired in the same conditions. Spectra for the same sample are averaged, then the background is subtracted, and finally the spectra are normalized to the Si substrate peak intensity.

To visualize the trend of the Raman activity with time, we present the aggregate Raman data as a function of time in Figure 7. Since the intensity of the bP peaks is not just related to aging, but as already discussed also to the initial thickness and lateral size of the aggregates, as a first indicator for the good preservation of the sample, we use the Raman activity in itself - namely the presence of the bP peaks in the Raman spectra, independent of their intensity. For each date, the number of Raman active aggregates with respect to the total number of measured aggregates is reported, attributing as an error bar the standard error of the distribution, evaluated as 1/√𝑁 , where N is the total number of measurements. Degradation is more pronounced for the pristine bP sample, but given the large variety in the aggregates' morphology, the two samples are still compatible within the error bars even after 12 months. However, in the last part of the experiment, and in particular for the last measurements performed after 20 months, the two distributions are well separated. This clearly shows that the Ni/2D bP sample degrades much more slowly than the pristine one. This result is even more striking when this information is complemented with the morphological information. From the beginning of the experiment, the Ni/2D bP aggregates are much smaller, as visible from the optical microscopy images in Fig. 4. In general, larger and thicker bP flakes degrade more slowly, because of their lower surfaceto-volume ratio and consequently their lower reactivity [6]. Here, however, we observe that the larger pristine  bP  aggregates  degrade  much  faster  than  the  smaller  Ni/2D  bP  aggregates.  This  is  an unequivocal proof of the enhanced environmental stability of the Ni/2D bP sample. At the end of the experiment, just (22 ± 15)% of the aggregates belonging to the non-functionalized sample were still Raman active, while of the Ni NPs-functionalized aggregates were (73 ± 15)%.

Figure 7 . Number of Raman active aggregates for the not functionalized sample (green bars) and for the  functionalized  sample  (red  bars).  Measurements  before  day  45  are  not  shown  to  enhance  the clarity of the graph, since they are all compatible within the error bar with 100% of Raman activity.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00237v1-figure-9cd077f9-9d91-4407-bba8-20336a4cdefe.jpg)

To quantify the decay in Raman activity, a plot of the A 2 g peak intensity, normalized to the intensity of the Si peak, is presented in Figure 8 as a function of time for (a) the pristine sample and (b) the functionalized one. This analysis has been done with the average spectra (with 5 days binning), to overcome as much as possible the scattering due to the morphological differences among flakes.

Figure 8 . Averaged intensity of the A 2 g peak normalized to the intensity of the Si peak, plotted as a function of time. The individual spectra are grouped with a 5 days binning and averaged. Then the background is subtracted and the intensity of the A 2 g peak is normalized to the Si peak. The trend is fitted with an exponential decay.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00237v1-figure-638e7ff4-5d61-4822-bec2-318dcd24f1bd.jpg)

The smaller intensity of the normalized A 2 g peak for the functionalized sample with respect to the pristine at the beginning of the experiment is consistent with the data presented in Figure 6.

The normalized A 2 g intensity data as a function of time have been fitted with an exponential decay ∝ 𝑒 -𝑡 𝜏 / ,  where 𝑡 is  time  and 𝜏 is  decay  time,  as  previously  reported  by  Favron et  al. [6].  The exponential decay captures well the trend of the data, as shown in Figure 8. We obtain a decay time 𝜏 ≈  90 𝑑𝑎𝑦𝑠 for the non-functionalized bP sample and 𝜏 ≈  300 𝑑𝑎𝑦𝑠 for the one functionalized with NPs, confirming the faster degradation of the non-functionalized sample. The decay for the Ni NPs-functionalized  sample  is  therefore  more  than  three  times  slower  than  for  pristine  bP,  in agreement with what observed in Figure 6. This clearly proofs that Ni NPs slow down the degradation of exfoliated bP.

## X-ray Photoemission and Absorption Spectroscopy study

To study quantitatively the role played by the presence of Ni NPs on the environmental stability of bP and to evaluate the amount and the identity of oxidized species formed on the surface, XPS were measured on a fresh sample of Ni/2D bP and after aging in ambient conditions for 18 days. In parallel, a sample of pristine bP kept in the same conditions was examined as a reference. The binding energy scale has been calibrated by measuring the Au 4f7/2 peak taken as reference [33]. As shown in Figure 9a, the XPS spectrum of a sample of fresh 2D bP exhibits two maxima at 130.0 and 130.85 eV, which are characteristic of the spin-orbit split phosphorus P 2p peak positions, namely P 2p3/2 and P 2p1/2, respectively [34]. The deconvolution of the XPS data was carried out using a fitting function with one spin-orbit doublet (P1) plus up to 3 Gaussian functions (P2, P3, P4) and the corresponding Shirley background [34, 35]. We name P1 the peak of not oxidized phosphorus, P4 the peak corresponding to phosphorous in its highest oxidation state (+5), and in some cases we added also two more peaks, P2 and P3, representing P-oxides with intermediate oxidation states +1 and +3. Figure 9b shows the XPS spectrum of a fresh sample of Ni/2D bP, which has the same features as pristine bP, plus a broad peak  at  134.2  eV  attributed  to  an  incipient  formation  of  phosphorus  oxide  species.  Indeed,  the presence of electronegative groups such as oxygen or  hydroxide on the surface of bP withdraws charge leading to more tightly bound phosphorus electrons. As a result, PxOy species experience a higher binding energy for P 2p core electrons in comparison to un-oxidized material [36].

Best Fit

138

Figure 9 . P 2p core level XPS spectra. a) XPS spectrum of a fresh sample of 2D bP dropcasted on Au/mica, b) XPS spectrum of a fresh Ni/2D bP sample on Au/mica. Both spectra were taken at h  = 415 eV.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00237v1-figure-7b1d0ca2-14b4-4f3c-8335-ce254f8fb09f.jpg)

Figure 10 shows a comparison between XPS spectra of non-functionalized bP (left side) and Ni/2D bP (right side) both aged in a dark environment for the same time (18 days). Intriguingly, while the peaks  related  to  P-P  bonds  remain  at  the  same  binding  energy  (B.E.)  (130.0  and  130.85  eV, respectively, for the P 2p doublet) in both samples, in the case of P-O bonds, the peak maximum is observed at different binding energies: 134.1 eV in Ni/2D bP and 134.9 eV in pristine bP, indicating different compositions of phosphorus oxide species. On the basis of calculations by D. Tománek et al . [37], 2D bP oxides with oxidation state of +1, +3, and +5 have core-level binding energies in the range of 132-136 eV. Thus, a shift of 0.8 eV to higher B.E. suggests a larger amount of P-oxides having  oxidation  state  +5, i.e. phosphorus  pentoxide,  P2O5.  Consequently,  a  higher  degree  of oxidation in the case of pristine bP with respect to Ni/2D bP is confirmed. To evaluate the thickness of the oxide layer, XPS measurements were carried out on both samples using three different energies for  the  incident  photons:  415  eV,  515  eV,  and  1000  eV,  in  order  to  increase  progressively  the photoelectron mean free path and thus the sampling depth of the experiment. Increasing the photon energy from 415 eV to 1000 eV, there is a gradual decrease of the intensity of the peak corresponding to P-O bonds. This is a clear signature that only a thin oxide layer is present atop of the bP flakes.

Figure 10 . Depth-dependent photoelectron spectroscopy of the P 2p core level of 2D bP (left side, panels a-c) and Ni/2D bP (right side, panels d-f) after 18 days of air exposure, taken at h  = 415 eV (a, d), 515 eV (b, e) and 1000 eV (c, f).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00237v1-figure-f25aeafd-bd47-4696-8d87-c11f53053c39.jpg)

The  thickness  of  the  phosphorus  oxide  layer  was  estimated  according  to  the  equation:  d  =  λ ln[(Ioxide/IbP)+1] where λ is the photoelectron inelastic mean free path [38], and Ioxide and IbP are the area in the XPS spectra under the component of P-oxides and elemental phosphorus, respectively.

Thus, we determine an average thickness of phosphorus oxide of (13.9 ± 0.4) Å for pristine bP and of (6.2 ± 0.4) Å for Ni/2D bP after aging. Being the thickness of a single layer of planar P2O5 equal to 5.5 Å according to literature [39, 31] this means that we have three layers of P2O5 on top of pristine bP and approximately one layer on top of Ni/2D bP, confirming that the corrosion caused by oxygen and water is much deeper in the case of pristine bP. A further confirmation of the protective role of Ni NPs comes from the comparison of the ratio between P1 and P4 areas as a function of the inelastic mean free path shown in Fig S9, where the signal corresponding to not oxidized P is dramatically damped down when Ni NPs are missing.

In the region where the XPS signature of Ni is expected, no significant signal could be acquired, because of the relatively low amount of Ni with respect to bP (see Experimental Section). Therefore, NEXAFS (Near Edge X-Ray Absorption Fine Structure) measurements were carried out, because this technique has a higher sensitivity. Figure S10 shows the NEXAFS spectrum collected on a sample of Ni/2D bP freshly prepared. The characteristic L3 and L2 peaks of metallic nickel at 852.7 eV and 870.0 eV, respectively, are detected while the peak at 854.0 eV due to Ni(II) is absent, which means that nickel nanoparticles are pure Ni(0), and there is no oxidation. Running NEXAFS on the aged sample of Ni/2D bP, we detected again the characteristic L3 signal of metallic nickel at 852.7 eV plus a shoulder at ~ 854.5 eV, suggesting the presence of Ni(II) due to the formation of a NiO layer as outer shell of the nanoparticles, being a consequence of the aging in ambient conditions, see Figure S11.

The different degradation rate in Ni/2D bP with respect to 2D bP can be undoubtedly related to the presence of Ni NPs, although they only cover roughly from 8% to 15% of the nanosheets' surface. Thus, we cannot ascribe to Ni NPs the role of a physical barrier against oxidation, we can invoke a role, mainly of chemical nature, played by Ni NPs which are expected to interact with defect sites of 2D bP surface, such as steps and edges, being the more reactive. In this way, they prevent those sites to be approached by oxygen and water molecules. The overall degradation process results markedly slowed down, meaning that also the reactivity of bare P atoms, i.e. not directly interacting with Ni NPs,  is  slowed  down.  This  can  be  interpreted  on  the  basis  of  our  recent  results  [20]  on  2D  bP functionalized with Pd NPs, where a coordination bond of covalent nature was found between P atoms and outer  shell  of  Pd  NPs.  On  this  ground,  a  covalent  bond  between  Ni  NPs  and  2D  bP  is  also expected, and since it is known [41] that covalent modifications of 2D bP may induce variations of electronic properties, and shifts of CBM (conduction band minimum) below the redox potential of O2-O2 , than it is easily explained the passivation observed in our functionalized bP nanosheets.

## Conclusions

The simultaneous influence of two factors, air and ambient humidity (up to 50%), has been evaluated in the stability of 2D bP functionalized with Ni NPs. The three techniques used, TEM, Raman, and high resolution XPS, are in very good agreement in showing an increased chemical stability of Ni/2D bP flakes, that maintain intact their size for one month. The presence of nickel nanoparticles on the surface of 2D bP makes the rate of the oxidation process more than three times slower as derived from  the  decay  of  A 2 g  Raman  signal,  and  a  thinner  layer  of  phosphorus  pentoxide  is  formed  as measured by XPS. Remarkably, in our experiments the surface coverage of bP flakes by Ni NPs is only partial, thus we cannot invoke a physical barrier of protection made by the metal nanoparticles but their chemical interaction with bP should be addressed. Thus, a large part of the bP surface is uncovered and remains available for further functionalization and applications in devices that are necessarily  kept  in  ambient  conditions  and  in  which  access  to  the  surface  is  crucial,  such  as  gas sensing.

## Acknowledgement

Thanks are expressed to EC for funding the project PHOSFUN ' Phosphorene  functionalization:  a  new  platform  for  advanced  multifunctional  materials '  (ERC ADVANCED GRANT N. 670173 to M.P.). F.T. thanks CNR-NANO for funding the SEED project SURPHOS ( Surface properties of few layer black phosphorus investigated by scanning tunneling microscopy. ) S.H. acknowledges support from Scuola Normale Superiore, project SNS16\_B\_HEUN - 004155. Mr. Carlo Bartoli (CNR ICCOM) is acknowledged for his skillfull technical assistance. A. V., A.G. and M.D.M. are partly supported by the NATO project G5140.

## Supporting information

Histogram relative to the size distribution of Ni NPs is given in S1. Graphs related to T (°C) and humidity to which samples were exposed along six months are shown in S2-S3. TEM images of a reference pristine bP flake is given in S4, graphs of the variation of the lateral dimensions versus time of selected Ni/2D bP flakes are given in S5-S7, crystallographic structures of bP and P2O5 are given in S8, the P1/P4 ratio versus photon energy in aged samples is given in S9, NEXAFS of Ni/2D bP are given in S10-S11.

## ORCID IDs

Caporali M http://orcid.org/0000.0001-6994-7313

Peruzzini M http://orcid.org/0000.0002-2708-3964

## Notes

The authors declare no competing financial interests.

## References

- [1]  Pumera  M,  Sofer  Z  and  Gusmão  R  2017  Black  Phosphorus  Rediscovered:  from  Bulk  Material  to Monolayers Angew. Chem. Int. Ed . 56 8052-72
- [2]  Peruzzini  M et  al. 2019  A  perspective  on  recent  advances  on  phosphorene  functionalization  and  its applications in devices Eur. J. Inorg. Chem . 1476-94
- [3] Das S, Zhang W, Demarteau M, Ho Mann A, Dubey M and Roelofs A 2014 Tunable transport gap in phosphorene Nano Lett . 14 5733-9
- [4] Castellanos-Gomez A 2015 Black Phosphorus: Narrow Gap, Wide Applications J. Phys. Chem. Lett . 6 4280-91
- [5] Long G, Maryenco D, Shen J, Xu S, Hou J, Wu Z, Wong W K, Han T, Lin J, Cai Y and Wang N 2016 Achieving Ultrahigh Carrier Mobility in Two-Dimensional Hole Gas Black Phosphorus Nano Lett. 16 776873
- [6] Favron A, Gaufrès E, Fossard F, Phaneuf-L'Heureux A-L, Tang N Y-W, Levesque P L, Loiseau A, Leonelli R, Francoeur S and Martel R 2015 Photooxidation and Quantum Confinement Effects in Exfoliated Black Phosphorus Nat. Mater . 14 826-32
- [7] Walia S et al . 2017 Ambient protection of few-layer black phosphorus via sequestration of reactive oxygen species Adv. Mater . 1700152
- [8] Zhang T, Wan Y, Xie H, Mu Y, Du P, Wang D, Wu X, H. Ji and Wan L 2018 Degradation Chemistry and Stabilization of Exfoliated Few-Layer Black Phosphorus in Water J. Am. Chem. Soc . 140 7561-67
- [9] Avsar  A  et  al.  2015  Air-stable  transport  in  graphene-contacted,  fully  encapsulated  ultrathin  black phosphorus-based field-effect transistors ACS Nano 94138-45
- [10] Doganov R A et al. 2015 Transport properties of pristine few-layer black phosphorus by van der Waals passivation in an inert atmosphere Nat. Commun . 66647
- [11]  Liu H, Neal A T, Si M, Du Y and Ye P D 2014 The effect of dielectric capping on few-layer phosphorene transistors: tuning the schottky barrier heights IEEE Electron Device Lett. 35 795-7
- [12] Abate Y, Akinwande D, Gamage S, Wang H, Snure M, Poudel N and Cronin S B 2018 Recent progress on stability and passivation of black phosphorus Adv. Mater . 1704749
- [13] Telesio F, Passaglia E, Cicogna F, Costantino F, Serrano-Ruiz M, Peruzzini M and Heun S, 2018 Hybrid 2D black phosphorus/polymer materials: new platforms for device fabrications Nanotechnology 29 295601 [14]  Passaglia  E et  al. 2018  Polymer-Based  Black  Phosphorus  (bP)  Hybrid  Materials  by  in  Situ  Radical Polymerization: An Effective Tool To Exfoliate bP and Stabilize bP Nanoflakes Chem. Mater . 30 2036-48

[15]  Pumera  M,  Sofer  Z  and  Gusmão  R  2018  Functional  Protection  of  Exfoliated  Black  Phosphorus  by Noncovalent Modification with Anthraquinone ACS Nano 12 5666-73

[16] Ryder C R, Wood J D, Wells S A, Yang Y, Jariwala D, Marks T J, Schatz G C, and Hersam M C 2016 Covalent Functionalization and Passivation of Exfoliated Black Phosphorus via Aryl Diazonium Chemistry Nature Chem . 8 597-602

[17] Koenig S P, Doganov R A, Seixas L, Carvalho A, Tan J Y, Watanabe K, Taniguchi T, Yakovlev N,  Castro Neto A H and Özyilmaz B 2016 Electron Doping of Ultrathin Black Phosphorus with Cu Adatoms Nano Lett . 16 2145-51

[18] Lei W, Zhang T, Liu P, Rodriguez J A, Liu G and Liu M 2016 Bandgap and Local Field-Dependent Photoactivity of Ag/Black Phosphorus Nanohybrid ACS Catal . 6 8009-20

[19] Cho S-Y, Koh H-J, Yoo H-W and Jung H-T, Tunable Chemical Sensing Performance of Black Phosphorus by Controlled Functionalization with Noble Metals. 2017 Chem. Mater . 29 7197-7205

[20] Vanni M, Serrano-Ruiz M, Telesio F, Heun S, Caporali M, Peruzzini M et al , Black Phosphorus/Palladium Nanohybrid: Unraveling the Nature of P-Pd Interaction and Application in Selective Hydrogenation 2019 Chem. Mater . accepted.

[21] Luo Z-Z, Zhang Y, Zhang C, Teng Tan H, Li Z, Abutaha A, Wu X-L, Xiong Q, Kor K A, Hippalgaonkar K, Xu J, Hng H H and Yan Q 2017 Multifuctional 0D-2D Ni2P nanocrystals-black phosphorus heterostructures Adv. Energy Mater. 7 1601285

[22] Wang H and Yu X-F 2018 Few-Layered Black Phosphorus: From Fabrication and Customization to Biomedical Applications Small 14 1702830

[23] Köpf M, Eckstein N, Pfister D, Grotz C, Krüger I, Greiwe M, Hansen T, Kohlmann H and Nilges T 2014 Access and in situ growth of phosphorene-precursor black phosphorus J. Crystal Growth 405 6-10

[24] Carenco S, Boissière C, Nicole L, Sanchez C, Le Floch P and Mézailles N 2010 Controlled Design of Size-Tunable Monodisperse Nickel Nanoparticles Chem. Mater . 22 1340-9

[25] Serrano-Ruiz M, Caporali M, Ienco A, Piazza V, Heun S and Peruzzini M 2016 The Role of Water in the Preparation and Stabilization of High-Quality Phosphorene Flakes Adv. Mater. Interfaces 3 1500441 [26] Caporali M, Serrano-Ruiz M, Telesio F, Heun S, Nicotra G, Spinella C and Peruzzini M 2017 Decoration of  exfoliated  black  phosphorus  with  nickel  nanoparticles  and  application  in  catalysis Chem. Commun . 53 10946-49

[27] Floreano L, Cossaro A, Gotter R, Verdini A, Bavdek G, Evangelista F, Ruocco A, Morgante A. and Cvetko D 2008 Periodic Arrays of Cu-Phthalocyanine Chains on Au(110) J. Phys. Chem. C 112 , 10794-802 [28] Hu Z, Li Q, Lei B, Zhou Q, Xiang D, Lyu Z, Hu F, Wang J, Ren Y, Guo R, Goki E, Wang L, Han C, Wang J and Chen W 2017 Water-Catalyzed Oxidation of Few-Layer Black Phosphorus in a Dark Environment Angew. Chem. Int. Ed . 56 9131-35

[29] Gamage S, Li Z, Yakovlev V S, Lewis C, Wang H, Cronin S B and Abate Y 2016 Nanoscopy of Black Phosphorus Degradation Adv. Mater. Interfaces 3 1600121

[30] Lange S, Schmidt P, Nilges T 2007 Au3SnP7@Black Phosphorus: An Easy Access to Black Phosphorus Inorg. Chem . 46 4028-35

[31] Stachel D, Svoboda I and Fuess H Phosphorus Pentoxide at 233 K Acta Cryst. C51 (6) 1049-50 [32] Ribeiro H B, Pimentab M A and de Matosa C J S 2018 Raman spectroscopy in black phosphorus J. Raman Spectrosc. 49 76-90

[33] Cossaro A, Floreano L, Verdini A, Casalis L and Morgante A 2009 Comment on ''Local Methylthiolate Adsorption Geometry on Au(111) from Photoemission Core-Level Shifts' Phys. Rev. Lett . 103 119601 [ 34] Edmonds M T, Tadich A, Carvalho A, Ziletti A, O'Donnell K M, Koenig S P, Coker D F, Özyilmaz B, Castro Neto A H and Fuhrer M S 2015 Creating a Stable Oxide at the Surface of Black Phosphorus ACS Appl. Mater. Interfaces 7 14557-62

[35] Taniguchi M, Suga S, Seki M, Sakamoto H, Kanzaki H, Akahama Y, Terada S, Endo S and Narita S 1983 Valence  Band  and  Core-Level  Photoemission  Spectra  of  Black  Phosphorus  Single  Crystals Solid  State Commun . 45 59-61

[36] Kuntz K L, Wells R A, Hu J, Yang T, Dong B, Guo H, Woomer A H, Druffel D L, Alabanza A, Tománek D and Warren S C 2017 Control of Surface and Edge Oxidation on Phosphorene ACS Appl. Mater. Interfaces 9 9126-35

[37] Yang T, Dong B, Wang J, Zhang Z, Guan J, Kuntz K, Warren S C and Tomanek D 2015 Interpreting core-level spectra of oxidizing phosphorene: theory and experiment Phys. Rev. B: Condens. Matter Mater. Phys . 92 125412

[38] The inelastic mean free path λ was calculated according to the following reference: Cumpson P J and Seah M P 1997 Surface and Interface Analysis 25 430-446 and resulted to be 9.14 Å, 11.2 Å, and 20.4 Å for photon energy 415 eV, 515 eV, and 1000 eV, respectively

[39]  Ziletti  A,  Carvalho  A,  Trevisanutto  P  E,  Campbell  D  K,  Coker  D  F  and  Castro  Neto  A  H  2015 Phosphorene oxides: band gap engineering of phosphorene by oxidation Phys. Rev. B 91 085407-416 [40] This is also observed in TEM measurements, see Figures S5-S7 in Supporting Info

[41]  Li  Q,  Zhou  Q,  Shi  L,  Chen  Q  and  Wang  J  2019  Recent  Advances  in  Oxidation  and  Degradation mechanism of ultrathin 2D materials under ambient conditions and their passivation strategies J. Mater. Chem. A 7 4291-4312

## Supplementary Information to 'Enhanced ambient stability of exfoliated black phosphorus by passivation with nickel nanoparticles'

Maria Caporali, a* Manuel Serrano-Ruiz,  Francesca Telesio,  Stefan Heun,  Alberto Verdini,  Albano a b b c Cossaro,  Matteo Dal Miglio,  Andrea Goldoni,  Maurizio Peruzzini c d d a*

- a CNR ICCOM, Via Madonna del Piano 10, 50019 Sesto Fiorentino, Italy.
- b NEST, Istituto Nanoscienze-CNR and Scuola Normale Superiore, Piazza S. Silvestro 12, 56127 Pisa, Italy.
- c TASC CNR IOM, Area Science Park, Basovizza, 34149 Trieste, Italy.
- d Elettra-Sincrotrone Trieste, Area Science Park, S.S. 14 km 163.5, I-34149 Trieste, Italy.

E-mail : maria.caporali@iccom.cnr.it, mperuzzini@iccom.cnr.it

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00237v1-figure-3a037864-d3a0-42a7-be3e-15f6a6c71b57.jpg)

Figure S1. Bright field TEM image of Ni/2D bP. Scale bar: 200 nm. The inset shows the relative size distribution histogram of the Ni NPs. The average size of the Ni NPs is (11.9 ±0.8) nm.

Figure S2. Ambient conditions, temperature (T) and humidity (h), to which TEM samples were exposed for six months. Average T = (23±3) °C, average h = (40±10)%.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00237v1-figure-8cd4083b-2ddc-4fe7-906c-a1640462fd9f.jpg)

Figure S3. Ambient conditions, temperature (T) and humidity (h), to which Raman samples were exposed for twenty months. Average T = (23±1.5) °C, average h = (30±8)%.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00237v1-figure-69830752-8635-4b09-b1ef-68321714b784.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00237v1-figure-9655b899-6415-43f3-8ed4-4ec6644dd901.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00237v1-figure-90db146a-37ce-49f9-bebf-46a9a2ed1acb.jpg)

(b)

(c)

Figure S4. TEM images of a pristine bP flake: a) freshly prepared, b) after two weeks, c) after two months under ambient conditions.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00237v1-figure-86736aff-2ae2-4566-a17c-2bd124ac01bc.jpg)

Figure S5. Variation of the lateral dimensions of flake 2, during six months of ambient exposure.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00237v1-figure-a6b63301-569c-4fc7-80f3-52b56933f202.jpg)

|   day |   lenght (nm) |   width (nm) |
|-------|---------------|--------------|
|     0 |          1107 |          768 |
|     7 |          1086 |          763 |
|    14 |          1070 |          760 |
|    24 |          1042 |          750 |
|    36 |          1036 |          764 |
|    58 |          1040 |          758 |
|    92 |          1290 |         1000 |
|   170 |          1276 |          983 |
|   192 |          1307 |          972 |

Figure S6. Variation of the lateral dimensions of flake 3, during six months of ambient exposure.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00237v1-figure-aeb51335-287f-4adc-b8fa-7ec11b5c7799.jpg)

|   day |   lenght (nm) |   width (nm) |
|-------|---------------|--------------|
|     0 |           704 |          165 |
|     7 |           718 |          173 |
|    14 |           689 |          161 |
|    24 |           686 |          160 |
|    36 |           688 |          168 |
|    58 |           742 |          166 |
|    92 |           910 |          260 |
|   170 |           970 |          265 |
|   192 |           971 |          266 |

Figure S7. Variation of the lateral dimensions of flake 4 during six months of ambient exposure.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00237v1-figure-68096c6c-7dba-4af4-9e10-a5ba3fb0a5b3.jpg)

| day   |   lenght (nm) |   width (nm) |
|-------|---------------|--------------|
|       |           737 |          381 |
|       |           757 |          392 |
| 14    |           756 |          377 |
| 24    |           753 |          373 |
| 36    |           754 |          371 |
| 58    |           760 |          394 |
| 92    |          1000 |          539 |
| 170   |          1011 |          540 |
| 192   |          1012 |          550 |

Figure S8. Comparison of the crystallographic structure of bP (left) and P2O5 (right).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00237v1-figure-e016e797-dc1b-47e0-be88-61427b1f7ea4.jpg)

Figure S9. The ratio between the areas underneath the peaks labeled P1 (related to P-P bonds) and P4 (related  to  P-oxides)  is  reported  as  a  function  of  the  photon  energy  (bottom  axis)  and  for  the calculated photoelectron inelastic mean free path (top axis).[1]

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00237v1-figure-7971432a-86d3-4d89-aeef-4bdbd51654f3.jpg)

Figure S10. NEXAFS of L2,3 edge of nickel in a fresh sample of Ni/2D bP.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00237v1-figure-be7be5f3-46eb-435f-bf60-b05b88bc5fa6.jpg)

Figure S11. NEXAFS of the L3 edge of nickel. a) Red line: fresh sample of Ni/2D bP; black line: Ni/bP aged 18 days.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00237v1-figure-c063d798-2979-4c27-8884-7bbaa15bd67e.jpg)

## References

[1] The inelastic mean free path λ was calculated according to the following reference: P. J. Cumpson, M. P. Seah, Surface and Interface Analysis, 25, 430-446 (1997), and resulted to be 9.14 Å, 11.2 Å, and 20.4 Å for photon energy 415 eV, 515 eV, and 1000 eV, respectively.
Maria Caporali Manuel Serrano-Ruiz Francesca Telesio Stefan Heun Alberto Verdini Albano Cossaro Matteo Dal Miglio Andrea Goldoni Maurizio Peruzzini
Enhanced ambient stability of exfoliated black phosphorus by passivation with nickel nanoparticles
Since its discovery, the environmental instability of exfoliated black
phosphorus (2D bP) has emerged as a challenge that hampers its wide application
in chemistry, physics, and materials science. Many studies have been carried
out to overcome this drawback. Here we show a relevant enhancement of ambient
stability in few-layer bP decorated with nickel nanoparticles as compared to
pristine bP. In detail, the behavior of the Ni-functionalized material exposed
to ambient conditions in the dark is accurately studied by TEM (Transmission
Electron Microscopy), Raman Spectroscopy, and high resolution X-ray
Photoemission and Absorption Spectroscopy. These techniques provide a
morphological and quantitative insight of the oxidation process taking place at
the surface of the bP flakes. In the presence of Ni NPs, the decay time of 2D
bP to phosphorus oxides is more than three time slower compared to pristine bP,
demonstrating an improved structural stability within twenty months of
observation.