Experimental Characterization of the Ultrafast, Tunable and Broadband Optical Kerr Nonlinearity in Graphene

Graphene’s giant nonlinear optical response along with its integrability has made it a vaunted material for on-chip photonics. Despite a multitude of studies confirming its strong nonlinearity, there is a lack of reports examining the fundamental processes that govern the response. Addressing this gap in knowledge we analyse the role of experimental parameters by systematically measuring the near-infrared spectral dependence, the sub-picosecond temporal evolution and pulse-width dependence of the effective Kerr coefficient (n2,eff) of graphene in hundreds of femtosecond regime. The spectral dependence measured using the Z-scan technique is corroborated by a density matrix quantum theory formulation to extract a n2,eff ∝ λ2 dependence. The temporal evolution obtained using the time-resolved Z-scan measurement shows the nonlinearity peaking at zero delay time and relaxing on a time-scale of carrier relaxation. The dependence of the n2,eff on pulse duration is obtained by expanding the input pulse using a prism-pair set-up. Our results provide an avenue for controllable tunability of the nonlinear response in graphene, which is limited in silicon photonics.

The Dirac band structure of graphene 1,2 has endowed this atomically thin two-dimensional material with exceptional optical and transport properties. As a scale-invariant host of chiral carriers 3 , graphene exhibits a universal optical response in absorption, with a single sheet capable of absorbing ~2.3% of normally incident light. Its reduced dimensionality together with the peculiar Dirac type dynamics of the quasiparticles enable graphene to have a relatively strong nonlinear optical response in the presence of intense electromagnetic illuminations 4,5 . Intrigued by its potential applications in photonics and optoelectronics several theoretical and experimental studies of the nonlinear optical properties of graphene have been carried out [3][4][5][6][7][8][9][10][11][12][13][14][15][16][17][18][19][20][21][22] . Due to the centrosymmetric crystalline structure of graphene, even-order nonlinearities are forbidden, and the first nonlinear contribution is a third order effect. In particular, the Kerr type nonlinearity is dominantly attributable to the third order nonlinear interactions in graphene. The refractive index, n, considering the first nonlinear term is intensity dependent and is given by the expression n = n 0 + n 2 I, where n 0 is the linear refractive index, n 2 is the Kerr coefficient, and I is incident intensity. This Kerr-type nonlinearity in graphene has been reported to be large, confirmed by studies characterising frequency mixing 4,6,7 , harmonic generation 8 , self-phase modulation 9,10 , and self-refraction 5,[11][12][13][14][15] .
Several theoretical works have attempted to provide a cohesive theory for the third order nonlinearity in graphene, addressing the magnitude 3,16-18 of the nonlinear refractive index (n 2 ), the spectral 3,4,16 and temporal 19,20 dependence, and effects of Fermi energy modulation 3,17,18,21,22 . Although topological anomalies in graphene hinder the adaption of perturbation theory in the treatment of the optical response of graphene 3,23 , fortunately, over the optical wavelength range, nonlinear response coefficients can still describe interactions of photons. Perturbative treatment of the nonlinear optical response of graphene yields explicit expressions for the higher orders conductivity tensors. Throughout this work we refer to the third order frequency mixing conductivity as σ ω ω ω ( , , ) p q r (3) where ω p , ω q and ω r are the frequencies of the photons undergoing nonlinear interaction through the graphene monolayer. The conductivity is proportional to the Kerr coefficient through the third order susceptibility. The studies show large tunability in n 2 through Fermi level modulation. Other factors such as the spectral and temporal dependence have seen relatively less explicit theoretical discussion.
The experimental studies probing the Kerr-nonlinearity in graphene have shown variation in n 2 spanning six order of magnitudes (10 −12 to 10 −7 cm 2 /W) 4,5,11,14 . A point of contention while comparing these results are the varying spectral and temporal properties of the excitation source, sample preparation techniques, and substrate material. The measurements are performed with the excitation wavelength ranging from the visible to the mid-infrared, while the pulse duration ranges from hundred femtoseconds to tens of picoseconds, under different sample preparation methods 24 . It has been recognised that the n 2 depends on the wavelength, pulse duration and power in semiconductor systems 25 . Particularly in graphene, the n 2 expected to become larger with longer wavelengths 4 . This trend has been observed and reported thereafter in the near 15 and mid-infrared 12 range, while spectral dependence in other regions remains experimentally unexplored. It has also been shown in graphene oxide, that the nonlinear parameters of the system change when excited with nanosecond and picosecond pulses 26 . The decrease of n 2 with increasing incident intensity due to saturation effects has been observed and reported previously in several studies 5,12,15,27 . Raman studies have shown p-type doping occurs when chemical vapour deposition (CVD) graphene is transferred onto quartz 24 . Considering these factors, it is more apt to consider the measured nonlinear parameters in conjunction with the experimental conditions. This makes it imperative to understand the dependence of these factors to gain a fundamental understanding of the governing processes.
In this study we report a systematic analysis of the near-infrared spectral dependence, the sub-picosecond temporal evolution and pulse-width dependence of the effective Kerr coefficient (n 2,eff ) of graphene in hundreds of femtosecond regime using our pump-probe integrated Z-scan (PPZS) set-up. We refer to the nonlinear refractive index (n 2 ) as n 2,eff because we find that the nonlinearity in graphene comprises of the conventional Kerr-type nonlinearity and a component stemming from graphene's unique structure. In addition to that, our set-up allowed for temporal measurements, unravelling the evolution of the nonlinearity in a pump-probe style measurement to reveal an ultrafast nonlinearity. Thus, this study provides a foundation for the development of predictable models imperative for device design, allowing for large controllable tunability of the nonlinear response of graphene.
Originally introduced by Sheik-Bahae et al. 28 , the Z-scan technique proved to be an experimentally facile yet sensitive method to extract the phase and magnitude of the Kerr coefficient. The optically induced self-refraction is quantified by relating the phase modulation of the traversing beam to the transmittance in the far-field in Closed Aperture (CA) configuration, while the Open Aperture (OA) configuration captures the effect of absorption. This technique was modified by Wang et al. 29 to extract the temporal evolution through the integration of a secondary time delayed beam to obtain the time-resolved Z-scan measurement. The PPZS set-up is schematically illustrated and further elucidated upon in the Methods and Materials section. The set-up can be operated in multiple modes of measurement. In the single-beam mode, the set-up is a standard Z-scan measurement. In dual-beam mode, with the addition of a lock-in amplifier, chopper and cross-polarisation filtering, the set-up can be used to perform a standard pump-probe measurement for the extraction of relaxation time constants. When used in conjunction with the Z-scan components, the dual-beam mode is used for the time-resolved Z-scan measurements.

Results and Discussion
spectral-dependence. In order to verify the quality and monolayer nature of the samples, the CVD fabricated graphene on quartz (GoQ) samples (1 × 1 cm 2 ) were characterised by Raman spectroscopy, shown in Fig. 1a. The G-band, 2D-band and D-band peaks appear at 1576 cm −1 , 2661 cm −1 and 1328 cm −1 , respectively. The relative intensities and widths of the G and 2D peaks confirm that the sample is single layer. We also perform a temporal correlation measurement on the sample, shown in Fig. 1b, to obtain a τ 1 relaxation time constant of about ~113 fs. The methodologies and results are further discussed in the Methods and Materials section. www.nature.com/scientificreports www.nature.com/scientificreports/ The spectral dependence of n 2,eff is obtained in the single-beam Z-scan mode and is measured by scanning the 100 fs excitation source from 800 to 1050 nm with 50 nm increments. The range of the excitation is limited by tunability of the laser (690-1050 nm). Since the measurement is sensitive to the beam quality, the beam is characterised using a knife edge beam profiler for quality and stability over the period of the measurement. To ensure that the Z-scan modulation observed in the far-field is originating from graphene, a bare quartz substrate is measured at 900 nm at high on-axis irradiance to show no effect (see Supporting Information Fig. S5). The on-axis irradiance ranges from 1.5 to 5 GW/cm 2 , limited by the available power of the laser at each wavelength. A standard acquired Z-scan trace showing the OA, CA and Fit is shown in Fig. 2a. Our study focuses on the nonlinear refraction in graphene and not the nonlinear or saturable absorption. The OA trace is simply acquired for the normalisation of the CA trace for absorption effects, and presented for illustrative purposes to explain the quantum theory as presented shortly. The effects of saturation resulting in decrease of n 2,eff with increasing irradiance is observed and shown at 950 nm in Fig. S6 in supporting information.
The cumulative results of the spectral dependence are presented in Fig. 2b along with the quantum theory based fit. A full set of Z-scan fits from which these results are compiled is also provided for reference ( Fig. S7 in supporting information). The values for n 2,eff in this regime are measured to be positive and range from 9.07 × 10 −9 to 1.76 × 10 −8 cm 2 /W in this excitation range which fits well with our quantum model with relaxation coefficients of γ 1 ≈ γ 2 ≈ 10 meV. The fitting parameters, γ 1/2 , are obtained from the time constants derived from the temporal cross-correlation measurement, shown in Fig. 1b. However, the pulse width of the pump and probes signals pose some uncertainties on the obtained decay constant of τ 1 ~ 113 fs (see Methods and Materials) as the correlated pulses are of similar pulse durations. This uncertainty is accounted for in our theoretical model and is illustrated by the shaded (yellow) region in Fig. 2b. Shorter wavelengths are observed to have a n 2,eff value that is lower as compared to longer wavelengths, with approximately the same on-axis irradiance. Our calculation suggests a quadratic dependence of n 2,eff on λ which agrees with the experimental results. It is noted that n 2,eff exhibits negligible dependence on the Fermi level for a low-doped graphene monolayer. In our case n 2,eff is dominated by contributions from interband transitions and since the λ 2 dependence is a direct consequence of the linear band diagram around the Dirac point, the most relevant transition occur at the zero detuning region where ω = 2|k|v F . Detuning is defined as where ω is the frequency of the excitation photon and cv  is the energy of the transition, see Fig. 3a. The theoretical model employed in our theory uses semiconductor Bloch equations (SBEs) to describe the cooperative intra-interband dynamics of the population difference  k t ( , ) (between the valence and conduction bands) and the polarisation (coherence)  k t ( , ) for the Bloch state k. The phenomenological relaxation coefficients, γ 1/2 , account for the collective broadening effects for the population and coherence decay, respectively. The theoretical derivation is succinctly provided in the Methods and Materials section. Under our theory, the electromagnetic coupling for normal illumination is defined by Fig. 3a. The nonlinear optical absorption has multiple origins namely bleaching effects culminating in saturation, and two photon absorption. Since the graphene sample used in our experiment is low-doped (i.e μ ≪ ℏω), absorption bleaching due to optically induced Pauli-blocking plays the leading role. This assumption is supported by observation of the absorption drops upon high optical irradiance, refer to the OA trace in Fig. 2a. Furthermore, since the Fermi energy level is small compared to the energy of the photons, it naturally follows that Pauli blocking around zero detuning region is widely demolished. Under high intensity illumination in a Kerr-type material, the charged carriers undergo ultrafast Rabi oscillations, much faster than their relaxation rate. As a consequence, the steady state population, st k  , is significantly modified by light, disturbing the distribution at equilibrium, eq k  . The origin of the Kerr-type nonlinearity is the optically induced change to the steady-state population difference. The relaxation dynamics as the population difference evolves is schematically shown in Fig. 3b. Before saturation takes www.nature.com/scientificreports www.nature.com/scientificreports/ place, the nonlinear contribution of the field to the population difference around the zero detuning region is a quadratic function of the field magnitude, The induced nonlinear current oscillating at the frequency ω is then given by accounts for the interband transitions. The nonlinear current J NL is the microscopic origin of the intensity dependence of the refractive index so that n 2 ∝ |J NL |/|E| 3 . For a small enough γ 2 , the Kerr-type nonlinear induced current is given by: where g s and g v are the spin and valley degeneracy factors, respectively, D = 1/4π 2 is the density of states, and β~π is a dimensionless quantity that arises from angular integration around the Dirac cone. Due to low doping,  ≈ 1 eq k over the zero detuning circle. Therefore, the frequency dependence of nonlinear current is dominated by the 1/k 2 term appearing in Eq. 1. Over the Δ k = 0 circle the Bloch wave-number k is linearly proportional to the frequency which in turn yields λ 2 dependence of the Kerr nonlinear coefficient.
It is worth pointing out that the quadratic wavelength dependence of the Kerr coefficient is a direct consequence of the linear energy-momentum dispersion of the Dirac quasi-particles. The calculations outlined above can be effectively applied to all two-level condensed matter systems by appropriately replacing interband coupling into Eq. 1. However, exclusive to graphene is the 1/k dependence of the interband coupling Φ k which results in the square wavelength dependence of the Kerr coefficient. temporal evolution. The temporal evolution of the nonlinearity is obtained using the dual-mode PPZS set-up. A single-beam (pump) Z-scan measurement is performed to locate the peak and valley positions of the sample. In order to locate the zero-delay position of the probe, the sample (GoQ) is placed at the peak position and the probe is scanned until a cross-correlation signal of the pump and probe pulses is obtained, shown in the inset of Fig. 4. The FWHM of the signal is 150 fs, which gives a pulse duration of 110 fs for the probe, considering a pulse duration of 102 fs for the pump. A similar scan is run at the valley position and both time-resolved data sets are used to extract the differential peak-valley transmittance, are the positions of the peak and valley, respectively, t d is the probe delay, and T OA is the OA transmittance for normalisation. The sign of ΔT pv is given by the sign of Z p − Z v . The pump power used for this measurement was 200 mW, equating to an intensity of ~3.5 GW/cm 2 . The pump to probe power ratio was kept at 20:1. In general, the far-field aperture should only allow about 1% of light transmittance to isolate the effect of wave front distortion due to phase modulation, however, due to the weak signal of the probe the aperture is opened to allow 10% of transmittance with additional averaging at each acquisition point. The temporal evolution of the nonlinearity is shown in Fig. 4b. The effect peaks at zero probe delay and relaxes on the time-scale of τ 1 , while at longer time scales, τ 2 , ΔT pv shows no discernible variation. The symmetric shape of the figure leads us to believe that in this measurement we are simply observing the probe following the pump signals rather than unravelling the nonlinear phenomena hidden at shorter time scales. This is also supported by the fact that the pulse duration of the pump and probe is similar to the measured temporal cross-correlation extracted in Fig. 1b. Therefore, we can conclude that the relaxation dynamics contributing to the observation of the nonlinear refraction are simply too fast to be measured by our laser pulses in this manner. The variation of ΔT P−V or the induced phase shift ΔΦ seen in Fig. 4 relates to the Gaussian power distribution within the pump pulse which reaches its maximum value (i.e. ~3.5 GW/cm 2 ) at the centre peak of the pulse. Therefore, according to =  www.nature.com/scientificreports www.nature.com/scientificreports/ As the intensity decreases moving away from the peak, the n 2,eff increasing in accordance with saturation. This measurement clearly shows that on the time scale greater than the heating and cooling times, the probe follows the pump and Kerr nonlinearity is tunable and controllable over this time scale. In this context, the heating and cooling times refers to the time-scale of carrier excitation to higher states and subsequent carrier relaxation to the ground state, respectively 30 . The modulation is understood through the relation n = n 0 + Δn, where Δn = n 2,eff I, and the Δn parameter is modified as the n 2,eff varies nonlinearly on the intensity prior to complete saturation 5 .
For the purposes of all-optical switching the on/off time of the nonlinearity is controlled by the pulse-duration and power of pulse. If the pulse-duration is longer than the effective relaxation times of the nonlinearity, then the evolution of the observed nonlinearity simply follows the Gaussian power distribution of the pulse. In this way through modulation of the Δn term, by changing the pulse properties, the switch can be controlled. However, there is saturation that takes place at relatively high intensities, leading to a deviation between the pump power within the pulse and the observed Δn.
pulse-width dependence. It is well understood that the relative timescales of the excitation pulse and the system response times determine the induced dynamics in the system. To investigate the effects of pulse duration on nonlinear refraction, we stretch the ~100 fs pulse upto ~475 fs at 900 nm to observe the effects in the long pulse regime; longer than the duration of the nonlinearity or relaxation time, τ 1 . Pulse stretching is achieved using a dispersion-based prism-pair apparatus, where the pulse can be temporally expanded by varying the separation distance of prisms 31 . A schematic of the set-up is shown in Fig. S8 in supplementary information. In our set-up we use N-SF11 flint glass prisms. The maximum allowed pulse expansion is limited by the space on the optical table. The autocorrelation curves for the expanded pulses at which the Z-scan is performed are given in Fig. 5a and the results of the pulse-width dependence are provided in Fig. 5b. A set of Z-scan fits for the presented data is also provided in Fig. S9 in supporting information.  www.nature.com/scientificreports www.nature.com/scientificreports/ The results show a clear dependence of n 2,eff on the excitation pulse duration, with n 2,eff becoming larger with increasing pulse-width. This trend was theoretically speculated by Vermeulen et al. 10 and shown here experimentally. The n 2,eff ranges from 1.02 × 10 −8 to 1.7 × 10 −8 cm 2 /W in this pulse duration regime with a minimum between 200-300 fs. The decay constant τ 1~1 13 fs measured in Fig. 1b falls near the beginning of the data set. However, doping can modify the decay constant values with p-and n-doping making the time constants long and shorter, respectively 32 . This data trend coincides with previously reported Z-scan studies on graphene performed with picosecond excitation where the n 2,eff is larger than what is reported in the femtosecond regime. In addition to this, a similar comparative analysis performed on carbon disulphide (CS 2 ), reference material used for calibration of Z-scan measurements, revealed a similar dependence of the n 2,eff on pulse duration 33 . In a recent publication 10 , this pulse-width dependence of n 2,eff is theoretically derived for a regime when the effective decay constant is larger than the pulse duration in Z-scan measurements. The interplay of relative carrier heating and cooling times is said to induce a nonlinear response that may not originate only from the conventional electronic Kerr-type nonlinearity but also from what they refer to as saturable photoexcited-carrier refraction (SPCR). As mentioned previously, the nonlinearity originates from the carrier population difference. Given a pulse of an arbitrary pulse duration that impinges on the system, the maximum excited carrier density is reached at the peak of the pulse. If the pulse duration is shorter than the system's relaxation constant i.e., τ 1~1 13 fs, then the maximum carrier density is reached on the time-scale of the pulse duration, the carriers then remain in the excited state even after the pulse passes, subsequently relaxing on a time-scale of τ 1 . However, if the pulse duration is longer than τ 1 then carriers continuously undergo excitation and relaxation as the pulse passes with the carrier population difference continuously evolving, and therefore the optical nonlinearity becomes dependent on the duration of the pulse. The saturability in graphene deviates from saturability in other 2D materials due to the presence of its unique gapless band structure which facilitates spontaneous saturation near the Dirac point even when there is no field 23 . Therefore, it is apt that we refer to the Kerr-type nonlinearity characterised using the Z-scan method as n 2,eff due to the presence of this contributing mechanism.

Conclusions
The systematic measurement and analysis of the dependence of n 2,eff on the spectral and temporal properties of the pulse, has allowed us to gain a more fundamental understanding of the parameters governing the observed nonlinear optical effect in graphene. This sheds significant light on the widespread debate in the field regarding the large variation in the n 2,eff value for graphene due to varying experimental conditions and sample preparation techniques. The dependence of n 2,eff on the exciting wavelength revealed a quadratic (n 2,eff ∝ λ 2 ) relationship observed experimentally and confirmed theoretically. The quantum theory developed can be used to calculate the nonlinear response under varying the relaxation coefficients and Fermi energy, and the results can be extended to other spectral regions. The time-resolved Z-scan measurement revealed that the heating and cooling dynamics within our graphene sample are simply too fast to be probed with ~100 fs pulses using this method. However, this method reveals a practical application in all-optical ultrafast switching of the nonlinearity, completely controlled by the pulse-duration and power of the impinging laser pulse.
Through the pulse-duration dependent measurement, we were able to confirm the theoretically predicted relationship between n 2,eff and the laser pulse-duration in the hundreds of femtosecond regime, with n 2,eff growing larger with longer pulse-duration. Our results also coincide with the trend that is observed over several papers 4,5,11,14 performing the measurement with different pulse durations, further reconciling our measurements, theory and experiments. In our measurements we used a prism-pair set-up to expand the pulses which is highly sensitive to beam alignment, posing a challenge if the pulse duration needs to be changed dynamically. Alternatively, transmission diffraction gratings can provide large pulse expansion ratios with high efficiencies and easier alignment 34 . Using the diffraction gratings a larger range of pulse durations can be scanned and a decisive relationship between n 2,eff and pulse duration can be established.
Throughout all our experiments, the value for n 2,eff remains positive. The sign of n 2,eff has been a constant point of disagreement in literature, with the sign alternating from positive to negative depending on the method of characterisation employed or the conditions of the experiment. Theoretical studies have revealed that n 2,eff is very sensitive to the Fermi energy of the sample 3,17,18,21,22 . Fermi energy modulation can be practically achieved with ease by electrically gating graphene samples in contrast to changing the spectral and temporal properties of the pulse to achieve the same tuning. Therefore, it becomes absolutely critical to perform an experimental characterisation of the dependence of n 2,eff on Fermi energy. By scanning this parameter and observing the nonlinear response, it is possible to observe if or when the change in sign takes place, which would resolve the supposedly contradictory reports. With this study we have gained critical fundamental understanding of the underlying processes governing the nonlinear optical phenomena and the effect of experimental parameters which can be used to largely tune the effect.

Methods and Materials
pump-probe integrated Z-scan measurement. The experimental set-up for this measurement is illustrated in Fig. 6. When operating in single beam mode, the set-up is a classic Z-scan measurement. The OA is used to normalise for absorption and laser fluctuations in the CA trace. In dual-beam mode, we use degenerate pump and probe at 900 nm, in a collinear configuration with a 20:1 pump/probe ratio. The pump and probe are orthogonally polarised using a λ/2 waveplate (λ/2), polarising beam-splitter (P-BS) and polariser (Pol1) before the sample. An analyser polariser (Pol2) placed after the sample is rotated to achieve extinction of the pump beam.
The normalised CA transmittance is fitted to the equation: , where x is the Rayleigh length (z R ) normalised position (z/z R ) and ΔΦ 0 is the nonlinearity induced phase shift. When a pulsed source is used with pulse-width comparable to the duration of the nonlineairty, the nonlinear refractive index is www.nature.com/scientificreports www.nature.com/scientificreports/ , where k 0 is the wave vector, I 0 is the on-axis irradiance and L eff is the effective length of the sample 35 . In our analysis the value for ΔΦ 0 is extracted using a nonlinear least-squares curve fitting routine in MATLAB. The routine takes in the aforementioned function along with a starting value for ΔΦ 0 , and a vector constraining the upper and lower bound for the value. The routine then works to find a minimum for the sum of squares of the input function in an iterative manner, with the allowed number of iterations being easily modified. In our measurements the value for ΔΦ 0 falls in the range of [−0.04, −0.15], which was the upper and lower bound input, while a starting point of 0 with 5000 iterations generally allowed for convergence. sample characterisation. In order to verify the quality of the samples, graphene on quartz (GoQ) samples (1 × 1 cm 2 ) were characterised by Raman spectroscopy, shown in Fig. 1a. The samples were purchased from an external supplier and were fabricated using chemical vapour deposition (CVD) onto a large copper sheet. The copper sheet was then cut into the respective sample size and graphene was transferred from the copper onto the quartz substrate. The Raman spectroscopy was performed using a home-built set-up, combining a 561 nm continuous wave (CW) laser and a high-resolution spectrometer. Prior to every measurement the laser input was read to normalise the spectral data to output the Raman shifts. The D band, G band and 2D band peaks appear at 1328 cm −1 , 1576 cm −1 and 2661 cm −1 respectively. The D band signifies disorder in the graphene lattice and is typically weak in high quality graphene, as seen in our sample measurement. A point to note as it pertains to the D band is that it is a resonant band and has very weak modes underlying it. The shape and position of this band can vary significantly depending on the excitation laser wavelength. The G band is a signature of an in-plane vibrational mode stemming from the stretching of the C-C bond, characteristic to all sp 2 carbon systems. The 2D band is the second overtone of the D band and is the result of a two-phonon lattice vibrational process. It is always strong in graphene and does not represent defects regardless of the presence of the D band. The peak at ~2460 cm −1 has been referred to as the G* band and is attributed to a double resonance 36 . The peak appearing at ~2325 cm −1 may arise from additional defects or resonances. The relative intensities of the 2D and G peaks confirm that the sample is single layer.
We also perform a temporal correlation measurement on the sample, shown in Fig. 1b. The measurement was performed in the dual-beam mode of the PPZS set-up using the translation stage and cross polarisers to obtain extinction of the pump. An extinction of ~10 −4 was achieved with this configuration. The pump/probe power ratio was kept at 20:1, with the pump power being 200 mW. Due to the short time-scale of the temporal correlation Figure 6. Pump-probe integrated Z-scan set-up (PPZS). The source is Ti:sapphire tunable laser capable of emitting from 690 to 1050 nm at a repetition rate of 80 MHz at ~100 fs. The laser has a variable high power attenuator on exit comprising of a half-waveplate (λ/2 (780) ) and a Glan-laser polariser (GP). The beam is directed by silver mirrors (M) into another half-waveplate (λ/2 (900) ) and polarising beamsplitter (P-BP) to prepare the orthogonal pump (transmitted) and probe (reflected) beams. Both beams are then directed through the dualfrequency chopper (Chp). The pump is directed into a small static retroreflector (Refl-S), while the probe is directed through a polariser (Pol1) and a large retroreflector (Refl-L) that is mounted on a translation stage that allows for one beam to be delayed relative to another. Upon return, both beams are combined at a pellicle beamsplitter (Pel) and guided to the Z-scan set-up. Both beams are then focused using an achromatic doublet lens (ADL) and impinge upong the graphene on quartz sample (GoQ), followed by another polariser (Pol2) which is oriented parallel to Pol1 for extinction of the pump beam. The beam is then bisected in the far-field by an non-polarising beamsplitter (NP-BP) with the reflected arm directed into the open aperture (OA) detector and the transmitted arm directed through an adjustable aperture (AP) into the closed aperture (CA) detector. For all Z-scan based measurements, single-beam or temporal, the power meter is utilised for data acquisition. In pump-probe mode the chopper and Lock-in amplifier are used for data acquisition.
www.nature.com/scientificreports www.nature.com/scientificreports/ of the two pulses, a stage capable of achieving a repeatable resolution of ~5 femtosecond was employed. Lower resolution could be reached but was seen to be unreliable for a mechanical stage. The data is obtained over a time-scale of 1.5 picoseconds, limited by speed of the stage and the oscilloscope sampling rate. The data is fit to a single exponential function of the form: exp(t/τ), to obtain a τ 1 relaxation time constant of about ~113 fs. theoretical modelling of the nonlinear optical response in graphene. The nonlinear response coefficients have been calculated via Semiconductor Bloch Equations (SBEs). In the continuous excitation regime, the optical response can be estimated based on the dynamics of (i) the microscopic population difference  t k ( , ) and (ii) the microscopic polarization t k ( , )  . In a moving frame {τ, k′} = {t, k − δk(t)}, where δk obeys δ and Γ is an phenomenological intraband relaxation coefficient, the dynamics of the charge carriers is governed by: ( ( , ) ) 2 ( , )Im{ ( , )}, is the matrix element of the external potential of the direct optical transition, and the unit vector φ k is defined as ϕ = ×ẑ k k/ k . The frequency E ϖ = 2 k k is the energy difference between the energy levels of the conduction and valence bands. At = t 0 before the electromagnetic field being applied, the population difference  is relaxed at its equilibrium eq F F where f is Fermi-Dirac distribution function. The moving frame accounts for intraband dynamics which has significant contribution in high harmonic generation via coherent Bloch oscillations.
Since the crystalline structure of graphene is centrosymmetric, even order nonlinearities are prohibited and therefore the first nonlinear term is the third order. We describe the effect of this third order nonlinearity by considering three complex fields with the time dependence of ω e i t p , ω e i t q and ω e i t r and their mixing through the third order conductivity of graphene. The third order optical response can be interpreted as a three-photon process and different terms contribute to the third order conductivity tensor namely pure intraband term, pure interband term and combination of the both. There are six distinct photon processes contributing to the third order optics of graphene. However, since intraband processes mostly concern harmonic generation, in optical range the interband process plays the leading role in Kerr type linearity. SBEs can be combined and recast in the form of a damp-driven harmonic oscillator equation. The distinctive nonlinear process can be then obtained using an iterative technique 23 .
The conductivity tensors for the six distinct photon processes contributing to the third order optical effect in graphene are presented in the supplemental document. The leading term (purely interband term) is: