Ultrafast nonlinear optical response in solution dispersions of black phosphorus

Ultrafast nonlinear optical response in solution dispersions of black phosphorus


Play all audios:


ABSTRACT We report the spatial self-phase modulation (SSPM) effect for solution dispersions of black phosphorus (BP). The experimental results suggest that this concentration-dependent


coherent light diffraction is due to the ultrafast and large third-order optical nonlinearity of BP. The third-order nonlinear susceptibility of BP has been simply obtained about 10−19 m2/V2


by analyzing the experimental results. The fast relaxation time during dynamic relaxation is obtained as 0.13 ps. Our experimental results imply novel potential application of BP in


ultrafast nonlinear phase modulation devices based on their nonlinear optical response. SIMILAR CONTENT BEING VIEWED BY OTHERS ULTRAFAST ENERGIZING THE PARITY-FORBIDDEN DARK EXCITON IN BLACK


PHOSPHORUS Article Open access 29 April 2025 NEGATIVE-MASS EXCITON POLARITONS INDUCED BY DISSIPATIVE LIGHT-MATTER COUPLING IN AN ATOMICALLY THIN SEMICONDUCTOR Article Open access 23


February 2023 FULL PHONON DISPERSION ALONG THE STACKING DIRECTION IN NANOSCALE VAN DER WAALS MATERIALS BY PICOSECOND ACOUSTICS Article Open access 13 June 2024 INTRODUCTION Low dimensional


optical materials with ultrafast response and strong nonlinearity have exhibited novel physics and tremendous application potential in optical communication, optical sensing, medical


diagnosis, etc1,2,3,4,5,6. Inspired by the emergence of two-dimensional atomic crystal graphene, the nonlinear optics of graphene and graphene-like materials have been explored intensively


in recent years endowed for their excellent nonlinear response and broadband applications7,8,9,10. Besides graphene, black phosphorus (BP), a newly rising graphene-like material with natural


band gap, has shown tremendous and intriguing physical, chemical, and electronic properties11,12,13. The BP has a thickness dependent direct band gap that changes from 0.3 eV in the bulk to


1.88 eV in a monolayer14. It can also be switched between insulating and conducting states, and it is still flat enough to confine electrons so that charge flows quickly, comparable with


those in single layer of molybdenum disulphide15. Sample mobility is also found to be thickness dependent, with the highest value up to ~1000 cm2/Vs obtained at thickness ~10 nm16.


Furthermore, the BP photodetector shows time response of 1 ms (rise) and 4 ms (fall), making few-layer BP a promising active 2D material in broadband and fast photodetectors across the


visible and near infrared region17. In addition, BP is made from a single element, in theory, the pure samples are easier to obtain. Compared to other graphene-like 2D materials, tunable


optical properties and direct bandgap for all thicknesses of BP may bring significant benefits to a variety of fascinating photonics applications. Besides its high electron mobility, BP has


unique optical properties, such as the strong in-plane anisotropic18. Considering the tunable band-gap, the electronic and optical performance of the materials can change dramatically


consequently. Beyond the linear optics regime, the optical properties, especially the ultrafast optical response and the nonlinear optical properties of BP in case of different thicknesses,


various doping conditions, and different polarization directions of the excitation light, have been paid more attention recently. The ultrafast relaxation dynamics of BP have been


investigated by the pump-probe technique19,20,21, and the fast recovery time in BP was determined to be ~24 fs at 1550 nm20, 16 and 32 fs for 800 and 2026 nm pulses21, which is much faster


than that of previous 2D crystal materials, such as graphene (~1.27 ps)22 and MoS2 (~2.1 ps)23. R. Suess _et al_. reported the anisotropic carrier dynamic with faster time to be ~180 ps at


780 nm in BP, which may be induced by hot carriers cooling19. In another aspect, the third-order nonlinear optical responses of BP have been investigated and characterized via Z-scan


technique, third-harmonic generation (THG) method, etc24,25,26. The open-aperture (OA) Z-scan experiments demonstrate that BP has perfect saturable absorption performance at both optical


communication band (1550 nm) and visible band (532 nm and 680 nm), making BP a preferred alternative material for saturable absorbers (SAs)20. The broadband nonlinear optics response of


few-layer BP towards mid-infrared has also been reported by S. Lu _et al_.24. In addition, X. Zheng _et al_. reported a transition from saturable absorption (SA) to reverse saturable


absorption (RSA) with the increase of laser intensity26. Moreover, both the passive Q-switching and the mode-locking operation of the erbium-doped fiber lasers with the BP-SA have been


experimentally demonstrated27. D. Li _et al_. studied the thickness and polarization dependent linear and nonlinear optical properties of BP thin films, and then utilized their nonlinear


absorption property to generate ultrafast and large-energy pulse with BP integrated fiber devices28. Besides Z-scan and THG methods, spatial self-phase modulation (SSPM), a phenomenon that


the intensity dependence of the refractive index in nonlinear optical media occurs, is a manifestation of the coherent response of the material. SSPM was originally observed in


liquid-crystal29, 30, then later in nanomaterials, such as carbon nanotubes (CNTs) solutions31, graphene32, graphene oxide33, topological insulators34, transition metal dichalcogenides


(TMDs)35, 36 and many other materials. While this manuscript was in preparation, a relevant work was published37. It reports the broadband spatial self-phase modulation of BP, and the


nonlinear refractive index of BP (~5 layers) is measured to be ~10−9 m2/W and the third order nonlinear susceptibility is χ(3) ~ 10−8 esu at multiple wavelengths (equals χ(3) to be 10−16 


m2/V2)37. Youngblood _et al_. report the χ(3) of multilayer (~29 layers) BP to be 10−19 m2/V2 via layer tunable THG25, while Zheng _et al_. reported a nonlinear refractive index value of


~10−13 m2/W (~50–100 layers) at 800 nm femtosecond pulsed Z-Scan measurement26. In addition, from other perspective, these methods have given the different nonlinear optical parameters of


BP25, 26, 37. It is exactly these disagreements that create the need for additional experimental and theoretical investigations to fully characterize the nonlinear susceptibility of BP


further. In this work, we studied the phenomenon of SSPM for solution dispersions of BP under different concentration and laser intensity. The ultrafast nonlinear optical response was


unambiguously observed and this coherent light diffraction is due to the remarkably large third-order nonlinearity. The third-order nonlinear refractive index of BP was then simply and


exactly obtained by analyzing the experimental results. RESULTS CHARACTERIZATIONS OF BLACK PHOSPHORUS Uniform few-layer BP dispersions have been prepared by the liquid phase exfoliation


(LPE) method, which is a simple and effective technique to prepare two dimensional (2D) materials from layered bulk crystals towards multiple layered structures. Figure 1 shows the


characterization of the BP sample. The scanning electron microscopy (SEM) image of BP crystal is shown in Fig. 1a,b, from which the high quality and flat layered surface of BPs can be


verified from both lateral (a) and top (b) view, respectively. As shown in Fig. 1c, there are three Raman peaks at 361.15 cm−1, 436.85 cm−1, and 464.68 cm−1 from the multi-layer phosphorus


corresponding to the \({{A}}_{{g}}^{1}\), \({B}_{2g}\), and \(\,{A}_{g}^{2\,}\) modes, respectively. Figure 1d shows the linear absorption spectrum of the BP solution dispersions in ethylene


glycol ranging from 300 to 1200 nm, measured by the spectrometer (Perkins Elmer Lambda750). The profile is nearly flat with a smooth absorption curve in the ultraviolet up to near-infrared


(UV-NIR) wavelength band, suggesting that multi-layer BP dispersions might be a promising broadband optical material. ULTRA-FAST SSPM RESPONSE Inspired by the unique electronic and optical


properties of BP, we have experimentally observed the SSPM phenomenon of BP. In the experiment, the ultrafast pulse laser was used to characterize the SSPM effects of BP solution


dispersions. The linearly polarized beam from a 1061 nm (repetition rate of 2.7 MHz and pulse duration of 220 ps) laser was firstly focused by a lens of focal length 150 mm, causing high


intensity illumination to the sample. A quartz cuvette of 5 mm thickness was used to contain the BP solution dispersions and placed before the beam focus. The distance between the focus and


the front surface of the cuvette was 20 mm. Then the transmitted spectrum was collected by a CCD camera (Coherent LaserCam HR), which is 200 mm away from the focal lens. The schematic


diagram was shown in Fig. 2. Owing to the nonlinear index change induced by the nonuniform Gaussian profile of the light beam, we observed obvious diffraction ring pattern in the far field.


Figure 3a is a typical pattern of diffraction rings caused by spatial self-phase modulation when the sample was irradiated by 1061 nm ultrafast laser. Figure 3b is the corresponding


intensity distribution of the experimental result, and Fig. 3c illustrates the corresponding far field intensity distribution simulated by means of the Fraunhofer approximation of the


Fresnel–Kirchhoff diffraction formula. From the figure, we conclude that the experimental result is in a good agreement with the theoretical one. Measurements of the pure ethylene glycol did


not show a nonlinear response for the intensity ranges that were investigated here and, therefore, confirm that the observed nonlinear behavior is from the presence of the BP (_see


supporting information_). To confirm the presence of the BP induces the nonlinear refractive index in the medium and causes the SSPM effect in laser beam propagation, we performed


concentration-dependence measurements. Figure 4a,b show the typical development of diffraction rings of SSPM when the input average power is 100 mW. They both emerged from the center of the


patterns with all their diameters enlarging gradually and they deformed with irradiation time evolution and finally become stable during one second. As shown in Fig. 4a, it takes about 0.95 


second for the development of diffraction rings, which includes emerging, enlarging and deforming. While varying the concentration of the BP solution dispersion, the number of diffraction


pattern increases, and the nonlinear response will become much faster with higher BP concentration. As the SSPM phenomena are proportional to the absorption of the nonlinear medium, pattern


formation time, the size and number of rings become larger with the increasing of the concentration. Figure 4c shows the final stable patterns for four different concentrations of solution


dispersions under the same input average power 140 mW, where the concentration decreases gradually from 15 vol.% to 7.5 vol.%, respectively. We can observe that both numbers and diameters of


the diffraction rings decrease with the concentration. These results verify that diffraction rings are mainly derived from the nonlinear optical properties of BP. We also performed


intensity-dependence measurements with the same experimental setup to convince the nonlinear response of the BP dispersion. The input laser power increases gradually from Fig. 5a–f as 60,


80, 100, 110, 130 and 140 mW, respectively. It can be observed that by increasing the laser power, the number of rings and the vertical asymmetry of the diffraction ring patterns increased.


Consequently, one can evidently confirm that when the input laser power becomes higher, the propagating beam suffers more serious distortions. The refractive index of many materials can be


described by the relation $$n={n}_{0}+{n}_{2}I$$ (1) where _n_ 0 is the linear refractive index and _I_ is the laser intensity and _n_ 2 is the coefficient of the intensity-dependent


refractive index38. The total refractive index increases with increasing optical intensity. A corresponding phase shift Δ_ψ_ of the beam traversing the nonlinear medium32: $${\rm{\Delta


}}\phi (r)=\frac{2\pi {n}_{0}}{\lambda }{\int }_{0}^{Leff}{n}_{2}I(r,z)dz$$ (2) where _r_ is the radial coordinate of the beam, _λ_ is the wavelength of incident laser, _L_ eff is the total


propagation length contributing to the SSPM. We assume, for simplicity, that $${\rm{\Delta }}\phi (r)={\rm{\Delta }}\phi (0)\exp (\frac{-2{r}^{2}}{{\omega }^{2}})$$ (3) where _ω_ is the beam


radius, as a constant. When \(r\in [0,+\infty )\), there is _r_ 1 and _r_ 2 possessing the same wave vector and can interfere, which can be expressed as29: $$\frac{d{\rm{\Delta }}\phi


({r}_{1})}{dr}=\frac{d{\rm{\Delta }}\phi ({r}_{2})}{dr}$$ (4) The constructive or destructive interference occurs and results appearance of diffraction rings. $${\rm{\Delta }}\phi


({r}_{1})-{\rm{\Delta }}\phi ({r}_{2})=m\pi $$ (5) where _m_ is an even or odd integer. The number of diffraction rings appeared can be estimated as29: $$N=\frac{{\rm{\Delta }}\phi


(0)-{\rm{\Delta }}\phi (r)}{2\pi }$$ (6) Assuming Δ_ψ_(_r_) = 0, thus, \(N={\rm{\Delta }}{\phi }_{0}/2\pi \), we can obtain the equation to solve _n_ 2 32: $${n}_{2}=(\frac{\lambda


}{2{n}_{0}{L}_{eff}})\frac{N}{I}$$ (7) The nonlinear refractive index and the third-order nonlinear susceptibility of BP at different concentrations have been summarized in Table 1. We


compared the rings number, and nonlinear refractive index _n_ 2 at different concentration and input power which corresponds to different input intensity. At lower concentration, _n_ 2


increases with input power or intensity, while the saturation of _n_ 2 occurs due to sample concentration increase, as shown in Fig. 6. In our experiment, _N_ is approximately proportional


to _I_ and the total refractive indexes increase with increasing optical intensity. The threshold for observing diffraction rings is about 0.58 × 105 W/cm2 for 15 vol.% and 0.97 × 105 W/cm2


for 7.5 vol.% BP solution dispersions. The nonlinearity _n_ 2 can be easily obtained by finding _N_/_I_ at the maximum laser power. We can estimate _n_ 2 as 4.35 × 10−12 m2/W for 15 vol.%


and 1.45 × 10−12 m2/W for 7.5 vol.% BP solution dispersions. The third order nonlinear susceptibility of BP was calculated to be 10−14 m2/ V2 (_see supporting information_). We use SI units


here to avoid confusion. We estimate that \({\chi }_{total}^{(3)}={M}^{2}{\chi }_{{single}\,{layer}}^{(3)}\) 32. The number of effective layers _M_ can be obtained as \(M=\frac{C\times


\,V\times \,{N}_{A}}{S/{a}^{2}}\), where _V_ is the volume of the solution of 1.5 mL. The total number of (BP) molecules in the solution is _C_ × _V_ × _N_ A, where _N_ A is the Avogadro’s


number, _S_ is the area of the cross-section of the cuvette ( = 1 cm × 3 cm), the total number of molecules in one effective layer is _S_/a2, where a (~4.5 Å) is the lattice constant of


BP39. In our experiment, _M_ is on the order of 300 to 600, thus the estimated \({\chi }_{single\,layer}^{(3)}\) is of the order of 10−19 m2/V2, which agrees with Z-scan measurements


(~50–100 layers)26, two order of magnitude larger than estimated from THG experiments (~29 layers)25, but more than two orders of magnitude smaller than that from ref. 37 (~5 layers). The


variation of nonlinear susceptibility by different methods can be attributed to the difference of BP thickness, as BP possesses a thickness-dependent energy band gap. Further efforts, from


both theoretical and experimental sides, are needed to fully understand the third-order nonlinear optical response of BP, e.g. the relationship between formation time and BP concentration.


The measurement uncertainty comes from several aspects: the calculation of incident pulse laser intensity, the measurement of distance between lens, sample, CCD camera, the concentration and


uniformity uncertainty of the BP sample during the formation of diffraction ring, and the counting of diffraction rings number, which is a key parameter in calculating the nonlinear


coefficient for BP. However, by well-designed experimental scheme and data analysis, SSPM phenomena can be regarded as an effective method to estimate the orders of magnitude of the


nonlinear refractive index. CARRIER DYNAMICS The ultrafast nonlinear response of BP was also investigated by time-resolved optical degenerate pump-probe transmission measurement. Femtosecond


(fs) laser (Coherent Libra-S), which gave an output with a pulse repetition rate of 1 kHz, a central wavelength of 800 nm and a pulse duration around 100 fs, was employed in the


time-resolved pump-probe experiment. The output laser beam was split into two different laser beams, one intense portion with a peak intensity of 75 GW/cm2 as the pump beam to generate the


photo-induced carriers, and the other weak portion with peak intensity of 0.3 GW/cm2 as the probe beam. Figure 7 presents the carriers dynamic relaxation process for the BP suspension. The


delay kinetics can be well fitted with a bi-exponential delay function: $$\frac{{\rm{\Delta }}T}{T}={A}_{1}{e}^{-t/{\tau }_{1}}+{A}_{2}{e}^{-t/{\tau }_{2}}$$ (8) where Δ_T_/_T_ is defined as


the relative change of the probe transmissivity caused by the pump, and τ1 is the decay time with the respective amplitude weights _A_ 1. The fast and slow relaxation time is fitted to be


0.13 ps and 1.15 ps, respectively. We attribute the SSPM effect a transient phenomenon. The interaction and diffraction can occur under the well-defined phase difference between photons. To


the ultrafast laser pulse, the carriers will relax in only about several pico-seconds, while every pulse is in the temporal slot of 0.37 μs (1061 nm, 2.7 MHz). Therefore, the incident laser


pulse is considered to be independent to each other, which confirms the transient effect of SSPM. The physical mechanism of the nonlinearity may be attributed to the reorientation and


alignment of the BP nanosheets induced by the electromagnetic field, which is similar to the case of liquid crystals29. Under the intense laser irradiation, the polarization state will


relocate to a new array which is parallel to each other. The reoriented sheets can help light beam become coherent and form SSPM rings. CONCLUSIONS In conclusion, the third-order nonlinear


susceptibility of BP dispersions was measured to be 10−19 m2/V2 and can be tuned via changing its concentration. Upon illumination, the BP solution dispersions show obvious SSPM response to


the excitation wavelengths, demonstrating its ultrafast nonlinear response. Experimental results of BP under pulse illumination are fitted with a good correspondence between experimental and


numerical results, which suggest that this coherent light scattering is due to the ultrafast, and large third-order optical nonlinearity of BP. The fast relaxation time during dynamic


relaxation is obtained as 0.13 ps. Our experimental results exhibit the practical potential of this promising material for various nonlinear and ultrafast optoelectronics applications (e.g.,


ultrafast lasers, optical switches and modulators). Moreover, our work may provide an inroad for measuring key parameters, such as molecular weights of the solvent and the nanomaterials


based on their novel nonlinear optical response. The technique can further be used to study the composition, purity and stability of other nanomaterials at broadband wavelength, including


those with potential for clinical applications. METHODS MATERIALS To begin, bulk BP crystal (purchased from Smart Elements) was grinded to BP powder. Then the powder was dispersed into


ethylene glycol, along with ultra-sonicated for 2 hours. After sonication, the phosphorene in ethylene glycol was centrifuged at a speed of 15000 rpm for 20 minutes. The upper ~80% which


contains few layers was decanted for analysis and the dispersions were settled for more than 20 hours to make sure that there is no large size sedimentations during our experiments. BP


dispersions were then diluted to 4 samples (15 vol.%, 12.5 vol.%, 10 vol.%, 7.5 vol.%, respectively) for SSPM experiments. REFERENCES * Tao, W. _et al_. Black Phosphorus Nanosheets as a


Robust Delivery Platform for Cancer Theranostics. _Adv. Mater._ 29, 1603276 (2017). Article  Google Scholar  * Zhang, Y.-H. _et al_. Improving gas sensing properties of graphene by


introducing dopants and defects: a first-principles study. _Nanotechnoloy_ 20, 185504 (2009). Article  ADS  Google Scholar  * Mueller, T., Xia, F. N. A. & Avouris, P. Graphene


photodetectors for high-speed optical communications. _Nat. Photonics_ 4, 297–301 (2010). Article  CAS  Google Scholar  * Song, Y. J., Qu, K. G., Zhao, C., Ren, J. S. & Qu, X. G.


Graphene Oxide: Intrinsic Peroxidase Catalytic Activity and Its Application to Glucose Detection. _Adv. Mater._ 22, 2206–2210 (2010). Article  CAS  PubMed  Google Scholar  * Geim, A. K.


& Novoselov, K. S. The rise of graphene. _Nat. Mater._ 6, 183–191 (2007). Article  ADS  CAS  PubMed  Google Scholar  * Wang, Q. H., Kalantar-Zadeh, K., Kis, A., Coleman, J. N. &


Strano, M. S. Electronics and optoelectronics of two-dimensional transition metal dichalcogenides. _Nat. Nanotechnol._ 7, 699–712 (2012). Article  ADS  CAS  PubMed  Google Scholar  * Lu, S.


B. _et al_. Third order nonlinear optical property of Bi2Se3. _Opt. Express_ 21, 2072–2082 (2013). Article  ADS  CAS  PubMed  Google Scholar  * Zhang, H. _et al_. Z-scan measurement of the


nonlinear refractive index of graphene. _Opt. Lett._ 37, 1856–1858 (2012). Article  ADS  CAS  PubMed  Google Scholar  * Woodward, R. I. _et al_. Wideband saturable absorption in few-layer


molybdenum diselenide (MoSe2) for Q-switching Yb-, Er- and Tm-doped fiber lasers. _Opt. Express_ 23, 20051–20061 (2015). Article  ADS  CAS  PubMed  Google Scholar  * Zhang, M. _et al_.


Solution processed MoS2-PVA composite for sub-bandgap mode-locking of a wideband tunable ultrafast Er:fiber laser. _Nano Res._ 8, 1522–1534 (2015). Article  CAS  Google Scholar  * Xia, F.


N., Wang, H. & Jia, Y. C. Rediscovering black phosphorus as an anisotropic layered material for optoelectronics and electronics. _Nat. Commun._ 5, 4458 (2014). ADS  CAS  PubMed  Google


Scholar  * Liu, H., Du, Y. C., Deng, Y. X. & Ye, P. D. Semiconducting black phosphorus: synthesis, transport properties and electronic applications. _Chem. Soc. Rev._ 44, 2732–2743


(2015). Article  CAS  PubMed  Google Scholar  * Wang, X. M. & Lan, S. F. Optical properties of black phosphorus. _Adv. Opt. Photonics_ 8, 618–655 (2016). Article  Google Scholar  * Guo,


Z. _et al_. From Black Phosphorus to Phosphorene: Basic Solvent Exfoliation, Evolution of Raman Scattering, and Applications to Ultrafast Photonics. _Adv. Funct. Mater._ 25, 6996–7002


(2015). Article  CAS  Google Scholar  * Liu, H. _et al_. Phosphorene: An Unexplored 2D Semiconductor with a High Hole Mobility. _ACS Nano_ 8, 4033–4041 (2014). Article  CAS  PubMed  Google


Scholar  * Li, L. _et al_. Black phosphorus field-effect transistors. _Nat. Nanotechnol._ 9, 372–377 (2014). Article  ADS  CAS  PubMed  Google Scholar  * Buscema, M. _et al_. Fast and


Broadband Photoresponse of Few-Layer Black Phosphorus Field-Effect Transistors. _Nano Lett._ 14, 3347–3352 (2014). Article  ADS  CAS  PubMed  Google Scholar  * Wang, X. M. _et al_. Highly


anisotropic and robust excitons in monolayer black phosphorus. _Nat. Nanotechnol._ 10, 517–521 (2015). Article  ADS  CAS  PubMed  Google Scholar  * Suess, R. J., Jadidi, M. M., Murphy, T. E.


& Mittendorff, M. Carrier dynamics and transient photobleaching in thin layers of black phosphorus. _Appl. Phys. Lett._ 107, 081103 (2015). Article  ADS  Google Scholar  * Wang, Y. _et


al_. Ultrafast recovery time and broadband saturable absorption properties of black phosphorus suspension. _Appl. Phys. Lett._ 107, 091905 (2015). Article  ADS  Google Scholar  * Wang, K.


_et al_. Ultrafast Nonlinear Excitation Dynamics of Black Phosphorus Nanosheets from Visible to Mid-Infrared. _ACS Nano_ 10, 6923–6932 (2016). Article  CAS  PubMed  Google Scholar  *


Nikolaenko, A. E. _et al_. Nonlinear graphene metamaterial. _Appl. Phys. Lett._ 100, 181109 (2012). Article  ADS  Google Scholar  * Wang, Q. _et al_. Valley Carrier Dynamics in Monolayer


Molybdenum Disulfide from Helicity-Resolved Ultrafast Pump–Probe Spectroscopy. _ACS Nano_ 7, 11087 (2013). Article  CAS  PubMed  Google Scholar  * Lu, S. B. _et al_. Broadband nonlinear


optical response in multi-layer black phosphorus: an emerging infrared and mid-infrared optical material. _Opt. Express_ 23, 11183–11194 (2015). Article  ADS  CAS  PubMed  Google Scholar  *


Youngblood, N., Peng, R., Nemilentsau, A., Low, T. & Li, M. Layer Tunable Third-Harmonic Generation in Multilayer Black Phosphorus. _ACS Photon._ 4, 8–14 (2016). Article  Google Scholar


  * Zheng, X. _et al_. Characterization of nonlinear properties of black phosphorus nanoplatelets with femtosecond pulsed Z-scan measurements. _Opt. Lett._ 40, 3480–3483 (2015). Article  ADS


  CAS  PubMed  Google Scholar  * Chen, Y. _et al_. Mechanically exfoliated black phosphorus as a new saturable absorber for both Q-switching and Mode-locking laser operation. _Opt. Express_


23, 12823–12833 (2015). Article  ADS  CAS  PubMed  Google Scholar  * Li, D. A. _et al_. Polarization and Thickness Dependent Absorption Properties of Black Phosphorus: New Saturable Absorber


for Ultrafast Pulse Generation. _Sci. Rep._ 5, 15899 (2015). Article  ADS  CAS  PubMed  PubMed Central  Google Scholar  * Durbin, S. D., Arakelian, S. M. & Shen, Y. R. Laser-Induced


Diffraction Rings from a Nematic-Liquid-Crystal Film. _Opt. Lett._ 6, 411–413 (1981). Article  ADS  CAS  PubMed  Google Scholar  * Khoo, I. C., Michael, R. R. & Yan, P. Y.


Optically-Induced Molecular-Reorientation in Nematic Liquid-Crystals and Nonlinear Optical Processes in the Nanoseconds Regime. _IEEE J. Quantum Elect._ 23, 267–272 (1987). Article  ADS 


Google Scholar  * Ji, W., Chen, W. Z., Lim, S. H., Lin, J. Y. & Guo, Z. X. Gravitation-dependent, thermally-induced self-diffraction in carbon nanotube solutions. _Opt. Express_ 14,


8958–8966 (2006). Article  ADS  CAS  PubMed  Google Scholar  * Wu, R. _et al_. Purely coherent nonlinear optical response in solution dispersions of graphene sheets. _Nano Lett._ 11,


5159–5164 (2011). Article  ADS  CAS  PubMed  Google Scholar  * Sadrolhosseini, A. R., Rashid, S. A., Shojanazeri, H., Noor, A. S. M. & Nezakati, H. Spatial self-phase modulation patterns


in graphene oxide and graphene oxide with silver and gold nanoparticles. _Opt. Quant. Electron._ 48, 1–13 (2016). Article  CAS  Google Scholar  * Shi, B. _et al_. Broadband ultrafast


spatial self-phase modulation for topological insulator Bi2Te3 dispersions. _Appl. Phys. Lett._ 107, 151101 (2015). Article  ADS  Google Scholar  * Wang, W., Wu, Y., Wu, Q., Hua, J. &


Zhao, J. Coherent Nonlinear Optical Response Spatial Self-Phase Modulation in MoSe2 Nano-Sheets. _Sci. Rep._ 6, 22072 (2016). Article  ADS  CAS  PubMed  PubMed Central  Google Scholar  *


Wang, G. _et al_. Tunable nonlinear refractive index of two-dimensional MoS2, WS2, and MoSe2 nanosheet dispersions [Invited]. _Photon. Res._ 3, A51 (2015). Article  Google Scholar  * Zhang,


J. D. _et al_. Broadband spatial self-phase modulation of black phosphorous. _Opt. Lett._ 41, 1704–1707 (2016). Article  ADS  PubMed  Google Scholar  * Boyd, R. W. Nonlinear Optics. 210


(Academic press, 2008,). * Qiao, J., Kong, X., Hu, Z.-X., Yang, F. & Ji, W. Few-layer black phosphorus: emerging direct band gap semiconductor with high carrier mobility. _Nat. Commun._


5, 4475 (2014). CAS  PubMed  PubMed Central  Google Scholar  Download references ACKNOWLEDGEMENTS This work was supported by National Natural Science Foundation of China (Grant No. 61475102)


and China Scholarship Council. AUTHOR INFORMATION AUTHORS AND AFFILIATIONS * Key Laboratory for Micro-/Nano- Optoelectronic Devices of Ministry of Education, School of Physics and


Electronics, Hunan University, Changsha, 410082, China Lili Miao, Jun Yi, Chujun Zhao & Shuangchun Wen * Femtosecond Optics Group, Department of Physics, Imperial College London, London,


SW7 2BW, UK Lili Miao * Wuhan Optics Valley Aerospace Sanjiang Laser Industrial Technology Research Institute Co. Ltd., Wuhan, 430000, China Bingxin Shi * Troop NO. 95829, Xiaogan, 432100,


China Yaqin Jiang Authors * Lili Miao View author publications You can also search for this author inPubMed Google Scholar * Bingxin Shi View author publications You can also search for this


author inPubMed Google Scholar * Jun Yi View author publications You can also search for this author inPubMed Google Scholar * Yaqin Jiang View author publications You can also search for


this author inPubMed Google Scholar * Chujun Zhao View author publications You can also search for this author inPubMed Google Scholar * Shuangchun Wen View author publications You can also


search for this author inPubMed Google Scholar CONTRIBUTIONS L.L.M., B.X.S., and C.J.Z. designed the research; L.L.M., B.X.S., and Y.Q.J. performed the experiments; J.Y. provided reagents


and technical support; L.L.M., B.X.S., and C.J.Z. analyzed the data and wrote the paper. S.C.W. supervised the project. All authors reviewed the manuscript. CORRESPONDING AUTHOR


Correspondence to Chujun Zhao. ETHICS DECLARATIONS COMPETING INTERESTS The authors declare that they have no competing interests. ADDITIONAL INFORMATION PUBLISHER'S NOTE: Springer


Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. ELECTRONIC SUPPLEMENTARY MATERIAL ULTRAFAST NONLINEAR OPTICAL RESPONSE IN


SOLUTION DISPERSIONS OF BLACK PHOSPHORUS RIGHTS AND PERMISSIONS OPEN ACCESS This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use,


sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative


Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated


otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds


the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/. Reprints and


permissions ABOUT THIS ARTICLE CITE THIS ARTICLE Miao, L., Shi, B., Yi, J. _et al._ Ultrafast nonlinear optical response in solution dispersions of black phosphorus. _Sci Rep_ 7, 3352


(2017). https://doi.org/10.1038/s41598-017-03667-z Download citation * Received: 07 March 2017 * Accepted: 03 May 2017 * Published: 13 June 2017 * DOI:


https://doi.org/10.1038/s41598-017-03667-z SHARE THIS ARTICLE Anyone you share the following link with will be able to read this content: Get shareable link Sorry, a shareable link is not


currently available for this article. Copy to clipboard Provided by the Springer Nature SharedIt content-sharing initiative