Introduction

The optical lenses are indispensable components in almost all aspects of science and technology including imaging, sensing, communications, medical diagnosis and treatment1,2,3. Recently, the rapid development in nano-optics and on-chip photonic systems has placed stringent demand for ultrathin flat lenses with three-dimensional (3D) subwavelength focusing capability4,5,6. Conventional optical lenses rely on sufficient phase accumulation when a beam propagates over optical paths much larger than the light wavelength to realize the desired wavefronts7. The phase modulation ϕ of a conventional lens can be expressed as ϕn × |k| × Δt, where Δn is the refractive index difference between the lens and the surrounding medium, k (|k|=2π/λ) is the wave vector of light and Δt is the thickness of the lens. Owing to the limited refractive index difference (Δn) of naturally available transparent dielectrics compared with the surrounding medium (mostly air), conventional geometrical lenses are fundamentally incapable of achieving subwavelength-scale (Δt<λ) ultrathin flat lenses with decent 3D diffraction-limited focusing8.

The recent breakthroughs in nanophotonics have led to the design and demonstration of a number of ultrathin flat lens concepts2,5,6,9,10,11,12,13,14, such as metamaterials12, metasurfaces13 or super-oscillations14, which have achieved subwavelength resolution in two dimensions. In particular, all-dielectric metasurface lenses show promising results in reducing the Ohmic losses15,16. However, the narrow operational bandwidth, the complex design and time-consuming multi-step manufacturing process have greatly limited their real-life applications especially when large-scale production and micro-devices integration are required. Therefore, it is highly desirable to develop ultrathin flat lenses with 3D subwavelength focusing capability, a large operational bandwidth, and a practically useful focusing efficiency (>20%) with a simple and low-cost fabrication method.

Here, we propose a new ultrathin flat lens concept that is able to effectively manipulate the phase and amplitude of an incident beam simultaneously. We report on a 200-nm-thick graphene oxide (GO) flat lens with 3D subwavelength focusing that is able to tightly focus broadband light from the visible to near infrared (VIS–NIR; 1,100 nm bandwidth) with an averaged absolute focusing efficiency of >32% over the entire band.

Results

Conceptual design of the GO lens

The conceptual design of the GO ultrathin flat lens and the corresponding theoretical and experimental results are illustrated in Fig. 1. The GO flat lens is made possible by the submicrometre concentric ring fabrication in a sprayed GO film by using the mask-free direct laser writing (DLW) method to covert the GO into reduced GO (rGO) via the photoreduction process, as illustrated in Fig. 1a. The laser power controlled removal of the oxygen groups in the exposed regions of the GO film leads to three continuously tunable local physical property variations: the reduction of film thickness, the increase of refractive index and the decrease of transmission/increase of extinction coefficient (Supplementary Figs 1 to 4)17. As a result, the three tunable local property variations provide an unprecedented flexibility in designing new ultrathin flat lens concept with capabilities of both amplitude and phase control (Fig. 1b), and thus make our GO flat lens concept fundamentally different from the conventional Fresnel lenses and other flat lenses, which mostly rely on either phase or amplitude modulation solely7,18.

Figure 1: The design of the GO lens.
figure 1

(a) Conceptual design and laser fabrication of the GO ultrathin lens. (b) Amplitude and phase modulations provided by the transmission and refractive index difference, respectively, between the GO and rGO zones. (c) Topographic profile of the GO lens measured with an optical profiler. (d) Left: Schematic of the wavefront manipulation by the GO lens converting the incident plane wave into a spherical wavefront. Right: Intensity distributions of the 3D focal spot predicted by the analytical model for a GO lens (a1=1.8 μm, N=3). (e,f) Theoretical focal intensity distributions in the lateral and axial directions. (g,h) Experimental focal intensity distributions along the lateral and axial directions.

To elucidate the light interaction with the GO flat lens, an analytical model based on the Rayleigh–Sommerfeld diffraction theory7 is developed to evaluate the focusing capability of the GO lens. When a uniform plane wave impinges on the GO lens, part of the beam is absorbed and refracted by the rGO zones, experiencing substantial amplitude as well as phase modulations. The other part of the beam propagating through the GO zones only experiences ignorable amplitude modulations. The 3D focusing is a result of the interferences of wavelets originated in the lens plane from different zones19, as illustrated in Fig. 1d. It should be emphasized that although the refractive index modulation from GO to rGO is giant (Δn0.8, one to two orders of magnitude larger than conventional refractive materials20), the resultant maximum phase change between adjacent GO and rGO zones is smaller than π, which is the required phase change for effective constructive interference in the focal region7. Under such a circumstance, the amplitude modulation in the rGO zones can make a substantial contribution to the final well-defined 3D focusing (Supplementary Figs 5 to 7 and Supplementary Note). The positions of each ring hold the key to allow both the amplitude and phase modulations positively contribute to the constructive interference in the focal region to generate a focal spot in the far field. To achieve constructive interference of all the GO and rGO zones in the focal region, the radius of the mth rGO zone, am, needs to satisfy the following condition (Supplementary Note 1):

where f is the designed focal length of the GO flat lens and Δϕ is the phase modulation between GO and rGO zones. It is thus evident that determines the positions of the rest rGO zones. As a result, we use a1 and the total rGO zone number (N) to define the geometry of our GO flat lenses in this paper.

To achieve the optimized diffraction efficiency in the focal region, Gaussian profiles for the phase and amplitude modulation functions (Fig. 1b,c), which can be naturally facilitated by the intensity profile of the fabrication laser (Supplementary Fig. 8), are employed to effectively direct the majority of the incident light to the first diffraction order21. In this way, the theoretical maximal focusing efficiency of 50% can be achieved. In contrast, the conventional rectangular profiles distribute the light energy to a series of diffraction orders that severely degrade the focusing efficiency and resolution of the main focal spot (Supplementary Fig. 9).

Experimental characterizations of the GO lenses

To validate the theoretical model, the dependence of the focusing intensity of the GO lenses on a1 is investigated. The calculated focusing intensity presents a nonlinear dependence on a1, reaching a maximum for a three-ring (N=3) GO lens when a1=1.8 μm. And the experimentally measured dependence of the normalized focusing intensities on a1 agrees remarkably well with the theoretical prediction (Supplementary Fig. 10). The calculated 3D focal intensity distributions for the optimized lens (a1=1.8 μm) is shown in Fig. 1d, with the cross-sectional plots in the lateral and axial directions shown in Fig 1e,f, respectively. Well-defined strong focus with ignorable side-lobes (<10% of the intensity of the central lobe) is observed. The full-width at half-maximum (FWHM) of the focal spot in the lateral direction (FWHMlateral) is 0.8λ and that in the axial direction (FWHMaxial) is 1.6λ, leading to a high-quality 3D focal spot with wavelength-scale resolution. The corresponding GO lens is experimentally fabricated with the surface profile shown in Fig. 1c, in which the Gaussian surface morphology is clearly evident. The measured focal intensity distributions in the lateral (Fig. 1g) and axial (Fig. 1h) directions under 700 nm light illumination almost reproduce the theoretical plots, with a measured FWHMlateral=567 nm (0.81λ) and a FWHMaxial=1.3 μm (1.8λ). The slight difference between the experiment and theory might be due to the imperfection of the fabrication and the surface roughness of the GO film.

Although 3D wavelength-level focusing has been demonstrated for the ultrathin GO flat lens, it is not the best focusing resolution that the GO lens can achieve. According to our theoretical predictions, increasing the ring number N (equivalent to increase the aperture of the lens) and/or reducing a1 can both improve the resolution (Supplementary Fig. 11). However, practically the maximum N and minimum a1 are both determined by the smallest feature size (l) that the laser can fabricate on the GO film. To separate two consecutive rings, the feature size has to be smaller than the distance between the two rings as l<am+1am. Meanwhile, the distance between the two outmost rings determines the maximum convergence angle β of the lens as . Therefore, the minimal feature size of each rGO zone is the key parameter practically limiting the achievable focusing resolution of the lens. However, laser fabrication based on thermal reduction in the literature can only achieve minimum feature size >1 μm (ref. 22), which severely restricts the lateral focusing resolution of the GO flat lens to be over 6.5λ according to our experiments using a MHz repetition rate laser (not shown). Recently, it has been reported that GO has a threshold photon energy for photochemical reduction at 3.2 eV (λ=390 nm)23. Therefore, for a high peak power femtosecond laser beam operating at 800 nm, it is potentially possible to harness the nonlinear photochemical mechanism to enhance the fabrication resolution24,25. As a result, to minimize the thermal effect and potentially make use of the high-resolution nonlinear reduction mechanism, we tightly focus low repetition rate, femtosecond pulsed, laser beam (100 fs, 10 kHz, 800 nm) with a high numerical aperture (1.4) objective to ensure the interval of the laser pulse to be sufficiently larger than the thermal diffusion constant of the GO22. As a result, a record small line width of 300 nm (0.38 λ) is achieved in the rGO laser fabrication (Supplementary Fig. 2), which represents a more than threefold resolution enhancement compared with the literature22. Consequently, maximum six rings with a minimum a1=730 nm can be fabricated. It should be noted that the notation a1 is the radius of the first ring of the lens, whereas the minimal feature size is the smallest fabricated line width l. The minimum a1 is determined by the minimum feature size l (300 nm).

Figure 2 presents the measured lateral and axial focusing intensity distributions of the experimentally fabricated six-ring ultrathin GO flat lenses with different a1 from 730 to 1,200 nm. It is obvious that well-defined 3D foci are achieved in all cases. The FWHM in both the lateral and axial directions of the foci reduce as the decrease of a1 following exactly the same trend predicted by our theory (Fig. 2c). The highest lateral resolution is 343 nm (0.49λ) and axial resolution is 584 nm (0.83λ) illuminated by a light beam of λ=700 nm, leading to a 3D diffraction-limited resolution of λ3/5. The demonstrated 3D subwavelength focal spot is to our knowledge the best far-field 3D resolution achieved so far for an ultrathin flat lens.

Figure 2: Dependence of focusing resolution on a1.
figure 2

Top: Measured cross-sectional intensity distributions in the lateral (a) and axial (b) directions, respectively. The lens plane is indicated as the white dash line in b. (c) The measured dependence of FWHMlateral (blue spheres) and FWHMaxial (red spheres) on a1. The error bars on both FWHMlateral and FWHMaxial are determined via resampling. The corresponding theoretical predictions (blue and red solid lines) are also presented and show excellent agreements with experimental results. Scale bar, 2 μm.

Focus intensity tuning and large-scale fabrication

The fact that the phase and amplitude modulations of rGO zones are highly dependent on the fabrication laser power offers a simple and unique approach to tune the focus of the ultrathin GO flat lens. To demonstrate this capability, GO lenses with an identical geometry (a1=1.6 μm, N=3) are fabricated with increasing laser powers from 1.5 to 8 μW. The normalized peak focusing intensity is found to increase with the fabrication laser power monotonically, as shown in Fig. 3a. To further illustrate the laser controllable tunability of the focus intensity and large-scale fabrication capability, a large size 200-nm-thick GO thin film is prepared on a glass substrate (Fig. 3b), and a contour of Australian map comprising of 162 GO lenses and a contour of a kangaroo comprising of 89 GO lenses are fabricated with 3 and 7 μW laser powers, respectively, in a 1.5 mm2 area (dashed circle in Fig. 3b), as shown in the microscopic image in Fig. 3c. By shining light on the fabricated region (Fig. 3d), the kangaroo contour shows a much larger brightness compared with the contour of Australian map because of the higher focusing intensity of the comprising lenses (Fig. 3e). The enlarged lateral cross-sectional intensity distributions of GO lenses taken from the kangaroo and Australian map are also shown in the insets of Fig. 3a, which highlight the focusing intensity differences. These results demonstrate not only the versatile laser addressable tunability of the GO flat lenses, but also the scalable and reliable laser fabrication of GO films ready for large-scale manufacturing and applications.

Figure 3: Laser controllable tunability of the focus intensity.
figure 3

(a) Normalized focusing intensities of GO lenses fabricated at increasing laser powers from 1.5 to 8 μW (black spheres). Theoretical prediction (red solid line) shows remarkable agreements with the experimental results. The error bar is based on the statistical average by sampling GO lenses with the identical fabrication power. Insets: Cross-sectional intensity distributions in the lateral direction of two GO lenses fabricated at 3 and 7 μW, respectively. Scale bar, 1 μm. (b) The prepared large size 200-nm-thick GO thin film on a glass substrate. (c) Two sets of GO lenses following the contours of an Australian map and a kangaroo are patterned in a 1.5 mm2 area (dashed circle in b). (d) Schematic diagram of focusing of the GO lenses under a plane wave illumination (white arrows). (e) Measured microscopic images of the contours of an Australian map and a kangaroo at the lens focusing plane showing different focusing intensities. Scale bar, 200 μm.

Broadband focusing

One of the critical challenges of the current ultrathin lenses, such as metalenses2, plasmonic lenses9 and super-oscillation lenses14, is the narrow bandwidth because of the constrains from both the dispersion properties of the constituting materials and the fundamental lens operation mechanisms26. The state-of-the-art bandwidth for an ultrathin flat lens is 280 nm (ref. 27), whereas most ultrathin flat lenses are only functional at a single wavelength2,10. In contrast, the wavefront shaping mechanism of our GO lens originates from the effective amplitude and phase modulations between the GO and rGO zones. Given the fact that the large modulations of both the refractive index and extinction coefficient between GO and rGO are almost maintained from visible (above 350 nm) to infrared regions (Supplementary Figs 3, 4 and 12), the GO lens provides a possible mechanism for achieving broadband focusing. To realize the broadband focusing, we design and fabricate a GO lens (a1=3.3 μm, N=3) with the optimal performance over the entire VIS–NIR wavelength range based on our theoretical guidance (Supplementary Fig. 13). The cross-sectional focusing intensity distributions in the lateral direction of this lens are experimentally characterized under different illumination wavelengths from 450 to 1,500 nm, as shown in the insets of Fig. 4a. The results clearly demonstrate that well-defined subwavelength foci have been achieved at all wavelengths. Moreover, the measured absolute focusing efficiencies, defined as the ratio of the power in the focal region to the total input power before the lens, are higher than 32% on average in the entire measured VIS–NIR band (Fig. 4a), representing a more than 30-time’s improvement compared with the state-of-the-art ultrathin flat lens13. It should be emphasized that such a broadband, high-resolution focusing over 1,050 nm bandwidth is achieved with a single GO lens. In addition, it is interesting to note that the measured FWHMlateral of the foci at different wavelengths shows an almost constant value of 0.8 μm, which is also predictable from the operating principle of our ultrathin lenses.

Figure 4: The broadband focusing capability from visible to infrared.
figure 4

(a) Characterization of the broadband focusing ability of a GO lens (a1=3.3 μm, N=3). Both FWHMlateral (black spheres) and focusing efficiency (red spheres) are measured over a broadband illumination of the wavelengths from 450 to 1,500 nm. The error bars on both FWHMlateral and efficiency are determined via resampling. Theoretical predictions of FWHMlateral (black solid line) are also presented matching well with the experimental results. Insets: Cross-sectional intensity distributions along the lateral direction captured by two CCDs operating at the visible and infrared, respectively. Different colours indicating different input wavelengths schematically. Scale bar, 1 μm. (bd) Focusing ability of a GO lens (a1=1.2 μm, N=3) under a white light illumination. Measured cross-sectional intensity distributions along the lateral (b) and axial (c) directions under a white light illumination. The intensity distributions along the red dashed line in b and blue dashed line in c are also plotted (red and blue spheres in d) and fitted (red and blue solid lines in d), respectively. Scale bar, 1 μm.

The broadband focusing of the GO lens allows us to minimize the chromatic aberration specifically in the visible regime, which is critical for multi-wavelength optical imaging and communications28. By carefully designing a1 of our GO lens, the shift of the focal length over all visible wavelengths can be smaller than the depth of focus, which suggests that our GO lens can be optimized to focus the incoherent white light within a common focal region. To achieve the minimal chromatic aberration in the visible regime, a GO lens (a1=1.2 μm, N=3) has been fabricated and characterized under white light illumination. The cross-sectional intensity distributions in both the lateral and axial directions are shown in Fig. 4b,c, exhibiting an obvious high-resolution 3D focal spot. By fitting the intensity distributions (Fig. 4d), the FWHM in the lateral and axial directions are found to be 577 nm and 1.56 μm, respectively, which even outperform most of the ultrathin lens operating under single wavelength illumination2,10. This is to our knowledge the first demonstration of a subwavelength-scale flat lens that can focus white light three-dimensionally with minimal chromatic aberration.

Mechanical robustness and wavefront engineering of GO films

Given its strong mechanical robustness, large-scale integratability and simultaneous patterning and beam manipulation abilities during the one-step mask-free DLW process, GO ultrathin films have the potential to revolutionize the next-generation integrated optical systems by making miniaturized and fully flexible photonics devices. To illustrate its mechanical flexibility and the versatility in wavefront manipulations, a large-scale GO thin film is integrated on a flexible polydimethylsiloxane substrate and then is patterned with complex logos as well as other functional diffractive optical elements, as shown in Fig. 5. After various bending and twisting (Fig. 5a–d), both the complex logo (Fig. 5e,f) and the 4 × 4 GO lens array (Fig. 5g,h) survive without any compromise of the morphology or optical performance, demonstrating the excellent mechanical strength of our GO films. Moreover, the wavefront shaping mechanism with our ultrathin GO film is readily applicable to other diffractive optical element design. Figure 5i presents a two-dimensional grating, in which rGO (darker regions) and GO offer periodic phase and amplitude modulations. The diffraction pattern under 405, 561 and 633 nm lasers simultaneous illumination clearly shows the high diffraction performance of the GO grating. Finally, as shown in Fig. 5j, a multi-foci GO lens (right column) can be realized by simply adjusting the position of the rGO rings in the single focal spot GO lens design (left column), which offers a great potential for super-resolution imaging, optical data storage and sensing applications29.

Figure 5: The flexibility and versatility in wavefront manipulations.
figure 5

(ad) Bending and twisting of the GO thin film on a polydimethylsiloxane substrate. (e,f) Microscopic images of Swinburne logo fabricated on the GO thin film before (e) and after (f) bending. Scale bar, 10 μm. (g,h) Microscopic images of the GO lens array at the focusing plane before (g) and after (h) bending. Scale bar, 30 μm. Insets: microscopic images of the GO lens array fabricated on the GO film. (i) Diffraction image of a two-dimensional grating, retrieved by using three beams at wavelength of 405, 561 and 633 nm simultaneously. Inset: microscopic image of the two-dimensional grating fabricated on the GO film. Scale bar, 20 μm. (j) Comparison of the single focal spot (left column) and multi-foci GO lenses (right column). The designs are shown schematically (top panel) and the focal intensity distributions along the axial directions are characterized experimentally (bottom panel). Scale bar, 2 μm.

Discussion

The new wavefront shaping concept with laser patterned GO ultrathin films provides new and variable solutions for ultralight weight, highly efficient, highly integratable and flexible optical systems, opening up new avenues for various multidisciplinary applications including non-invasive 3D biomedical imaging, laser tweezing, all-optical broadband photonic chips, light harvesting, aerospace photonics, optical microelectromechanical systems and lab-on-chip devices.

Methods

GO film formation

The high-quality GO solutions are synthesized by the chemical reduction of graphite via a modified Hummers Method30. Then homogenous GO thin films of controllable thicknesses can be prepared by the spraying method with varied spraying speed and solution concentration.

Femtosecond laser fabrication

A homemade DLW system is used to fabricate the GO lens on the as-prepared GO film directly31. To minimize the thermal effect, a low repetition rate, femtosecond pulsed, laser beam (100 fs pulse, 10 kHz, 800 nm) is used to reduce the GO film through the DLW method. The sample is mounted on a 3D nanometric piezo stage. A computer-controlled system is used to create 3D arbitrary structures.

Optical characterization of GO films

A scanning electron microscope (RAITH150-TWO), an optical profiler (Bruker ContourGT InMotion), an ellipsometer (J.A. Woollam M-2000) and a spectrometer (PerkinElmer UV/Vis Spectrophotometer) have been used to characterize the morphology and optical parameters of the GO and rGO films.

Imaging system and broadband measurement

A Nikon N-STORM microscope with 700 nm plane wave illumination is used for the visible lens characterization. The cross-sectional distributions of the generated focal spots of the GO lenses are captured with a numerical aperture=1.4, × 100 objective into a charge-coupled device (CCD) camera. By gradually adjusting the distance between the objective and the GO lens in a step of 100 nm, we are able to obtain the optical intensity distributions at different axial positions, and the 3D images of the focal spot can be reconstructed. By normalizing the sensitivity and exposure time of the CCD camera, the lateral cross-sectional intensity distributions can be captured and the peak focusing intensities can be compared directly.

For the broadband focusing characterization, a super-continuum laser with illumination wavelengths from 450 to 1,500 nm is used as the light source. Two CCD cameras, operating at visible (Watec 902H3 SUPREME) and infrared (Xenics Xeva-1.7-320) regions, respectively, are used to capture the focal plane images of the GO lens at different wavelengths.

Additional information

How to cite this article: Zheng, X. et al. Highly efficient and ultra-broadband graphene oxide ultrathin lenses with three-dimensional subwavelength focusing. Nat. Commun. 6:8433 doi: 10.1038/ncomms9433 (2015).