Processing math: 83%
Research article

Asymmetric Effects on Risks of Virtual Financial Assets (VFAs) in different regimes: A Case of Bitcoin

  • Received: 12 November 2018 Accepted: 20 November 2018 Published: 22 November 2018
  • JEL Codes: C32, C44, C51

  • The rapid development of VFAs allows investors to diversify their choices of investment products. In this paper, we measure the return risk of VFAs based on GARCH-type model. By establishing a Markov regime-switching Regression (MSR) Model, we explore the asymmetric effects of speculation, investor attention, and market interoperability on return risks in different risk regimes of VFAs. The results show that the influences of speculation and investor attention on the risks of VFAs are significantly positive at all regimes, while market interoperability only admits a positive impact on risk under high risk regime. All of the three factors exert asymmetric effects on risks in different regimes. Further study presents that the risk regime-switching also shows asymmetric characteristic but the medium risk regime is more stable than any others. Therefore, transactions of investors and arbitrageurs are monitored by certain policies, such as limiting the number of transactions or restricting the trading amount at high risk regime. However, when return risk is low, it will return to a medium level if we encourage investors to access.

    Citation: Zhenghui Li, Hao Dong, Zhehao Huang, Pierre Failler. Asymmetric Effects on Risks of Virtual Financial Assets (VFAs) in different regimes: A Case of Bitcoin[J]. Quantitative Finance and Economics, 2018, 2(4): 860-883. doi: 10.3934/QFE.2018.4.860

    Related Papers:

    [1] Bilal Akbar, Khuram Pervez Amber, Anila Kousar, Muhammad Waqar Aslam, Muhammad Anser Bashir, Muhammad Sajid Khan . Data-driven predictive models for daily electricity consumption of academic buildings. AIMS Energy, 2020, 8(5): 783-801. doi: 10.3934/energy.2020.5.783
    [2] Zhongjiao Ma, Zichun Yan, Mingfei He, Haikuan Zhao, Jialin Song . A review of the influencing factors of building energy consumption and the prediction and optimization of energy consumption. AIMS Energy, 2025, 13(1): 35-85. doi: 10.3934/energy.2025003
    [3] Kaovinath Appalasamy, R Mamat, Sudhakar Kumarasamy . Smart thermal management of photovoltaic systems: Innovative strategies. AIMS Energy, 2025, 13(2): 309-353. doi: 10.3934/energy.2025013
    [4] Feng Tian . A wind power prediction model with meteorological time delay and power characteristic. AIMS Energy, 2025, 13(3): 517-539. doi: 10.3934/energy.2025020
    [5] Xiaojing Tian, Weiqi Ye, Liang Xu, Anjian Yang, Langming Huang, Shenglong Jin . Optimization research on laminated cooling structure for gas turbines: A review. AIMS Energy, 2025, 13(2): 354-401. doi: 10.3934/energy.2025014
    [6] Emrah Gulay . Forecasting the total electricity production in South Africa: Comparative analysis to improve the predictive modelling accuracy. AIMS Energy, 2019, 7(1): 88-110. doi: 10.3934/energy.2019.1.88
    [7] Rahaman Abu, John Amakor, Rasaq Kazeem, Temilola Olugasa, Olusegun Ajide, Nosa Idusuyi, Tien-Chien Jen, Esther Akinlabi . Modeling influence of weather variables on energy consumption in an agricultural research institute in Ibadan, Nigeria. AIMS Energy, 2024, 12(1): 256-270. doi: 10.3934/energy.2024012
    [8] Yan Li, Yaheng Su, Qixin Zhao, Bala Wuda, Kaibo Qu, Lei Tang . An electricity price optimization model considering time-of-use and active distribution network efficiency improvements. AIMS Energy, 2025, 13(1): 13-34. doi: 10.3934/energy.2025002
    [9] Anbazhagan Geetha, S. Usha, J. Santhakumar, Surender Reddy Salkuti . Forecasting of energy consumption rate and battery stress under real-world traffic conditions using ANN model with different learning algorithms. AIMS Energy, 2025, 13(1): 125-146. doi: 10.3934/energy.2025005
    [10] Ghafi Kondi Akara, Benoit Hingray, Adama Diawara, Arona Diedhiou . Effect of weather on monthly electricity consumption in three coastal cities in West Africa. AIMS Energy, 2021, 9(3): 446-464. doi: 10.3934/energy.2021022
  • The rapid development of VFAs allows investors to diversify their choices of investment products. In this paper, we measure the return risk of VFAs based on GARCH-type model. By establishing a Markov regime-switching Regression (MSR) Model, we explore the asymmetric effects of speculation, investor attention, and market interoperability on return risks in different risk regimes of VFAs. The results show that the influences of speculation and investor attention on the risks of VFAs are significantly positive at all regimes, while market interoperability only admits a positive impact on risk under high risk regime. All of the three factors exert asymmetric effects on risks in different regimes. Further study presents that the risk regime-switching also shows asymmetric characteristic but the medium risk regime is more stable than any others. Therefore, transactions of investors and arbitrageurs are monitored by certain policies, such as limiting the number of transactions or restricting the trading amount at high risk regime. However, when return risk is low, it will return to a medium level if we encourage investors to access.


    In the present paper, we study the following coupled nonlinear wave equations (see [1])

    uttdiv(ρ(|u|2)u)+|ut|m1ut+f1(u,v)=0, (1.1)
    vttdiv(ρ(|v|2)v)+|vt|n1vt+f2(u,v)=0, (1.2)

    where ut=ut,vt=vt, m,n>1, is the gradient operator, div is the divergence operator and fi(.,.):R2R,i=1,2 are known functions. For arbitrary solutions of (1.1) and (1.2), the function ρ is supposed to satisfy one or the other of the two conditions:

    A1). 0<ρ2(q2)q2m1˜ρ(q2),

    or

    A2). 0<1K1<ρ1(q2)K2˜ρ(q2)1q2,

    where ˜ρ(q2)=q20ρ(ζ)dζ, q2=|u|2,m1,K1,K2>0.

    Quintanilla [2] imposed condition A1 and obtained the growth or decay estimates of the solution to the type Ⅲ heat conduction. The condition A1 can be satisfied easily, e.g., ρ(q2)=11+q2 or ρ(q2)=(1+bpq2)p1,b>0,0<p1. The similar condition as A2 was also considered by many papers (see [3,4]). ρ(q2)=1+q2 satisfies the condition A2.

    In addition, we introduce a function F(u,v) which is defined as Fu=f1(u,v), Fv=f2(u,v), where F(0,0)=0.

    The wave equations have attracted many attentions of scholars due to their wide application, and a large number of achievements have been made in the existence of solutions (see [1,5,6,7,8,9,10,11,12]). The fuzzy inference method is used to solve this problem. The algebraic formulation of fuzzy relation is studied in [13,14]. In this paper, we study the Phragmén-Lindelöf type alternative property of solutions of wave equations (1.1)–(1.2). It is proved that the solution of the equations either grows exponentially (polynomially) or decays exponentially (polynomially) when the space variable tends to infinite. In the case of decay, people usually expect a fast decay rate. The Phragmén-Lindelöf type alternative research on partial differential equations has lasted for a long time (see [2,15,16,17,18,19,20,21,22,23]).

    It is worth emphasizing that Quintanilla [2] considered an exterior or cone-like region. Under some appropriate conditions, the growth/decay estimates of some parabolic problems are obtained. Inspired by [2], we extends the research of to the nonlinear wave model in this paper. However, different from [2], in addition to condition A1 and condition A2, we also consider a special condition of ρ. The appropriate energy function is established, and the nonlinear differential inequality about the energy function is derived. By solving this differential inequality, the Phragmén-Lindelöf type alternative results of the solution are obtained. Our model is much more complex than [2], so the methods used in [2] can not be applied to our model directly. Finally, a nonlinear system of viscoelastic type is also studied when the system is defined in an exterior or cone-like regions and the growth or decay rates are also obtained.

    Letting that Ω(r) denotes a cone-like region, i.e.,

    Ω(r)={x||x|2r2,rR0>0},

    and that B(r) denotes the exterior surface to the sphere, i.e.,

    B(r)={x||x|2=r2,rR0>0},

    Equations (1.1) and (1.2) also have the following initial-boundary conditions

    u(x,0)=v(x,0)=0, in Ω, (2.1)
    u(x,t)=g1(x,t), v(x,t)=g2(x,t), in B(R0)×(0,τ), (2.2)

    where g1 and g2 are positive known functions, x=(x1,x2,x3) and τ>0.

    Now, we establish an energy function

    E(r,t)=t0B(r)eωηρ(|u|2)uxruηdSdη+t0B(r)eωηρ(|v|2)vxrvηdSdηE1(r,t)+E2(r,t). (2.3)

    Let r0 be a positive constant which satisfies r>r0R0. Integrating E(z,t) from r0 to r, using the divergence theorem and Eqs (1.1) and (1.2), (2.1) and (2.2), we have

    E(r,t)E(r0,t)=t0rr0B(ξ)eωη[ρ(|u|2)uuη]dsdξdη+t0rr0B(ξ)eωη[ρ(|v|2)vvη]dsdξdη+t0rr0B(ξ)eωη[uηη+|uη|m+1+f1(u,v)]uηdsdξdη+t0rr0B(ξ)eωη[vηη+|vη|n+1+f2(u,v)]vηdsdξdη+12t0rr0B(ξ)eωηη˜ρ(|u|2)dsdξdη+12t0rr0B(ξ)eωηη˜ρ(|v|2)dsdξdη=12eωtrr0B(ξ)[|ut|2+|vt|2+˜ρ(|u|2)+˜ρ(|v|2)+2F(u,v)]dsdξ+12ωt0rr0B(ξ)eωη[|uη|2+|vη|2+˜ρ(|u|2)+˜ρ(|v|2)+2F(u,v)]dSdξdη+t0rr0B(ξ)eωη[|uη|m+1+|vη|n+1]dsdξdη, (2.4)

    from which it follows that

    rE(r,t)=12eωtB(r)[|ut|2+|vt|2+˜ρ(|u|2)+˜ρ(|v|2)+2F(u,v)]ds+12ωt0B(r)eωη[|uη|2+|vη|2+˜ρ(|u|2)+˜ρ(|v|2)+2F(u,v)]dSdη+t0B(r)eωη[|uη|m+1+|vη|n+1]dsdη, (2.5)

    where ω is positive constant.

    Now, we show how to bound E(r,t) by rE(r,t). We use the Hölder inequality, the Young inequality and the condition A1 to have

    |E1(r,t)|[t0B(r)eωηρ2(|u|2)|u|2dsdηt0B(r)eωηu2ηdsdη]12m1[t0B(r)eωη˜ρ(|u|2)dsdηt0B(r)eωηu2ηdsdη]12m12[t0B(r)eωη˜ρ(|u|2)dsdη+t0B(r)eωηu2ηdsdη], (2.6)

    and

    |E2(r,t)|m12[t0B(r)eωη˜ρ(|v|2)dsdη+t0B(r)eωηv2ηdsdη]. (2.7)

    Inserting (2.6) and (2.7) into (2.3) and combining (2.5), we have

    |E(r,t)|m1ω[rE(r,t)]. (2.8)

    We consider inequality (2.8) for two cases.

    I. If r0>R0 such that E(r0,t)0. From (2.5), we obtain E(r,t)E(r0,t)0,rr0. Therefore, from (2.8) we have

    E(r,t)m1ω[rE(r,t)],rr0. (2.9)

    Integrating (2.9) from r0 to r, we have

    E(r,t)[E(r0,t)]eωm1(rr0),rr0. (2.10)

    Combing (2.4) and (2.10), we have

    limr{eωm1r[12eωtrr0B(ξ)[|ut|2+|vt|2+˜ρ(|u|2)+˜ρ(|v|2)+2F(u,v)]dsdξ+12ωt0rr0B(ξ)eωη[|uη|2+|vη|2+˜ρ(|u|2)+˜ρ(|v|2)+2F(u,v)]dsdξdη+t0rr0B(ξ)eωη[|uη|m+1+|vη|n+1]dsdξdη]}E(R0,t)eωm1R0. (2.11)

    II. If r>R0 such that E(r,t)<0. The inequality (2.8) can be rewritten as

    E(r,t)m1ω[rE(r,t)],rR0. (2.12)

    Integrating (2.12) from r0 to r, we have

    E(r,t)[E(R0,t)]eωm1(rR0),rR0. (2.13)

    Inequality (2.13) shows that limr[E(r,t)]=0. Integrating (2.5) from r to and combining (2.13), we obtain

    12eωtrB(ξ)[|ut|2+|vt|2+˜ρ(|u|2)+˜ρ(|v|2)+2F(u,v)]dsdξ+12ωt0rB(ξ)eωη[|uη|2+|vη|2+˜ρ(|u|2)+˜ρ(|v|2)+2F(u,v)]dsdξdη+t0rB(ξ)eωη[|uη|m+1+|vη|n+1]dsdξdη[E(r0,t)]eωm1(rR0). (2.14)

    We summarize the above result as the following theorem.

    Theorem 2.1. Let (u,v) be solution of the (1.1), (1.2), (2.1), (2.2) in Ω(R0), and ρ satisfies condition A1. Then for fixed t, (u,v) either grows exponentially or decays exponentially. Specifically, either (2.11) holds or (2.14) holds.

    If the function ρ satisfies the condition A2, we recompute the inequality (2.6) and (2.7). Therefore

    |E1(r,t)|[t0B(r)eωηρ2(|u|2)|u|2dsdηt0B(r)eωηu2ηdsdη]12K1K1[t0B(r)eωηρ1(|u|2)dsdηt0B(r)eωηu2ηdsdη]12K1K1K22[t0B(r)eωη˜ρ(|u|2)dsdη+t0B(r)eωηu2ηdsdη], (3.1)

    and

    |E2(r,t)|K1K1K22[t0B(r)eωη˜ρ(|v|2)dsdη+t0B(r)eωηv2ηdsdη]. (3.2)

    Inserting (3.1) and (3.2) into (2.3) and combining (2.5), we have

    |E(r,t)|K1K1K2ω[rE(r,t)]. (3.3)

    By following a similar method to that used in Section 2, we can obtain the Phragmén-Lindelöf type alternative result.

    Theorem 3.1. Let (u,v) be solution of the (1.1), (1.2), (2.1), (2.2) in Ω(R0), and ρ satisfies condition A1. Then for fixed t, (u,v) either grows exponentially or decays exponentially.

    Remark 3.1. Clearly, the rates of growth or decay obtained in Theorems 2.1 and 3.1 depend on ω. Because ω can be chosen large enough, the rates of growth or decay of the solutions can become large as we want.

    Remark 3.2. The analysis in Sections 2 and 3 can be adapted to the single-wave equation

    uttdiv(ρ(|u|2)u)+|ut|m1ut+f(u)=0 (3.4)

    and the heat conduction at low temperature

    autt+butcΔu+Δut=0, (3.5)

    where a,b,c>0.

    In this section, we suppose that ρ satisfies ρ(q2)=b1+b2q2β, where b1,b2 and β are positive constants. Clearly, ρ(q2)=b1+b2q2β can not satisfy A1 or A2. In this case, we define an "energy" function

    F(r,t)=t0B(r)eωη(b1+b2|u|2β)uxruηdsdη+t0B(r)eωη(b1+b2|v|2β)vxrvηdsdηF1(r,t)+F2(r,t). (4.1)

    Computing as that in (2.4) and (2.5), we can get

    F(r,t)=F(r0,t)+12eωtrr0B(ξ)[|ut|2+b1|u|2+1β+1b2|u|2(β+1)]dsdξ+12eωtrr0B(ξ)[|vt|2+b1|v|2+1β+1b2|v|2(β+1)+F(u,v)]dsdξ+12ωt0rr0B(ξ)eωη[|uη|2+b1|u|2+1β+1b2|u|2(β+1)]dsdξdη+12ωt0rr0B(ξ)eωη[|vη|2+b1|v|2+1β+1b2|v|2(β+1)+F(u,v)]dsdξdη+t0rr0B(ξ)eωη|uη|m+1dsdξdη+t0B(r)eωη|vη|n+1dsdξdη, (4.2)

    and

    rF(r,t)=12eωtB(r)[|ut|2+b1|u|2+1β+1b2|u|2(β+1)]ds+12eωtB(r)[|vt|2+b1|v|2+1β+1b2|v|2(β+1)+F(u,v)]ds+12ωt0B(r)eωη[|uη|2+b1|u|2+1β+1b2|u|2(β+1)]dsdη+12ωt0B(r)eωη[|vη|2+b1|v|2+1β+1b2|v|2(β+1)+F(u,v)]dsdη+t0B(r)eωη|uη|m+1dsdη+t0B(r)eωη|vη|n+1dsdη. (4.3)

    Using the Hölder inequality and Young's inequality, we have

    |F1(r,t)|b1[t0B(r)eωη|u|2dsdηt0B(r)eωηu2ηdsdη]12+b2[t0B(r)eωη|u|2(β+1)dsdη]2β+12(β+1)t0B(r)eωη|uη|m+1dsdη]1m+1|t2πr|12(β+1)1m+1b1R0ω[b12rωt0B(r)eωη|u|2dsdη+12rωt0B(r)eωηu2ηdsdη]+b2|2tπ|12(β+1)1m+1Rββ+12m+10[2β+12(β+1)r2β+12(β+1)+1m+1t0B(r)eωη|u|2(β+1)dsdη+1m+1r2β+12(β+1)+1m+1t0B(r)eωη|uη|m+1dsdη]2β+12(β+1)+1m+1, (4.4)

    and

    |F2(r,t)|b1R0ω[b12rωt0B(r)eωη|v|2dsdη+12rωt0B(r)eωηv2ηdsdη]+b2|2tπ|12(β+1)1n+1Rββ+12n+10[2β+12(β+1)r2β+12(β+1)+1n+1t0B(r)eωη|v|2(β+1)dsdη+1n+1r2β+12(β+1)+1n+1t0B(r)eωη|vη|n+1dsdη]2β+12(β+1)+1n+1, (4.5)

    where we have chosen 12(β+1)>1n+1. Inserting (4.4) and (4.5) into (4.1), we obtain

    |F(r,t)|c1[rrF(r,t)]+c2[rrF(r,t)]2β+12(β+1)+1m+1+c3[rrF(r,t)]2β+12(β+1)+1n+1, (4.6)

    where c1=b1R0ω,c2=b2|2tπ|12(β+1)1n+1Rββ+12m+10,c3=b2|2tπ|12(β+1)1n+1Rββ+12n+10.

    Next, we will analyze Eq (4.6) in two cases

    I. If r0R0 such that F(r0,t)0, then F(r,t)F(r0,t)0,rr0. Therefore, (4.6) can be rewritten as

    F(r,t)c1[rrF(r,t)]+c2[rrF(r,t)]2β+12(β+1)+1m+1+c3[rrF(r,t)]2β+12(β+1)+1n+1,rr0. (4.7)

    Using Young's inequality, we have

    [rrF(r,t)]2β+12(β+1)+1m+1(4β+34(β+1)+12(m+1))[rrF(r,t)]12+(14(β+1)12(m+1))[rrF(r,t)], (4.8)
    [rrF(r,t)]2β+12(β+1)+1n+1(4β+34(β+1)+12(n+1))[rrF(r,t)]12+(14(β+1)12(n+1))[rrF(r,t)]. (4.9)

    Inserting (4.8) and (4.9) into (4.7), we have

    F(r,t)c4[rrF(r,t)]12+c5[rrF(r,t)],rr0. (4.10)

    where c4=c2(4β+34(β+1)+12(m+1))+c3(4β+34(β+1)+12(n+1)) and c5=c1+c2(14(β+1)12(m+1))+c3(14(β+1)12(n+1)). From (4.10) we have

    F(r,t)c5[rrF(r,t)+c42c5]2c244c25,rr0

    or

    rF(r,t)[F(r,t)c5+c244c25c42c5]21r,rr0 (4.11)

    Integrating (4.11) from r0 to r, we get

    2c5[ln(F(r,t)c5+c244c25c42c5)ln(F(r0,t)c5+c244c25c42c5)]c4{[F(r,t)c5+c244c25c42c5]1[F(r0,t)c5+c244c25c42c5]1}ln(rr0),rr0 (4.12)

    Dropping the third term on the left of (4.12), we have

    F(r,t)c5+c244c25c42c5Q1(r0,t)(rr0)12c5, (4.13)

    where Q1(r0,t)=(F(r0,t)c5+c244c25c42c5)exp{c42c5[F(r0,t)c5+c244c25c42c5]1}.

    In view of

    F(r,t)c5+c244c25F(r,t)c5+c42c5,

    we have from (4.13)

    F(r,t)c5Q21(r0,t)(rr0)1c5. (4.14)

    Combining (4.2) and (4.14), we have

    limr{r1c5[12eωtrr0B(ξ)[|ut|2+b1|u|2+1β+1b2|u|2(β+1)]dsdξ+12eωtrr0B(ξ)[|vt|2+b1|v|2+1β+1b2|v|2(β+1)+F(u,v)]dsdξ+12ωt0rr0B(ξ)eωη[|uη|2+b1|u|2+1β+1b2|u|2(β+1)]dsdξdη+12ωt0rr0B(ξ)eωη[|vη|2+b1|v|2+1β+1b2|v|2(β+1)+F(u,v)]dsdξdη+t0rr0B(ξ)eωη|uη|m+1dsdξdη+t0B(r)eωη|vη|n+1dsdξdη]}c5Q21(r0,t)r1c50. (4.15)

    II. If rR0 such that F(r,t)<0, then (4.6) can be rewritten as

    F(r,t)c1[rrF(r,t)]+c2[rrF(r,t)]2β+12(β+1)+1m+1+c3[rrF(r,t)]2β+12(β+1)+1n+1,rR0. (4.16)

    Without losing generality, we suppose that m>n>1.

    [rrF(r,t)]2β+12(β+1)+1n+11n+11m+12β+12(β+1)1m+1[rrF(r,t)]+2β+12(β+1)1m+12β+12(β+1)1m+1[rrF(r,t)]2β+12(β+1)+1m+1, (4.17)
    [rrF(r,t)]2β(m+3)+4(2β+1)(m+2)+1[rrF(r,t)]2β+12(β+1)+1m+1+m+12(β+1)(2β+1)(m+2)+1[rrF(r,t)]2[2β+12(β+1)+1m+1]. (4.18)

    Inserting (4.17) and (4.18) into (4.16), we get

    F(r,t)c6[rrF(r,t)]2β+12(β+1)+1m+1+c7[rrF(r,t)]2[2β+12(β+1)+1m+1], rR0, (4.19)

    where c6=[c1+c31n+11m+12β+12(β+1)1m+1]2β(m+3)+4(2β+1)(m+2)+1+c32β+12(β+1)1m+12β+12(β+1)1m+1,c7=c3m+12(β+1)(2β+1)(m+2)+1. From (4.19) we obtain

    rrF(r,t)[F(r,t)c7+c264c27c62c7]12β+12(β+1)+1m+1, rR0

    or

    2c7F(r,t)c7+c264c27[F(r,t)c7+c264c27c62c7]12β+12(β+1)+1m+1d{F(r,t)c7+c264c27c62c7}1r, rR0 (4.20)

    Integrating (4.20) from R0 to r, we obtain

    2c7(2β+12(β+1)+1m+1)β2(β+1)+2m+1[F(r,t)c7+c264c27c62c7]β2(β+1)+2m+12β+12(β+1)+1m+12c7(2β+12(β+1)+1m+1)β2(β+1)+2m+1[F(R0,t)c7+c264c27c62c7]β2(β+1)+2m+12β+12(β+1)+1m+1c6(2β+12(β+1)+1m+1)12(β+1)1m+1[F(r,t)c7+c264c27c62c7]12(β+1)1m+12β+12(β+1)+1m+1+c6(2β+12(β+1)+1m+1)12(β+1)1m+1[F(R0,t)c7+c264c27c62c7]12(β+1)1m+12β+12(β+1)+1m+1ln(R0r). (4.21)

    Dropping the first and fourth terms on the left of (4.21), we obtain

    F(r,t)c7+c264c27[c8ln(rR0)Q2(R0,t)]2β+12(β+1)+1m+112(β+1)1m+1+c62c7, (4.22)

    where Q2(R0,t)=2c7(12(β+1)1m+1)c6(β2(β+1)+2m+1)[F(R0,t)c7+c264c27c62c7]β2(β+1)+2m+12β+12(β+1)+1m+1 and c8=12(β+1)1m+1c6(2β+12(β+1)+1m+1).

    Squaring (4.22) we have

    F(r,t)c7[c8ln(rR0)Q2(R0,t)]2β+1β+1+2m+112(β+1)1m+1+c6[c8ln(rR0)Q2(R0,t)]2β+12(β+1)+1m+112(β+1)1m+1. (4.23)

    From (4.15) and (4.23) we can obtain the following theorem.

    Theorem 4.1. Let (u,v) be the solution of (1.1), (1.2), (2.1) and (2.2) with ρ(q2)=b1+b2q2β, where 12(β+1)>max{1m+1,1n+1}. Then for fixed t, when r, (u,v) either grows algebraically or decays logarithmically. The growth rate is at least as fast as z1c5 and the decay rate is at least as fast as (lnr)2β+12(β+1)+1m+112(β+1)1m+1.

    Remark 4.1. Obviously, in this case of ρ(q2)=b1+b2q2β, the decay rate obtained by Theorem 4.1 is slower than that obtained by Theorem 2.1 and Theorem 3.1.

    In this section, we concern with a system of two coupled viscoelastic equations

    uttΔu+t0h1(tη)Δu(η)dη+f1(u,v)=0, (5.1)
    vttΔv+t0h2(tη)Δv(η)dη+f2(u,v)=0, (5.2)

    which describes the interaction between two different fields arising in viscoelasticity. In (5.1) and (5.2), 0<t<T and h1,h2 are differentiable functions satisfying h1(0),h2(0)>0 and

    2(T0h21(Tη)dη), 4Th1(0)[T0(h1(Tη)ωh1(Tη))2dη]h1(0), (5.3)
    2(T0h22(Tη)dη), 4Th2(0)[T0(h2(Tη)ωh2(Tη))2dη]h2(0). (5.4)

    Messaoudi and Tatar [24] considered the system (5.1) and (5.2) in a bounded domain and proved the uniform decay for the solution when t. For more special cases, one can refer to [25,26,27]. They mainly concerned the well-posedness of the solutions and proved that the solutions decayed uniformly under some suitable conditions. However, the present paper extends the previous results to Eqs (5.1) and (5.2) in an exterior region. We consider Eqs (5.1) and (5.2) with the initial-boundary conditions (2.1) and (2.2) in Ω.

    We define two functions

    G1(r,t)=t0B(r)eωηuxruηdsdηt0B(r)eωη(η0h1(ηs)uds)xruηdsdηI1+I2, (5.5)
    G2(r,t)=t0B(r)eωηvxrvηdsdηt0B(r)eωη(η0h2(ηs)vds)xrvηdsdηJ1+J2. (5.6)

    Integrating (5.5) from r0 to r and using (5.1), (5.2), (2.1), (2.2) and the divergence theorem, we have

    G1(r,t)=G1(r0,t)+12rr0B(ξ)eωt[|ut|2+|u|2]dsdξ+12ωt0rr0B(ξ)eωη|uη|2dsdξdη+12ωt0rr0B(ξ)eωη|u|2dsdξdη+t0rr0B(ξ)eωηh1(0)|u|2dsdξdηrr0B(ξ)eωt(t0h1(tτ)udτ)udsdξ+t0rr0B(ξ)eωη[η0(h1(ητ)ωh1(ητ))udτ]udsdξdη+t0rr0B(ξ)eωηf1(u,v)uηdsdξdη. (5.7)

    From (5.7) it follows that

    zG1(r,t)=12B(r)eωt[|ut|2+|u|2]ds+12ωt0B(r)eωη|uη|2dsdη+12ωt0B(r)eωη|u|2dsdη+t0B(r)eωηh1(0)|u|2dsdηB(r)eωt(t0h1(tτ)udτ)uds+t0B(r)eωη[η0(h1(ητ)ωh1(ητ))udτ]udsdη+t0B(r)eωηf1(u,v)uηdsdη. (5.8)

    By the Young inequality and the Hölder inequality, we have

    \begin{align} \Big|&-\int_{B(r)}e^{-\omega t}\Big(\int_0^t h_1(\eta-\tau)\nabla ud\tau\Big)\cdot\nabla uds\Big|\\ &\leq\int_{B(r)}e^{-\omega t}\Big[\Big(\int_0^t h_1(t-s)\nabla u(s)ds\Big)^2+\frac{1}{4}|\nabla u|^2\Big]ds\\ &\leq\Big(\int_0^t h_1^2(t-\tau)d\tau\Big)\int_0^t\int_{B(r)}e^{-\omega \eta}|\nabla u|^2dsd\eta+\frac{1}{4}\int_{B(r)}e^{-\omega t}|\nabla u|^2ds, \end{align} (5.9)

    and

    \begin{align} \Big|&\int_0^t\int_{B(r)}e^{-\omega \eta}\Big[\int_0^\eta \Big(h_1'(\eta-\tau)-\omega h_1(\eta-\tau)\Big)\nabla ud\tau\Big]\cdot\nabla uds d\eta\Big|\\ &\leq \frac{t}{h_1(0)}\Big[\int_0^t \Big(h_1'(\eta-\tau)-\omega h_1(\eta-\tau)\Big)^2d\tau\Big]\int_0^t\int_{B(r)}e^{-\omega \eta}|\nabla u|^2dsd\eta\\ &+\frac{1}{4}\int_0^t\int_{B(r)}e^{-\omega \eta}h_1(0)|\nabla u|^2dsd\eta. \end{align} (5.10)

    Inserting (5.9) and (5.10) into (5.8) and using (5.3), we have

    \begin{align} \frac{\partial}{\partial z}G_1(r, t)&\geq\frac{1}{2}\int_{B(r)}e^{-\omega t}\Big[|u_t|^2+\frac{1}{2}|\nabla u|^2\Big]ds+\frac{1}{2}\omega\int_0^t\int_{B(r)}e^{-\omega \eta}|u_\eta|^2ds d\eta\\ &+\frac{1}{2}\omega\int_0^t\int_{B(r)}e^{-\omega \eta}|\nabla u|^2dsd\eta+\int_0^t\int_{B(r)}e^{-\omega \eta}f_1(u, v)u_\eta ds d\eta. \end{align} (5.11)

    Similar to (5.11), we also have for G_2(r, t)

    \begin{align} \frac{\partial}{\partial z}G_2(r, t)&\geq\frac{1}{2}\int_{B(r)}e^{-\omega t}\Big[|v_t|^2+\frac{1}{2}|\nabla v|^2\Big]ds+\frac{1}{2}\omega\int_0^t\int_{B(r)}e^{-\omega \eta}|v_\eta|^2dsd\eta\\ &+\frac{1}{2}\omega\int_0^t\int_{B(r)}e^{-\omega \eta}|\nabla v|^2dsd\eta+\int_0^t\int_{B(r)}e^{-\omega \eta}f_2(u, v)v_\eta dsd\eta. \end{align} (5.12)

    If we define

    G(r, t) = G_1(r, t)+G_2(r, t),

    then by (5.11) and (5.12) we have

    \begin{align} \frac{\partial}{\partial z}G(r, t)&\geq\frac{1}{2}\int_{B(r)}e^{-\omega t}\Big[|u_t|^2+|v_t|^2+\frac{1}{2}|\nabla u|^2+\frac{1}{2}|\nabla v|^2+2F(u, v)\Big]ds \\ &+\frac{1}{2}\omega\int_0^t\int_{B(r)}e^{-\omega \eta}\Big[|u_\eta|^2+|v_\eta|^2+|\nabla u|^2+|\nabla v|^2+2F(u, v)\Big]dsd\eta. \end{align} (5.13)

    On the other hand, we bound G(r, t) by \frac{\partial}{\partial r}G(r, t) . Using the Hölder inequality, the AG mean inequality, (5.3) and combining (5.13), we have

    \begin{align} \Big|I_1\Big|+\Big|J_1\Big|&\leq\Big(\int_0^t\int_{B(r)}e^{-\omega\eta}|\nabla u|^2 dAd\eta\cdot\int_0^t\int_{B(r)}e^{-\omega\eta}|u_\eta|^2 dsd\eta\Big)^\frac{1}{2} \\ &+\Big(\int_0^t\int_{B(r)}e^{-\omega\eta}|\nabla v|^2 dAd\eta\cdot\int_0^t\int_{B(r)}e^{-\omega\eta}|v_\eta|^2 dsd\eta\Big)^\frac{1}{2} \\ &\leq\frac{1}{\omega}\Big[\frac{\partial}{\partial r}G(r, t)\Big], \end{align} (5.14)

    and

    \begin{align} \Big|I_2\Big|&\leq\Big(\int_0^t\int_{B(r)}e^{-\omega\eta}\Big|\int_0^\eta h_1(\eta-\tau)\nabla ud\tau\Big|^2 dsd\eta\cdot\int_0^t\int_{B(r)}e^{-\omega\eta}|u_\eta|^2dsd\eta\Big)^\frac{1}{2} \\ &\leq\Big(\int_0^t\int_{B(r)}e^{-\omega\eta}\Big(\int_0^\eta h^2_1(\eta-\tau)d\tau\Big)\Big(\int_0^\eta |\nabla u|^2d\tau\Big) dsd\eta\cdot\int_0^t\int_{B(r)}e^{-\omega\eta}|u_\eta|^2dsd\eta\Big)^\frac{1}{2} \\ &\leq t\Big(\int_0^t h^2_1(\eta-\tau)d\tau\Big)^\frac{1}{2}\Big(\int_0^t\int_{B(r)}e^{-\omega\eta} |\nabla u|^2 dsd\eta\cdot\int_0^t\int_{B(r)}e^{-\omega\eta}|u_\eta|^2dsd\eta\Big)^\frac{1}{2} \\ &\leq\frac{t}{2}\Big(\int_0^t h^2_1(\eta-\tau)d\tau\Big)^\frac{1}{2}\Big[\int_0^t\int_{B(r)}e^{-\omega\eta} |\nabla u|^2 dsd\eta+\int_0^t\int_{B(r)}e^{-\omega\eta}|u_\eta|^2dsd\eta\Big]\\ &\leq\frac{T}{2\omega}\sqrt{h_1(0)}\Big[\frac{\partial}{\partial r}G(r, t)\Big]. \end{align} (5.15)

    Similar to (5.15), we have

    \begin{align} \Big|J_2\Big|\leq\frac{T}{2\omega}\sqrt{h_2(0)}\Big[\frac{\partial}{\partial r}G(r, t)\Big]. \end{align} (5.16)

    Inserting (5.14)–(5.16) into (5.5) and (5.6), we have

    \begin{align} \Big|G(r, t)\Big|\leq c_9\frac{1}{\omega}\Big[\frac{\partial}{\partial r}G(r, t)\Big], \end{align} (5.17)

    where c_9 = \frac{T}{2}(\sqrt{h_1(0)}+\sqrt{h_2(0)})+2 .

    We can follow the similar arguments given in the previous sections to obtain the following theorem.

    Theorem 5.1. Let (u, v) be the solution of (5.1), (5.2), (2.1) and (2.2) in \Omega , and (5.3) and (5.4) hold. For fixed t ,

    (1) If \exists\ R_0\geq r_0 , G(R_0, t)\geq0 , then

    G(r, t)\geq G(R_0, t)e^{\frac{\omega}{b_9}(r-R_0)}.

    (2) If \forall\ r\geq r_0 , G(r, t) < 0 , then

    -G(r, t)\leq \Big[-G(r_0, t)\Big]e^{-\frac{\omega}{b_9}(r-r_0)}.

    Again, the rate of growth or decay obtained in this case is arbitrarily large

    Remark 5.1. It is clear that the above analysis can be adapted without difficulties to the equation (see [28,29])

    u_{tt}-k_0\triangle u+\int_0^tdiv[a(x)h(t-s)\nabla u(s)]ds+b(x)g(u_t)+f(u) = 0

    and the equation (see [30])

    |u_{t}|^\sigma u_{tt}-k_0\triangle u-\Delta u_{tt}+\int_0^th(t-s)\Delta u(s)]ds-\gamma\Delta u_t = 0

    with some suitable g and a(x)+b(x)\geq b_{10} > 0 and k_0, \sigma, \gamma > 0.

    In this paper, we have considered several situations where the solutions of Eqs (1.1) and (1.2) either grow or decay exponentially or polynomially. We emphasize that the Poincaré inequality on the cross sections is not used in this paper. Thus, our results also hold for the two-dimensional case. On the other hand, there are some deeper problems to be studied in this paper. We can continue to study the continuous dependence of coefficients in the equation as that in [31]. These are the issues we will continue to study in the future.

    The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.

    The authors express their heartfelt thanks to the editors and referees who have provided some important suggestions. This work was supported by the Guangzhou Huashang College Tutorial Scientific Research Project (2022HSDS09), National Natural Science Foundation of China (11371175) and the Research Team Project of Guangzhou Huashang College (2021HSKT01).

    The authors declare that there is no conflict of interest.

    [1] Andersen TG, Bollerslev T (1997) Intraday periodicity and volatility persistence in financial markets. J Empir Financ 4: 115–158. doi: 10.1016/S0927-5398(97)00004-2
    [2] Balcilar M, Elie B, Rangan G, et al. (2017). Can volume predict Bitcoin returns and volatility? A quantiles-based approach. Econ Model 64: 74–81.
    [3] Blau BM (2017) Price dynamics and speculative trading in bitcoin. Res in Int Business and Financ 41: 493–499. doi: 10.1016/j.ribaf.2017.05.010
    [4] Baur DG, Hong K, Lee AD (2017) Medium of exchange or speculative Assets? Working paper, SSRN.
    [5] Bollerslev T (1986) Generalized autoregressive conditional heteroscedasticity. J Econom 3:307–327.
    [6] Bouri E, Jalkh N, Molnar P, et al. (2017a) Bitcoin for energy commodities before and after the December 2013 crash: diversifier, hedge or safe haven? Applied Economic 49: 5063–5073.
    [7] Bouri E, Rangan G, Aviral KT, et al. (2017b) Does Bitcoin hedge global uncertainty? Evidence from wavelet-based quantile-in-quantile regressions. Financ Res Lett 23: 87–95.
    [8] Carrick J (2016) Bitcoin as a Complement to Emerging Market Currencies. Emerg Mark Financ Trade 52: 2321–2334. doi: 10.1080/1540496X.2016.1193002
    [9] Cabedoa JD, Moyab I (2003) Estimating oil price 'Value at Risk' using the historical simulation approach. Energy Econ 25: 239–253.
    [10] Chen Y (2018) Blockchain tokens and the potential democratization of entrepreneurship and innovation. Bus Horiz 61: 567–575. doi: 10.1016/j.bushor.2018.03.006
    [11] Cocco L, Marchesi M (2016) Modeling and Simulation of the Economics of Mining in the Bitcoin Market. PloS ONE 11: e0164603. doi: 10.1371/journal.pone.0164603
    [12] Davidson S, Filippi PD, Potts J (2018) Blockchains and the economic institutions of capitalism. J Inst Econ 14: 639–658.
    [13] De Bondt WFM, Thaler RH (1995) Financial decision-making in markets and firms: A behavioral perspective. Handbooks in operations res and management sci 9: 385–410. doi: 10.1016/S0927-0507(05)80057-X
    [14] Dyhrberg AH, Foley S, Svec J (2018) How investible is Bitcoin? Analyzing the liquidity and transaction costs of Bitcoin markets. Econ Lett 171: 140–143.
    [15] Feng WJ, Wang YM, Zhang ZJ (2018) Informed trading in the Bitcoin market. Financ Res Lett 26: 63–70. doi: 10.1016/j.frl.2017.11.009
    [16] Garcia D, Tessone C, Mavrodiev P, et al. (2014) The digital traces of bubbles: feedback cycles between socio-economic signals in the Bitcoin economy. J of the Royal Society Interface 11: 1–8.
    [17] Gkillas K, Katsiampa P (2018). An application of extreme value theory to cryptocurrencies. Econ Lett, 164: 109–111. doi: 10.1016/j.econlet.2018.01.020
    [18] Glaser F, Zimmarmann K, Haferhorn M, et al. (2014) Bitcoin-Asset or currency? Revealing users' hidden intentions. ECIS 2014, Tel Aviv.
    [19] Grinberg R (2012) Bitcoin: An Innovative Alternative Digital Currency. Hastings Sci & Technology Law J 4: 159–207.
    [20] Groshoff D (2014) Kickstarter My Heart: Extraordinary Popular Delusions and the Madness of Crowdfunding Constraints and Bitcoin Bubbles. William & Mary Business Law Rev 5: 489–557.
    [21] Gunter FR (2017) Corruption, costs, and family: Chinese capital flight, 1984–2014. China Econ Rev 43: 105–117. doi: 10.1016/j.chieco.2017.01.010
    [22] Hamilton JD (2008) Regime switching models. In: Durlauf, N., & Blume, L.E. (Eds.), The New Palgrave Dictionary of Economics, 2nd edition. London: Palgrave Macmillan.
    [23] Hayes AS (2016) Cryptocurrency value formation: an empirical study leading to a cost of production model for valuing bitcoin. Telematics and Informatics 34: 1308–1321.
    [24] Jacobs E (2011) Bitcoin: A Bit Too Far? J of IntBanking and Commerc 16: 1–4.
    [25] Kristoufek L (2015) What are the main drivers of the bitcoin price? evidence from wavelet coherence analysis. PLOS ONE 10: e0123923.
    [26] Kristoufek L (2013) Bitcoin meets Google trends and wikipedia: quantifying the relationship between phenomena of the Internet era. Scientific Reports 3:3415. doi: 10.1038/srep03415
    [27] Lan J, Timothy L, Tu Z (2016) Capital Flight and Bitcoin Regulation. Int Rev Financ 16: 445–455. doi: 10.1111/irfi.12072
    [28] Llorente G, Michaely R, Saar G, et al. (2002) Dynamic volume-Return relation of individual stocks. Rev Financ Stud 15: 1005–1047. doi: 10.1093/rfs/15.4.1005
    [29] Li Z, Chen S, Chen S (2017) Statistical Measure of Validity of Financial Resources Allocation. EURASIA J of Mathematics, Science and Technology Education 13: 7731–7741.
    [30] Lukáš P, Taisei K (2017) Volatility Analysis of Bitcoin Price Time Series. Quantitative Financ and Econ 1: 474–485. doi: 10.3934/QFE.2017.4.474
    [31] Luther WJ (2016) Bitcoin and the Future of Digital Payments. The Independent Rev 20: 397–404.
    [32] Mertzanis C (2018) Complexity, big data and financial stability. Quantitative Financ and Econ 2: 637–660. doi: 10.3934/QFE.2018.3.637
    [33] Murad Z, Sefton M, Starmer C (2016) How do risk attitudes affect measured confidence? J Risk Uncertain 52: 21–46. doi: 10.1007/s11166-016-9231-1
    [34] Nelson DB (1991) Conditional Heteroskedasticity in Asset Returns: A New Approach. Econometrica 59: 347–370. doi: 10.2307/2938260
    [35] Osterrieder J, Lorenz J (2017) A statistical risk assessment of Bitcoin and its extreme tail behavior. Annals of Financ Econ 12: UNSP 1750003.
    [36] Pilkington M (2013) Bitcoin and Complexity Theory: Some Methodological Implications. Working Paper, University of Burgundy, France.
    [37] Plassaras NA (2013) Regulating Digital Currencies: Bringing Bitcoin within the Reach of the IMF. Chicago J of Int Law 14: 377–407.
    [38] Rockafellar RT, Uryasev S (2002) Conditional value-at-risk for general loss distributions. J of Banking & Financ 26: 1443–1471.
    [39] Sapuric S, Kokkinaki A (2014) Bitcoin Is Volatile! Isn't that Right?. In: Abramowicz W., Kokkinaki A. (eds) Business Information Systems Workshops. BIS 2014. Lecture Notes in Business Information Processing, vol 183. Springer, Cham.
    [40] Sornette D, Cauwels P, Smilyanov G (2018) Can we use volatility to diagnose financial bubbles? lessons from 40 historical bubbles. Quantitative Financ and Econ 2: 1–105.
    [41] Yermack D (2013) Is Bitcoin A Real Currency? An Economic Appraisal. Working paper, National Bureau of Economic Research.
    [42] Yermack D (2015) Bitcoin, innovation, financial instruments, and big data. Handbook of Digital Currency 31–43.
  • This article has been cited by:

    1. S.A. Pedro, J.M. Tchuenche, HIV/AIDS dynamics: Impact of economic classes with transmission from poor clinical settings, 2010, 267, 00225193, 471, 10.1016/j.jtbi.2010.09.019
    2. C. P. Bhunu, S. Mushayabasa, H. Kojouharov, J. M. Tchuenche, Mathematical Analysis of an HIV/AIDS Model: Impact of Educational Programs and Abstinence in Sub-Saharan Africa, 2011, 10, 1570-1166, 31, 10.1007/s10852-010-9134-0
    3. Sarada S. Panchanathan, Diana B. Petitti, Douglas B. Fridsma, The development and validation of a simulation tool for health policy decision making, 2010, 43, 15320464, 602, 10.1016/j.jbi.2010.03.013
    4. Cristiana J. Silva, Delfim F.M. Torres, A SICA compartmental model in epidemiology with application to HIV/AIDS in Cape Verde, 2017, 30, 1476945X, 70, 10.1016/j.ecocom.2016.12.001
    5. S. D. Hove-Musekwa, F. Nyabadza, H. Mambili-Mamboundou, Modelling Hospitalization, Home-Based Care, and Individual Withdrawal for People Living with HIV/AIDS in High Prevalence Settings, 2011, 73, 0092-8240, 2888, 10.1007/s11538-011-9651-7
    6. Emmanuelina L. Kateme, Jean M. Tchuenche, Senelani D. Hove-Musekwa, HIV/AIDS Dynamics with Three Control Strategies: The Role of Incidence Function, 2012, 2012, 2090-5572, 1, 10.5402/2012/864795
    7. F. Nyabadza, Z. Mukandavire, S.D. Hove-Musekwa, Modelling the HIV/AIDS epidemic trends in South Africa: Insights from a simple mathematical model, 2011, 12, 14681218, 2091, 10.1016/j.nonrwa.2010.12.024
    8. Xiaodong Wang, Chunxia Wang, Kai Wang, Extinction and persistence of a stochastic SICA epidemic model with standard incidence rate for HIV transmission, 2021, 2021, 1687-1847, 10.1186/s13662-021-03392-y
  • Reader Comments
  • © 2018 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(7718) PDF downloads(953) Cited by(55)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog