Stable mid-infrared polarization imaging based on quasi-2D tellurium at room temperature

Next-generation polarized mid-infrared imaging systems generally requires miniaturization, integration, flexibility, good workability at room temperature and in severe environments, etc. Emerging two-dimensional materials provide another route to meet these demands, due to the ease of integrating on complex structures, their native in-plane anisotropy crystal structure for high polarization photosensitivity, and strong quantum confinement for excellent photodetecting performances at room temperature. However, polarized infrared imaging under scattering based on 2D materials has yet to be realized. Here we report the systematic investigation of polarized infrared imaging for a designed target obscured by scattering media using an anisotropic tellurium photodetector. Broadband sensitive photoresponse is realized at room temperature, with excellent stability without degradation under ambient atmospheric conditions. Significantly, a large anisotropic ratio of tellurium ensures polarized imaging in a scattering environment, with the degree of linear polarization over 0.8, opening up possibilities for developing next-generation polarized mid-infrared imaging technology.

P olarized infrared imaging systems (PIRIS) are of great importance for the development of astronomy, biological diagnosis, national security, and industrial applications [1][2][3][4] , which require precise recognition of a target object in the presence of clutter. At present, the trend of modernizing polarized infrared imaging applications increasingly requires convenience, miniaturization, good workability at room temperature and severe environments, and larger pixel density for future devices [5][6][7] . Recently, the emerging researches into two-dimensional (2D) materials have potential to provide another routine for nextgeneration photodetectors. The inherent characteristics of strong light-matter interaction and low dark current in 2D materials can provide high photoresponsivity and low noise at room temperature [8][9][10][11][12][13][14] . The lack of dangling bonds on the surface of 2D materials can reduce surface-induced performance degradation, facilitating easy integration with other materials or substrates 15 . More importantly, some 2D materials with intrinsic in-plane anisotropy structure exhibit the potential to realize polarization photodetection without polarization filters [16][17][18] . Hence, it is a feasible approach to develop state-of-the-art PIRIS devices by implementing in-plane anisotropic 2D materials. Although many 2D materials have been used to successfully fabricate polarization photodetectors, there are still some limitations of these 2D materials which hinder their practical applications, such as instability under ambient conditions, narrow photoresponse range, and insufficient polarized sensitivity 9,[19][20][21][22][23] . Besides, there are reports of 2D material mid-infrared (MIR) polarization photodetectors 24,25 , but the works about polarization imaging under scattering based on 2D materials have not yet been realized, addressing the feasibility of 2D materials-based photodetectors to further achieve contrast-enhanced polarized infrared imaging.
Here, we have systematically investigated polarized infrared imaging applications based on quasi-2D tellurium (Te) photodetector. The narrow bandgap of quasi-2D Te provides strong light-matter interactions in the spectral range from visible to MIR (∼3.0 μm), and the high carrier mobility and small effective carrier mass provide the potential to achieve high photodetecting performance [26][27][28] . Furthermore, the anisotropic crystal structure of quasi-2D Te gives rise to strong anisotropic photoresponse, explained by the energy band structure in theory 27,[29][30][31] . Detailed Raman spectrum characterizations have been performed to verify the high quality of quasi-2D Te, carrier mobility measured by field-effect transistor (FET)-based devices can reach up to ∼9 × 10 2 cm 2 V −1 s −1 at room temperature, which is at the forefront in the field of 2D materials. The wide spectral photoresponse of the Te photodetector is confirmed from 0.52 to 3.0 μm, and the most sensitive photoresponsivity (R) under 1.06 μm illumination is ∼1.36 × 10 3 AW −1 with corresponding measured detectivity (D M *) of ∼1.15 × 10 10 Jones. The photoresponsivity also remains higher than ∼3.53 × 10 2 AW −1 and the measured detectivity higher than ∼3.01 × 10 9 Jones, when illuminated by the incident light with a wavelength of 3.0 μm. Due to the unique asymmetric crystal structure, the polarized photoresponse is also demonstrated, with the photoresponsivity anisotropic ratio of ∼8 under the incident light with a wavelength of 2.3 μm, which is sufficient for potential polarized imaging. Importantly, polarized infrared imaging tests under scattering environments have also been carried out for the Te-based photodetectors, without the integration of complex polarization filters. Under strongly scattering environment, achieving target photoelectric imaging is hard for 2H-phase MoTe 2 -based photodetectors which have no polarized photoresponse, however, the quasi-2D Te-based photodetector reported here can realize imaging for the target with the degree of linear polarization (DoLP) over 0.8 at the wavelength of 2.3 μm. Furthermore, we have also implemented polarized imaging for incident light over a wide wavelength range in more detail. Significantly, the feasibility of implementing Te for polarized imaging proved in this work is at the first stage, for the next step, the integration for authentic facilities in the future requires efforts upon achieving the design and fabrication of quasi-2D Te-based area-array, which worth more investigations.

Results
Anisotropic characterization of Te. The quasi-2D Te crystals (Supplementary Note 1, Supplementary Fig. 1a) are synthesized through a solution method 27,32 , the crystal quality is verified through Raman spectrum tests ( Supplementary Fig. 1b-h). The thickness of Te flake measured by atomic force microscope (AFM) is ∼27.5 nm ( Supplementary Fig. 1i-j), and its thicker thickness is beneficial to reach stronger absorption, leading to higher photocurrent, but the thickness should not be too thick, otherwise the overlarge dark current will lead to poor photodetectivity and other performance. Besides, the thicker Te flake possesses narrower bandgap 8 , which is beneficial to extend the detectable bandwidth to longer infrared wavelength. Quasi-2D Te consists of a helical chain in the [0001] direction, with adjacent helical chains stacked together through a non-covalent weak bonding [26][27][28]32 , forming an anisotropic crystal structure, as shown in Fig. 1a. To investigate its electrical properties in more detail, we have fabricated a Te-based FET device with five parallel electrodes marked by numbers from 1 to 5 as shown in the schematic of the device in Fig. 1b, resulting in different channel lengths between electrodes, and the channel length between two electrodes is measured through optical images as shown in Fig. 1c. Figure 1d shows the transfer curves for the devices with different channel lengths at room temperature, with the drain bias fixed at 1.0 V and the gate bias is varied from −10.0 V to 10.0 V. The on-off ratio is slightly decreased in shorter channel length, which is due to slightly weaker gate control over the drain current. High mobility of ∼9 × 10 2 cm 2 V −1 s −1 is calculated from the transfer curves, which verifies the potential to achieve high photodetecting performance (more detailed analysis in Supplementary Notes 2 and 3, Supplementary Table 1, Supplementary Figs. 2 and 3, 25). The anisotropic crystal structure of Te causes a certain degree of intrinsic anisotropy for the electrical performances of Te-based FET (optical image is shown in Supplementary Fig. 1k), with the electrical conductivity anisotropic ratio (along x-and y-axis) of ∼2.24 under gate bias of 0 V (inset of Fig. 1e). From the transfer curves (Fig. 1e), we have calculated the on-off ratio and fieldeffect mobility for both x-axis and y-axis directions respectively, as shown in Fig. 1f, similarly, the field-effect mobility anisotropic ratio is ∼1.69 (Fig. 1g). The polarized Raman spectra of Te is measured as shown in Fig. 1h-k. Based on the results, the intensities of both E 1 -TO peak and A 1 peak show two-lobe shape, while the intensity of E 2 peak shows a four-lobe shape ( Fig. 1i-k). To more intuitively present the high anisotropy in Raman spectra, we have defined Raman anisotropic ratio as the maximum intensity of Raman peak over the minimum, producing Raman anisotropic ratios of ∼13, ∼4, and ∼12.5 for E 1 -TO, A 1 , and E 2 mode, respectively, in experiments. In theory, quasi-2D Te crystal belongs to D 3 4 space group, the anisotropic Raman intensity can be fitted by Raman tensor 33,34 : absorption edge extended to more than 3.0 μm. We have also measured the polarized absorption spectra with polarized illumination along the y-and x-axis direction, respectively. From the results, the absorption along the x-axis (green curve) is higher than that along the y-axis (blue curve) for the measured wavelength range, which is relevant to the anisotropic band structure of the Te Brillouin zone ( Supplementary Fig. 5). The anisotropic ratio of absorption increases with the increase of absorption wavelength ( Supplementary Fig. 6) 27,29 . Notably, in the range of absorption wavelength from 2.5 to 3.0 μm, the y-axis direction absorption is very low. Therefore, the absorption anisotropic ratio can reach up to ∼10 for the absorption range from 2.3 to 3.0 μm, promising excellent potential in polarized MIR photodetectors.
Broadband stable optoelectronic performance of Te. The photoresponse performances are measured between the electrode 2 and 4 in Fig. 1c, with the device channel length of ∼6 μm and channel width of ∼2.5 μm, located at the center of the Te flake.
The gate bias is fixed at 0 V, the drain bias is varied from −1.0 to 1.0 V, and the incident lights are unpolarized with the output power from 0 to 6.0 mW. As shown in Fig Supplementary Fig. 7k-n). The drain bias is fixed at 1.0 V, and the gate bias is fixed at 0 V, the frequency of the noise current is chosen to be 1 kHz, the corresponding noise current density is 4:55 10 À11 A ffiffiffiffiffiffiffiffiffiffi ffi Hz À1 p at room temperature, and the calculation details are described in the Supplementary Notes 5. The lower illumination power always leads to higher R and D M *, and the evolution of R and D M * under fixed laser power is in accordance with absorption results. Due to effective absorption of Te at 3.0 μm, high R of ∼3.54 × 10 2 AW −1 and D M * of ∼3.01 × 10 9 Jones are realized under the illumination at the wavelength of 3.0 μm, which is competitive among various 2D materialsbased photodetectors (Supplementary Table 2 Fig. 1i-j) also contributes to such high performance due to its stronger absorption character, as the transmitted light to the substrate is limited. The high performance is also related to a large photoconductive gain, which is assisted by the fast carrier transition time due to the high mobility of Te and the short channel length of ∼6.0 μm. In addition, the defects located inside the Te crystal can contribute to the photoconductive gain through shallow states of photo-induced carriers at the band tail 44,45 , which is a fast process to provide a large gain. The existence of photogain is revealed from the EQE values in Supplementary Fig. 8a-c. Besides, we have also measured the transfer curves under 1.06 μm illumination and 1.0 V drain bias at room temperature, as shown in Fig. 2g, the dark state is shown in Fig. 2h, and the net photocurrent is summarized in Fig. 2i, the maximum net photocurrent is realized at the neutral point with ∼5.9 V gate bias, also suggesting that the photoresponses of the device are not dominated by the trap-induced photogating effect. To further show the origin of photogain, we have measured the photocurrent distribution in the Te channel with another device (Fig. 2j), the photocurrent is generated at the channel/electrode interface instead of the whole channel, excluding the photogating as a dominant factor for the high photogain. The response time is measured at 1.0 V drain bias, 0 V gate bias, and 6.0 mW laser output power for several selected wavelengths, as summarized in Supplementary Fig. 9, the rise time is at the range of 48.7-62.7 μs, and the fall time is at the range of 62.7-78.0 μs. The carrier lifetime is extracted from transient photoluminescence decay curves of Te, as shown in Supplementary Fig. 10. Based on carrier lifetime and transit time, a photogain (G) at the level of ∼2.3 × 10 3 can be calculated, indicating the high performance in the Tebased device. Besides, to make it more convenient to compare the device performances with previous works, the calculated detectivity (D C *) is also summarized in Supplementary Fig. 11. It is still worth noting that the calculated detectivity (D C *) can always overestimate the actual detectivity of 2D materials-based devices, as the dark current is lower than the actual noise current. Moreover, we have also investigated the device performances under low-temperature conditions, to more completely analyze its photodetecting performances. The drain bias is also varied from −1.0 to 1.0 V, and the gate bias is fixed at 0 V. At lower temperature (77 K), the photodetecting performances can be further enhanced due to lower dark current, as shown in Fig. 2c, d (all photocurrent curves are summarized in Supplementary  Fig. 12a-j). The dependences of R and D M * on the wavelength and laser power of the incident light, are summarized in Fig. 2k, l, respectively (scatter diagrams are included in Supplementary  Fig. 12k-n, external quantum efficiency (EQE) is shown in Supplementary Fig. 8d-f). Here the results show similar trends to those measured at 300 K (Fig. 2e, f). At 77 K, a higher R can reach up to ∼8.36 × 10 2 AW −1 and D M * is at a higher level of ∼9.04 × 10 9 Jones, under the illumination wavelength of 3.0 μm. Moreover, device stability is very important for 2D materials-based MIR photodetectors. Although some 2D materials-based photodetectors are showing high MIR photoresponse, they cannot realize long-lasting stable photodetecting under ambient environment, such as easily oxidized black phosphorus and black phosphorus-arsenic. Here, the Te photodetector reveals high stability under ambient conditions, and the device has been measured after 3 months, exhibiting unnoticeable degradation of performance ( Fig. 2m and Supplementary Fig. 13). The corresponding R and D M * in Fig. 2n, o are calculated from the shaded rectangle region in Fig. 2m, which also indicate the stable high performance of the device.
Polarized photoresponses of Te. The polarized photoresponses of the device have also been investigated under linear polarized illuminations. The drain bias is varied from −1.0 to 1.0 V, the gate bias is fixed at 0 V, the laser power is fixed at 6.0 mW, the laser polarization direction is fixed and initially parallel to the longer side of the device (y-axis which is defined as 0°), and the device is rotated with a step of 10°to measure its polarized photoresponse. The polarized photocurrents under the illuminations with the wavelengths of 0.52, 0.637, 0.785, 1.55, and 2.3 μm are presented in Supplementary Fig. 14, and the corresponding line shape curves are included in the Supplementary  Figs. 15-19. To better illustrate the photoresponse anisotropy, here we define the net polarization current ΔI p : The net polarization photocurrents under the illuminations at 0.52, 1.55, and 2.3 μm are shown in Fig. 3a-c, respectively, and the results measured with other wavelengths are shown in Supplementary Fig. 20. Strong polarization photodetection is clearly observed to demonstrate that the Te photodetector can identify the slight differences between different polarized illuminations. The photocurrents under 1.0 V drain bias for different incident polarization angles are shown in Fig. 3d-f, corresponding to the incident wavelengths of 0.52, 1.55, and 2.3 μm, respectively. Even at 0°and 180°polarization direction where the absorption is relatively weaker, obvious photoresponse is still observed, which also ensures the high performances of the device for polarization applications. The visual polarized photoresponsivities are presented in polar diagrams, as shown in   26,27 . Moreover, although black phosphorus and black phosphorus-arsenic may realize similar or higher anisotropic ratio 49,50 , stable Te with a strong anisotropic ratio is more appealing for future polarization applications when taking the easy oxidization nature of these materials at ambient conditions into account.   Supplementary Fig. 21b-c, when the device is under different orientation, to confirm the strongly polarized photosensitivity of the device. Based on Supplementary  Fig. 21, the polarization extinction ratio of Te-based device under 2.3 μm illumination is reaching ∼8, which is acceptable in the field of 2D materials. Based on absorption results, it can reach even higher extinction ratio at a longer wavelength, making it a potential material for polarization applications in the MIR range. Under stronger optical scattering conditions, such as foggy and cloudy weather, colloid and emulsion media, the incident light is severely scattered to obstruct imaging based on non-polarizing facilities 26,52 . To test the performance in scattering conditions, we have carried out imaging tests for the Te photodetector, without the integration of a complicated polarization filter, and the scattering media is placed on the surface of the target. Since no polarization filter is needed, the fabricated device can be used directly to get polarized images, while the traditional method by rotating the polarization filter is not applicable. According to the imaging method as shown the schematic in the inset of Fig. 4a, the scattering signal illuminated onto the photodetector is approximately linear polarized (light scattering mechanism in Supplementary Note 6), the polarized photocurrents along different orientations are collected by rotating the photodetector. The signal processing procedure calculates the DoLP (Fig. 4b), defined by Eqs. (8)-(11) where I θ x; y ð Þ refers to the polarized intensity detected at θ angle (more details in Supplementary Note 7) 6,51,52 . For the conventional division-of-focal-plane polarimeter (DoFP) structures, it requires at least four pixels to acquire the first three components of the Stokes vector to calculate one DoLP data ( Supplementary  Fig. 22) 5,6,52 , which means that to get a m × n DoLP polarization image, it requires a 2m × 2n device matrix. On the contrary, only a single-pixel Te photodetector is adequate to calculate a DoLP data, which means that it requires only a m × n device matrix to get a m × n DoLP polarization image. It suggests that the Tebased device needs fewer pixels to realize the same resolution for DoLP, making it easier to further promote the imaging resolution. To better demonstrate the advantage of polarized imaging based on anisotropic Te flakes, a non-polarization 2H-phase MoTe 2 photodetector is used to image the HUST and SITP target sheltered by the scattering media. Figure 4c shows the imaging results for the MoTe 2 photodetector with hardly distinguishable DoLP. However, the imaging results for the Te photodetector show much better contrast for the target after calculating the DoLP, as shown in Fig. 4d, e. At 2.3 μm illumination, the larger absorption anisotropic ratio can realize clearer polarization imaging, comparing the results with 1.55 μm illumination (Fig. 4f), demonstrating that the Te photodetector is promising for polarization imaging applications.

Discussion
Imaging applications based on the DoFP structure are applicable to all regions from visible to LWIR, which has attracted significant attention. At present, most DoFP structures require a layer of polarization filter and a layer of photo-sensing array 5 , making it requires four single-point devices with different polarization filtering directions to detect polarized signals along four different directions for the calculation of one pixel in DoLP. Theoretically, wire grid polarizers with matured fabrication technology can also reach a high polarization extinction ratio as only light polarized parallel to the grid can pass through the polarizer. As for 2D materials including Te, their intrinsic polarized photosensitivity can provide another route to realize polarization functionalities. For example, the intrinsic polarized photosensitivity of Te promotes polarization imaging without polarization filters. The photocurrent in one single-point device along different polarization orientations can be collected by rotating the Te device, for the calculation of DoLP. However, it is worth noting that the Te-based devices and many other 2D materials-based polarization devices are still at the infancy with some existing challenges at present. Currently, the fabrication of large-scale device arrays is a great challenge in the field of 2D materials. With the development of matured methods for producing large films, the Te photodetector array could be one promising option to realize the practical polarization imaging. This work serves as the first step toward the development of polarization imaging arrays based on 2D materials, which is of great significance.
It is also worth noting that a large photoresponse anisotropy of ∼8 is realized under 2.3 μm illumination, and the anisotropic ratio can be further enhanced for longer wavelength according to the absorption results, promising Te as an ideal material for NIR and MIR polarized imaging devices to reach a higher DoLP. As one of the 2D materials, Te can be integrated with on-chip structures easily, such as waveguides, photonic crystals, and metasurfaces, to further enhance optoelectronic responses and increase the photoresponse anisotropy, which is worth more investigation. 2D Te is also applicable for flexible devices with polarization sensitivity, promising the high level of integration and flexibility 53 . In addition, most DoFP structure lacks the capability to detect circularly polarized signals, which is the fourth element of a Stokes vector: RH means right-handed light and LH means left-handed light, thus, to integrate Te with chiral meta-structures and make full use of its strong spin-orbit coupling, may be a promising path to fully understand the imaging of circularly polarized signals [54][55][56][57] .
In conclusion, we have systematically studied the photodetecting properties of quasi-2D tellurium photodetectors and used the photodetector to successfully realize polarized infrared imaging under scattering. High photoresponsivity and detectivity were achieved for illumination with the wavelengths from 0.52 to 3.0 μm. Due to its anisotropic crystal structure, the Te photodetector shows a photoresponsivity anisotropic ratio reaching up to ∼8 under 2.3 μm illumination and can be even higher than ∼10 under the MIR range illumination. Based on the anisotropic photoresponses, polarized imaging capability is preserved for quasi-2D Te-based photodetector under scattering environments, with the DoLP reaches over 0.8 at the wavelength of 2.3 μm. As a rising star in the field of 2D materials, after overcoming the challenge of controllable grown of large-scale tellurium film, Te photodetector is promising for the next-generation imaging facilities, which is significant for the development of national security, remote sensing, biological diagnosis, and optical engineering.

Methods
Solution synthesis method. The crystals are fabricated through a solution synthesis method based on ref. 27 .
Raman and AFM measurements. Raman spectra is measured by LABRAM ARAMIS system (Horiba), the incident laser is 532 nm Nd:YAG laser, the scattering signal is reflected to the detector with 1/10 attenuation.
The polarized Raman spectra are measured by Alpha 300R confocal Raman imaging system (WITec), with 532 nm laser. The polarization directions of the incident laser and the scattering signal are modulated to be parallel, then the sample is rotated in a step of 10°to measure the polarized Raman signal. AFM data are measured by using Digital TNSTRUments/Veeco Dimension 3000 system in the tapping mode of 0.5 Hz. indicating that MoTe 2 device shows no polarization detection ability, which is not applicable for polarized imaging under scattering environments. d, e Imaging DoLP results for the Te device under 1.55 and 2.3 μm illumination, respectively. The imaging contrast is strong for the target due to its intrinsic excellent polarization detection ability. Scale bars in (c)-(e), 1 mm. f Normalized imaging contrast of S 0 and DoLP. For MoTe 2 , the contrast becomes lower due to its isotropic photoresponse along any direction, and for Te device, the contrast is higher after calculating DoLP, and DoLP enhancement under 2.3 μm illumination is higher than 1.55 μm illumination due to stronger absorption anisotropy.
intensity can be determined to be the x-direction. Based on the calculation results, in Fig. 1d, we have defined the x-direction to be 0°and rotated the pattern.
Absorption measurements. The absorption spectra is measured at room temperature with confocal microscopes. The broadband laser is focused on the sample. At first, the total reflection r total on silver mirror is measured to confirm the profile of the laser source. Then the reflected signal on a bare silicon wafer r SiO 2 is collected, finally, the reflected signal on the tellurium sample r Te is collected. Based on the following equations: where α SiO 2 is the absorption of the silicon wafer, α Te is the absorption of the Te sample, and the absorption rate of the Te sample can be calculated: For polarized absorption measurements, a half-wave plate is added before the light source, the polarization direction is modulated to be along (y-direction, 0°) and perpendicular (x-direction, 90°) to the longer side of tellurium sample, the polarized absorption spectra is calculated by using the same method mentioned above. Electronic and optoelectronic measurements. The unpolarized current data are measured by using a probe station system (SEMISHARE SE-6), a semiconductor characterization system (KEYSIGHT B1500A) and a home-made photocurrent mapping station at atmospheric pressure conditions. The photocurrent is measured in the range of 0.52-3.0 μm. In the range of 0.52-1.55 μm, the selected singlewavelength laser sources include 0.52 μm (Thorlabs, LP520-SF15), 0.637 μm (Thorlabs, LP637-SF70), 0.785 μm (Thorlabs, LP785-SF20), 0.83 μm (Thorlabs, LP830-SF30), 0.94 μm (Thorlabs, LP940-SF30), 1.06 μm (Thorlabs, LP1060-FC), 1.31 μm (Thorlabs, LPSC1310-FC), and 1.55 μm (Thorlabs, LPSC1550-FC). For 2.3 and 3.0 μm illumination, the output laser wavelength is acquired by filtering a supercontinuum broadband white laser source. The photocurrent from 0.52 to 0.94 μm are collected through the semiconductor characterization system, the laser is directly illuminated on the device without focusing, the diameter of the laser spot is determined through optical images, which is 1000 μm in radius. The photocurrent from 1.06 to 3.0 μm are collected through the home-made photocurrent mapping station, the laser is focused on the sample by a ×20 objective, so the radius is much smaller at 50 μm. In addition, the output laser power illuminated on the device throughout the mapping station is also decreased due to losses through the light path and lenses. The realistic light power illuminated on the device location is calibrated by using power meters (Thorlabs, S130C, S132C, S175C), the light power density is similar for photocurrent measurements (0.32, 0.64, 0.95, 1.27, 1.59, and 1.91 nW μm −2 for 1.0, 2.0, 3.0, 4.0, 5.0, and 6.0 mW output power of the laser source, respectively).
Polarized photocurrent measurements. The polarized photocurrent is measured by using a home-made photocurrent mapping station, we have selected five wavelengths in our measurements based on the operating wavelength range of the half-wave plate, 0.52 μm (Thorlabs, LP520-SF15), 0.637 μm (Thorlabs, LP637-SF70), 0.785 μm (Thorlabs, LP785-SF20), 1.55 μm (Thorlabs, LPSC1550-FC), and 2.3 μm (filtered through a supercontinuum broadband white laser source). The laser spot is focused on the sample by a ×20 objective with the laser spot of 50 μm in radius from optical images collected through CCD camera, a half-wave plate is placed between the objective and our sample, the laser polarization direction is fixed, and initially parallel to the y-axis of our device, then the polarized photocurrent is acquired by rotating the device with a step of 10°.
Response speed measurements. The chopped laser is controlled by the laser diode/temperature controller (Thorlabs ITC4020). The photocurrent is measured by an oscilloscope for response speed analysis.
Noise spectral density measurements. The noise spectral density is measured by using a spectrum analyzer (SR770).
Imaging application measurements. The imaging measurements are carried out through a home-made imaging system. The light source is Fianium/NKT Photonics LLFT Contrast supercontinuum white laser source, the output wavelength is between 1.0 and 2.4 μm. For unpolarized imaging tests, we have selected four wavelengths of 1.06, 1.550, 2.0, and 2.3 μm. For polarized imaging tests, we have measured 1.55 and 2.3 μm by adding a half-wave plate between the target and the device to acquire linear polarized illumination, the laser polarization direction is fixed, and initially parallel to the y-axis of our device, the polarized imaging results are collected by rotating the device with a step of 10°. For scattering tests, the halfwave plate is removed from the light path, photocurrent at different orientations is also measured by rotating the device. One resistance is cascaded with the device, and the photocurrent changes are measured through changes of bias for the series resistance. The device location is controlled by two x-axis and y-axis step motor, by changing the device location, photocurrent signal is acquired for each pixel. Details for the measurements in Fig. 4 are also provided in the Supplementary Information.

Data availability
The datasets for this study are available from the corresponding author on reasonable request.

Code availability
The simulation codes for this study are available from the corresponding authors on reasonable request.