1. Introduction
The Hyper-Rayleigh scattering (HRS) technique have been employed for several years to investigate the nonlinear optical properties of molecules in liquid suspensions as well as colloidal metal particles in aqueous solutions, nanocrystals in liquid suspensions and various other materials. It entails directing a laser beam onto a sample and detecting the scattered light, which has double the frequency of the incident light. The intensity of the scattered light is contingent on the first hyperpolarizability of the sample, a measure of its ability to alter the polarization of light. In principle, HRS technique provides a direct assessment of the magnitude of the first hyperpolarizability of molecules and nanoparticles that describes the second-order nonlinear optical (NLO) response. Substantial NLO activity is a crucial prerequisite for emerging technological photonic applications such as photonics and all-optical switching and optical parametric oscillators, that have intensified the quest for NLO-active materials. For example, previous research has underscored the remarkable NLO response in solutions of metallic nanometer-sized particles in solutions. The origin of this response lies in the partial one- or two-photon resonance of the incident laser light used in the experiments, coupled with the particles’ strong surface plasmon absorption. Graphene quantum dots (GQDs) are nanometer-size graphene segments that are small enough to exhibit exciton confinement and quantum size effects. Unlike graphene, they present non-zero energy bandgap that is responsible for their characteristic electrical and optical properties.
Generally, the optical absorption spectra of GQDs present a strong peak in the ultraviolet (
due to the
excitation of the
bonds of aromatic C=C and a weaker peak at
due to
transitions of C=O bonds [
1,
2]. Introduction of functional groups in the GQDs may lead to new absorption features and may change the GQDs optical properties [
3,
4]. Moreover, since they present molecule-like structure its bandgap can be tuned by changing the GQD size and by surface chemistry. Therefore, the optical properties of GQDs are dependent on the preparation method as well as the functional groups added to the graphene segments [
3,
4,
5,
6,
7,
8].
When compared to heavy metal-based semiconductor quantum dots that causes concern for in-vivo bioimaging, GQDs due to their low toxicity [
9,
10] and strong two-photon induced fluorescence have been exploited for biological applications by many authors [
3,
11,
12]. Likewise, there are reported applications of GQDs for energy conversion, environmental monitoring, and electronic sensors for humidity detection, among other applications [
3,
4,
7,
13,
14]. Recently, Qi
et al. reported second harmonic generation (SHG) in boron doped-graphene quantum dots (B-GQD) and their application in stem cell imaging and tracking wound healing. These applications were achieved because the doping promoted the break of symmetry in the B-GQD allowing the SHG signal [
15].
Although optical absorption, linear luminescence, second harmonic generation and the third-order nonlinearity of GQDs have been studied [
5,
6,
7,
8,
16,
17,
18,
19], the second-order nonlinearity, associated to the first hyperpolarizability,
, did not receive much attention. It is true that the experiments to be done are not very easy to set up, such as measuring nonlinear refraction in a 4f setup (see for example [
20,
21]). Computational results were reported in [
22,
23,
24] for some GQDs but no experimental result is available. Here we want to estimate the order of magnitude of these coefficients knowing that the manufacturing methods relative to new materials, sometimes different from one lab to another, could influence slightly the result. The absorption spectra of certain DIY items, assumed to be identical, exhibit partial mismatches owing to the presence of significant residues (impurities) within the final solution. Hence, it proves beneficial to examine readily accessible, commercially available samples. Preparations employing hybrid manufacturing methods may yield slightly divergent results but should not alter the core conclusion in this context.
In the present work the hyper-Rayleigh scattering (HRS) technique [
25,
26] was employed for the first time to determine the effective first hyperpolarizability,
, of GQDs suspended in water. The HRS is an incoherent nonlinear process in which two photons with frequency
are instantaneously combined to generate new photon with frequency
due to the interaction with a material. Unlike the two-photon induced luminescence process the generation of second harmonic in the HRS process does not require on-resonance excitation of the medium to high energy levels. The HRS technique has been routinely used to study molecules [
25,
26,
27,
28,
29], metal nanoparticles in liquid suspensions [
30,
31,
32], and dielectric nanocrystals [
33,
34,
35,
36,
37]. The first hyperpolarizabilities of molecules and nanoparticles can be determined as well as the origin of the incoherent second harmonic generation (SHG) signal detected.
In the present paper we report values for the first hyperpolarizabilities of two different GQDs [
38] investigated by measuring the intensity of the incoherent SHG signal upon excitation using a laser off-resonance with absorption transitions of the GQDs. Values of
were determined applying the HRS technique using for calibration para-nitroaniline (
p-NA) which is the usual reference standard for HRS measurements. To determine the static hyperpolarizability per nanosheet,
, the classical two-level model [
37,
39] was applied in combination with the results for
and considering that each GQD has five nanosheets with thickness of 0.3 nm each. Studying the nonlinear optical properties of these materials can provide insights into their structure and functionality involving the interaction of light with nanoparticles at a nanoscale level.
3. Results and Discussion
Figure 2a shows the absorbance spectra of both GQDs samples. The absorption bands in the ultraviolet are associated to the
and
transitions. The arrows in
Figure 2a indicate the laser and the HRS wavelengths. Notice that for both samples the laser and its second harmonic are off resonance with the absorption bands; therefore, we do not expect a PL signal superimposed on the HRS signal.
Figure 2b shows the normalized PL spectra obtained by excitation at 386 nm for the Acqua-Cyan GQDs and 486 nm for the Acqua-Green GQDs using a commercial spectrofluorometer (Fluoromax HORIBA). The temporal decay of the PL signal with maximum at 550 nm (Acqua-Green) and 480 nm (Acqua-Cyan), excited at 850 nm (via two-photon absorption), was fit by a single exponential with decay time of 26 ns (Acqua-Green) and 19 ns (Acqua-Cyan).
Figure 3 shows the HRS spectra and the temporal behavior of the HRS signal (obtained for a single laser shot). The solid red lines in
Figure 3 are the Gaussian fit to the experimental data and the solid black lines on the inset of
Figure 3b show the temporal evolution of the excitation laser. The duration for both GQDs was determined to be 4.6 ns from a Gaussian fit. Notice in
Figure 3a that the spectra are centered at 532 nm, the second harmonic of the incident laser beam, with the linewidth limited by the spectrometer resolution. The temporal evolution of the HRS signal exhibited in the inset of
Figure 3b follows the laser pulse corroborating the above statement that no long-lived luminescence is overlapping with the HRS signal. The HRS is a parametric phenomenon, and the strong correlation between the temporal behavior of the signal and the excitation source supports this assertion. These measurements were performed with laser peak intensity of
, and the minimum intensity that allowed detection of the HRS signal was
.
The HRS signal is due to the nonlinear scattering by the randomly oriented GQDs and water molecules. The scattered intensity, at the second harmonic frequency of the incident laser beam, is described by
, where
is the laser intensity,
is the orientational average of the first hyperpolarizability and the sub-index represents the constituents of the suspension. The symbol 〈 〉 indicates orientational average. The factor
depends on the scattering geometry and contains information on the transformation of coordinates from the GQD to the laboratory reference system [
33,
34,
35,
36,
37].
Figure 4a shows plots of the HRS signal,
versus the incident laser intensity where a quadratic dependence with
was verified.
Figure 4b shows linear dependence of
versus the GQDs concentration,
that confirms the contribution of isolated GQDs to the HRS signal. Moreover, we recall that aggregation of GQDs would redshift the linear absorption spectrum, a behavior that was not observed for the various concentration used.
Since the samples consist of GQDs with different dimensions, the experiments provided effective values for the GQD hyperpolarizability,
. The external reference method [
41], using
p-NA dissolved in methanol as the reference standard, was applied to measure
and the numerical values were obtained using the equation:
where
refers to methanol,
[
41],
molecules/cm
3,
[24,25]],
molecules/cm
3,
[
25],
molecules/cm
3. The concentration of GQDs was
for both samples that were suspended in water. Considering the results obtained for
(see
Figure 4a), the effective values of
were
and
for Acqua-Cyan and Acqua-Green, respectively. The quantity
is an effective value contributed by GQDs with different number of nanosheets. Consequently, we considered the following approach to estimate the hyperpolarizability per nanosheet,
. First we considered the distribution of GQDs sizes shown in the
Appendix A; then, assuming the nanosheet thickness 1.5 nm, and a Gaussian distribution of GQDs sizes represented by
, where
is the number of nanosheets in a GQD, we obtain
for Acqua-Cyan and
for Acqua Green. The value of
is obtained from the AFM measurements given in the SI. The value of
is determined by
. Using these results, we could estimate the effective first hyperpolarizability per nanosheet,
, considering that each GQDs have five nanosheets with 0.3 nm of thickness.
The two-level model [
39,
42] was used to relate
with the first static hyperpolarizability per nanosheet,
in each GQD
First, it is necessary to estimate the static hyperpolarizability,
, for the samples. Then, considering the amount of graphene sheets per GQD,
is estimated. The following expression was obtained from [
39,
42]
where
= 22222 cm
-1 (Acqua-Cyan) and
= 19608 cm
-1 (Acqua-Green). The frequency
correspond to the smallest band frequency in the absorption spectrum.
The results for are and for Acqua-Cyan and Acqua-Green, respectively. Notice that the value of is approximately equal for both samples but the values of are different because of the different laser frequency detuning with respect to the first excited resonance frequency .
We recall that imaging measurements have been performed using GQDs exploiting their large PL, but SHG microscopy presents an important advantage with respect to PL experiments because lasers operating in a large infrared wavelength range can be used since no resonance with excited states is required. This is a benefit because one may select the laser wavelength according to the transparency window of the biological system of interest. Obviously, care must be taken not to damage living cells and the underlying nanoparticles with too high intensities, as can often happen with thin films irradiated with pulsed lasers in the picosecond regime [
43]. To characterize the origin of the HRS signal, polarization-resolved experiments were performed following the usual procedure [
30,
31,
32]. The laser beam was linearly polarized and the angle
between the optical field and the vertical direction, was varied from 0 to 2
. The polarization of the HRS signal either in the vertical direction or parallel to the beam propagation axis, was measured. At this point, due to the insertion of half-wave plate and a polarizer (see
Figure 1), the laser intensity was
0.7 GW/cm
2. The results are shown in
Figure 5 and to fit the data it was used the expression:
where
is the hyper – Rayleigh scattering intensity and Γ refers to the horizontal (H) or vertical (V) polarization. The coefficients
,
,
,
,
, were obtained from the fit. The coefficient of multipolarity,
, were also determined and their values are indicated in
Table 1.
The origin of the second-order nonlinear effect arising from nanoparticles (NPs) with crystalline structure presenting inversion symmetry was discussed in [
44,
45,
46,
47]. Two important cases to consider are: NPs with noncentrosymmetrical shape and NPs with centrosymmetric shape [
47]. In the former, SHG is due to the induced electric dipole in the NPs volume and in the interfaces. In the latter, the induced electric dipole does not occur because of the inversion symmetry in the crystalline structure of the NPs. Then, besides the contributions from the interfaces an important mechanism of SHG is related to the field’s gradients occurring within the NPs bulk [
46,
47]. However, it is difficult to synthesize NPs with perfect shape and therefore, in real NPs there is a competition between the two main cases above mentioned. Thereby, the coefficient of multipolarity,
, is applied to distinguish both contributions [
47]. The case for pure centrosymmetrical shape requires
and the case for the pure noncentrosymmetrical shapes requires
. As one can see in
Table 1, the values of
for both samples suggest the dominant mechanism for harmonic generation is the field retardation in the NPs. This fact should be expected in our samples since there is no guarantee that the nanosheets are perfectly compacted and discordances should also occur between the sheets. As a result of such discordances, discontinuities in the electric field should occur in the NPs volume and the quadrupole contributions become dominant.
Ongoing research is focused on assessing the toxicity of graphene-family nanoparticles, with specific attention to GQDs [
48]. Factors such as particle size, synthesis methods, and chemical doping play a role in determining their toxicity, both in vivo and in terms of cytotoxicity [
49]. Some researchers argue that GQDs, composed primarily of organic materials, exhibit low toxicity and high biocompatibility, providing an advantage over semiconductor quantum dots. In vitro studies using cell cultures have shown minimal bad effects of GQDs on human cell viability [
50]. Fluorescence imaging, a non-radioactive method, allows the visualization of morphological details in various biological specimens, from living cells to animals. Notably, GQDs, unlike many semiconductor quantum dots containing heavy metals, consist of graphene lattices containing light elements. Because they are primarily made of carbon, the most abundant element in biological systems, GQDs are generally considered biocompatible.