• Opto-Electronic Science
  • Vol. 2, Issue 12, 230024 (2023)
Xinyu Wang1、†, Xuke Qiu2、†, Mulong Liu3、†, Feng Liu2, Mengmeng Li2, Linpei Xue2, Bohan Chen2, Mingran Zhang4, and Peng Xie2、5、*
Author Affiliations
  • 1School of Future Technology, University of Chinese Academy of Sciences, Beijing 100049, China
  • 2Department of Engineering Science, University of Oxford, Parks Road, Oxford OX1 3PJ, UK
  • 3School of Science, Northwest A&F University, Yangling 712100, China
  • 4Department of Physics, City University of Hong Kong, Kowloon, Hong Kong 999077, SAR, China
  • 5School of Electrical and Electronic Engineering, Nanyang Technological University, Singapore 639798, Singapore
  • show less
    DOI: 10.29026/oes.2023.230024 Cite this Article
    Xinyu Wang, Xuke Qiu, Mulong Liu, Feng Liu, Mengmeng Li, Linpei Xue, Bohan Chen, Mingran Zhang, Peng Xie. Flat soliton microcomb source[J]. Opto-Electronic Science, 2023, 2(12): 230024 Copy Citation Text show less

    Abstract

    Mode-locked microcombs with flat spectral profiles provide the high signal-to-noise ratio and are in high demand for wavelength division multiplexing (WDM)-based applications, particularly in future high-capacity communication and parallel optical computing. Here, we present two solutions to generate local relatively flat spectral profiles. One microcavity with ultra-flat integrated dispersion is pumped to generate one relatively flat single soliton source spanning over 150 nm. Besides, one extraordinary soliton crystal with single vacancy demonstrates the local relatively flat microcomb lines when the inner soliton spacings are slightly irregular. Our work paves a new way for soliton-based applications owing to the relatively flat spectral characteristics.

    Introduction

    Optical microcavity provides the abundant nonlinear investigations owing to the strong field enhancement. The emergence of optical frequency comb in the high-quality microcavity brings to the optical ruler to chip-scale integration and mode-locking pulse with ultrahigh repetition rate. Dissipative Kerr solitons (DKSs) emerged from the balances of cavity dissipation and gain, as well as nonlinearity and dispersion maintain their temporal waveform during propagating in the Kerr medias1-6. This chip-scale optical soliton source is envisioned to play the essential roles in integrated microwave photonics7, high-speed accurate ranging measurements8-10, dual-comb spectroscopy11-13, and optical information processing assisted by wavelength division multiplexing techniques14-18. To date, various techniques have been introduced to generate the coherent soliton microcomb sources, including fast frequency tuning19, power-kicking20, self-injecting solitons21, auxiliary modes22, auxiliary lasers heating the cavity23, photorefractive effects24, and extraordinary soliton states25, 26.

    Optical amplifier is one fundamental path of adapting the light signal power to the demands of engineering and scientific applications. Signal-to-noise ratio (SNR) is one important criterion for exploring the soliton-based applications when external amplifiers are introduced to enhance the microcomb power. The amplifier spontaneous emission (ASE) can significantly impact the noise level of the microcombs and microcomb lines with lower power are challenging to be utilized. One feasible method is focused on increasing the soliton transformation efficiency, e.g. platicon or dark pulse at the normal dispersion regime25, periodic or quasi-periodic arranged soliton sequences termed as soliton crystals (SCs)27, and solitons in coupled dual microcavities28-30. However, these strategies still result in the power attenuation when the microcomb lines are far away from the pump laser, all the microcomb lines are still not efficiently used. An alternative approach is constructing one specifical solion which is not sensitive to the cavity dispersion or soliton source with large-scale flat spectral profiles. The widespread optical source with flat spectral profile is widely prepared in electro-optics modulation combs31, dispersion-less fiber cavities assisted by dispersion optimization32, microcomb coherent combining33, artificial intelligence-driven meta-microcavities34, 35 and micro-cavity with near-zero dispersion characteristics36, 37. Microcavity with dispersion-less characteristics at a larger wavelength domain could facilitate the soliton with a large-scale, relatively flat spectral profile.

    In this paper, we present two schemes for forming a soliton microcomb (SMC) source with the relatively flat spectral envelope. One special microcavity with an ultra-flat integrated dispersion spanning from 1500 nm to 1650 nm is pumped by the counter-propagating dual-color lasers. Single soliton at this dispersion-less regime demonstrates the relatively flat spectral envelope. Besides, microcavity with pure quadratic dispersion also ensures the implementation of local relatively flat spectrum via forming one-defected SC with irregular inter-soliton spacings. Our work not only presents novel soliton dynamics but also offers a high SNR microcomb source, which is crucial for exploring soliton-based applications.

    Numerical explorations and experimental results

    A numerical model suited to describe the microcomb formation can be modelled by the Lugiato-Lefever equation (LLE)38:

    A(Φ,t)t=(γ2+iδ)A(Φ,t)+iD222A(Φ,t)Φ2+ig|A(Φ,t|)2A(Φ,t)+kexPinωp,

    where Φ is the azimuthal domain along with the circumference of the cavity, t is the evolution time of the field envelope. A(Φ,t) is the propagating field envelope circulating in the microcavity. γ, D2/2π and Pin denote as the total dissipation rate, second-order dispersion and on-chip power, respectively. g is single-photon-induced Kerr frequency shift. kex is the external coupling rate between microcavity and waveguide. δ is the cavity detune between the degenerate mode eigenfrequency ω0 and the pump frequency ωp. Figure 1 shows the simulated optical spectrum of single SMC under different dispersion conditions. The simulation parameters for the condition of anomalous dispersion regime are set as: input wavelength λp=1560 nm, second-order dispersion D2/2π=200kHz, soliton repetition rate D1/2π=200GHz, single-photon-induced Kerr frequency shift g=1.2, total quality Q=1.5 million and external quality Qext=6 million, on-chip pump power 150 mW. When the micro-cavity is located at the anomalous dispersion regime (red line in Fig. 1(b)), a bright soliton with a sech² function profile exhibits a prominent power attenuation with a power variation of approximately 10 dB in the 160 nm range, as shown in Fig. 1(a). As shown in Fig. 1(b), by adjusting the local dispersion to 0 (from the −50th resonance mode to 50th resonance mode), a single soliton microcomb with a local relatively flat spectrum (blue line in Fig. 1(a)) is achieved, demonstrating ~2.5 dB power variation at the same 160 nm wavelength range. Furthermore, we also investigate the spectral envelope of single soliton at the normal dispersion regime in Fig. 1(c). The second-order dispersion (D2/2π) is −200 kHz and one frequency shift of −200 MHz on the pump mode is introduced to achieve the local anomalous dispersion. A typical dark platicon is observed with high local power levels (red curve in Fig. 1(c)) but still exhibits the significant power variations in one larger frequency domain. When local integrated dispersion is tuned to zero (from −100th resonance mode to 100th resonance mode as shown in Fig. 1(d)), the stable single soliton (blue line in Fig. 1(c)) has one relatively flat domain with about 320 nm spectrum width. At this region, the maximum power variation is ~6 dB. Meanwhile, temporal waveform of single soliton demonstrates the typical bright pulsed waveform owing to the implementations of large-scale weak anomalous dispersion. The presence of large-scale ultra-flat integrated dispersion promotes the formation of flat solitons, offering the feasibility to achieve the high-SNR microcomb sources when the external amplifier is introduced.

    Optical spectrums of single SMC source and dispersion curves. (a, c) and (b, d) are the simulated optical spectrums of single soliton and dispersion curves, respectively. Blue line and red line in dispersion curves correspond to the conditions of local flat integrated dispersion and initial dispersion. Large-scale flat dispersion reduces the microcomb power variations at the ultra-weak dispersion regime.

    Figure 1.Optical spectrums of single SMC source and dispersion curves. (a, c) and (b, d) are the simulated optical spectrums of single soliton and dispersion curves, respectively. Blue line and red line in dispersion curves correspond to the conditions of local flat integrated dispersion and initial dispersion. Large-scale flat dispersion reduces the microcomb power variations at the ultra-weak dispersion regime.

    To overcome the strong thermo-optic effect, an auxiliary laser is constructed by using a commercial acoustic optical modulator and an external radio-frequency (RF) signal. Both the pump laser and auxiliary laser are amplified by using a polarization-maintaining amplifier and counter-coupled in the micro-resonator. The output signals are monitored by using an optical spectrum analyzer (OSA), an electrical spectrum analyzer (ESA), and an oscilloscope (OSC) (Fig. 2(a)). The core device responsible for microcomb generation is a CMOS-compatible all-pass micro-ring resonator (MRR), fabricated on high-index doped silica glass platform, as is shown in Fig. 2(b). The micro-resonator is coupled with a commercial standard fiber-array with a coupling loss of 3 dB per facet. One thermoelectric cooler (TEC) chip is used for temperature stabilization. The MRR has a radius of 148.1 μm and the free spectral range (FSR) is 190 GHz. The dispersion curve is shown in Fig. 2(c). A specific integrated dispersion curve (green line in Fig. 2(c)) serves as one reference, and over 40 modes are closed to the green reference line, indicating an ultra-flat or ultra-weak integrated dispersion characteristic ranging from 1500 nm to 1650 nm. This ultra-flat integrated dispersion covers the S, C, and L bands, enabling relatively flat SMC generation across a wide spectral range. The pump mode is marked by a red dot in Fig. 2(c) and has a Q-factor of 1.18 million as shown in Fig. 2(d).

    Experimental setup. (a) Experimental setup for the robust single SMC formation. An auxiliary laser is introduced for the stable soliton formation. ECDL, External cavity diode laser; EDFA, Erbium-doped fiber amplifier; OSA, Optical spectrum analyzer; OSC, Oscilloscope; ESA, Electrical spectrum analyzer; MRR, Micro-ring resonator; PD, Photodetector; Cir, Circulator. (b) Microscope image of the high-index doped silica glass micro-ring resonator with a radius of 148.1 μm (lower panel). Butterfly-packaged device (upper panel). (c) Dispersion characteristic of the MRR. The green line (Dint=0) is one referenced integrated dispersion curve. The micro-cavity demonstrates the ultra-flat dispersion characteristic. The red-dot is the soliton mode. (d) The transmission spectra of the soliton mode.

    Figure 2.Experimental setup. (a) Experimental setup for the robust single SMC formation. An auxiliary laser is introduced for the stable soliton formation. ECDL, External cavity diode laser; EDFA, Erbium-doped fiber amplifier; OSA, Optical spectrum analyzer; OSC, Oscilloscope; ESA, Electrical spectrum analyzer; MRR, Micro-ring resonator; PD, Photodetector; Cir, Circulator. (b) Microscope image of the high-index doped silica glass micro-ring resonator with a radius of 148.1 μm (lower panel). Butterfly-packaged device (upper panel). (c) Dispersion characteristic of the MRR. The green line (Dint=0) is one referenced integrated dispersion curve. The micro-cavity demonstrates the ultra-flat dispersion characteristic. The red-dot is the soliton mode. (d) The transmission spectra of the soliton mode.

    To verify the dynamics of the single SMC evolution in the flat integrated dispersion regime and ensure the implementation of thermal balance, the initial modulation frequency on the acoustic optical modulator (AOM) is set to 100 MHz, and the pump laser wavelength is 1564.2 nm. The pump laser and auxiliary laser maintain the identical polarization. Figure 3(a) shows the power trace from the continuous wave (CW) state to the single SMC state. Thermal tuning of the operation temperature of microcavity to make the pump laser entry into the resonance mode and the accumulated intra-cavity power gradually meets the threshold of four-wave mixing (FWM) as shown in Fig. 3(c). The first parametric sidebands (state I) are 25-FSR away from the pump. Continuously decrease the operation temperature of the MRR, soliton steps are observed in Fig. 3(a). Once the soliton existence regime is reached, a typical bi-directional tuning technique is employed: the modulation frequency of the AOM is increased, and the pump laser frequency is decreased by using piezoelectric (PZT) control. Typical soliton switching as shown in Fig. 3(a) leads to the generation of a single SMC (Ⅲ) from the dual-soliton microcomb (Ⅱ) with typical dual-pulse interference characteristics in Fig. 3(c). After obtaining a stable single SMC, the modulation frequency on the AOM is backward tuned to 100 MHz to enhance the stability of micro-soliton, and the beat signal is demonstrated in Fig. 3(b). The ideal generated single SMC exhibits a rectangular spectrum, marked by the green line in Fig. 3(c), and the experimental spectrum shows the trapezoidal profile. The power of the spectrum shows some fluctuations due to irregular mode shifts39. These fluctuations can be mitigated through better optimization methods in the dispersion curve and loss control, resulting in flatter microcombs. Long-term stability of soliton is of vital importance for exploring soliton-based applications. Figure 3(d) shows the long-term power traces by the feedback on the laser frequency and the single SMC has higher stability toward the environment.

    Experimental results at the flat dispersion regime. (a) Power trace of microcomb evolution from the continuous wave (CW) state to single soliton state. MI: Modulation instability. I: Turning pattern. II: Dual-soliton microcomb. III: Single soliton microcomb. (b) Radio frequency spectrum of single soliton state. (c) Optical spectrums of Turning pattern (I), dual-soliton microcomb (II) and single SMC source (III). Via introducing the perfect rectangular spectrum as one standard reference, the single SMC demonstrates the trapezoidal spectrum with weak power-varied spectral profile. (d) Long-term power trace of single SMC.

    Figure 3.Experimental results at the flat dispersion regime. (a) Power trace of microcomb evolution from the continuous wave (CW) state to single soliton state. MI: Modulation instability. I: Turning pattern. II: Dual-soliton microcomb. III: Single soliton microcomb. (b) Radio frequency spectrum of single soliton state. (c) Optical spectrums of Turning pattern (I), dual-soliton microcomb (II) and single SMC source (III). Via introducing the perfect rectangular spectrum as one standard reference, the single SMC demonstrates the trapezoidal spectrum with weak power-varied spectral profile. (d) Long-term power trace of single SMC.

    Experimental results for SC formation with local relatively flat spectrum. (a) Dispersion characteristic of microcavity. (b) Soliton crystal with irregular inter-soliton spacings. Zoom-in local spectrum. At this frequency region, the microcomb demonstrates the relatively flat spectral profile. (c) Temporal evolution of SC with local relatively flat spectral lines.

    Figure 4.Experimental results for SC formation with local relatively flat spectrum. (a) Dispersion characteristic of microcavity. (b) Soliton crystal with irregular inter-soliton spacings. Zoom-in local spectrum. At this frequency region, the microcomb demonstrates the relatively flat spectral profile. (c) Temporal evolution of SC with local relatively flat spectral lines.

    Soliton crystal (SC) is an extraordinary soliton state characterized by collectively ordered soliton sequences, which is known for enhancing the microcomb power due to their high conversion efficiency. The formation of SCs relies on the interference between the modulated background wave and the chaotic modulation instability (MI) field. Intracavity modulated background wave could be constructed by perturbation on the local dispersion. Figure 4 shows the experimental results for SC with local relatively spectral lines. The core device is a high-index doped silica glass MRR, which has a FSR of 48.9 GHz for the fundamental TM mode. The microcavity is designed with anomalous dispersion and has the cross section of 2 μm × 3 μm. Figure 4(a) shows the dispersion curve of TM00 mode and the second-order dispersion (D2/2π) is 141.35 kHz. Typical mode-crossing effect is around 1514.4 nm. In our experiments, the wavelength of the pump laser is 1556.4 nm. The spacing between pump laser and mode-crossing mode is 105 FSR. The on-chip power is about 850 mW. Via controlling the operation temperature of microcavity, various SCs could be formed. SCs exhibit various temporal structures, including periodic arranged solitons termed as perfect soliton crystals, vacancies, superstructure and other complex arranged structure. Among these SC states, a specific defected soliton crystal with slight irregular inter-soliton spacings provides the local flat spectral lines, as shown in Fig. 4(b). In this case, 34 solitons stably exist in the micro-cavity as shown in Fig. 4(c), and the local flat spectral range covers approximately 25 nm. The local perturbations of microcomb power in this type of SC are mainly induced by different dissipation rates and Raman-induced frequency shifts40. This specific SC can be converted from a soliton crystal with one defect by slightly tuning the operation temperature or laser frequency41.

    Conclusion

    In conclusion, we have demonstrated two methods for generating the coherent mode-locked microcomb source with the local relatively flat spectral profile in the micro-ring resonator. One micro-cavity with large-range ultra-flat integrated dispersion is verified and single SMC at this ultra-flat dispersion regime demonstrates relatively flat spectral lines from 1480 nm to 1650 nm, effectively covering the S, C, and L bands. Besides, we explore the generation of local flat microcomb lines using a single-defected soliton crystal with slightly irregular soliton spacing. This WDM-compatible flat soliton microcomb is an excellent candidate as multi-wavelength source for parallel high-capacity informational processing and parallel photonic neuromorphic computing.

    Acknowledgements

    We appreciate the funding support from DreamX International Innovation Team and acknowledge the support from the startup grant from Nanyang Technological University (022527-00001)

    XY Wang, XK Qiu and P Xie conceived this concept. All authors contributed to the manuscript. P Xie supervised the whole project.

    The authors declare no competing financial interests.

    References

    [1] N Akhmediev, A Ankiewicz. Dissipative Solitons(2005).

    [2] L Chang, ST Liu, JE Bowers. Integrated optical frequency comb technologies. Nat Photonics, 16, 95-108(2022).

    [3] AL Gaeta, M Lipson, TJ Kippenberg. Photonic-chip-based frequency combs. Nat Photonics, 13, 158-169(2019).

    [4] TJ Kippenberg, AL Gaeta, M Lipson, ML Gorodetsky. Dissipative Kerr solitons in optical microresonators. Science, 361, eaan8083(2018).

    [5] AFJ Runge, DD Hudson, KKK Tam, CM de Sterke, A Blanco-Redondo. The pure-quartic soliton laser. Nat Photonics, 14, 492-497(2020).

    [6] P Grelu, N Akhmediev. Dissipative solitons for mode-locked lasers. Nat Photonics, 6, 84-92(2012).

    [7] D Marpaung, JP Yao, J Capmany. Integrated microwave photonics. Nat Photonics, 13, 80-90(2019).

    [8] J Riemensberger, A Lukashchuk, M Karpov, WL Weng, E Lucas et al. Massively parallel coherent laser ranging using a soliton microcomb. Nature, 581, 164-170(2020).

    [9] MG Suh, KJ Vahala. Soliton microcomb range measurement. Science, 359, 884-887(2018).

    [10] P Trocha, M Karpov, D Ganin, MHP Pfeiffer, A Kordts et al. Ultrafast optical ranging using microresonator soliton frequency combs. Science, 359, 887-891(2018).

    [11] MG Suh, QF Yang, KY Yang, X Yi, KJ Vahala. Microresonator soliton dual-comb spectroscopy. Science, 354, 600-603(2016).

    [12] A Dutt, C Joshi, XC Ji, J Cardenas, Y Okawachi et al. On-chip dual comb source for spectroscopy. Sci Adv, 4, e1701858(2018).

    [13] BX Xu, XY Fan, S Wang, ZY He. Sub-femtometer-resolution absolute spectroscopy with sweeping electro-optic combs. Opto-Electron Adv, 5, 210023(2022).

    [14] AS Raja, S Lange, M Karpov, K Shi, X Fu et al. Ultrafast optical circuit switching for data centers using integrated soliton microcombs. Nat Commun, 12, 5867(2021).

    [15] J Feldmann, N Youngblood, M Karpov, H Gehring, X Li et al. Parallel convolutional processing using an integrated photonic tensor core. Nature, 589, 52-58(2021).

    [16] XY Xu, MX Tan, B Corcoran, JY Wu, A Boes et al. 11 TOPS photonic convolutional accelerator for optical neural networks. Nature, 589, 44-51(2021).

    [17] XY Wang, P Xie, BH Chen, XC Zhang. Chip-based high-dimensional optical neural network. Nano-Micro Lett, 14, 221(2022).

    [18] A Rizzo, A Novick, V Gopal, BY Kim, XC Ji et al. Massively scalable Kerr comb-driven silicon photonic link. Nat Photonics, 17, 781-790(2023).

    [19] T Herr, V Brasch, JD Jost, CY Wang, NM Kondratiev et al. Temporal solitons in optical microresonators. Nat Photonics, 8, 145-152(2014).

    [20] X Yi, QF Yang, KY Yang, MG Suh, K Vahala. Soliton frequency comb at microwave rates in a high-Q silica microresonator. Optica, 2, 1078-1085(2015).

    [21] BQ Shen, L Chang, JQ Liu, HM Wang, QF Yang et al. Integrated turnkey soliton microcombs. Nature, 582, 365-369(2020).

    [22] Q Li, TC Briles, DA Westly, TE Drake, JR Stone et al. Stably accessing octave-spanning microresonator frequency combs in the soliton regime. Optica, 4, 193-203(2017).

    [23] XY Wang, P Xie, WQ Wang, Y Wang, ZZ Lu et al. Program-controlled single soliton microcomb source. Photonics Res, 9, 66-72(2021).

    [24] Y He, QF Yang, JW Ling, R Luo, HX Liang et al. Self-starting bi-chromatic LiNbO3 soliton microcomb. Optica, 6, 1138-1144(2019).

    [25] XX Xue, Y Xuan, Y Liu, PH Wang, S Chen et al. Mode-locked dark pulse Kerr combs in normal-dispersion microresonators. Nat Photonics, 9, 594-600(2015).

    [26] AW Bruch, XW Liu, Z Gong, JB Surya, M Li et al. Pockels soliton microcomb. Nat Photonics, 15, 21-27(2021).

    [27] DC Cole, ES Lamb, P Del’Haye, SA Diddams, SB Papp. Soliton crystals in Kerr resonators. Nat Photonics, 11, 671-676(2017).

    [28] JMC Boggio, D Bodenmüller, S Ahmed, S Wabnitz, D Modotto et al. Efficient Kerr soliton comb generation in micro-resonator with interferometric back-coupling. Nat Commun, 13, 1292(2022).

    [29] XX Xue, XP Zheng, BK Zhou. Super-efficient temporal solitons in mutually coupled optical cavities. Nat Photonics, 13, 616-622(2019).

    [30] ÓB Helgason, M Girardi, ZC Ye, FC Lei, J Schröder et al. Surpassing the nonlinear conversion efficiency of soliton microcombs. Nat Photonics, 17, 992-999(2023).

    [31] A Parriaux, K Hammani, G Millot. Electro-optic frequency combs. Adv Opt Photonics, 12, 223-287(2020).

    [32] XX Xue, P Grelu, BF Yang, M Wang, SY Li et al. Dispersion-less Kerr solitons in spectrally confined optical cavities. Light Sci Appl, 12, 19(2023).

    [33] BY Kim, Y Okawachi, JK Jang, XC Ji, M Lipson et al. Coherent combining for high-power Kerr combs. Laser Photonics Rev, 17, 2200607(2023).

    [34] Y Meng, YZ Chen, LH Lu, YM Ding, A Cusano et al. Optical meta-waveguides for integrated photonics and beyond. Light Sci Appl, 10, 235(2021).

    [35] E Lucas, SP Yu, TC Briles, DR Carlson, SB Papp. Tailoring microcombs with inverse-designed, meta-dispersion microresonators. Nat Photonics, 17, 943-950(2023).

    [36] ZY Xiao, TY Li, ML Cai, HY Zhang, Y Huang et al. Near-zero-dispersion soliton and broadband modulational instability Kerr microcombs in anomalous dispersion. Light Sci Appl, 12, 33(2023).

    [37] MH Anderson, WL Weng, G Lihachev, A Tikan, JQ Liu et al. Zero dispersion Kerr solitons in optical microresonators. Nat Commun, 13, 4764(2022).

    [38] C Godey, IV Balakireva, A Coillet, YK Chembo. Stability analysis of the spatiotemporal Lugiato–Lefever model for Kerr optical frequency combs in the anomalous and normal dispersion regimes. Phys Rev A, 89, 063814(2014).

    [39] T Herr, V Brasch, JD Jost, I Mirgorodskiy, G Lihachev et al. Mode spectrum and temporal soliton formation in optical microresonators. Phys Rev Lett, 113, 123901(2014).

    [40] ZZ Lu, WQ Wang, WF Zhang, ML Liu, LR Wang et al. Raman self-frequency-shift of soliton crystal in a high index doped silica micro-ring resonator [Invited]. Opt Mater Express, 8, 2662-2669(2018).

    [41] WQ Wang, ZZ Lu, WF Zhang, ST Chu, BE Little et al. Robust soliton crystals in a thermally controlled microresonator. Opt Lett, 43, 2002-2005(2018).

    Xinyu Wang, Xuke Qiu, Mulong Liu, Feng Liu, Mengmeng Li, Linpei Xue, Bohan Chen, Mingran Zhang, Peng Xie. Flat soliton microcomb source[J]. Opto-Electronic Science, 2023, 2(12): 230024
    Download Citation