This work focuses on the optimization of coupling coefficient (k) of the inductive link for the wireless power transfer (WPT) system to be used in implantable medical devices (IMDs) of centimeter size. The analytic expression of k is presented. Simulations are conducted by using the high-frequency structure simulator (HFSS). Analytic results are verified with simulations. The receiving (Rx) coil is implanted in the body and set as a circular coil with a radius of 5 millimeters for reducing the risk of tissue inflammation. The inductive link under misalignment scenarios is optimized to improve k. When the distance between the transmitting (Tx) and Rx coils is fixed at 20 mm, it is found that, to maximize k, the Tx coil in a planar spiral configuration with an average radius of 20 mm is preferred, and the Rx coil in a solenoid configuration with a wire pitch of 0.7 mm is recommended. Based on these optimization results, an inductive link WPT system is proposed; the coupling coefficient k, the power transfer efficiency (PTE), and the maximum power delivered to the load (MPDL) of the system are obtained with both simulation and experiment. Different media of air, muscle, and bone separating the Tx and Rx coils are tested. For the muscle (bone) medium, PTE is 44.14% (43.07%) and MPDL is 145.38 mW (128.13 mW), respectively.

1. Introduction

Implantable medical devices (IMDs) have been widely used in the clinical society, such as cardiac pacemakers, cochlear implants, left vagus nerve stimulators, spinal cord stimulators, retinal implants, and deep brain stimulators [19]. At present, most IMDs are powered by batteries for the convenience and comfort of patients. The IMD has to be replaced with a new one by surgery when the battery is out of energy, which may be quite painful for the patients. Wireless power transfer (WPT) has emerged as a promising solution for these drawbacks, enable minimally invasive long-term implantation of those devices, and eliminate the need for another surgery for IMD replacement.

There are several techniques to realize WPT, such as ultrasound, optical, microwave radiation, electric field coupling, and inductive coupling [1018]. Among these methods, inductive coupling is the most preferred for IMDs, where power is transmitted from one coil and received by the other (usually integrated with the IMD inside a human body) through electromagnetic induction. To enhance the power transfer efficiency (PTE), the inductive coupling coefficient (k) between the transmitting (Tx) and receiving (Rx) coils needs to be maximized. As a key parameter of the WPT system, k can be obtained with the relationship , where L1 and L2 are the self-inductances of the Tx and Rx coils, respectively, and M is the mutual inductance, which typically depends on the geometries of the two coils, their distance, and the material separating them.

Existing works focused on the mutual inductance for the coils in different geometries, such as a solenoid, thin disk, filamentary circular, and spiral rectangular coils [1930]; as well as in different relative positions, such as coaxial, noncoaxial, and inclined axial [3133]. In addition, some works focused on calculating the self-inductance of the coils in different geometries, such as square, hexagonal, octagonal, and circular planar [34,35]. Analytical expressions for self-inductance were obtained in [36,37]. The self-inductance and mutual inductance of the printed square-shaped spiral coils were explored, which were used to design an inductive WPT system [38].

To avoid or eliminate the side effects of large-sized IMDs, such as inflammation, cell death, and users’ pain, the size of the IMD has been continuously scaling down [17,18]. With the rapid development in integrated circuit (IC) technologies, IMDs for cardiac pacemakers or nerve stimulators can likely be miniaturized to a centimeter (cm) scale [42]. Therefore, the dimension of the Rx coil is limited, and PTE and the maximum power delivered to the load (MPDL) are reduced. New WPT systems have to be designed for cm-sized IMDs to enhance the coupling coefficient, which plays a critical role in improving the PTE and MPDL.

An LCC-C compensated WPT system for cardiac pacemakers was designed in [39], and the size of the designed system to be embedded was quite large, 31.4 mm × 47.4 mm × 8.1 mm, which was against the spirit of minimally invasive applications. The distance between the Tx and Rx coils was only 8 mm, which may be too short for real-world usage. A sandwiched WPT system was designed for recharging the battery of a micromedical robot [40]. The system is quite complicated and is inconvenient for the patient to use.

For midrange (∼2 m) inductive link power transmission, it is necessary to transmit power at resonance to reduce the dissipated energy in the objects between the Tx and Rx coils and improve the transmission efficiency [41]. For the short-range (∼2 cm) power transmission, power losses in the Tx and Rx coils become significant. Resonance is generally employed to reduce power and tissue losses and improve the PTE [42, 43].

In this paper, we aim to optimize the coupling coefficient of an inductive link, where the operating frequency, the configurations, and sizes of the Tx and Rx coils are selected to maximize k. After that, we propose a WPT system with the optimized link for cm-sized IMDs. We consider two types of candidate Tx coils, namely, single-turn circular coil and planar spiral circular coil. We examine the mutual inductances of Tx-Rx, and in turn, the coupling coefficients, with the Tx coil in various sizes through extensive simulations, under different Tx-Rx coil distances, operating frequency, wire pitch of the Rx coil, and turn number of the Tx coil. According to the optimized results, for the specific applications such as a cardiac pacemaker, the planar spiral coil is selected as the Tx coil for the WPT system, along with other optimized parameters for the Tx and Rx coil configurations. We proceed to investigating the PTE and MPDL of the planar spiral Tx coil when the Tx and Rx coils are separated with different media, including air, muscle, and bone. At last, we carry out experiments to corroborate the simulation results and validate the efficiency of the proposed WPT system in real-world applications. Such WPT systems can largely extend the lifespan of cm-sized IMDs and greatly benefit the patients.

We focus on the power transmission for IMDs, and data communication is out of the scope of this work and will be an interesting research direction for future works.

The novelties of the present work are summarized as follows:(1)The optimizations of the coupling coefficient of an inductive link WPT system for cm-sized IMDs are presented, with the misalignment scenarios of the inductive link coils included.(2)Based on the optimizations of k, an inductive link WPT system with the Rx coil in cm size is proposed, which can be applied to recharge the button cell in a cm-sized IMD.(3)Simulation results of k, PTE, and MPDL of the proposed inductive link WPT system are verified against experiments for different media separating the Tx and Rx coils, for example, air, muscle, and bone.(4)The deembedding calibration method is adopted in experiments. When the operation frequency is 25 MHz, for the muscle (bone) medium, a PTE of 44.14% (43.07%) and an MPDL of 145.38 mW (128.13 mW) are demonstrated with measurements.

The rest of the paper is arranged as follows: in Section 2, the theoretical expressions of k, PTE, and MPDL for an inductive link WPT system are presented. In Section 3, extensive simulations are provided. Different inductive link models are optimized by simulations, and an inductive link WPT system is proposed. In Section 4, experimental systems of the WPT links are set up, and the measured results are obtained by the using deembedding calibration method. Discussions are provided in Section 5, and concluding remarks are presented in Section 6.

2. Analytical Model

2.1. Single-Turn Circular Tx Coil

For the Tx coil, a single-turn coil can be used to reduce the specific absorption rate (SAR) and increase MPDL [17]. As a multiturn coil can be equaled to a set of single-turn coils, we only deal with a single-turn Tx coil in this work. Figure 1(a) shows the model of a coaxial WPT system with a single-turn circular Tx coil. The inductive link consists of a Tx coil and an Rx coil. For an implantable WPT system, the Rx coil should be small sized as aforementioned, so R2 is limited to 5 mm in this paper. The Tx coil is placed outside the body to be excited to generate magnetic fluxes. A portion of the fluxes is picked up by the implanted Rx coil. In this way, power can be transmitted wirelessly, overcoming the obstacles of the human body. The tangential component of the vector magnetic potential at the location of the Rx coil can be expressed as [31, 44]where h is the spacing between the two coils, R1 and R2 are the radius of Tx and Rx coils, respectively, and I1 is the excitation current in the Tx coil.is the second kind of complete elliptic integral,is the first kind of complete elliptic integral, andis the modulus of the elliptic integral. µ0 is the permeability in the vacuum.

The magnetic flux excited on the Rx coil by the Tx coil can be obtained as [31]

Then, the mutual inductance, M, can be written as [45, 46]

The self-inductance of the filamentary circular coil, Lfcc, is given by [47]where R is the radius of the coil, for the Tx and Rx coils, it is R1 and R2, respectively. r is that of the wire, and the constant p depends on the distribution of the current at the cross-section of the wire: p = 1 when the current is uniformly distributed, and p = 0 when there is an obvious electric current skin effect. Under high operating frequencies (>10 MHz), the skin effect will be large and p = 0. O(…) stands for a small term, when , the term μ0O(r2/R) is negligible, then we have

With equations (6) and (8), we can obtain the coupling coefficient of the single-turn circular Tx and Rx coils case as follows:where

2.2. Planar Spiral Circular Tx Coil

Usually, a planar spiral Tx coil can be used to enhance the coupling coefficient [38]; therefore, we deal with this type of configuration for the Tx coil. Figure 1(b) shows the WPT system established in HFSS, and the Tx coil is configured as a planar spiral. In the model, Rin and Rout are the inner and outer radii of the Tx coil, respectively. The average radius Rav equals to 0.5(Rout + Rin), b equals to 0.5(Rout − Rin), and n1 is the turn number of the coil. Then, the mutual inductance of the Tx coil with a single-turn circular Rx coil, Mp, can be obtained as [31,33]where n = 1, 2, 3, … is a positive integer, and a larger value of n leads to a more precise result of Mp.where

The self-inductance of a planar spiral coil, Lpsc, can be written as [36]where . The k of the planar spiral circular Tx coil with the aforementioned Rx coil kps can be obtained aswhere Mp can be obtained with equation (11), and L1 and L2 can be obtained with equations (8) and (15), respectively.

2.3. PTE and MPDL

Figure 2 shows the illustration of the circuit model for an inductive link WPT system. RTx, RRx, and RL are the resistances of the Tx coil, the Rx coil, and the load, respectively. C1 and C2 are the capacitors. It should be noted that the values of RTx and RRx are frequency-dependent, owing to the electric current skin effect under high frequency. The PTE, denoted with η, is defined as the ratio of the power received by the load, PL, over the power input to the Tx coil, Pin, and at resonance η can be expressed as [18, 48]where , , and are the quality factors of the Tx, the Rx, and the load, respectively. ω is the angular frequency.

The PDL can be obtained as follows [18]:where Rref is the reflected impedance to account for the effect of the Rx coil on the Tx coil at resonance and is expressed as follows:

It should be noted that PDL increases with the square of I1, but on the other hand, the SAR increases also with the square of I1. The local SAR limit for a human body is 2 W/kg [49]. Typically, energy absorption occurs on the tissue surface, and the peak value of SAR, denoted with SARp, is assumed to be the maximum radiation value when the Tx coil is driven by a constant current, I1 = 1 A in this work. Therefore, the maximum current in the Tx coil within the SAR limit can be expressed as [50]

Then, the MPDL is

With equations (17) and (21), PTE and MPDL can be obtained.

3. Simulation Results

3.1. Single-Turn Circular Tx Coil

It is said that, for the cm-sized IMDs, the operating frequency of a few tens of MHz is good for both PTE and SAR-constrained PDL. If the frequency is raised, the PTE would be enhanced; but on the other hand, the SAR-constrained PDL would be depressed, and vice versa [18]. Therefore, the operating frequency range in this work is from 10 MHz to 60 MHz.

Figure 3(a) shows M versus R1 of the Tx coil with different values of h, and R2 is assumed to be 5 mm. It is noted in the figure that when the value of h is smaller, which means Tx is closer to Rx, Rx will receive more magnetic fluxes emitted from the Tx coil, leading to a larger value of M. Figure 3(b) illustrates k versus R1 at the operation frequency fp = 30 MHz. It is observed that, as the value of h increases, the maximized k can be obtained by selecting a proper value of R1. Moreover, when h is large (>12 mm), the effect of the variation of R1 on k will be small. If the distance between the Tx and Rx coils is set as 20 mm, the maximum k of 0.013 can be obtained when R1 = 20 mm.

The effect of the Tx coil deformed from a circular circle to two ellipses is investigated. For an ellipse, the aspect ratio is defined as the ratio of the major axis length over minor axis length, and the aspect ratios of the two ellipses are set to be 1.5 and 2, respectively. The radius of the circular coil is 20 mm, and the circumference of every coil is 126 mm. The distance between the Tx and Rx coils is 20 mm. Figure 4 shows the simulation results of value k between one of the three shapes of a Tx coil and an Rx coil, respectively. It is observed that, k will decrease with the increase of the aspect ratio. For instance, when fp = 25 MHz, k is 0.0121 for the circular circle of the Tx coil, k drops to 0.0120 and 0.0113 for the elliptic coils of the aspect ratio of 1.5 and 2, respectively. Therefore, a deformed circular Tx coil leads to a smaller value of k.

Figure 5 shows the simulation results of k with multiturn Rx coils. The turn number of the Rx coil, n2, is swept from 2 to 6, and the wire pitch of the coil, s, is swept from 0.1 mm to 1.1 mm. It is revealed that the coupling coefficient k increases with n2. On the other hand, as s increases, k first increases and then decreases, and peaks around s = 0.7 mm. The Tx coil is a single-turn circular coil with R1 = 20 mm, and r1 = 0.287 mm.

3.2. Planar Spiral Circular Tx Coil

For the planar spiral circular Tx coil case, we can obtain the mutual inductance between the Tx and Rx coils with equation (11), and the self-inductance of the Tx coil with equation(15). After that, with the use of equation(16), we can obtain k.

Figure 6 depicts the k value versus the turn number of the Tx coil, n1, with different b at fp = 30 MHz. Simulations show that, k increases slightly with n1, and/or with b. The reason is that, when Rav is fixed, as b increases, Rout will be increased while Rin will be reduced; then a few more magnetic fluxes emitted from the Tx coil will be collected by the Rx coil, as such, k increases slightly with b. Therefore, with a fixed value of Rav, a large b of the Tx coil is preferred.

3.3. Optimizations of Tx and Rx Coils for the Misalignment Scenarios

The above researches are based on the premise that the Tx and Rx coils are perfectly aligned, and they are coaxial. However, in practice, misalignment may happen in two ways as illustrated in Figure 7: (1) horizontal distance deviation, where the two coils are not coaxial, and the distance between the two axes of the coils is d0, and (2) vertical angle deviation (θ) of the Rx coil, where the two coils are not parallel to each other. When d0 = 0 and θ = 0, the two coils are perfectly aligned. We use HFSS to simulate the inductive link of the two misalignment cases. The parameters optimized in Section 3.2 are adopted, i.e., Rav = 20 mm, n1 = 7, b = 6 mm, r1 = 0.287 mm, R2 = 5 mm, r2 = 0.04 mm, h = 20 mm, and fp = 30 MHz.

Figure 8 shows the simulated k of the inductive links when the horizontal deviation distance (d0) is swept from 0 mm to 20 mm, and the vertical angle deviation (θ) of the Rx coil is swept from 0° to 90°. It can be seen that the value of k will be greatly reduced by increasing either d0 or θ. However, when there is a horizontal distance deviation, k can be improved by introducing a proper vertical angle deviation intentionally. The reason is that, when the vertical angle deviation is introduced, the two coils are adjusted to be coaxial, and the Rx coil will pick up more magnetic fluxes than before. Similarly, if the Rx coil has a vertical angle deviation, k can be improved by introducing a proper horizontal distance intentionally.

Table 1 shows the simulated values of PTE and MPDL versus vertical angle deviation (θ), when there is no horizontal distance deviation, d0 = 0 mm. It can be seen that the values of PTE and MPDL will decrease with θ. When θ = 60°, PTE will decrease from 12.29% to 6.49%, and MPDL will drop from 27.72 mW to 4.33 mW.

It is noted that observed from Figure Rong et al. [51] have investigated the optimization design of resonance coils for drone wireless charging system, where the horizontal misalignment between the Tx and Rx coils may happen. For the first time, the authors proposed an unequally spaced Tx coil and a varying radius size per turn Rx coil to well address the misalignment problem. They obtained a maximum PTE improvement of 56.23% over traditional coils. In our work, the application is for IMDs, the horizontal and angular misalignments between the Tx and Rx coils may also happen, and the problem can be solved by moving the Tx coil outside the human body.

3.4. A Proposed Inductive Link for the WPT System

Based on the simulation results, an inductive link for the WPT system is proposed, which may be applied to recharging the battery of a cm-sized IMD. With the concern of coupling coefficient, the Tx coil in a planar spiral circular configuration is used. To increase the MPDL, the Rx coil with 5 turns in a solenoid configuration is adopted, and from the results of Figure 5, the wire pitch is set to be 0.7 mm. Typically, the optimal load resistance of an IMD for cardiac implantable electric devices is 1 kΩ [52], and this load resistance value is assumed in this work. The geometric parameters are listed in Table 2.

Figure 9 demonstrates the maximum current allowed in the Tx coil (solid line with circles) and the peak value of SAR at the surface of the tissue (dashed line with squares), which are presented with the left and the right vertical axes, respectively. For the peak value of SAR, the excitation current in the Tx coil is set to be 1 A, while for the maximum current allowed in the Tx coil, the peak value of SAR is assumed to be 2 W/kg. As expected, the peak value of SAR increases with the frequency, while the maximum current allowed in the Tx coil decreases with frequency.

Figure 10 presents the simulation results of k, PTE, and MPDL of the proposed inductive link WPT system, with different media of air, muscle, and bone separating the two coils. The frequency-dependent dielectric properties of bone and muscle are adopted in simulations [53]. Typically, for the frequency of 25 MHz, the relative permittivity, the conductivity, and dielectric loss tangent of muscle and bone are listed in Table 3.

Figure 10(a) shows that, the values of k are almost the same with all three different media at low frequencies (<30 MHz). In the frequency range of 30 MHz < fp < 60 MHz, there are some differences of values of k, the largest is about 3% which occurs between the air and the muscle for fp = 60 MHz. The reason is that, k is associated with the permeability of the medium, which are almost the same for the three media at low frequencies. As frequency increases, so does permeability. The fastest increase of permeability is the muscle, followed by the bone. Therefore, in the frequency range of 30 MHz < fp < 60 MHz, there are some differences of values of k.

It is observed from Figure 10(b) that, in a tissue (muscle or bone) environment, the PTE first increases and then decreases with frequency, and it reaches a maximum value at fp = 25 MHz. For the frequency range of interest, the electromagnetic field loss in the air is almost frequency independent, and as k in the air increases with frequency, the PTE will also increase. However, for the tissue media, the loss of electromagnetic field will increase with frequency, and as a result, the PTE in tissue increases first (owing to the increase of k) and then decreases (owing to the increase of field loss).

Figure 10(c) shows the SAR-constrained PDL. Owing to the fact that SAR is proportional to the square of the frequency, the PDL within SAR limit thus decreases with frequency monotonically. For the air environment, there is no physical meaning for SAR; therefore, the MPDL in air is not available.

4. Experimental Verification

4.1. Measurement Methodology

To verify the simulation results, experiments are carried out, and the detailed view of the Tx and Rx coils is shown in Figure 11. r1 = 0.287 mm (#23 AWG) and r2 = 0.04 mm (#40 AWG) are chosen as the wire radius of the Tx and Rx coils, respectively. The experimental setup is similar to that of other works [17,18]. The Tx and Rx coils are supported by a nonconducting Plexiglas frame so that the distance between the two coils can be adjusted. The two coils are connected to a vector network analyzer (VNA) via a pair of subminiature version A (SMA) connectors. Since the size of the Rx coil is close to that of the SMA connector, the magnetic interference caused by the connector will affect the measured results. To overcome this drawback, the deembedding calibration method is adopted [54]. The method subtracted the parasitic components generated by the measuring fixture from the original device-under-test (DUT) measurements. The original scattering parameters, S, are written as Sdut and can be converted into admittance parameter, Y, following the relationship:where G0 is the unit matrix, Z0 is the characteristic impedance matrix of the ports and is expressed as 50G0. Similar toYdut, Yshort and Yopen can be obtained by using Sshort and Sopen, where Sshort and Sopen are the scattering parameters obtained when the circuit is short and open connected, respectively. After removing the parasitic capacitive coupling and self-capacitance, the Y parameter of the system can be expressed as [55]

The resultant Y parameters are converted into Z parameters, Z = Y−1, to obtain [18]

4.2. Measured Results of the Proposed Inductive Link WPT System

According to the simulation results, the PTE will reach the maximum value in the muscle around the operation frequency of 25 MHz; therefore, fp = 25 MHz is chosen in the experiments. Beef is used to mimic the muscle tissue, and pig’s rib bone is used to mimic the bone tissue. The parameters listed in Table 2 are adopted in the experiments.

The Z parameters are obtained by the calculations with S and Y parameters. Having obtained the Z parameters, the self-inductance, the quality factor, and the coupling coefficient of the Tx and Rx coils can be achieved with equations 2426. The SAR constrained maximum current allowed in the Tx coil, I1,max, can be obtained from Figure 9. At the frequency of 25 MHz, it is 0.384 A and 0.349 A for muscle and bone tissues, respectively. With these parameters, the PTE and the SAR-constrained maximum PDL can be calculated with equations (17) and (21). All the results are demonstrated in Table 4.

In Table 4, the measurement shows that the proposed inductive link WPT system can achieve a PTE of 52.35%, 44.14%, and 43.07% in air, muscle, and bone environments, respectively. The MPDL of the system can reach 145.38 mW in the muscle tissue and 128.13 mW in the bone tissue. Agreements between the simulation and measurement are observed.

In order to show the significant contributions of our proposed inductive link WPT system for cm-sized IMDs, we compare our work with others in Table 5. It is easy to observe that the system proposed in this work has advantages in performance.

5. Discussions

This work researches the optimizations of coupling coefficient of the inductive link in the WPT system used for IMDs. Simulations and experiments are carried out to test the proposed inductive link. Over 40% of PTE and 120 mW MPDL of the WPT system can be obtained in air, muscle, and bone environments, demonstrating the feasibility of the system. However, in real-world applications, there are several issues that need to be addressed.

5.1. Biocompatibility

All implantable devices should be sealed in a biocompatible material to be protected from the harsh environment in the human body. The Rx coil of the WPT system can be embedded in or wrapped around the biomedical implant. For biosafety reasons, if the coil is embedded in an implant, the implant cannot be sealed with metal, as the metal will shield the electromagnetic field from the Rx coil. In this case, the implant should be encapsulated with a polymeric biocompatible material such as polydimethylsiloxane or silicone. If the Rx coil is wrapped around an implant, the coil should be covered with the aforementioned biocompatible material.

5.2. Temperature Rise in the Human Body

According to human physiology, some people will go into convulsions at 41°C of body temperature, and this value is usually used as the basic temperature restriction [44]. The energy obtained in the receiving loop can be used to recharge the button cell embedded in the implant. With the low power consumption design in IC technology, the power consumed in the Rx circuit should be kept as low as possible such that the heat generated in the circuit is small, and the recharge process time can be increased.

5.3. The Recharge Cycle of the Button Cell

Generally, the power for a modern cardiac pacemaker is 6 μW (according to the reference manual of Nanostim TM, St. Jude Medical) [52]. Assuming that the inductive link WPT system provides a PDL of 120 mW, the recharge efficiency to the button cell is 50%, and the recharge time is 20 minutes; then, the recharged energy can last 3.9 years. Therefore, the button cell is only needed to be recharged once in every 3 years.

5.4. The SAR Safety

In this work, the IMD is intended to be used for left vagus nerve stimulation or cardiac pacemaker, which can be embedded in the chest. Therefore, the local SAR value of 2 W/kg is adopted. The Tx coil is set to be 2 mm away from the skin. The Tx coil can be embedded in wearable things, such as cloth, and the 2-mm separation space should be guaranteed. In practice, the local value of SAR should be measured with equipment such that the security is absolutely assured.

In the practical applications, the aforementioned issues should be carefully explored. Furthermore, the magnetic field distribution results for all verification cases on Rx movement and the changes in tissue characteristics are of interests. Future researches may be focused on all of these issues.

6. Conclusions

The analytical expression for the coupling coefficient of the inductive link is presented, and the optimizations of the coupling coefficient are provided. Based on the optimization results, an inductive link for the WPT system is proposed, and simulation and measurement results are obtained. Measurement demonstrates that, in a muscle environment, the system can achieve a PTE of 44.14%, with the SAR-constrained PDL of 145.38 mW.

With the rapid progress in IC technologies, the IMDs can be miniaturized to mm-scale and can be implanted for deep brain stimulation [5]. The advancement in biocompatible and novel nanomaterial electrode interfaces has made it possible to attach the IMDs directly to the neural tissues [42]. However, for some specific applications, such as a cardiac pacemaker or nerve stimulator, IMDs in the cm scale are preferred. With the advancement in modern techniques, this target will surely be implemented in the near future. For both convenience and comfort, these IMDs are usually powered by a button cell, as such, the scale of the IMD is restricted by the size of the cell. The inductive link WPT system proposed in this work is to recharge the cell; with this technology, surgery of replacing the IMD is avoided and the patients can be benefited both physiologically and economically.

Data Availability

The data are available from the corresponding author upon request ([email protected]).

Conflicts of Interest

The authors declare that they have no conflicts of interest.


This work was supported in part by the Chinese Postdoctoral Science Foundation (Grant no. 2020M681168), in part by the National Natural Science Foundation of China (Grant no. 62101135), and in part by the Shanghai Pujiang Program (Grant no. 20PJ1400900).