IR PHOTODETECTOR WITH GRAPHENE AND PHASE CHANGE LAYERS AND RELATED METHODS
20230082228 · 2023-03-16
Inventors
Cpc classification
H01L31/1013
ELECTRICITY
G01J5/0853
PHYSICS
International classification
Abstract
An IR photodetector includes an electrically conductive layer, a first dielectric layer over the electrically conductive layer, and a phase change material layer over the first dielectric layer. The IR photodetector further includes first and second electrically conductive contacts coupled to the phase change material layer, and a graphene layer over the phase change material layer and having a perforated pattern therein. The IR photodetector includes circuitry configured to apply a bias voltage between the first and second electrically conductive contacts, and detect a sensing current in the phase change material layer caused by IR radiation received by the graphene layer, the IR radiation having a frequency range based upon the bias voltage.
Claims
1. An infrared (IR) photodetector comprising: an electrically conductive layer; a first dielectric layer over the electrically conductive layer; a phase change material layer over the first dielectric layer; first and second electrically conductive contacts coupled to the phase change material layer; a graphene layer over the phase change material layer and having a perforated pattern therein; and circuitry configured to apply a bias voltage between the first and second electrically conductive contacts, and detect a sensing current in the phase change material layer caused by IR radiation received by the graphene layer, the IR radiation having a frequency range based upon the bias voltage.
2. The IR photodetector of claim 1 wherein the phase change material layer has a thickness gradient.
3. The IR photodetector of claim 1 wherein the bias voltage comprises a pulse train.
4. The IR photodetector of claim 3 wherein the pulse train comprises a multi-level pulse train.
5. The IR photodetector of claim 1 wherein the phase change material layer comprises vanadium oxide.
6. The IR photodetector of claim 1 further comprising: a transparent electrically conductive layer over the phase change material layer; a second dielectric layer over the phase change material layer; and a protective layer over the graphene layer.
7. The IR photodetector of claim 6 wherein the first dielectric layer comprises a polymer layer; and wherein the second dielectric layer comprises silicon nitride.
8. The IR photodetector of claim 1 wherein the perforated pattern comprises an array of elliptical holes.
9. The IR photodetector of claim 1 wherein the graphene layer is configured to receive at least one of mid-wavelength IR (MWIR) radiation and long wavelength IR (LWIR).
10. The IR photodetector of claim 1 wherein the electrically conductive layer comprises at least one of gold, silver, and platinum.
11. An infrared (IR) photodetector comprising: an electrically conductive layer; a first dielectric layer over the electrically conductive layer; a phase change material layer over the first dielectric layer and having a thickness gradient; first and second electrically conductive contacts coupled to the phase change material layer; a graphene layer over the phase change material layer and having a perforated pattern therein; and circuitry configured to apply a pulse train bias voltage between the first and second electrically conductive contacts, and detect a sensing current in the phase change material layer caused by IR radiation received by the graphene layer, the IR radiation having a frequency range based upon the pulse train bias voltage.
12. The IR photodetector of claim 11 wherein the pulse train bias voltage comprises a multi-level pulse train.
13. The IR photodetector of claim 11 wherein the phase change material layer comprises vanadium oxide.
14. The IR photodetector of claim 11 further comprising: a transparent electrically conductive layer over the phase change material layer; a second dielectric layer over the phase change material layer; and a protective layer over the graphene layer.
15. The IR photodetector of claim 14 wherein the first dielectric layer comprises a polymer layer; and wherein the second dielectric layer comprises silicon nitride.
16. The IR photodetector of claim 11 wherein the perforated pattern comprises an array of elliptical holes.
17. The IR photodetector of claim 11 wherein the graphene layer is configured to receive at least one of mid-wavelength IR (MWIR) radiation and long wavelength IR (LWIR).
18. A method of making an infrared (IR) photodetector, the method comprising: forming a first dielectric layer over an electrically conductive layer; forming a phase change material layer over the first dielectric layer; forming first and second electrically conductive contacts coupled to the phase change material layer; forming a graphene layer over the phase change material layer and having a perforated pattern therein; and coupling circuitry to apply a bias voltage between the first and second electrically conductive contacts, and detect a sensing current in the phase change material layer caused by IR radiation received by the graphene layer, the IR radiation having a frequency range based upon the bias voltage.
19. The method of claim 18 wherein the phase change material layer has a thickness gradient.
20. The method of claim 18 wherein the phase change material layer comprises vanadium oxide.
21. The method of claim 18 further comprising: forming a transparent electrically conductive layer over the phase change material layer; forming a second dielectric layer over the phase change material layer; and forming a protective layer over the graphene layer.
Description
BRIEF DESCRIPTION OF THE DRAWINGS
[0009]
[0010]
[0011]
[0012]
[0013]
[0014]
[0015]
[0016]
[0017]
[0018]
[0019]
[0020]
[0021]
[0022]
DETAILED DESCRIPTION
[0023] The present disclosure will now be described more fully hereinafter with reference to the accompanying drawings, in which several embodiments of the invention are shown. This present disclosure may, however, be embodied in many different forms and should not be construed as limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the present disclosure to those skilled in the art. Like numbers refer to like elements throughout, and base 100 reference numerals are used to indicate similar elements in alternative embodiments.
[0024]
[0025]
[0026]
[0027]
[0028]
[0029]
[0030] Vanadium oxide (VO.sub.2) is considered as one of the standard materials for building microbolometers with broadband MWIR photodetection because it features a reversible IMT when heated above the phase transition temperature T.sub.c, which is slightly above and close to room temperature [2]. Bulk VO.sub.2 undergoes a phase transition from an insulating state with monoclinic crystal structure below 68° C. to a metallic state with rutile crystal structure above T, =68° C. (=341 K) [3-5]. This phase transition is fully reversible with a hysteresis loop, occurs on a sub picosecond timescale [6, 7], and can be initiated either thermally, electrically [8], or optically [9]. While for bulk VO.sub.2 optically induced IMT can only be achieved by pumping above the band gap of E.sub.g=670 meV, it is possible to induce IMT in thin films at energies of 200 meV and above (corresponding to wavelengths of Λ=6.2 μm and below) due to electronic defects inside the band gap [10]. This is the reason why VO.sub.2 thin films containing defects can be used for IR detection in a wide IR range, in particular in the 3-5 μm range. Interestingly, when VO.sub.2 is in the form of a thin film, its transition temperature T.sub.c depends strongly on the thickness d of the film, i.e. T.sub.c decreases from 65° C. (=338 K) for d=25 nm down to 52° C. (=325 K) for d=3 nm [11]. While IR radiation with wavelengths above about Λ=1 μm cannot detect the change in thickness of around 20 nm, it certainly distinguishes between the insulating and the metallic phase of VO.sub.2. All these properties make VO.sub.2 the ideal material for developing MWIR photodetectors based on the IMT effect.
[0031] However, photodetection of MWIR light with wavelength above about 6 μm is inefficient with VO.sub.2 bulk or thin films because of the relatively weak interaction between the incident photons and the optical phonons in VO.sub.2. This interaction is so weak that the IMT cannot be achieved. That is why in the wavelength regime of 8 to 12 μm the photodetection is based purely on the bolometric effect in the semiconducting phase of VO.sub.2 [13, 14]. Advantageously, by adding a single layer of nanopatterned graphene (NPG) on top of a layer of VO.sub.2, an IR photodetector 100, 200 is disclosed that not only greatly enhances the absorption of MWIR light energy in the longer wavelength regime from Λ=6 μm and exceeding 12 μm but also narrows the absorption bandwidth to 0.1 μm within the MWIR range of 3 to 12 μm, thereby enabling plasmonically enhanced spectrally selective absorption of MWIR light for the IMT effect in a heterostructure made of NPG and VO.sub.2.
[0032] Described herein is a model of a MWIR microbolometer that includes a hexagonal boron nitride (h-BN) coated NPG, silicon nitride (Si.sub.3N.sub.4), indium tin oxide (ITO), VO.sub.2, polymer, and gold (Au) mirror, as shown in the IR photodetector 100 in
[0033] For the electronic response of the graphene sheet and the VO.sub.2 layer to the incident MWIR photons, the intraband conductivity of graphene and the dielectric function of VO.sub.2 in the insulating and metallic regimes is considered. Using the linear dispersion relation, the intraband optical conductivity of graphene is [12, 18]
which in the case of E.sub.F »kBT is reduced to
where τ is determined by impurity scattering and electron-phonon interaction τ.sup.−1=τ.sub.imp.sup.−1+τ.sub.e−ph.sup.−1. Using the mobility μ of the NPG sheet, it can be presented in the form τ.sup.−1=ev.sup.2F/(μE.sub.F), where V.sub.F=10.sup.6 m/s is the Fermi velocity in graphene.
ω.sub.P=√{square root over (e.sup.2E.sub.F/2ε.sub.m)}
is the bulk graphene plasma frequency.
Following Refs. 19 and 13, the dielectric function of VO.sub.2 can be modeled by means of a Drude Lorentz model for the insulating phase,
ϵ.sub.i=ϵ∞ (3)
where j=1, . . . , N denotes the N phonon modes, ω.sub.j are the phonon vibration frequencies, y.sub.j are the scattering rates, and S.sub.j are the oscillation strengths. The value of these parameters can be found in Ref. 13. The metallic phase of VO.sub.2 can be modeled by means of the Drude formula,
where Ω.sub.p=N.sub.f e.sup.2/.sub.0m*=8000 cm.sup.−1 is the plasma frequency, Γ=e/m*,μ.sub.vo2=10000 cm.sup.−1 is the collision frequency, with μ.sub.vo2=2 cm.sup.2/Vs being the mobility, N.sub.f=1.3×10.sup.22 cm.sup.−3 the free-carrier concentration, and m*=2m.sub.e the effective mass of the charge carriers. m.sub.e is the free electron mass.
[0034] The finite-difference time domain method (FDTD) is used to calculate the absorbance of the hybrid IR photodetector 100, 200 as shown in diagram 1000 of
[0035] The absorbance resonance peak as a function of wavelength can be tuned by means of the Fermi energy of NPG, as shown in diagram 1010 of
[0036] For modeling the operation of the IR photodetector 100, 200, COMSOL is used, and the following theory for the thermoelectric properties of VO.sub.2 close to the phase transition temperature T.sub.c. The VO.sub.2 layer is operated around the IMT phase transition temperature T.sub.c. The performance of the bolometric detection can be analyzed by means of the heat equation and a hysteresis model [21]. The heat equation reads
[0037] where C is the heat capacity, a is the absorbance, P is the power of the incident radiation, I is the time-independent bias current, R(T) is the temperature-dependent resistance, G is the thermal conductivity of the heat sink, and T.sub.h is the time independent temperature of the heat sink. The hysteretic behavior of R(T) for VO.sub.2 layer can be calculated by
where the semiconductor volume fraction is given by
[0038] where w is the width of the hysteresis, β is a function of dg/dT at T.sub.c, P (x) is an arbitrary monotonically decreasing function, and δ=sign (dT/dt). The proximity temperature is given by
Equations (5)-(8) describe the hysteretic behavior of the plasmonically driven bolometric photodetector. The incident power P.sub.inc is then given by the energy pumped into the plasmonic nanostructure.
[0039] Using this thermoelectric theory and combining it with our FDTD results, a photothermoelectric theory of the NPG-VO.sub.2 heterostructure for the IR photodetector 100, 200 is provided. The geometry of the IR photodetector 100, 200 is optimized in terms of performance. A channel width of 5 μm (distance between the source and drain contacts) and a channel length of 5 μm are chosen. The larger the channel length, the smaller is the resistance of the VO.sub.2 channel, resulting in reduced Joule heating. This is the method used to reduce the resistance instead of increasing the thickness of the VO.sub.2 layer, which would increase the absorption of the VO.sub.2 layer in the metallic phase. For the initial study, a thickness of 3 nm is used. Note that the current flows only through VO.sub.2, not through NPG. NPG is used only as a photothermoelectric heating element. When the incident MWIR light is off, a very weak dark current I on the scale of μA is flowing when a bias voltage V.sub.b is applied in the range from 0.1 to 1.7 V. In stark contrast, when the incident MWIR light is on, a much larger light current I on the scale of mA is flowing with the same applied bias voltage V.sub.b. This effect is due to the phase transition of VO.sub.2 between insulating and metallic phases. Using the advantage of this effect, the model of an ultrasensitive photodetector based on the NPG-VO.sub.2 heterostructure is developed.
[0040] After modeling the heating and cooling of the NPG-VO.sub.2 heterostructure as a function of time, the optimum photodetection process may be identified. Diagrams 1040, 1050 in
[0041] Diagram 1050 shows an enlarged portion of diagram 1080, which comprises a magnification between times t=0 ms and t=3 ms in order to demonstrate that the temperature increase of VO.sub.2 with pristine graphene or without graphene is not sufficient to trigger an IMT in VO.sub.2. This means that NPG may be necessary for triggering the IMT in VO.sub.2. The photocurrent I.sub.ph through the VO.sub.2 layer as a function of temperature T for a constant applied bias voltage V.sub.b=1.7 V exhibits the expected hysteresis loop as shown in diagram 1080 of
[0042] Since it is desirable to realize a linear dependence of the photocurrent I.sub.ph as a function of input power Pi. of the MWIR light, a gradient in the thickness of the VO.sub.2 layer is added, as shown in the IR photodetector 200 of
[0043] The responsivity R of the IR photodetector 100, 200 can be calculated by means of the formula [16]
where R is the resistance of VO.sub.2 in the metallic phase, I.sub.light is the light current when the incident light is on, I.sub.dark is the dark current when the incident light is off, and P.sub.inc is the power of the incident light. The responsivity as a function of Fermi energy E.sub.F of NPG is shown in diagram 1170 of
where
v.sub.n=√{square root over (v.sub.t.sup.2+v.sub.b.sup.2+v.sub.f.sup.2)}, (11)
is the root-mean-square of the total noise voltage, which consists of the sum over all possible noise voltages, such as the thermal Johnson-Nyquist noise v.sub.t, due to thermal motion of the charge carriers and independent of the bias voltage V.sub.b, the shot noise V.sub.b, due to the discrete nature of uncorrelated charge carriers, and the 1/f noise v.sub.f, also called flicker noise, due to random resistance fluctuations. The Johnson noise is given by [26]
v.sub.t+√{square root over (4k.sub.BTR)}, (12)
where k.sub.B is the Boltzmann constant, T is the temperature, and R is the resistance. The shot noise is given by [26]
v.sub.b+√{square root over (2eI.sub.dR.sup.2)}, (13)
where e is the elementary charge and Id is the dark current. Since the dark current is very low and the IR photodetector 100, 200 operates close to room temperature, the shot noise is much smaller than the Johnson noise. Therefore, it is possible to safely neglect the shot noise. At a modulation frequency of V.sub.b of around 1 kHz, it is possible to also neglect the 1=f noise. Using the NEP, the detectivity of the IR photodetector 100, 200 is calculated by means of the formula [16]
where A is the area of the photodetector. The results of these figures of merit are shown in Table I.
TABLE-US-00001 TABLE I Figures of Merit of the NPG-VO.sub.2 gradient photodetector at a Fermi energy of E.sub.F = 1.0 E.sub.v ΔT [K] P [μW] T.sub.0 [K] [V/W] NEP [fW/√{square root over (Hz)}] D* [Jones] 5.0 18.0 315.0 1.1 × 10.sup.4 347 0.50 × 10.sup.10 3.0 14.0 317.0 1.4 × 10.sup.4 273 0.64 × 10.sup.10 1.2 10.6 318.8 1.9 × 10.sup.4 203 0.85 × 10.sup.10
[0044] Table I shows that the detectivity D* of the IR photodetector 100, 200 operating close to room temperature is close to D* of cryogenically cooled HgCdTe photodetectors. The IR photodetector 100, 200 reaches a sensitivity close to VO.sub.2 microbolometers while exhibiting a shorter detection time of around 1 ms and being able to detect photons also in the LWIR regime, which is impossible for VO.sub.2 microbolometers.
[0045] In conclusion, the model of an ultrasensitive MWIR photodetector based on a heterostructure made of NPG and VO.sub.2 is provided, thereby extending the responsivity of a VO.sub.2 microbolometer to the LWIR domain. Moreover, this hybrid IR photodetector 100, 200 has a narrowband absorption in the MWIR and LWIR that can be tuned by means of a gate voltage. Our results show that the IR photodetector 100, 200 can reach a large responsivity R˜10.sup.4 V/W, a detectivity D˜10.sup.10 Jones, and a sensitivity in terms of NEP, NEP˜10 fW/sqrt(Hz) close to room temperature by taking advantage of the phase change of a thin VO.sub.2 layer. The NPG sheet achieves an absorption of nearly 100% due to LSPs around the patterned circular holes in a hexagonal lattice symmetry. The electrostatic gate potential can be used to tune the wavelength peak in the MWIR and LWIR regimes between 3 and 12 microns, thereby overcoming the intrinsic upper limit of 6 microns for microbolometers based on VO.sub.2. COMSOL simulations show that the IR photodetector 100, 200 is able to operate on a time scale of 1 ms, much shorter than the response times of current microbolometers based on VO.sub.2 alone. The IR photodetector 100, 200 reaches detectivities of cryogenically cooled HgCdTe photodetectors and sensitivities close to and field of view similar to VO.sub.2 microbolometers while operating close to room temperature.
[0046] Referring now to
[0047] The IR photodetector 100 illustratively includes first and second electrically conductive contacts 107a-107b coupled to the phase change material layer 105. Each of the first and second electrically conductive contacts 107a-107b may comprise one or more of gold, silver, platinum, copper, and aluminum.
[0048] The phase change material layer 105 may comprise vanadium oxide (VO.sub.X), for example. Of course, other phase change materials may be used, but the vanadium oxide embodiment is helpful due to the room temperature phase change property. In some embodiments, the phase change material layer 105 may have a thickness gradient, which is shown with the dashed line in
[0049] The IR photodetector 100 illustratively includes a transparent electrically conductive layer 106 over the phase change material layer 105. The transparent electrically conductive layer 106 may comprise indium tin oxide, for example. In some embodiments, such as shown in
[0050] The IR photodetector 100 includes a second dielectric layer 110 over the phase change material layer. The second dielectric layer 110 may comprise silicon nitride, for example.
[0051] The IR photodetector 100 includes a graphene layer 111 over the phase change material layer 105 and having a perforated pattern 112 therein. In the illustrated embodiment, the perforated pattern comprises an array of elliptical holes (e.g. 5×8 array of circle-shaped holes (
[0052] The IR photodetector 100 illustratively comprises circuitry 102 configured to apply a bias voltage (V.sub.g) between the first and second electrically conductive contacts 107a-107b. The bias voltage may comprise a pulse train or a multi-level pulse train. As discussed hereinabove, the multi-level pulse train may comprise a pulse train as depicted in
[0053] The circuitry 102 is configured to detect a sensing current in the phase change material layer 105 caused by IR radiation received by the graphene layer 111. In particular, the circuitry 102 is configured to apply a bias voltage (V.sub.b) to the phase change material layer 105, and detect a current through the phase change material layer due to phase changes.
[0054] As will be appreciated, the spectral sensitivity of the IR photodetector 100 is tuned by the pulse train. In other words, the sensed IR radiation has a frequency range based upon the bias voltage V.sub.g. In some embodiments, the graphene layer 111 may be configured to receive at least one of MWIR radiation and LWIR.
[0055] Another aspect is directed to a method of making an IR photodetector 100. The method comprises forming a first dielectric layer 104 over an electrically conductive layer 103, and forming a phase change material layer 105 over the first dielectric layer. The method further comprises forming first and second electrically conductive contacts 107a-107b coupled to the phase change material layer 105, and forming a graphene layer 111 over the phase change material layer and having a perforated pattern 112 therein. The method also includes coupling circuitry 102 to apply a bias voltage between the first and second electrically conductive contacts 107a-107b, and detect a sensing current in the phase change material layer 105 caused by IR radiation received by the graphene layer 111, the IR radiation having a frequency range based upon the bias voltage.
[0056] Referring now additionally to
[0057] Many modifications and other embodiments of the present disclosure will come to the mind of one skilled in the art having the benefit of the teachings presented in the foregoing descriptions and the associated drawings. Therefore, it is understood that the present disclosure is not to be limited to the specific embodiments disclosed, and that modifications and embodiments are intended to be included within the scope of the appended claims.
REFERENCES (INCORPORATED BY REFERENCE IN THEIR ENTIRETY)
[0058] [1] A. Rogalski, Infrared Physics & Technology 43, 187 (2002). [0059] [2] A. Rogalski, Infrared Physics & Technology 54, 136 (2011). [0060] [3] J. Goodenou, Journal of Solid State Chemistry 3, 490 (1971). [0061] [4] F. J. Morin, Physical Review Letters 3, 34 (1959). [0062] [5] S. Chandra, D. Franklin, J. Cozart, A. Safaei, and D. Chanda, ACS Photonics 5, 4513 (2018). [0063] [6] K. Appavoo and J. Haglund, Richard F., Nano Letters 11, 1025 (2011). [0064] [7] M. J. Dicken, K. Aydin, I. M. Pryce, L. A. Sweatlock, E. M. Boyd, S. Walavalkar, J. Ma, and H. A. Atwater, Optics Express 17, 18330 (2009). [0065] [8] G. Stefanovich, A. Pergament, and D. Stefanovich, Journal of Physics-Condensed Matter 12, 8837 (2000). [0066] [9] A. Cavalleri, C. Toth, C. W. Siders, J. A. Squier, F. Raksi, P. Forget, and J. C. Kieer, Physical Review Letters 87, 237401 (2001). [0067] [10] M. Rini, Z. Hao, R. W. Schoenlein, C. Giannetti, F. Parmigiani, S. Fourmaux, J. C. Kieer, A. Fujimori, M. Onoda, S. Wall, and A. Cavalleri, Applied Physics Letters 92, 181904 (2008). [0068] [11] G. Xu, P. Jin, M. Tazawa, and K. Yoshimura, Applied Surface Science 244, 449 (2005). [0069] [12] A. Safaei, S. Chandra, A. VŽzquez-Guardado, J. Calderon, D. Franklin, L. Tetard, L. Zhai, M. N. Leuenberger, and D. Chanda, Physical Review B 96, 165431 (2017). [0070] [13] A. S. Barker, H. W. Verleur, and H. J. Guggenheim, Phys. Rev. Lett. 17, 1286 (1966). [0071] [14] C. Chen, X. Yi, J. Zhang, and B. Xiong, International Journal of Infrared and Millimeter Waves 22, 53 (2001). [0072] [15] A. Safaei, S. Chandra, M. N. Leuenberger, and 8 D. Chanda, Acs Nano 13, 421 (2019). [0073] [16] A. Safaei, S. Chandra, M. W. Shabbir, M. N. Leuenberger, and D. Chanda, Nature Communications 10, 3498 (2019). [0074] [17] M. W. Shabbir and M. N. Leuenberger, Scientic Reports 10, 17540 (2020). [0075] [18] H. P. Paudel, A. Safaei, and M. N. Leuenberger, in Nanoplasmonics—Fundamentals and Applications, edited by G. Barbillon (Intech, London, 2017) Chap. 3, p. 1142. [0076] [19] Z. Wang, P. Li, Y. Chen, J. Liu, F. Qi, H. Tian, B. Zheng, and J. Zhou, Applied Surface Science 307, 712 (2014). [0077] [20] P. U. Jepsen, B. M. Fischer, A. Thoman, H. Helm, J. Y. Suh, R. Lopez, and R. F. Haglund, Phys. Rev. B 74, 205103 (2006). [0078] [21] L. A. L. de Almeida, G. S. Deep, A. M. N. Lima, I. A. Khrebtov, V. G. Malyarov, and H. Ne, Applied Physics Letters 85, 3605 (2004). [0079] [22] M. M. Qazilbash, M. Brehm, G. O. Andreev, A. Frenzel, P. C. Ho, B.-G. Chae, B.-J. Kim, S. J. Yun, H.-T. Kim, A. V. Balatsky, O. G. Shpyrko, M. B. Maple, F. Keilmann, and D. N. Basov, Physical Review B 79, 075107 (2009). [0080] [23] X. Zhong, X. Zhang, A. Gupta, and P. LeClair, Journal of Applied Physics 110, 084516 (2011). [0081] [24] S. Samanta, A. K. Raychaudhuri, X. Zhong, and A. Gupta, Physical Review B 92, 195125 (2015). [0082] [25] J. Ordonez-Miranda, Y. Ezzahri, K. Joulain, J. Drevillon, and J. J. Alvarado-Gil, Physical Review B 98, 075144 (2018). [0083] [26] W. Guo, Z. Dong, Y. Xu, C. Liu, D. Wei, L. Zhang, X. Shi, C. Guo, H. Xu, G. Chen, L. Wang, K. Zhang, X. Chen, and W. Lu, Advanced Science 7, 1902699 (2020).