Skip to main content

ORIGINAL RESEARCH article

Front. Phys., 09 March 2023
Sec. Interdisciplinary Physics
This article is part of the Research Topic Nonlocal Integrable System and Nonlinear Waves View all 8 articles

The shifted parity and delayed time-reversal symmetry-breaking solutions for the (1+1)-dimensional Alice–Bob Boussinesq equation

Peng DongPeng Dong1Zheng-Yi Ma
Zheng-Yi Ma2*Jin-Xi FeiJin-Xi Fei3Wei-Ping CaoWei-Ping Cao3
  • 1Department of Mathematics, Zhejiang Sci-Tech University, Hangzhou, China
  • 2Department of Mathematics, Lishui University, Lishui, China
  • 3Department of Photoelectric Engineering, Lishui University, Lishui, China

The integrable Alice–Bob system with the shifted parity and delayed time reversal is presented through the Lax pair for the (1 + 1)-dimensional Boussinesq equation. After introducing an extended Bäcklund transformation, this system shows abundant exact solutions with the auxiliary functions consisting of hyperbolic functions or rational functions. The corresponding soliton structures contain line solitons, breathers, and lumps, all which satisfied the shifted parity and delayed time-reversal symmetry for the states of Alice A and Bob B. In particular, some lower-order circumstances can be expressed through their explicit solutions and their dynamic structures.

1 Introduction

For one (1 + 1)-dimensional model, except for identity transformation, there are the shifted parity P̂s and delayed time-reversal T̂d transformations for the spatial variable x and time variable t, respectively. In other words,

x=x+x0P̂sx,t=t+t0T̂dt,(1)

with x0 and t0 being arbitrary constants [1, 2]. However, the Alice–Bob system, which can be successively used to describe two-place physical problems, may be entangled with each other through the following relation:

B=f̂A=Af̂,(2)

with the state of Alice AA(x, t) and the Bob’s state BB(x′, t′); f̂ is a suitable operator (such as f̂2=1 or f̂Θ{1,P̂s,T̂d,P̂sT̂d}) [37]. Usually, this intrinsic Alice–Bob system is non-local since {x′, t′} is far away from {x, t}. However, for one Alice–Bob system, through the P̂sT̂d transformation, there indeed exist some types of multiple soliton structures with symmetry-breaking solutions according to Lou’s research. In other words, by applying the operator P̂sT̂d on one solution S, one can find P̂sT̂dSS.

For an illustrated model, the (1 + 1)-dimensional Boussinesq equation has the following form:

wtt+wxxw2xx13wxxxx=0,(3)

or

wtt=wxx+2wx2+2wwxx+13wxxxx,(4)

where wtt2t2w,wixixiw. Eq. 3 is an integrable equation as it has the following Lax pair:

ψxxx=34Iwtdxwxψ+3432wψx,(5)
ψt=I44w1ψ+Iψxx,(6)

and its adjoint version is as follows:

ψxxx=34Iwtdxwxψ+3432wψx,(7)
ψt=I44w1ψIψxx,(8)

where I is an imaginary unit.

In non-linear science, this equation is one of the important prototypic models. It can be used to study the dynamics of thin inviscid layers with a free surface, the non-linear string, and the propagation of waves in elastic rods and in the continuum limit of lattice dynamics or coupled electrical circuits. Multiple complex soliton solutions through multiple exponential function schemes, interactions between solitons and cnoidal periodic waves using the truncated Painlevé expansion method, and soliton solutions by the extended Kudryashov’s approach can all be presented [810].

Except for the former works where the physical quantity A+B2 is taken directly, we derive the Alice–Bob system of Eq. 3 through its Lax pair and the dark parameterization approach [1113].

After adopting the f̂=P̂sxT̂dt symmetry principle through B=f̂A=Af̂, the Boussinesq Eq. 3 can be induced into the following Alice–Bob system:

Att=ABA+32A212B2xx+13Axxxx,(9)
Btt=ABB+32B212A2xx+13Bxxxx.(10)

The corresponding Lax pairs of Eqs 5, 6areas follows:

Ψxxx=34IWtdxWxΨ+34E32WΨx,(11)
Ψt=I44WEΨ+IΨxx,(12)

with

Ψ=ψ1ψ2,W=wσvvw,E=1001,(13)

and we can obtain the following coupled equations:

wtt+wxxw2xx13wxxxx=0,(14)
vtt+vxx2wvxx13vxxxx=0,(15)

when σ = 0. After letting w = A+ B and v = AB, the non-local systems (9) and (10) are a direct result from Eqs 14, 15.

Another method to derive the non-local systems (9) and (10) is the dark parameterization approach [1417]. For the coupling Boussinesq system,

wi,tt=wi,xx+2j=0iwj,xwij,x+wjwij,xx+13wi,xxxx,i=0,1,2,,n.(16)

w0 is the usual solution of Eq. 3; when taking w = u + (α is a dark parameter) and n = 1, the coupled equations are as follows:

utt+uxxu2xx13uxxxx=0,(17)

and

vtt+vxx2uvxx13vxxxx=0.(18)

These equations can directly derive the non-local systems (9) and (10) through u = A+ B and v = AB.

The rest of this paper is organized as follows: in Section 2, after introducing an extended Bäcklund transformation, the Hirota bilinear form is presented through an undetermined function f, which may contain some soliton solutions for the Alice–Bob systems (9) and (10) of the (1 + 1)-dimensional Boussinesq Eq. 3. Then, the hyperbolic function solution and the rational solution for this system are shown subsequently. In order to illustrate more clearly, three kinds of explicit solutions and their corresponding soliton structures are given for the lower-order circumstances. All of these results satisfy the symmetry of B=P̂sT̂dA. A short summary is given in Section 3.

2 The symmetry-breaking solutions

We first introduce an extended Bäcklund transformation:

A=lnf+blnfx+clnftxx,B=lnfblnfxclnftxx,(19)

with b and c being two constants; ff (x, t) is an undetermined function and satisfies the following conditions:

f=f̂f=ff̂.(20)

By substituting Eq. 19 into the non-local systems (9) and (10), we obtain the following bilinear form:

Dt2+Dx213Dx4ff=0,(21)

where Hirota’s bilinear derivative operators Dt2,Dx2 and Dx4 are written as [18, 19]

DxmDtnfg=xxmttnfx,tgx,t|x=x,t=t.(22)

Eq. 21 also has the following explicit expression:

ftt=13f3ft2+3fx23ffxx+3fxx24fxfxxx+ffxxxx.(23)

2.1 The hyperbolic function solutions and their soliton structures

Through the bilinear form (21) or Eq. 23, the function f exists in the following form of the hyperbolic function for the Boussinesq Eq. 3 [46].

f=fN=νKνcoshi=1Nνiξi,ξi=kix+ωit,ωi=δiki312ki29,δi2=1,(24)

where {ν} = {νi = ±1} and ki(i = 1, 2, ⋅⋅⋅, N) are arbitrary constants, and

Kν=i<jNaij,aij=3δiδj4ki234kj2324ki236νiνjkikj24kj2+9,(25)

where (i, j = 1, 2, ⋅⋅⋅, N, ij).

When N = 1, Eq. 24 has the following simple form:

f=coshξ1,ξ1=k1x+ω1t,ω1=δ1k1312k129,δ12=1.(26)

After substituting this form into the Bäcklund transformation (19), the non-local solution of Eqs 9, 10 can be derived as follows:

A=2k12bk1+cω1tanhξ11sech2ξ1,(27)
B=2k12bk1+cω1tanhξ1+1sech2ξ1.(28)

This single-soliton solution satisfies the condition of P̂sT̂d symmetry B=f̂A=A(x,t). By introducing Eq. 26, the Alice–Bob system is the coupling form of two solitary waves. The two solitary waves move along the X-axis at the speed 12k1293; the direction is determined by δ1, and the amplitude and wave width are determined by k1, b, and c, which is also confirmed by Eqs 27, 28. Figure 1 shows this structure when the related parameters are taken as follows:

b=1,c=1,k1=32,δ1=1,ω1=1.(29)

FIGURE 1
www.frontiersin.org

FIGURE 1. (A,B) are single soliton structures of Eqs 27, 28, respectively, when the parameters taken as Eq. 29.

When N = 2, Eq. 24 becomes as follows:

f=K0coshk1+k2x+ω1+ω2t+K1coshk1k2x+ω1ω2t,(30)

with

K0=3δ1δ24k1234k22324k12+36k1k224k22+9,(31)
K1=3δ1δ24k1234k22324k1236k1k224k22+9.(32)

The corresponding two-soliton solution can be obtained by substituting Eq. 30 with Eqs 31, 32 into Eq. 19. From the perspective of algebra, it is natural to consider the simplification of the function of Eq. 30 by quantifying the double variables of the hyperbolic function into single variables, that is, k1 = ±k2, ω1 = ±ω2. These four cases may produce the corresponding soliton or breather solutions for the Alice–Bob system, respectively. Two typical cases are presented here for N = 2. Figure 2 presents this structure when the related parameters are taken through the real constants as follows:

b=1,c=1,k1=1,k2=1,δ1=1,δ2=1,ω1=33,ω2=33.(33)

Therefore,

K0=6I,K1=78I,f=6Icosh2x+78Icosh233t.(34)

FIGURE 2
www.frontiersin.org

FIGURE 2. (A,B) are the two-soliton structures of the Alice-Bob systems (9) and (10), respectively through Eq. 19 after selecting the conditions are Eqs 33, 34.

It is not difficult to find that Eq. 34 is coupled by two hyperbolic functions similar to Eq. 26, and its image also shows this phenomenon.

On the other hand, by restricting the parameters k1, k2 to the assumed units on the two-soliton solution, the breather can be obtained. For example, by setting the following parameters

b=1,c=1,k1=I,k2=I,δ1=1,δ2=1,ω1=213,ω2=213,(35)

the following equation can be derived:

K0=42,K1=114,f=42cos2x+114cosh2213t.(36)

Here, the cosine part of Eq. 36 makes the Alice–Bob system periodic, and the corresponding breather structure is obtained, as shown in Figure 3.

FIGURE 3
www.frontiersin.org

FIGURE 3. Breather structures of the Alice–Bob system (9) and (10) through Eq. 19 after selecting the conditions are Eqs 35, 36. (C,D) are the corresponding sectional plots of (A,B) at t =0, respectively.

When N = 3, Eq. 24 has the following more complicated situation:

f=K0coshξ1+ξ2+ξ3+K1coshξ1ξ2ξ3+K2coshξ1ξ2+ξ3+K3coshξ1+ξ2ξ3,(37)

with

K0=3δ1δ24k1234k22324k12+36k1k224k22+9,
3δ1δ34k1234k32324k12+36k1k324k32+93δ2δ34k2234k32324k22+36k2k324k32+9,
K1=3δ1δ24k1234k22324k1236k1k224k22+9
3δ1δ34k1234k32324k1236k1k324k32+93δ2δ34k2234k32324k22+36k2k324k32+9,
K2=3δ1δ24k1234k22324k1236k1k224k22+9,
3δ1δ34k1234k32324k12+36k1k324k32+93δ2δ34k2234k32324k2236k2k324k32+9,
K3=3δ1δ24k1234k22324k12+36k1k224k22+9,
3δ1δ34k1234k32324k1236k1k324k32+93δ2δ34k2234k32324k2236k2k324k32+9,

and

ξi=kix+ωit,ωi=δiki312ki29,δi2=1i=1,2,3.

Based on the selecting parameters b, c, k1, k2, δ1, and δ2 of Eqs 33, 35, two kinds of interactions for the solitons can be constructed by considering the following equation:

k3=1110,δ3=1,ω3=11138150,(38)

and

k3=32,δ3=1,ω3=322.(39)

For this time, K{i}(i = 0, 1, 2, 3) are expressed as follows:

K0=9146I 25,K1=17296578231660046I25,
K2=17296578+231660046I25,K3=636594I25,(40)

and

K0=321233+1369721233+13697I,K1=31144153642,
K2=3114415+3642,K3=321233+13697+21233+13697I.(41)

The corresponding functions of Eq. 37 are expressed as follows:

f=K0cosh3110x+11138150t+K1cosh1110x+23311138150t
+K2cosh1110x+233+11138150tK3cosh910x+11138150t,(42)

and

f=621233+13697sin2xsinh32x+322t+621233+13697cos2xcosh32x+322t
+3114415+3642cosh32x+3222213t+31144153642cosh32x+3222213t.(43)

Figure 4 and Figure 5 show two interaction structures of the Alice–Bob systems (9) and (10) through Eqs 42, 43.

FIGURE 4
www.frontiersin.org

FIGURE 4. (A,B) are the interaction structures of three solitons for the Alice–Bob systems (9) and (10) through Eq. 42. (C,D) are the corresponding density plots of (A,B), respectively.

FIGURE 5
www.frontiersin.org

FIGURE 5. (A,B) are the interaction structures between the breather and one-soliton structures for the Alice–Bob systems (9) and (10) through Eq. 43. (C,D) are the corresponding density plots of (A,B), respectively.

2.2 The rational solutions and lump structures

The Alice–Bob systems (9) and (10) have a series of rational solutions and hence contain the corresponding lump structures. For this purpose, we introduce the following polynomial function:

f=fN=m=012NN+1j=0maj,mx2mjt2j,(44)

where aj,m(j=0,1,2,,m,m=0,1,2,,12N(N+1),N=1,2,) are constants determined by the powers of the variables x and t [20, 21].

When N = 1, Eq. 44 has the following simple form:

f=a0,0+a0,1x2+a1,1t2.(45)

When

b=c=a0,0=a0,1=a1,1=1,(46)

the lump solution of the Alice–Bob systems (9) and (10) has the following rational form:

A=216x2t+2t2+2x3+6x2t6xt22t3x4+t41+x2+t23,(47)
B=21+6x+2t+2t22x36x2t+6xt2+2t3x4+t41+x2+t23,(48)

which is obtained through the Bäcklund transformation (19) (Figure 6).

FIGURE 6
www.frontiersin.org

FIGURE 6. (A,B) are two lump structures of A and B from Eqs 47, 48. (C,D) are the corresponding density plots of (A,B), respectively.

When N = 2, Eq. 44 has the following function form:

f=a0,0+a0,1x2+a0,2x4+a0,3x6+a1,1t2+a1,2x2t2+a1,3x4t2+a2,2t4+a2,3x2t4+a3,3t6.(49)

A pair of lumps of A and B for the Alice–Bob systems (9) and (10) can be shown after the constants are taken as b = c = a0,0 = 1, just as Eq. 46, while

a0,1=15,a0,2=325,a0,3=9625,a1,1=1925,a1,2=54125,a1,3=27625,
a2,2=57625,a2,3=27625,a3,3=9625.(50)

These two pairs of lumps for the Alice–Bob systems (9) and (10) through Eq. 19, Eq. 49, and Eq. 50 are shown in Figure 7.

FIGURE 7
www.frontiersin.org

FIGURE 7. (A,B) are two pairs of lumps for the Alice–Bob systems (9) and (10) through Eqs 19, 49, 50. (C,D) are the corresponding density plots of (A,B), respectively.

When N = 3, Eq. 44 has the more complicated function form:

f=a0,0+a0,1x2+a0,2x4+a0,3x6+a0,4x8+a0,5x10+a0,6x12+a1,1t2+a1,2x2t2+a1,3x4t2+a1,4x6t2+a1,5x8t2+a1,6x10t2+a2,2t4+a2,3x2t4+a2,4x4t4+a2,5x6t4+a2,6x8t4+a3,3t6+a3,4x2t6+a3,5x4t6+a3,6x6t6+a4,4t8+a4,5x2t8+a4,6x4t8+a5,5t10+a5,6x2t10+a6,6t12.(51)

The lumps of A and B for the Alice–Bob systems (9) and (10) can also be constructed after the constants b = c = a0,0 = 1; therefore, the following equations are obtained:

a0,1=1811,a0,2=135847,a0,3=32446585,a0,4=21873587045,a0,5=437417935225,a0,6=6561878826025,
a1,1=2610847,a1,2=4869317,a1,3=43740717409,a1,4=554043587045,a1,5=301806175765205,a1,6=39366878826025,
a2,2=361287717409,a2,3=2916717409,a2,4=1560065021863,a2,5=87482282665,a2,6=19683175765205,
a3,3=184906825109315,a3,4=670682282665,a3,5=638604175765205,a3,6=26244175765205,
a4,4=632043175765205,a4,5=249318175765205,a4,6=19683175765205,
a5,5=126846878826025,a5,6=39366878826025,a6,6=6561878826025.(52)

These lump structures of the Alice–Bob systems (9) and (10) obtained through Eqs 19, 51, and 52 are shown in Figure 8.

FIGURE 8
www.frontiersin.org

FIGURE 8. (A,B) are lump structures of A and B from Eqs 9, 10 through Eqs 19, 51, and 52. (C,D) are the corresponding density plots of (A,B), respectively.

3 Summary

In this paper, according to the (1 + 1)-dimensional Boussinesq Eq. 3, the Alice–Bob systems (9) and (10) for this equation are first derived through the Lax pair and the dark parameterization approach. This non-local system owns the bilinear form and may exist in the explicit solution. Therefore, the N-soliton solutions of the Alice–Bob systems (9) and (10) are presented with the aid of an undetermined function f after introducing an extended Bäcklund transformation. Typically, the auxiliary function can be taken as the hyperbolic function or rational function. These two kinds of functions induce the system having solutions that satisfy B=P̂sT̂dA. The lower-order circumstances for N = 1, 2, 3 are presented through their auxiliary functions, and the symmetry-breaking solutions can be constructed. With the special parameters, the antisymmetric local structures are depicted, which contain line solitons, breathers, and lumps. Whether the induced Alice–Bob systems (9) and (10) of the (1 + 1)-dimensional Boussinesq Eq. 3 or the derived results through the hyperbolic and rational functions satisfy the symmetry of B=P̂sT̂dA is first shown here for our understanding. We believe that this approach is important to solve the Alice–Bob system for one integrable equation, which may possess rich local structures.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material; further inquiries can be directed to the corresponding author.

Author contributions

PD: writing—original draft and software. Z-YM: supervision and funding acquisition. J-XF: conceptualization and software. W-PC: formal analysis and investigation.

Funding

This work was supported by the National Natural Science Foundation of China (11775104) and the Zhejiang Province Natural Science Foundation of China (2022SJGYZC01).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

1. Bender CM, Boettcher S. Real spectra in non-hermitian Hamiltonians HavingPTSymmetry. Phys Rev Lett (1998) 80:5243–6. doi:10.1103/physrevlett.80.5243

CrossRef Full Text | Google Scholar

2. Ablowitz MJ, Musslimani ZH. Integrable nonlocal nonlinear Schrödinger equation. Phys Rev Lett (2013) 110:064105–5. doi:10.1103/physrevlett.110.064105

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Lou SY. Alice-Bob systems, P̂-T̂-Ĉ symmetry invariant and symmetry breaking soliton solutions. J Math Phys 2018;59:083507–20. doi:10.1063/1.5051989

CrossRef Full Text | Google Scholar

4. Lou SY. Prohibitions caused by nonlocality for nonlocal Boussinesq-KdV type systems. Stud Appl Math (2019) 142:1–16. doi:10.1111/sapm.12265

CrossRef Full Text | Google Scholar

5. Li H, Lou SY. Multiple soliton solutions of Alice-Bob Boussinesq equations. Chin Phys Lett (2019) 36:050501–5. doi:10.1088/0256-307x/36/5/050501

CrossRef Full Text | Google Scholar

6. Zhao QL, Lou SY, Jia M. Solitons and soliton molecules in two nonlocal Alice-Bob Sawada-Kotera systems. Commun Theor Phys (2020) 72:085005–7. doi:10.1088/1572-9494/ab8a0e

CrossRef Full Text | Google Scholar

7. Jia M, Lou SY. Exact PsTd invariant and PsTd symmetric breaking solutions, symmetry reductions and Bäcklund transformations for an AB-KdV system. Phys Lett A (2018) 382:1157–66. doi:10.1016/j.physleta.2018.02.036

CrossRef Full Text | Google Scholar

8. Wazwaz AM. Multiple soliton solutions and multiple complex soliton solutions for two distinct Boussinesq equations. Nonlinear Dyn (2016) 85:731–7. doi:10.1007/s11071-016-2718-0

CrossRef Full Text | Google Scholar

9. Yang D, Lou SY, Yu WF. Interactions between solitons and cnoidal periodic waves of the Boussinesq equation. Commun Theor Phys (2013) 60:387–90. doi:10.1088/0253-6102/60/4/01

CrossRef Full Text | Google Scholar

10. Adem AR, Yildirim Y, Yasar E. Soliton solutions to the non-local Boussinesq equation by multiple exp-function scheme and extended Kudryashov’s approach. Pramana-j Phys (2019) 92:24–11. doi:10.1007/s12043-018-1679-x

CrossRef Full Text | Google Scholar

11. Fei JX, Ma ZY, Cao WP. Controllable symmetry breaking solutions for a nonlocal Boussinesq system. Sci Rep-uk (2019) 9:19667–14. doi:10.1038/s41598-019-56093-8

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Cao WP, Fei JX, Li JY. Symmetry breaking solutions to nonlocal Alice-Bob Kadomtsev-Petviashivili system. Chaos Soliton Fract (2021) 144:110653–8. doi:10.1016/j.chaos.2021.110653

CrossRef Full Text | Google Scholar

13. Wu WB, Lou SY. Exact solutions of an Alice-Bob KP equation. Commun Theor Phys (2019) 71:629–32. doi:10.1088/0253-6102/71/6/629

CrossRef Full Text | Google Scholar

14. Lou SY. Dark parameterizations, equivalent partner fields and integrable systems. Commun Theor Phys (2011) 55:743–6. doi:10.1088/0253-6102/55/5/02

CrossRef Full Text | Google Scholar

15. Peng L, Yang XD, Lou SY. Integrable KP coupling and its exact solution. Commun Theor Phys (2012) 58:1–4. doi:10.1088/0253-6102/58/1/01

CrossRef Full Text | Google Scholar

16. Cao WP, Fei JX, Ma ZY, Liu Q. Bosonization and new interaction solutions for the coupled Korteweg-de Vries system. Wave Random Complex (2020) 30:130–41. doi:10.1080/17455030.2018.1489565

CrossRef Full Text | Google Scholar

17. Fei JX, Cao WP, Ma ZY. Dark parameterization approach to the Benjamin Ono equation. Int J Mod Phys B (2020) 34:2050247–8. doi:10.1142/s0217979220502471

CrossRef Full Text | Google Scholar

18. Hirota R. Exact solution of the Korteweg-de Vries equation for multiple collisions of solitons. Phys Rev Lett (1971) 27:1192–4. doi:10.1103/physrevlett.27.1192

CrossRef Full Text | Google Scholar

19. Hirota R. The direct method in soliton theory. Cambridge: Cambridge University Press (2004).

Google Scholar

20. Clarkson PA, Dowie E. Rational solutions of the Boussinesq equation and applications to rogue waves. Trans Math Appl (2017) 1:1–26. doi:10.1093/imatrm/tnx003

CrossRef Full Text | Google Scholar

21. Ma YL, Li BQ. Analytic rogue wave solutions for a generalized fourth-order Boussinesq equation in fluid mechanics. Math Meth Appl Sci (2019) 42:39–48. doi:10.1002/mma.5320

CrossRef Full Text | Google Scholar

Keywords: (1+1)-dimensional Alice–Bob Boussinesq equation, Lax pair, Bäcklund transformation, PsTd symmetry-breaking solution, hybrid structure

Citation: Dong P, Ma Z-Y, Fei J-X and Cao W-P (2023) The shifted parity and delayed time-reversal symmetry-breaking solutions for the (1+1)-dimensional Alice–Bob Boussinesq equation. Front. Phys. 11:1137999. doi: 10.3389/fphy.2023.1137999

Received: 05 January 2023; Accepted: 16 February 2023;
Published: 09 March 2023.

Edited by:

Yunqing Yang, Zhejiang Ocean University, China

Reviewed by:

Sakkaravarthi Karuppaiya, Asia Pacific Center for Theoretical Physics (APCTP), Republic of Korea
Alejandro Aceves, Southern Methodist University, United States

Copyright © 2023 Dong, Ma, Fei and Cao. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Zheng-Yi Ma, ma-zhengyi@163.com

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.