Noncovalent Functionalization of 2D Black Phosphorus with

May 29, 2019 - We studied the chemical–physical nature of interactions involved in the ...... AFM microscopy of the oxidized adducts and oxidized 2D...
1 downloads 0 Views 2MB Size
This is an open access article published under an ACS AuthorChoice License, which permits copying and redistribution of the article or any adaptations for non-commercial purposes.

Research Article Cite This: ACS Appl. Mater. Interfaces 2019, 11, 22637−22647

www.acsami.org

Noncovalent Functionalization of 2D Black Phosphorus with Fluorescent Boronic Derivatives of Pyrene for Probing and Modulating the Interaction with Molecular Oxygen

Downloaded via 178.57.65.87 on August 3, 2019 at 08:18:01 (UTC). See https://pubs.acs.org/sharingguidelines for options on how to legitimately share published articles.

Margherita Bolognesi,† Salvatore Moschetto,† Mariachiara Trapani,‡ Federico Prescimone,† Claudia Ferroni,§ Gabriele Manca,∥ Andrea Ienco,∥ Silvia Borsacchi,⊥ Maria Caporali,∥ Michele Muccini,† Maurizio Peruzzini,∥ Manuel Serrano-Ruiz,∥ Lucia Calucci,⊥ Maria Angela Castriciano,*,‡ and Stefano Toffanin*,† †

Istituto per lo Studio dei Materiali Nanostrutturati (ISMN)Consiglio Nazionale delle Ricerche (CNR), Via P. Gobetti 101, 40129 Bologna, Italy ‡ Istituto per lo Studio dei Materiali Nanostrutturati (ISMN)Consiglio Nazionale delle Ricerche (CNR), c/o Dipartimento di Scienze Chimiche, Biologiche, Farmaceutiche e Ambientali, University of Messina, V.le F. Stagno d’Alcontres 31, 98166 Messina, Italy § Istituto per la Sintesi Organica e la Fotoreattività (ISOF)Consiglio Nazionale delle Ricerche (CNR), Via P. Gobetti 101, 40129 Bologna, Italy ∥ Istituto di Chimica dei Composti Organometallici (ICCOM)Consiglio Nazionale delle Ricerche (CNR), Via Madonna del Piano 10, 50019 Sesto Fiorentino, Florence, Italy ⊥ Istituto di Chimica dei Composti Organometallici (ICCOM)Consiglio Nazionale delle Ricerche (CNR), via G. Moruzzi 1, 56124 Pisa, Italy S Supporting Information *

ABSTRACT: We studied the chemical−physical nature of interactions involved in the formation of adducts of twodimensional black phosphorus (2D BP) with organoboron derivatives of a conjugated fluorescent molecule (pyrene). Timeresolved fluorescence spectroscopy showed a stabilization effect of 2D BP on all derivatives, in particular for the adducts endowed with the boronic functionalities. Also, a stronger modulation of the fluorescence decay with oxygen was registered for one of the adducts compared to the corresponding organoboron derivative alone. Nuclear magnetic resonance experiments in suspension and density functional theory simulations confirmed that only noncovalent interactions were involved in the formation of the adducts. The energetic gain in their formation arises from the interaction of P atoms with both C atoms of the pyrene core and the B atom of the boronic functionalities, with a stronger contribution from the ester with respect to the acid one. The interaction results in the lowering of the band gap of 2D BP by around 0.10 eV. Furthermore, we demonstrated through Raman spectroscopy an increased stability toward oxidation in air of 2D BP in the adducts in the solid state (more than 6 months). The modification of the electronic structure at the interface between 2D BP and a conjugated organic molecule through noncovalent stabilizing interactions mediated by the B atom is particularly appealing in view of creating heterojunctions for optoelectronic, photonic, and chemical sensing applications. KEYWORDS: phosphorene, black phosphorus, 2D materials, non-covalent functionalization, boron



INTRODUCTION

the discovery of other semiconducting 2D materials, such as 2D black phosphorus (BP) whose single layer is named phosphorene.2 The most appealing property of 2D BP is surely the direct thickness-dependent band gap, which changes from 0.33 eV for

The ongoing technological race toward low-dimensional and high-performing nanodevices has opened the way for the study of nanostructured materials. In particular, the discovery of graphene1 in 2004, which resulted in the Nobel Prize in Physics in 2010, has drawn the attention of materials scientists toward two-dimensional (2D) materials. The exploration of new 2D materials starting from other elements of the periodic table led to © 2019 American Chemical Society

Received: March 11, 2019 Accepted: May 29, 2019 Published: May 29, 2019 22637

DOI: 10.1021/acsami.9b04344 ACS Appl. Mater. Interfaces 2019, 11, 22637−22647

ACS Applied Materials & Interfaces



the bulk to 2.26 eV for the single layer.3 Owing to the high optical absorbance efficiency, combined with a high anisotropy, and the direct band gap, this material is a promising candidate for optoelectronic and photonic applications in the near-mid infrared range.4 One of the major drawbacks of 2D BP is its high instability due to oxidation under ambient conditions.5 The covalent and noncovalent functionalization of 2D BP with organic molecular materials is a very common approach that is used to protect its surface from air degradation6,7 and also to tailor its electrical properties.8−10 Till date, there are some reports in the literature dealing with the covalent functionalization of 2D BP, such as the covalent bonding of phenyl groups,11 of fluorine atoms12 and of C60 molecules13 or the reaction with metal-halide reagents,14 organic radicalic species,15 and conjugated polymers.16 Alternatively, the protection of 2D BP can be achieved through noncovalent interactions, which are particularly relevant to phosphorene, for which covalent bonds alter strongly the electronic structure. For example, literature reports on the adsorption of 7,7,8,8-tetracyano-p-quinodimethane and various perylene derivatives on 2D BP flakes.8,9,17−19 This strategy allowed to improve the performance of 2D BP-based devices17,20−22 and to implement heterojunctions in biophotonics17 or optoelectronics.23 Indeed, hybrid inorganic semiconductor (2D BP)/organic (molecular) systems combine together the high carrier mobility and high excitation densities of inorganic semiconductors and the strong light−matter coupling featured by the organic compounds. Molecular phosphines easily react with the suitable boron precursor leading to either phosphinoboranes (containing a covalent B−P double bond) or borylphosphines (having a covalent single B−P bond):24 thus, 2D BP is expected to show a high affinity and/or a specific chemical reactivity with compounds containing B units. Moreover, 2D BP flakes display relevant interactions with planar and conjugated organic molecules on its surface by chemisorption. However, a detailed study on the functionalization of 2D BP by means of conjugated organo-boron compounds is still lacking. In this work, this twofold approach is defined, implemented, and validated by using boronic derivatives (acid or ester) of pyrene (Py, PBA, and PBE, Figure 1). Pyrene derivatives have been selected because they

Research Article

EXPERIMENTAL SECTION

General. Pyrene (Sigma-Aldrich, 98% purity), pyren-1-boronic acid (Sigma-Aldrich, 95% purity), and all reagents and solvents employed for the synthesis and preparations reported below were used as received, without further purification. Synthesis of BP. BP was synthesized following a procedure described in the literature.7,32 Red phosphorus was used as a starting material, and Sn and SnI4 were used as catalysts. The solids were heated in a quartz vial at a temperature of 650 °C for 3 days and then cooled down to room temperature very slowly to obtain high quality crystals of BP. Synthesis of 4,4,5,5-Tetramethyl-2-pyren-1-yl-1,3,2-dioxaborolane. 4,4,5,5-Tetramethyl-2-pyren-1-yl-1,3,2-dioxaborolane (PBE) was synthesized as reported in the literature.33 Briefly, to a well degassed solution of 1-bromopyrene (180 mg, 0.64 mmol) in anhydrous 1,4-dioxane (1.6 mL) were added bis(pinacolato) diboron (245 mg, 0.96 mmol), potassium acetate (KOAc, 138 mg, 1.41 mmol), and 7 mg of [1,1′-Bis(diphenylphosphino)ferrocene]dichloropalladium(II) (Pd(dppf)Cl2). The resulting mixture was stirred at 90 °C for 5 h under N2. After cooling, the solvent was removed under pressure, and CH2Cl2 was added. The organic layers were washed with H2O, dried over Na2SO4, filtered, and concentrated. The crude was purified on a silica gel column, eluting with 0−10% CH2Cl2 in cyclohexane, affording a yellowish solid (148 mg, 0.45 mmol, yield = 70%). The 1H and 13C NMR spectra were recorded on Varian spectrometers (500 MHz for 1H and 125 MHz for 13C) at 298 K. Deuterated chloroform was used as the solvent for NMR experiments. 1 H and 13C chemical shift values (δ) are reported in ppm with the residual solvent signal as the internal standard (chloroform at 7.26 and 77.00 ppm, respectively). Abbreviation for signal coupling are as follows: s, singlet; d, doublet; t, triplet; q, quartet; quin., quintet; m, multiplet; br., broad. Flash chromatography was performed on Teledyne Isco CombiFlash Rf 200 using RediSep Normal-phase Silica Flash Columns (230−400 mesh). Thin-layer chromatography was performed on silica gel 60 F254 plastic sheets. Purity of PBE was determined by high-performance liquid chromatography (HPLC)ultraviolet (UV) analysis (Waters 600 HPLC instrument connected to a photodiode array detector 996). 1 H NMR (500 MHz, CDC13): δ 1.53 (s, 12H), 8.01−8.25 (m, 7H), 8.60 (d, J = 9.5 Hz, 1H), 9.14 (d, J = 12.0 Hz, 1H); 13C NMR (125 MHz, CDC13): δ 25.03, 83.85, 124.04, 124.36, 124.58, 125.15, 125.30, 125.65, 127.44, 127.74, 128.00, 128.49, 130.74, 131.09, 133.43, 133.84, 136.42. BP Exfoliation in Solution. Five milligrams of BP crystals were suspended in a vial with 5.0 mL of dry dimethylsulfoxide, and deoxygenated distilled water (3−5 μL) was added. The vial was sealed under an inert atmosphere and sonicated in an ultrasonic bath (200 W, 37 KHz) for 120 hours.34 The Raman spectrum of the 2D BP flakes drop-casted onto a freshly cleaved silicon substrate is reported in the Supporting Information (Figure S1). Preparation of the Adducts in Suspension. The 2D BP suspension was dried and then dispersed again in anhydrous tetrahydrofuran (THF) (1 mg/mL), sonicated in an ultrasonic bath for 10 min, and centrifuged at 1000 rpm for 1 h. The supernatant fraction was separated from the solid and centrifuged at 2000 rpm for 1 h. The supernatant fraction was again separated from the solid and mixed with either Py, PBA, or PBE in anhydrous THF to obtain a suspension with a total dye concentration of 3 mM to have an excess of 2D BP. The resulting suspension was sonicated in an ultrasonic bath for 10 min and diluted to a total dye concentration of 30 μM. All solutions and suspensions were handled in an inert argon atmosphere. Preparation of the Adducts on Solid Substrates. BP was mechanically exfoliated with a blue foil tape on freshly cleaved quartz substrates to produce supported 2D BP. The samples were functionalized by dip coating in anhydrous THF solutions of either Py, PBA, or PBE 10 mM. Both exfoliation and functionalization steps were performed in a glovebox with an inert N2 environment (O2 < 0.1 ppm, H2O < 0.1 ppm).

Figure 1. Chemical structure of Py, PBA, and PBE.

are a class of largely investigated solvatochromic, conjugated, planar fluorophores, which find wide use in optoelectronic and sensing applications.24,25 They are characterized by a long excited-state lifetime, the values of which are affected by the chemical environment, such as the nature of the solvent and the presence of oxygen in solution.26−28 Furthermore, in concentrated solution and/or in the solid state, pyrene derivatives can lead to the formation of characteristic excimers.29−31 Thus, we used computational modelling, NMR and Raman spectroscopies, fluorescence emission, and microscopy to study the interactions of pyrene and/or its derivatives with 2D BP. 22638

DOI: 10.1021/acsami.9b04344 ACS Appl. Mater. Interfaces 2019, 11, 22637−22647

Research Article

ACS Applied Materials & Interfaces

Figure 2. Optimized structure of the adducts of phosphorene with Py (a), PBA (b), and PBE (c). Scanning Electron Microscopy. The morphology of the adducts prepared on quartz, as described above, was studied by scanning electron microscopy (SEM) using a dual beam instrument, Gaia 3 TESCAN. The elemental analysis on the surface of the samples was carried out by an energy-dispersive X-ray (EDAX) system interfaced with SEM (AMETEK, Mahwah NJ, USA, software TEAM EDS Basic Suite). Atomic Force Microscopy. The morphology of the adducts supported on solid quartz substrates, as described above, was studied also by atomic force microscopy (AFM) in the PeakForce Tapping mode on a Bruker MultiMode 8 microscope with a NanoScope V controller and a J scanner. Raman Spectroscopy. The Raman spectra on liquid exfoliated samples deposited on silicon substrates were studied using a Renishaw inVia confocal Raman microscope equipped with a 532 nm laser and a motorized stage for 2D mapping of samples. A laser spot size of approximately 2 μm in diameter was used. Confocal micro-Raman spectroscopy of the samples prepared with mechanical exfoliation and functionalized by dip coating was performed on single flakes of approximately the same thickness, for better comparison, by using a Renishaw micro-Raman 1000 system equipped with a 488 nm laser and a motorized stage for 2D mapping of samples. The laser beam was focused through a 20× objective. Ultraviolet−Visible and Fluorescence Spectroscopies. Ultraviolet−visible (UV−vis) spectra were obtained with a Hewlett-Packard mod. 8453 diode array spectrophotometer. Static and time-resolved fluorescence emission measurements were performed on a Jobin YvonSpex Fluoromax 4 spectrofluorimeter using time-correlated singlephoton counting technique. A NanoLED (λ = 390 nm) has been used as the excitation source. Emission spectra were not corrected for the absorbance of the samples. The fluorescence decay profiles were analyzed using iterative fitting algorithm (Levenberg−Marquardt, simplex), achieving an instrumental resolution (corresponding to the minimum measurable time value) of about few tens of picoseconds. Spectroscopy was done on adduct suspensions, in anhydrous and degassed THF, prepared as described above. Fluorescence and Optical Microscopies. For these experiments, the adducts supported on solid quartz substrates were used. The samples were prepared as described above. The images were recorded with a Nikon TE2000 optical microscope, equipped with a 100× oil immersion objective with a 1.20 numerical aperture and a CCD camera from Hamamatsu as the detector. Excitation was done through a mercury arc lamp with 300 W of power equipped with a low-pass filter for UV excitation and with no filters for visible excitation. Emission was recorded in the whole UV−vis spectrum. NMR. NMR spectra were recorded on a Varian InfinityPlus spectrometer, working at a Larmor frequency of 400.35, 162.07, and 128.45 MHz for 1H, 31P, and 11B, respectively, using a direct excitation pulse sequence under high-power decoupling from 1H nuclei. The suspensions of BP nanoflakes (4.5 mg in 0.5 mL of THF-d5) or 2D BP/ PBA (prepared following the above described procedure for the adducts in suspension, but using 0.5 mL of THF-d5) were directly put

into zirconia rotors for solid-state NMR, with a diameter of 7.5 mm, equipped with a caps kit for liquid samples. Spectra were recorded both under static conditions and under magic angle spinning (MAS) at a variable frequency up to 3 kHz. For 31P spectra, the duration of the excitation pulse was 5 μs, a recycle delay of 30 s was used, and 100−400 transients were accumulated for the suspension of pristine BP nanoflakes and that of 2D BP/PBA. The chemical shift scale was referred to the signal of H3PO4 (80%) at 0 ppm. For 11B spectra, the duration of the excitation pulse was 2.5 μs, a recycle delay of 30 s was used, and 200−800 transients were accumulated under static and MAS conditions. The chemical shift scale was referred to the signal of NaBH4 at −42.06 ppm. Density Functional Theory Simulations. The optimized geometries and energetics of the bare and functionalized phosphorene surface have been obtained at the PBE-D3 density functional theory (DFT) level of theory,35,36 with the inclusion of the dispersion contribution by using the CRYSTAL17 software package.37 A selected supercell of 64 phosphorus atoms has been adopted for performing the optimization of both the atomic positions and lattice parameters. To take into account the effect of oxidation of the surface of 2D BP on the absorption of the organic moiety, the surface has been optimized together with the adduct with PBA. The optimized structure of the adduct is reported in Figure S2. The TZVP basis set has been used for all atomic species.38 Preliminary calculations on molecular models with hydrogen atoms at the boundary sides have been performed by using Gaussian 16 software.39 The complete list of all optimized structural and energy parameters is available in the Supporting Information.



RESULTS AND DISCUSSION DFT Calculations. To understand the nature of the interaction between 2D BP and pyrene (Py) or its boronic derivatives (PBA or PBE), we carried out a detailed solid-state DFT computational analysis. All efforts to optimize a covalent adduct involving the boron center and phosphorene failed because of electronic and steric reasons, recently deeply illustrated from a computational viewpoint.40,41Hence, we focused our attention on the noncovalent possible interaction between the organic moiety and the material.17 The formation of the first adduct, phosphorene/Py in Figure 2a, was estimated to be favored by −2.2 eV. The P−C distances between the organic moiety and the 2D material are within 3.20−3.40 Å, confirming the noncovalent nature of the involved interactions. Although noncovalently bonded, the structures of both phosphorene and pyrene are rearranged in the phosphorene/Py adduct. Indeed, the C−C distances of pyrene are slightly lengthened, being in between 1.35 and 1.42 Å in the pristine pyrene moiety and in between 1.37 and 1.44 Å in the adduct (Figure S3a in the Supporting Information). On the opposite, 22639

DOI: 10.1021/acsami.9b04344 ACS Appl. Mater. Interfaces 2019, 11, 22637−22647

Research Article

ACS Applied Materials & Interfaces

NMR Spectroscopy. To confirm the DFT predictions, the interaction between 2D BP and PBA was investigated by means of multinuclear NMR spectroscopy directly on the 2D BP/PBA suspension. To the best of our knowledge, this is the first time that NMR has been applied on suspensions of exfoliated BP, allowing chemical features and interactions of BP nanoflakes to be characterized directly in the dispersed state. Figure 4 reports the 31P NMR spectra of pristine 2D BP in suspension, both under static and MAS conditions.

the P−P nonbonding distances of phosphorene are shortened, especially the ones directly interacting with the organic moiety, passing from an estimated medium value of 3.5 Å for pristine phosphorene to 3.2 Å for the first and fourth channels of phosphorene in the adduct (indicated as first ch. and fourth ch. in Figure 2). The favourable noncovalent P−B interaction gives a relevant contribution to the further stabilization of the other two adducts, phosphorene/PBA in Figure 2b and phosphorene/PBE in Figure 2c. Indeed, the interaction energy between phosphorene and the organic molecules increases by −0.30 eV in the case of PBA and by −0.50 eV in the case of PBE, compared to Py. This occurs even though a noncovalent P−B bond is formed, the calculated P−B distance being as long as 3.00 Å. This is in perfect agreement with NMR results (see the following section). However, the main contribution to the stabilization energy of the moieties of the adduct has still to be ascribed to the noncovalent interaction between the conjugated core of pyrene and the phosphorus atoms. The slight increase (−0.1 eV) in the stabilization energy of phosphorene/PBE compared to phosphorene/PBA could be attributed to the dependence of the dispersion interaction on the number of involved atoms. Therefore, the phosphorene/PBE adduct is the most stable within all, even though the optimized P···B and H···P distances in this adduct are longer than in phosphorene/PBA likely due to the higher steric hindrance of the ester around the B atom in PBE with respect to the OH groups in PBA. In addition, for phosphorene/PBA and phosphorene/PBE adducts, the interaction leads to a geometrical rearrangement of the two organic moieties, evidenced by a variation of the interatomic distances (Figure S3b,c in the Supporting Information). From an electronic viewpoint, the band gap calculated for the phosphorene/PBA and phosphorene/PBE adducts with the PBE-D3 functional decreases by 0.10 and 0.08 eV with respect to bare phosphorene. The corresponding density of state (DOS) plots of the adduct (Figure 3) reveal that the band gap lowering

Figure 4. 31P NMR spectra at 9.4 T of the suspension of pristine BP nanoflakes in THF-d5 under (a) static and (b) MAS (frequency of 3 kHz) conditions. (c) 31P NMR spectrum at 9.4 T of the suspension of 2D BP/PBA in THF-d5 under MAS (frequency of 3 kHz) conditions. Asterisks indicate spinning sidebands.

The static spectrum (Figure 4a) is constituted by a broad signal, similar to that observed for pristine solid bulk BP,42,43 which arises from the anisotropic spin interactions of 31P nuclei. The anisotropy is partially averaged out by MAS, and indeed in the 31P MAS spectrum (Figure 4b), it is possible to observe an intense signal of 2D BP at an isotropic chemical shift of about 18.5 ppm, similar to that observed for bulk BP crystals and 2D BP, mechanically exfoliated or embedded in polymer hybrids,6,15,32,42,44,45 together with spinning sidebands arising from the residual anisotropy of the spin interactions. In the 31P NMR spectrum of the 2D BP/PBA suspension (Figure 4c), the signal of 2D BP does not show any modification, suggesting that substantial chemical functionalization of 2D BP with PBA did not occur. Two weak narrow peaks, ascribable to minor amounts of oxidized (phosphate and phosphite) species, appear between 0 and 10 ppm.6,46 11 B NMR spectra of the 2D BP/PBA suspension (Figure 5), under static and MAS conditions, confirm the substantial absence of covalent functionalization of 2D BP. The static spectrum (Figure 5a) shows only one approximately Lorentzian peak at a chemical shift of about 30.5 ppm typical of boronic acids in solution.47 As expected, this peak slightly narrows in the MAS spectrum (Figure 5b), in which, on the other hand, no

Figure 3. Band structure and DOS for the adduct between phosphorene and PBA. In the DOS plot, the contribution highlighted in blue is due to the PBA moiety.

is due to a contribution of the organic molecules in the DOS region next to the top of the valence band. This contribution arises from a partial mixing of the π* of the pyrene moiety with the 2p orbitals of the lone pairs of phosphorene, as observed in the DOS plots (Figure 3). This is in accordance with the structural data of the optimized adducts, which, as described above, feature an elongation of the C−C bonds of the pyrene moiety and a shortening of the P−P distances within the channels of phosphorene. 22640

DOI: 10.1021/acsami.9b04344 ACS Appl. Mater. Interfaces 2019, 11, 22637−22647

Research Article

ACS Applied Materials & Interfaces

Figure 6. (a) UV−vis extinction spectra of pristine 2D BP in suspension (green line) and of the adducts of 2D BP with Py (black line), PBA (red line), and PBE (blue line) in a suspension of THF, 3 × 10−5 M. (b) Fluorescence spectra of the pristine fluorophores Py (black dotted line), PBA (red dotted line), and PBE (blue dotted line) and their adducts with 2D BP in THF suspensions (same colors, continuous lines); λexc = 260 nm.

Figure 5. 11B NMR spectra at 9.4 T of the suspension of 2D BP/PBA in THF-d5 under (a) static and (b) MAS (frequency of 3 kHz) conditions.

additional signals appear. Analogous results were obtained from 13 C and 1H NMR spectra (Figure S4 in the Supporting Information). Photophysical Investigation. Photophysical investigation has been carried out, thanks to the presence of luminescent probes. The interactions between 2D BP and pyrene-based molecules have been explored in a medium polarity solvent (THF) to have a good compromise between the solubility of the fluorophores and a good dispersion of 2D BP. For comparison with NMR, the dye concentration for photophysical investigation has been fixed at 3 × 10−5 M. In any case, control experiments have been performed on more dilute samples (3 × 10−6 M), and no difference in photophysical behavior has been observed. The extinction spectrum of pristine 2D BP in THF suspension is broad and featureless, except for a peak at 275 nm, and covers the whole UV−vis region. It is worth noting that, working with suspensions, in Figure 6, we report the “extinction spectra” measured by a conventional spectrophotometer as the sum of the absorption and scattering contributes, whereas absorbance is typically reported for the electronic spectra of the pure fluorophores in solution (see Figure S5 in the Supporting Information). The three pristine pyrene derivatives show, in the electronic spectra, absorption bands with a characteristic vibronic fine structure between 300 and 350 nm. According to the literature, the absorption maxima of PBA and PBE present a slight red shift compared to Py.48,49 Thus, the optical band gaps (Eopt g ) of the three samples, calculated from the λ onsets of the absorption spectra, vary in the trend Py > PBA > PBE. The electronic spectra of the three adducts present profiles corresponding to the sum of the spectra of the single components (Figures 6a and S5 in the Supporting Information). This excludes the formation of absorbing charge transfer (CT) electronic states. According to the literature, the emission spectra of the three pure fluorophores (Figure 6b), taken as a reference, show structured bands spanning from 380 to 500 nm typical for this class of compounds.29,48,50 For Py, the most intense vibronic band is III, whereas the presence of the boronic functionalities in PBA and PBE causes a significant enhance-

ment of the higher energy vibronic bands (I) at the expense of the lower ones. The II vibronic band is slightly detectable. The emission spectra of the adducts were taken by exciting the suspensions in the UV region, where both 2D BP and the investigated fluorophores absorb. However, no emission in the red−NIR range, characteristic of 2D BP, was observed. In the presence of 2D BP, the emission spectra of all pyrenes show a slight bathochromic shift of the bands (by +2 nm) and an evident change in the intensity ratio of the first to third vibronic emission bands (I/III ratio), more evident for PBA and PBE. Since for this class of compound, the I/III intensities ratio is very sensitive to the solvent polarity,24,29,48,49,51−56 the different emission profiles of the fluorophores in the adducts with respect to free ones can be ascribed to a change of the local environment polarity because of the presence of 2D BP. Even if the emission spectra are not corrected for the absorbance of the samples, an evident increase of the emission intensity has been observed. To well define this aspect and better investigate the nature of the adducts, time-resolved fluorescence emission experiments have been carried out. Lifetimes and relative amplitude values are collected in Table 1. The fluorescence emission decay of pristine Py presents a biexponential behavior with a long lifetime value of about 7.1 ns and a shorter one of about 1.2 ns with almost equal amplitudes (Table 1). Differently, the PBA and PBE fluorescence decay profiles are well-fitted by triexponential curves (Table 1). Decays of both PBA and PBE present a short (2.1 and 1.6 ns) and long (8.8 and 5.0 ns) component, similar to those of Py; but differently from Py, they also present a third and longer lifetime value. In detail, pristine PBA has a τ3 of 53.3 ns with 0.19 of amplitude, and PBE has a τ3 of 40.1 ns with an amplitude as high as 0.61 (Table 1). A similar lifetime (around 50−60 ns) is reported in the literature for the fluorescence decay of molecular Py in solution.29,57 As already reported, the apparent discrepancy in lifetime values can be due to the nature of the solvent and the presence of excimers. To confirm the absence of excimers in solution, emission spectra at various sample concentrations and excitation spectra at different emission 22641

DOI: 10.1021/acsami.9b04344 ACS Appl. Mater. Interfaces 2019, 11, 22637−22647

Research Article

ACS Applied Materials & Interfaces

To note, the alignment of the energetic levels of 2D BP with Py, PBA, and PBE would favor a quenching of the emission of the organic molecules by energy or CT process to 2D BP.59,60 However, this did not occur. Pyrene derivatives and their nanohybrids are used in the literature as simple, cheap biocompatible fluorescent compounds for the detection of molecular oxygen.61 Moreover, 2D BP has attractive physical properties for high-performance chemical sensing applications, such as a high molecular adsorption energy and surface-to-volume ratio.62−64 In addition, by exploiting the remarkable attitude to fast oxidation in the presence of oxygen, 2D BP can be used as an efficient getter like other molecular compounds containing phosphorus. In this scenario, we investigated the possibility to use the asdesigned heterosystem as an efficient fluorescent chemosensor for molecular oxygen. As a feasibility test, we performed timeresolved fluorescence spectroscopy on the pyrenic derivatives and the corresponding 2D BP-based heterosystem in solution with and without exposition to air. We first carried out timeresolved fluorescence experiments on the PBA and PBE solutions exposed to air. In both cases, fluorescence decays show after exposure to air overnight a biexponential behavior with lifetime values, which are shorter compared to PBA and PBE in the degassed solution. According to the literature, the decrease of the fluorescence lifetime is ascribable to the presence of oxygen dissolved in the system.48 Moreover, the fluorescence of the 2D BP/PBA and 2D BP/PBE adducts in solution exposed to air under the same conditions presents shorter lifetime components in the case of a triexponential decay with respect to the corresponding adducts in an inert atmosphere. For a quantitative comparison between the pyrenic compounds alone and on 2D BP-based adduct, the weighted average lifetimes (τav) are calculated as

Table 1. Fitting Parameters of the Time-Resolved Fluorescence Decays; λexc = 390 nm λem = 450 nm sample Py PBA

PBE

PBA + O2

PBE + O2

τi (ai)a 1.2 ± 0.1 (0.55) 7.1 ± 0.2 (0.45) 2.1 ± 0.1 (0.42) 8.8 ± 0.3 (0.39) 53.3 ± 0.6 (0.19)b 1.6 ± 0.1 (0.18) 5.0 ± 0.2 (0.21) 40.1 ± 0.5 (0.61) 2.0 ± 0.1 (0.17) 9.6 ± 0.3 (0.83)b 2.3 ± 0.1 (0.04) 9.6 ± 0.2 (0.96)

sample 2D BP/Py 2D BP/PBA

2D BP/PBE

2D BP/PBA + O2

2D BP/PBE + O2

τi (ai)a 3.1 ± 0.2 (0.64) 12.0 ± 0.4 (0.36) 2.2 ± 0.1 (0.48) 8.7 ± 0.3 (0.36) 61.3 ± 0.1 (0.16)b 2.0 ± 0.2 (0.21) 7.0 ± 0.2 (0.17) 62.2 ± 0.1 (0.62) 1.0 ± 0.1 (0.22) 4.0 ± 0.1 (0.43) 15.0 ± 0.3 (0.25)b 0.7 ± 0.2 (0.24) 3.3 ± 0.2 (0.52) 15.2 ± 0.1 (0.24)

Data were fitted through bi- or triexponential decays where τi and ai are the resulting lifetimes (in ns) and relative amplitudes, respectively. For each fitting, χ2 = 1. bWeighted average lifetimes (τav) are τav(PBA) = 14.4 ns; τav(PBA + O2) = 8.3 ns; τav(2D BP/PBA) = 14.0 ns; τav(2D BP/PBA + O2) = 5.7 ns; τav(PBE) = 25.7 ns; τav(PBE + O2) = 9.3 ns; τav(2D BP/PBE) = 40.2 ns; τav(2D BP/PBE + O2) = 5.5 ns. a

wavelengths have been performed. In any case, no evidence of the characteristic excimer bands has been observed. Considering these evidences, the bi- or triexponential decay fitting curves are clearly indicative of the presence in solution of two or three different species that we can justify by only speculating on the monomeric, dimeric, or oligomeric nature of the fluorophore, as previously reported for pyrene sulphonate derivatives.58 The fluorescence decays of the three adducts were fitted with biexponential decay for Py and triexponential decays for PBA and PBE, respectively. In all three cases, the presence of 2D BP leads to an increase in the fluorescence lifetime values. In particular, for Py, both τi values increase (τ1 passing from 1.2 to 3.1 ns and τ2 passing from 7.1 to 12.0 ns), whereas the relative amplitudes are not very much affected. For PBA and, even more, for PBE, the stabilization effect is evidenced as an increase in the long lifetime component, τ3, which increases from 53.3 to 61.3 ns for PBA and from 40.1 to 62.2 ns for PBE, whereas τ1 and τ2 and the relative amplitudes remain almost constant. This stabilization, which becomes stronger with the trend Py < PBA < PBE, is in perfect accordance with theoretical calculations, with optical and fluorescence microscopy, and with Raman spectroscopy experiments (see following sections). Overall, the interaction of the pyrene derivatives with 2D BP has a stabilization effect on their emissive electronic excited states. This could be due to a decrease in the nonradiative constant of deactivation of the S1 excited state through the suppression of some vibrational modes (see the Raman spectroscopy section). In detail, the stabilization effect becomes stronger when passing from Py to PBA and to PBE. This could be due either to the different energetic displacement and geometry of the vibronic levels of Py, PBA, and PBE or to a different coupling of these vibronic states with those of 2D BP. Indeed, while in the 2D BP/Py adduct, the coupling occurs only through the pyrene conjugated core, in 2D BP/PBA and 2D BP/ PBE adducts, the coupling is also mediated by the boronic functionality and, in particular, is stronger through the ester group of PBE than through the acid group of PBA.

n

τav =

∑ τi ·ai i=0

In detail, τav(PBA + O2) is shorter by 1.7 times compared to τav(PBA), whereas τav(2D BP/PBA + O2) is shortened by 2.6 times compared to τav(2D BP/PBA) (Table 1). Therefore, the signal variation due to the flux of oxygen, passing from a ratio of 1.7 to 2.6, is amplified by 1.5 times in the system containing 2D BP. Such amplification is even higher (2.6) for the 2D BP/PBE system with respect to PBE alone. By degassing again the suspensions of the adducts, the fluorescence decays do not recover their initial state. Moreover, the fluorescence decays of the suspensions of the adducts, prepared with a 2D BP suspension previously exposed to air, show the same behavior as the pristine suspensions of the adducts then exposed to air. Thus, the observed variation of the lifetimes of the 2D BP/PBA and 2D BP/PBE adducts in the presence or absence of oxygen is predominantly due to the interaction of PBA and PBE directly with the oxidized sites of 2D BP (POx). This indicates that, in suspension, oxygen is able to reach the 2D BP surface and oxidize it irreversibly creating POx species. DFT modelling on the 2D BP/PBA system confirmed that the interaction energy between PBA and oxidized 2D BP increases by −0.22 eV compared to pristine 2D BP. This is due to the formation of hydrogen bonds between the acid group of PBA and the POx sites on the surface of 2D BP, as also evidenced by the much shorter O······H distances compared to the P···H ones (2.88 vs 1.95 Å, see DFT section and Figure S2 in the Supporting Information). Analogous considerations can be 22642

DOI: 10.1021/acsami.9b04344 ACS Appl. Mater. Interfaces 2019, 11, 22637−22647

Research Article

ACS Applied Materials & Interfaces

BP flakes is substantially homogeneous, suggesting a homogeneous distribution of the Py, PBA, and PBE molecules on the 2D BP flakes. This was further confirmed by both SEM analysis, correlated with EDAX analysis on the same samples (Figure S6 in the Supporting Information), and by AFM analysis. AFM gave very low root-mean-square roughness values, comprised between 1.20 and 1.50 nm, in the three cases (Figure S7 in the Supporting Information), which is nonetheless slightly higher than the roughness of bare 2D BP (0.2 nm) that is atomically flat.11 Raman Spectroscopy and Aging Studies. Simulations describe that the coupling between 2D BP and pyrene derivatives is energetically favored and mediated by the boronic functionalities, with a stabilization energy that follows the trend 2D BP/Py < 2D BP/PBA < 2D BP/PBE. In addition, the spectroscopic data above discussed evidence that the excited state of the pyrene derivatives is stabilized by 2D BP with the same trend. In this scenario, Raman spectroscopy can shed light on the stabilization of 2D BP when terminated by a layer of Py, PBA, or PBE. The Raman spectra of 2D BP/Py, 2D BP/PBA, and 2D BP/ PBE samples show the three typical peaks of 2D BP at 360, 435, and 460 cm−1 attributed to the A1g, B2g, and A2g vibrational modes, respectively (Figure 8).65 For all three samples, the relative

done for the oxidized 2D BP/PBE system, where stabilization probably occurs through favorable boron−oxygen interactions. The properties of 2D BP as an oxygen getter, the tunable interaction of PBA and PBE with either pristine and oxidized 2D BP, and the fluorescent properties of the pyrenic derivatives allow to envisage the potential use of the proposed heterostructures as fluorescent active materials in the chemosensor for oxygen detection with amplified dynamic range and possibly enhanced sensitivity with respect to the pyrenic derivatives alone. Optical and Fluorescence Microscopies. To get insight into the quality of the 2D BP flakes and their surface and to analyze the amount and coverage of Py, PBA, and PBE molecules on the surface of the functionalized samples, we performed optical and fluorescence microscopies on mechanically exfoliated 2D BP functionalized with either Py, PBA, or PBE by dip coating. The optical microscopy images of pristine and functionalized 2D BP collected under backscattering conditions have very similar contrast (Figure 7a,c,e,g).

Figure 7. Optical (a,c,e,g) and fluorescence (b,d,f,h) microscopy images of the pristine and functionalized 2D BP flakes mechanically exfoliated on quartz: pristine 2D BP flakes (a,b), 2D BP/Py (c,d), 2D BP/PBA (e,f), and 2D BP/PBE (g,h). The scale bar is 50 μm.

Figure 8. Micro-Raman confocal spectroscopy on pristine 2D BP (blue line); 2D BP after 3 months of air exposure (green line); fresh 2D BP/ Py (brown line) or after 3 months (red line) or 6 months (orange line) of air exposure; fresh 2D BP/PBA (violet line) or after 3 months (dark pink line) or 6 months (light pink line) or air exposure; and fresh 2D BP/PBE (black line) or after 3 months (dark gray line) or 6 months (light gray line) of air exposure. A silicon sample (sky-blue line) with the peak at 521 cm−1 is also reported as a reference.

The excitation light is efficiently reflected by the surface of the 2D BP flakes made of atomically flat superimposed sheets, whereas it is transmitted through the transparent substrate (quartz). The same areas of the samples, excited with UV light, give fluorescence images (Figure 7b,d,f,h). The fluorescence microscopy image of pristine 2D BP flakes is almost completely dark (Figure 7b), with the only contrast on the flakes given by the reflection of the residual excitation light passing through the long-pass filter. The fluorescence images of 2D BP flakes functionalized with either Py, PBA, or PBE show almost the same contrast with respect to the background when the images are collected under the same experimental conditions (Figure 7d,f,h). The fluorescence emission of the heterostructures is weak either due to the low-thickness functional layer, which is probably limited to a single or few molecular layers, or to the lower fluorescence quantum yield of the pyrene-based adducts in the solid state compared to that in solution. However, the increase in the collected-light intensity is evaluated in a factor 1.5−2 going from bare to functionalized 2D BP flakes. To note, the fluorescence signal from the surface of the functionalized 2D

intensity of the A1g peak, which corresponds to out-of-plane atomic vibrations, decreases with functionalization, as compared to the relative intensity of the B2g and A2g peaks. Interestingly, the A1g peak decreases with the trend 2D BP > 2D BP/Py > 2D BP/ PBA > 2D BP/PBE. This confirms that the interaction forces involved in the formation of the adducts increase in strength following the trend 2D BP/Py < 2D BP/PBA < 2D BP/PBE, in parallel with the gradual suppression of the vibrational modes of both the organic moiety (as evidenced by time-resolved fluorescence emission) and the inorganic moiety (as evidenced by Raman spectroscopy). Notably, the position and relative intensity of the peaks in the Raman spectra of freshly prepared and aged samples of functionalized 2D BP remain constant in time even after 6 months of exposure to air. On the contrary, the spectrum of a 3 month-aged sample of bare 2D BP is flat and featureless. Thus, 22643

DOI: 10.1021/acsami.9b04344 ACS Appl. Mater. Interfaces 2019, 11, 22637−22647

Research Article

ACS Applied Materials & Interfaces

based heterostructures indicated that the 2D BP/PBA and 2D BP/PBE systems can be implemented as active materials in a fluorescent chemosensor for oxygen detection. Indeed, the two systems show relevant modulation of their fluorescence decay time constants with respect to the pyrenic derivatives alone. Finally, Raman spectroscopy and fluorescence microscopy showed that 2D BP also is stabilized in the adducts in solid state. Specifically, 2D BP when terminated with a layer of pyrene derivatives forms heterostructures, which are stable in air for more than 6 months. The presented study represents an innovative contribution to the understanding of the chemistry at the interface between 2D BP and organoboron compounds. The band gap tailoring and the modulation of the reactivity of 2D BP by means of pyrenic derivatives are interesting tools for designing heterosystems with multifunctional applications, from passivation and protection of 2D semiconductors to the sensing of analytes of interest such as molecular oxygen.

while pristine 2D BP undergoes severe degradation in air, its crystalline structure in the adducts is preserved from oxidation in air for over 6 months.66 Interestingly, the optical images of the 6 month-aged sample reveal that while many of the functionalized 2D BP flakes are preserved, some of them seem to have suffered degradation from the edges. Edges are indeed presumably less covered by the organic layer and thus more prone to air oxidation (Figure S8 in the Supporting Information). This indicates that the organic layer, acting as a physical barrier for water/oxygen and as an energy stabilizing agent, has a strong protecting functionality. To note, in suspension, the intimate interaction between 2D BP flakes and organic compounds, which is at the basis of the protecting functionality, is reasonably much weaker compared to the solid state because of scavenging solvent molecules. In solid state, the protecting functionality seems equal for all three investigated samples. However, the optical and fluorescence microscopy images of the adducts, taken over a period of 1 week, reveal the condensation of droplets of atmospheric water on the surface of the functionalized flakes. The images also reveal that the pyrene derivatives partially solubilize into droplets (Figure S9 in the Supporting Information). In detail, a more rapid condensation of atmospheric water is observed on the surface of the 2D BP/Py and 2D BP/PBA adducts (1−2 days) as compared to 2D BP/PBE (1 week). The condensation of atmospheric water, in the case of bare 2D BP, accelerates the degradation processes through the dissociation of oxidized P atoms from the surface of 2D BP (Figure S10 in the Supporting Information).5,46,67 Therefore, the higher wettability/solubility of the PBA and Py layers with/in atmospheric water, compared to the PBE layer, combined with the slightly lower stabilization energy of their adducts with 2D BP, could cause a lower long-term stability of the corresponding 2D BP/Py and 2D BP/PBA samples compared to the 2D BP/PBE samples. However, the stability in air for at least 6 months is demonstrated for all functionalized 2D BP samples by Raman spectroscopy, which is in good accordance with the most relevant literature up to date.68



ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acsami.9b04344. Raman spectrum of liquid exfoliated 2D BP, absorption spectra of pyrene derivatives in solution, optimized geometry and energy parameters and bond lengths from DFT simulations, 13C and 1H NMR spectra of the adducts, SEM and EDAX analysis of the adducts, and optical, fluorescence and AFM microscopy of the oxidized adducts and oxidized 2D BP (PDF)



AUTHOR INFORMATION

Corresponding Authors

*E-mail: [email protected] (M.A.C.). *E-mail: s.toff[email protected] (S.T.). ORCID

Margherita Bolognesi: 0000-0003-0080-3279 Claudia Ferroni: 0000-0002-7386-1624 Gabriele Manca: 0000-0003-2068-1731 Andrea Ienco: 0000-0002-2586-4943 Maria Caporali: 0000-0001-6994-7313 Lucia Calucci: 0000-0002-3080-8807 Maria Angela Castriciano: 0000-0002-1514-8820 Stefano Toffanin: 0000-0003-4099-8664



CONCLUSIONS The aim of this work was to study the chemical−physical nature of the interactions involved in the formation of adducts of 2D BP and boronic derivatives of a conjugated organic fluorescent molecule (pyrene). DFT simulations predicted a favorable coupling between 2D BP and pyrene derivatives. The main contribution to the stabilization energy arises from the noncovalent interactions between the conjugated core of pyrene and the phosphorus atoms of 2D BP. A smaller but still relevant contribution (up to 0.5 eV) arises from the noncovalent P−B interaction, which is stronger in the boronic ester than in the boronic acid group of the pyrene derivatives. NMR spectra in suspension confirmed the noncovalent nature of the interactions involved in the formation of the adducts. Although noncovalently bonded, modelling showed that both the atomic and electronic structures of 2D BP and of the pyrene derivatives are rearranged in the adducts, resulting in a lowering of the band gap of 2D BP by around 0.10 eV. Time-resolved fluorescence emission data evidenced that, in the adducts, the pyrene derivatives are stabilized in their excited state by 2D BP, with a stabilization effect that is again stronger for the pyrene derivative bearing the boronic ester rather than the boronic acid group, thus following the theoretically predicted trend. Moreover, further investigation in the excitonic dynamics of the 2D BP-

Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS The computing resources and the related technical support used for this work have been provided by CRESCO/ENEAGRID High Performance Computing infrastructure and its staff.69 CRESCO/ENEAGRID High Performance Computing infrastructure is funded by ENEA, the Italian National Agency for New Technologies, Energy and Sustainable Economic Development and by Italian and European research programmes; see http://www.cresco.enea.it/english for information. We also acknowledge the ISCRA-CINECA HP grants “HP10CFMSCC and HP10C2QI78” and CREA (Centro Ricerche Energia e Ambiente) of Colle Val D’Elsa (Siena, Italy) for computational resources. The authors thank the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation program (grant agreement no. 670173) for funding 22644

DOI: 10.1021/acsami.9b04344 ACS Appl. Mater. Interfaces 2019, 11, 22637−22647

Research Article

ACS Applied Materials & Interfaces

(17) Abellán, G.; Lloret, V.; Mundloch, U.; Marcia, M.; Neiss, C.; Gö rling, A.; Varela, M.; Hauke, F.; Hirsch, A. Noncovalent Functionalization of Black Phosphorus. Angew. Chem., Int. Ed. 2016, 55, 14557−14562. (18) Wild, S.; Lloret, V.; Vega-Mayoral, V.; Vella, D.; Nuin, E.; Siebert, M.; Koleśnik-Gray, M.; Löffler, M.; Mayrhofer, K. J. J.; Gadermaier, C.; Krstić, V.; Hauke, F.; Abellán, G.; Hirsch, A. Monolayer Black Phosphorus by Sequential Wet-Chemical Surface Oxidation. RSC Adv. 2019, 9, 3570−3576. (19) Artel, V.; Guo, Q.; Cohen, H.; Gasper, R.; Ramasubramaniam, A.; Xia, F.; Naveh, D. Protective Molecular Passivation of Black Phosphorous. npj 2D Mater. Appl. 2017, 1, 6. (20) Ryder, C. R.; Wood, J. D.; Wells, S. A.; Hersam, M. C. Chemically Tailoring Semiconducting Two-Dimensional Transition Metal Dichalcogenides and Black Phosphorus. ACS Nano 2016, 10, 3900−3917. (21) Du, Y.; Yang, L.; Zhou, H.; Ye, P. D. Performance Enhancement of Black Phosphorus Field-Effect Transistors by Chemical Doping. IEEE Electron Device Lett. 2016, 37, 429−432. (22) Telesio, F.; Passaglia, E.; Cicogna, F.; Costantino, F.; SerranoRuiz, M.; Peruzzini, M.; Heun, S. Hybrid Nanocomposites of 2D Black Phosphorus Nanosheets Encapsulated in PMMA Polymer Material: New Platforms for Advanced Device Fabrication. Nanotechnology 2018, 29, 295601. (23) Liu, F.; Zhu, C.; You, L.; Liang, S.-J.; Zheng, S.; Zhou, J.; Fu, Q.; He, Y.; Zeng, Q.; Fan, H. J.; Ang, L. K.; Wang, J.; Liu, Z. 2D Black Phosphorus/SrTiO3 -Based Programmable Photoconductive Switch. Adv. Mater. 2016, 28, 7768−7773. (24) Figueira-Duarte, T. M.; Müllen, K. Pyrene-Based Materials for Organic Electronics. Chem. Rev. 2011, 111, 7260−7314. (25) Singh, M.; Holzinger, M.; Biloivan, O.; Cosnier, S. 3DNanostructured Scaffold Electrodes Based on Single-Walled Carbon Nanotubes and Nanodiamonds for High Performance Biosensors. Carbon 2013, 61, 349−356. (26) Craig, B. B.; Kirk, J.; Rodgers, M. A. J. Fluorescence Decay Parameters for Pyrene (S1) in Micellar and Homogeneous Liquid Dispersions. Chem. Phys. Lett. 1977, 49, 437−440. (27) van Bünau, G.; Birks, J. B. Photophysics of Aromatic Molecules; Wiley-Interscience: London, 1970. (28) Murov, S. L.; Carmichael, I.; Hug, G. L. Handbook of Photochemistry, 2nd ed.; Marcel Dekker, 1993. (29) Kalyanasundaram, K.; Thomas, J. K. Environmental Effects on Vibronic Band Intensities in Pyrene Monomer Fluorescence and Their Application in Studies of Micellar Systems. J. Am. Chem. Soc. 1977, 99, 2039−2044. (30) Mohr, A.; Talbiersky, P.; Korth, H.-G.; Sustmann, R.; Boese, R.; Bläser, D.; Rehage, H. A New Pyrene-Based Fluorescent Probe for the Determination of Critical Micelle Concentrations. J. Phys. Chem. B 2007, 111, 12985−12992. (31) Winnik, F. M. Photophysics of Preassociated Pyrenes in Aqueous Polymer Solutions and in Other Organized Media. Chem. Rev. 1993, 93, 587−614. (32) Lange, S.; Schmidt, P.; Nilges, T. Au3SnP7@Black Phosphorus: An Easy Access to Black Phosphorus. Inorg. Chem. 2007, 46, 4028− 4035. (33) Seifert, S.; Schmidt, D.; Shoyama, K.; Würthner, F. Base-Selective Five- versus Six-Membered Ring Annulation in Palladium-Catalyzed CC Coupling Cascade Reactions: New Access to Electron-Poor Polycyclic Aromatic Dicarboximides. Angew. Chem., Int. Ed. 2017, 56, 7595−7600. (34) Serrano-Ruiz, M.; Caporali, M.; Ienco, A.; Piazza, V.; Heun, S.; Peruzzini, M. The Role of Water in the Preparation and Stabilization of High-Quality Phosphorene Flakes. Adv. Mater. Interfaces 2016, 3, 1500441. (35) Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865−3868. (36) Grimme, S.; Antony, J.; Ehrlich, S.; Krieg, H. A Consistent and Accurate Ab Initio Parametrization of Density Functional Dispersion Correction (DFT-D) for the 94 Elements H-Pu. J. Chem. Phys. 2010, 132, 154104.

the project PHOSFUN by an ERC Advanced Grant to M.P. and “MatISSE” Materiali Innovativi e Sostenibili per la Salute e l’Energia, CUP G77B17000180009, grant to M.T. Thanks are also expressed to E. Berretti and A. Lavacchi who contributed to the SEM and EDAX analyses of the solid samples.



REFERENCES

(1) Novoselov, K. S.; Geim, A. K.; Morozov, S. V.; Jiang, D.; Zhang, Y.; Dubonos, S. . V.; Grigorieva, I. V.; Firsov, A. A. Electric Field Effect in Atomically Thin Carbon Films. Science 2004, 306, 666−669. (2) Peruzzini, M.; Bini, R.; Bolognesi, M.; Caporali, M.; Ceppatelli, M.; Cicogna, F.; Coiai, S.; Heun, S.; Ienco, A.; Benito, I. I.; Kumar, A.; Manca, G.; Passaglia, E.; Scelta, D.; Serrano-Ruiz, M.; Telesio, F.; Toffanin, S.; Vanni, M. A Perspective on Recent Advances in Phosphorene Functionalization and Its Applications in Devices. Eur. J. Inorg. Chem. 2019, 1476−1494. (3) Liu, H.; Neal, A. T.; Zhu, Z.; Luo, Z.; Xu, X.; Tománek, D.; Ye, P. D. Phosphorene: An Unexplored 2D Semiconductor with a High Hole Mobility. ACS Nano 2014, 8, 4033−4041. (4) Liu, H.; Hu, K.; Yan, D.; Chen, R.; Zou, Y.; Liu, H.; Wang, S. Recent Advances on Black Phosphorus for Energy Storage, Catalysis, and Sensor Applications. Adv. Mater. 2018, 30, 1800295. (5) Island, J. O.; Steele, G. A.; van der Zant, H. S. J.; CastellanosGomez, A. Environmental Instability of Few-Layer Black Phosphorus. 2D Mater 2015, 2, 011002. (6) Passaglia, E.; Cicogna, F.; Costantino, F.; Coiai, S.; Legnaioli, S.; Lorenzetti, G.; Borsacchi, S.; Geppi, M.; Telesio, F.; Heun, S.; Ienco, A.; Serrano-Ruiz, M.; Peruzzini, M. Polymer-Based Black Phosphorus (BP) Hybrid Materials by in Situ Radical Polymerization: An Effective Tool to Exfoliate BP and Stabilize BP Nanoflakes. Chem. Mater. 2018, 30, 2036−2048. (7) Caporali, M.; Serrano-Ruiz, M.; Telesio, F.; Heun, S.; Nicotra, G.; Spinella, C.; Peruzzini, M. Decoration of Exfoliated Black Phosphorus with Nickel Nanoparticles and Its Application in Catalysis. Chem. Commun. 2017, 53, 10946−10949. (8) Wang, C.; Niu, D.; Liu, B.; Wang, S.; Wei, X.; Liu, Y.; Xie, H.; Gao, Y. Charge Transfer at the PTCDA/Black Phosphorus Interface. J. Phys. Chem. C 2017, 121, 18084−18094. (9) He, Y.; Xia, F.; Shao, Z.; Zhao, J.; Jie, J. Surface Charge Transfer Doping of Monolayer Phosphorene via Molecular Adsorption. J. Phys. Chem. Lett. 2015, 6, 4701−4710. (10) Jing, Y.; Tang, Q.; He, P.; Zhou, Z.; Shen, P. Small Molecules Make Big Differences: Molecular Doping Effects on Electronic and Optical Properties of Phosphorene. Nanotechnology 2015, 26, 095201. (11) Ryder, C. R.; Wood, J. D.; Wells, S. a.; Yang, Y.; Jariwala, D.; Marks, T. J.; Schatz, G. C.; Hersam, M. C. Covalent Functionalization and Passivation of Exfoliated Black Phosphorus via Aryl Diazonium Chemistry. Nat. Chem. 2016, 8, 597−602. (12) Tang, X.; Liang, W.; Zhao, J.; Li, Z.; Qiu, M.; Fan, T.; Luo, C. S.; Zhou, Y.; Li, Y.; Guo, Z.; Fan, D.; Zhang, H. Fluorinated Phosphorene: Electrochemical Synthesis, Atomistic Fluorination, and Enhanced Stability. Small 2017, 13, 1702739. (13) Zhu, X.; Zhang, T.; Jiang, D.; Duan, H.; Sun, Z.; Zhang, M.; Jin, H.; Guan, R.; Liu, Y.; Chen, M.; Ji, H.; Du, P.; Yan, W.; Wei, S.; Lu, Y.; Yang, S. Stabilizing Black Phosphorus Nanosheets via Edge-Selective Bonding of Sacrificial C60 Molecules. Nat. Commun. 2018, 9, 4177. (14) Sofer, Z.; Luxa, J.; Bouša, D.; Sedmidubský, D.; Lazar, P.; Hartman, T.; Hardtdegen, H.; Pumera, M. The Covalent Functionalization of Layered Black Phosphorus by Nucleophilic Reagents. Angew. Chem., Int. Ed. 2017, 56, 9891−9896. (15) Hu, H.; Gao, H.; Gao, L.; Li, F.; Xu, N.; Long, X.; Hu, Y.; Jin, J.; Ma, J. Covalent Functionalization of Black Phosphorus Nanoflakes by Carbon Free Radicals for Durable Air and Water Stability. Nanoscale 2018, 10, 5834−5839. (16) Cao, Y.; Tian, X.; Gu, J.; Liu, B.; Zhang, B.; Song, S.; Fan, F.; Chen, Y. Covalent Functionalization of Black Phosphorus with Conjugated Polymer for Information Storage. Angew. Chem., Int. Ed. 2018, 57, 4543−4548. 22645

DOI: 10.1021/acsami.9b04344 ACS Appl. Mater. Interfaces 2019, 11, 22637−22647

Research Article

ACS Applied Materials & Interfaces (37) Dovesi, R.; Saunders, V. R.; Roetti, C.; Orlando, R.; ZicovichWilson, C. M.; Pascale, F.; Civalleri, B.; Doll, K.; Harrison, N. M.; Bush, I. J.; D’Arco, P.; Llunell, M.; Causà, M.; Noël, Y.; Maschio, L.; Erba, A.; Rerat, M.; Casassa, S. CRYSTAL17 User’s Manual; University of Torino: Torino, 2017. (38) Peintinger, M. F.; Oliveira, D. V.; Bredow, T. Consistent Gaussian Basis Sets of Triple-Zeta Valence with Polarization Quality for Solid-State Calculations. J. Comput. Chem. 2013, 34, 451−459. (39) Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Scalmani, G.; Barone, V.; Petersson, G. A.; Nakatsuji, H.; Li, X.; Caricato, M.; Marenich, A. V.; Bloino, J.; Janesko, B. G.; Gomperts, R.; Mennucci, B.; Hratchian, H. P.; Ortiz, J. V.; Izmaylov, A. F.; Sonnenberg, J. L.; Williams-Young, D.; Ding, F.; Lipparini, F.; Egidi, F.; Goings, J.; Peng, B.; Petrone, A.; Henderson, T.; Ranasinghe, D.; Zakrzewski, V. G.; Gao, J.; Rega, N.; Zheng, G.; Liang, W.; Hada, M.; Ehara, M.; Toyota, K.; Fukuda, R.; Hasegawa, J.; Ishida, M.; Nakajima, T.; Honda, Y.; Kitao, O.; Nakai, H.; Vreven, T.; Throssell, K.; Montgomery, J. A., Jr.; Peralta, J. E.; Ogliaro, F.; Bearpark, M. J.; Heyd, J. J.; Brothers, E. N.; Kudin, K. N.; Staroverov, V. N.; Keith, T. A.; Kobayashi, R.; Normand, J.; Raghavachari, K.; Rendell, A. P.; Burant, J. C.; Iyengar, S. S.; Tomasi, J.; Cossi, M.; Millam, J. M.; Klene, M.; Adamo, C.; Cammi, R.; Ochterski, J. W.; Martin, R. L.; Morokuma, K.; Farkas, O.; Foresman, J. B.; Fox, D. J. Gaussian 16, Revision B.01; Gaussian, Inc.: Wallingford CT, 2016. (40) Ienco, A.; Manca, G.; Peruzzini, M.; Mealli, C. Modelling Strategies for the Covalent Functionalization of 2D Phosphorene. Dalton Trans. 2018, 47, 17243−17256. (41) Mealli, C.; Ienco, A.; Peruzzini, M.; Manca, G. The Atomic Level Mechanism of White Phosphorous Demolition by Di-Iodine. Dalton Trans. 2018, 47, 394−408. (42) Wang, Y.; Yang, B.; Wan, B.; Xi, X.; Zeng, Z.; Liu, E.; Wu, G.; Liu, Z.; Wang, W. Degradation of Black Phosphorus: A Real-Time 31 P NMR Study. 2D Mater 2016, 3, 035025. (43) Márkus, B. G.; Simon, F.; Nagy, K.; Fehér, T.; Wild, S.; Abellán, G.; Chacón-Torres, J. C.; Hirsch, A.; Hauke, F. Electronic and Magnetic Properties of Black Phosphorus. Phys. Status Solidi B 2017, 254, 1700232. (44) Qiu, P.; Xu, C.; Zhou, N.; Chen, H.; Jiang, F. Metal-free black phosphorus nanosheets-decorated graphitic carbon nitride nanosheets with C P bonds for excellent photocatalytic nitrogen fixation. Appl. Catal., B 2018, 221, 27−35. (45) Zhu, M.; Kim, S.; Mao, L.; Fujitsuka, M.; Zhang, J.; Wang, X.; Majima, T. Metal-Free Photocatalyst for H2 Evolution in Visible to Near-Infrared Region: Black Phosphorus/Graphitic Carbon Nitride. J. Am. Chem. Soc. 2017, 139, 13234−13242. (46) Plutnar, J.; Sofer, Z.; Pumera, M. Products of Degradation of Black Phosphorus in Protic Solvents. ACS Nano 2018, 12, 8390−8396. (47) Weiss, J. W. E.; Bryce, D. L. A Solid-State 11B NMR and Computational Study of Boron Electric Field Gradient and Chemical Shift Tensors in Boronic Acids and Boronic Esters. J. Phys. Chem. A 2010, 114, 5119−5131. (48) Wang, T.; Zhang, N.; Zhang, K.; Dai, J.; Bai, W.; Bai, R. Pyrene Boronic Acid Cyclic Ester: A New Fast Self-Recovering Mechanoluminescent Material at Room Temperature. Chem. Commun. 2016, 52, 9679−9682. (49) Feng, X.; Qi, C.; Feng, H.-T.; Zhao, Z.; Sung, H. H. Y.; Williams, I. D.; Kwok, R. T. K.; Lam, J. W. Y.; Qin, A.; Tang, B. Z. Dual Fluorescence of Tetraphenylethylene-Substituted Pyrenes with Aggregation-Induced Emission Characteristics for White-Light Emission. Chem. Sci. 2018, 9, 5679−5687. (50) Zhang, H.; Ren, T.; Ji, Y.; Han, L.; Wu, Y.; Song, H.; Bai, L.; Ba, X. Selective Modification of Halloysite Nanotubes with 1-Pyrenylboronic Acid: A Novel Fluorescence Probe with Highly Selective and Sensitive Response to Hyperoxide. ACS Appl. Mater. Interfaces 2015, 7, 23805− 23811. (51) Feng, X.; Hu, J.-Y.; Redshaw, C.; Yamato, T. Functionalization of Pyrene To Prepare Luminescent Materials-Typical Examples of Synthetic Methodology. Chem.Eur. J. 2016, 22, 11898−11916.

(52) Lerner, M. B.; Resczenski, J. M.; Amin, A.; Johnson, R. R.; Goldsmith, J. I.; Johnson, A. T. C. Toward Quantifying the Electrostatic Transduction Mechanism in Carbon Nanotube Molecular Sensors. J. Am. Chem. Soc. 2012, 134, 14318−14321. (53) Porcu, P.; Vonlanthen, M.; Ruiu, A.; González-Méndez, I.; Rivera, E. Energy Transfer in Dendritic Systems Having Pyrene Peripheral Groups as Donors and Different Acceptor Groups. Polymers 2018, 10, 1062. (54) Wang, Q.; Li, Z.; Tao, D.-D.; Zhang, Q.; Zhang, P.; Guo, D.-P.; Jiang, Y.-B. Supramolecular Aggregates as Sensory Ensembles. Chem. Commun. 2016, 52, 12929−12939. (55) Zuo, M.; Qian, W.; Xu, Z.; Shao, W.; Hu, X.-Y.; Zhang, D.; Jiang, J.; Sun, X.; Wang, L. Multiresponsive Supramolecular Theranostic Nanoplatform Based on Pillar[5]arene and Diphenylboronic Acid Derivatives for Integrated Glucose Sensing and Insulin Delivery. Small 2018, 14, 1801942. (56) Micali, N.; Mallamace, F.; Castriciano, M.; Romeo, A.; Monsú Scolaro, L. Separation of Scattering and Absorption Contributions in UV/Visible Spectra of Resonant Systems. Anal. Chem. 2001, 73, 4958− 4963. (57) Somerharju, P. Pyrene-Labeled Lipids as Tools in Membrane Biophysics and Cell Biology. Chem. Phys. Lipids 2002, 116, 57−74. (58) Castriciano, M. A.; Cardiano, P.; Fazio, E.; Mineo, P. G.; Nicosia, A.; Zagami, R.; Trapani, M.; Monsù Scolaro, L.; Lo Schiavo, S. Novel Luminescent Ionic Adducts Based on Pyrene-1-Sulfonate. ACS Omega 2018, 3, 18811−18820. (59) Zöphel, L.; Enkelmann, V.; Müllen, K. Tuning the HOMOLUMO Gap of Pyrene Effectively via Donor-Acceptor Substitution: Positions 4,5 versus 9,10. Org. Lett. 2013, 15, 804−807. (60) Hu, W.; Yang, J. Two-Dimensional van Der Waals Heterojunctions for Functional Materials and Devices. J. Mater. Chem. C 2017, 5, 12289−12297. (61) Amelia, M.; Lavie-Cambot, A.; McClenaghan, N. D.; Credi, A. A ratiometric luminescent oxygen sensor based on a chemically functionalized quantum dot. Chem. Commun. 2011, 47, 325−327. (62) Cho, S.-Y.; Lee, Y.; Koh, H.-J.; Jung, H.; Kim, J.-S.; Yoo, H.-W.; Kim, J.; Jung, H.-T. Superior Chemical Sensing Performance of Black Phosphorus: Comparison with MoS2 and Graphene. Adv. Mater. 2016, 28, 7020−7028. (63) Cho, S.-Y.; Koh, H.-J.; Yoo, H.-W.; Jung, H.-T. Tunable Chemical Sensing Performance of Black Phosphorus by Controlled Functionalization with Noble Metals. Chem. Mater. 2017, 29, 7197− 7205. (64) Jain, R.; Singh, Y.; Cho, S.-Y.; Sasikala, S. P.; Koo, S. H.; Narayan, R.; Jung, H.-T.; Jung, Y.; Kim, S. O. Ambient Stabilization of Few Layer Phosphorene via Noncovalent Functionalization with Surfactants: Systematic 2D NMR Characterization in Aqueous Dispersion. Chem. Mater. 2019, 31, 2786−2794. (65) Sugai, S.; Shirotani, I. Raman and Infrared Reflection Spectroscopy in Black Phosphorus. Solid State Commun. 1984, 53, 753−755. (66) Abellán, G.; Wild, S.; Lloret, V.; Scheuschner, N.; Gillen, R.; Mundloch, U.; Maultzsch, J.; Varela, M.; Hauke, F.; Hirsch, A. Fundamental Insights into the Degradation and Stabilization of Thin Layer Black Phosphorus. J. Am. Chem. Soc. 2017, 139, 10432−10440. (67) Zhang, T.; Wan, Y.; Xie, H.; Mu, Y.; Du, P.; Wang, D.; Wu, X.; Ji, H.; Wan, L. Degradation Chemistry and Stabilization of Exfoliated Few-Layer Black Phosphorus in Water. J. Am. Chem. Soc. 2018, 140, 7561−7567. (68) Han, X.; Han, J.; Liu, C.; Sun, J. Promise and Challenge of Phosphorus in Science, Technology, and Application. Adv. Funct. Mater. 2018, 28, 1803471. (69) Ponti, G.; Palombi, F.; Abate, D.; Ambrosino, F.; Aprea, G.; Bastianelli, T.; Beone, F.; Bertini, R.; Bracco, G.; Caporicci, M.; Calosso, B.; Chinnici, M.; Colavincenzo, A.; Cucurullo, A.; Dangelo, P.; De Rosa, M.; De Michele, P.; Funel, A.; Furini, G.; Giammattei, D.; Giusepponi, S.; Guadagni, R.; Guarnieri, G.; Italiano, A.; Magagnino, S.; Mariano, A.; Mencuccini, G.; Mercuri, C.; Migliori, S.; Ornelli, P.; Pecoraro, S.; Perozziello, A.; Pierattini, S.; Podda, S.; Poggi, F.; 22646

DOI: 10.1021/acsami.9b04344 ACS Appl. Mater. Interfaces 2019, 11, 22637−22647

Research Article

ACS Applied Materials & Interfaces Quintiliani, A.; Rocchi, A.; Sciò, C.; Simoni, F.; Vita, A. Proceedings of the 2014 International Conference on High Performance Computing and Simulation in The Role of Medium Size Facilities in the HPC Ecosystem: The Case of the New CRESCO4 Cluster Integrated in the ENEAGRID Infrastructure, 2014; pp 1030−1033.

22647

DOI: 10.1021/acsami.9b04344 ACS Appl. Mater. Interfaces 2019, 11, 22637−22647