Skip to main content

Upper bounds for the decay rate in a nonlocal p-Laplacian evolution problem

Abstract

We obtain upper bounds for the decay rate for solutions to the nonlocal problem t u(x,t)= R n J(x,y) | u ( y , t ) u ( x , t ) | p 2 (u(y,t)u(x,t))dy with an initial condition u 0 L 1 ( R n ) L ( R n ) and a fixed p>2. We assume that the kernel J is symmetric, bounded (and therefore there is no regularizing effect) but with polynomial tails, that is, we assume a lower bounds of the form J(x,y) c 1 | x y | ( n + 2 σ ) , for |xy|> c 2 and J(x,y) c 1 , for |xy| c 2 . We prove that u ( , t ) L q ( R n ) C t n ( p 2 ) n + 2 σ ( 1 1 q ) for q1 and t large.

MSC: 35K05, 45P05, 35B40.

1 Introduction

In this paper we deal with nonlocal Cauchy problems of the form

t u(x,t)= R n J(x,y)|u(y,t)u(x,t) | p 2 ( u ( y , t ) u ( x , t ) ) dy
(1.1)

for t R + and x R n with n2, a fixed p>2 and an initial condition u(x,0)= u 0 (x) satisfying u 0 L 1 ( R n ) L ( R n ). As regards the kernel J, we will always assume that it is a bounded and symmetric function defined for (x,y) R n × R n together with the integrability condition J(,y) L 1 ( R n ) for all y R n . Under these hypotheses existence and uniqueness of a solution follow from a fixed point argument as in [1].

Nonlocal problems have been recently widely used to model diffusion processes (see [2] and [3] for a general nonlocal vector calculus). Problem (1.1) and its stationary version have been considered recently in connection with real applications, for example to peridynamics or a recent model for elasticity. We quote for instance [48] and the recent book [1].

Our main goal here is to obtain upper bounds for the asymptotic behavior of the solution of (1.1) as t+. It is expected that the diffusive nature of the equation implies that the solution goes to zero when t+.

To obtain our results the key assumptions are the following lower bounds for J:

J ( x , y ) c 1 | x y | ( n + 2 σ ) for  | x y | > c 2 and J ( x , y ) c 1 for  | x y | c 2
(1.2)

for certain constants c 1 , c 2 >0, and σ(0,1). For simplicity we will assume c 2 =1.

The main result of this paper reads as follows.

Theorem 1.1Letn2, q[1,+)andσ(0,1). LetJbe a kernel satisfying (1.2). Then the solution of (1.1) associated to an initial condition u 0 L 1 ( R n ) L ( R n )decays in L q ( R n )with the upper bound

u ( , t ) L q ( R n ) C t n ( p 2 ) n + 2 σ ( 1 1 q ) ,
(1.3)

where the constantCdepends on u 0 , q, σ, andn.

Let us end the introduction with some comments on the previous bibliography. For the linear case, p=2, and for smooth kernels J with compact support, it is proven in [9] that the solution u of (1.1) has the decay estimate

u ( , t ) L q ( R n ) C t n 2 ( 1 1 q )

for any q[1,). Note that this decay rate is the same as the one that holds for solutions of the classical Heat equation. In the case of an equation in convolution form, that is, when J(x,y)=K(xy) with K a nonnegative radial function, not necessarily compactly supported, it is proven in [10] that the solutions of equations with the form (1.1) have the decay estimate

u ( , t ) L q ( R n ) C t n 2 σ ( 1 1 q ) ,

provided the function K has a Fourier transform satisfying the expansion K ˆ (ξ)=1A | ξ | 2 σ +o( | ξ | 2 σ ) for ξ0, where A>0 is a constant. In this case the decay estimate is analogous to the one for the σ-order fractional heat equation, v t = ( Δ ) σ v, with σ(0,1). We also note that the convolution form of the equation allows the use of Fourier analysis to obtain this result. However, the use of Fourier analysis is not helpful here due to fact that our operator is not in convolution form. In spite of this difficulty, energy methods can be applied; see [9, 11]. We also mention the recent reference [12], where similar estimates can be found for nonlocal equations with unbounded kernels. We borrow ideas and techniques from these references. In particular we use Proposition 3.2 of [11] (whose proof is included here for completeness). However, we have to point out that in [12] and [11] only the linear case, that is, p=2, was treated, while here we deal with (1.1) for any p2. For examples of kernels with exponential decay bounds we refer to [13] and [14]. Finally, we point out that our results are also valid for unbounded kernels (in the spirit of [12]) since only lower bounds are assumed in (1.2).

The case 1p<2 remains open as well as the corresponding estimate for the L -norm.

2 Basic facts and preliminaries

First, we need to introduce fractional Sobolev spaces and their seminorms, we refer to [15] for details. For σ(0,1) and r[1,), W σ , r ( R n ) is the fractional Sobolev space of all L r ( R n ) functions with finite fractional seminorm [ v ] σ , r , given by

[ v ] σ , r r = R 2 n | v ( x + z ) v ( x ) | r | z | n + r σ dxdz.
(2.1)

Under these definitions, we have the following fractional Sobolev-type inequality: there exists a constant C>0 such that, for each v W σ , r ( R n ) with σr<n, we have

v L s ( R n ) r C [ v ] σ , r r ,
(2.2)

where s=nr/(nσr) (see [15]).

First, we consider a positive smooth function ψ: R n R with the following properties:

supp(ψ) B 1 and R n ψ(x)dx=1.
(2.3)

With the aid of this function, we split a function u into two parts. We will denote the ‘smooth’ part of u as v and the remaining as w. We let

v(x,t):= R n ψ(xz)u(z,t)dz;w(x,t):=u(x,t)v(x,t).
(2.4)

Sometimes, for simplicity in the notation and where the context is clear, we will write u, v, and w as functions depending only of x.

As a first property of this decomposition we find that each L r norm of the functions v and w is controlled by the corresponding norm of u.

Lemma 2.1Letvandwbe given by (2.4). For eachr(1,+), we have

v L r ( R n ) u L r ( R n ) and w L r ( R n ) u L r ( R n ) .

Proof We start with v. Denoting r =r/(r1) the Hölder conjugate of r and using the definition of v, we have

R n | v ( x ) | r d x = R n | R n ψ ( x y ) u ( y ) d y | r d x = R n | R n ψ ( x y ) 1 / r ψ ( x y ) 1 / r u ( y ) d y | r d x R n [ ( R n ψ ( x y ) d y ) 1 / r ( R n ψ ( x y ) | u ( y ) | r d y ) 1 / r ] r d x = R n | u ( y ) | r R n ψ ( x y ) d x d y R n | u ( y ) | r d y .

The inequality for w easily follows immediately from the triangular inequality in L r . □

Now we state a key result to get the desired estimate on the decay rate.

Proposition 2.2Letn2and letJ: R n × R n R + be a kernel satisfying (1.2), ψsatisfying (2.3), β[σ,1), andr>max{1,2β}. Then there exists a constantC>0such that, for allu L r ( R n )andv, wdefined in (2.4), we have

[ v ] 2 β r 1 , r r + w L r ( R n ) r C R 2 n J(x,y)|u(x)u(y) | r dxdy.
(2.5)

The constantCdepends onψ, β, r, andn.

Proof For the estimate concerning w, we have

R n | w ( x ) | r d x = R n | u ( x ) v ( x ) | r d x = R n | u ( x ) R n ψ ( x z ) u ( z ) d z | r d x = R n | R n ψ ( x z ) ( u ( x ) u ( z ) ) d z | r d x = R n | R n ψ ( x z ) 1 / r ψ ( x z ) 1 / r ( u ( x ) u ( z ) ) d z | r d x .

Applying Hölder’s inequality, we get

R n | w ( x ) | r d x R n ( R n ψ ( x z ) d z ) r / r ( R n ψ ( x z ) | u ( x ) u ( z ) | r d z ) d x R n R n ψ ( x z ) | u ( x ) u ( z ) | r d z d x ,

where r =r/(r1).

Since ψ is supported in B 1 , we have ψ(xz)J(x,z) for all |xz|1, and, since J verifies J(x,z) c 1 for |xz|1, there exists a constant C depending only on | ψ | such that ψ(xz)CJ(x,z). Then

w L r ( R n ) r C R 2 n J(x,y)|u(x)u(y) | r dxdy.

Now we deal with the term with v. We split the fractional seminorm as

[ v ] 2 β r 1 , r r = | x y | > 1 | v ( x ) v ( y ) | r | x y | n + 2 β d x d y + | x y | 1 | v ( x ) v ( y ) | r | x y | n + 2 β d x d y = : I ext + I int

and look at these integrals separately. For I ext , using the definition of v we have

I ext = | x y | > 1 | R n ( u ( x z ) u ( y z ) ) ψ(z)dz | r | x y | ( n + 2 β ) dxdy.

Now, we can think of the measure μ(dz)=ψ(z)dz as a probability measure (because of (2.3)), and since the function t | t | r is convex in , we can apply Jensen’s inequality on the dz-integral in right-hand side of the last expression to obtain

I ext | x y | > 1 R n |u(xz)u(yz) | r ψ(z)dz | x y | ( n + 2 β ) dxdy,

which, after an application of Fubini’s theorem, gives

I ext R n ψ(z) ( | x y | > 1 | u ( x z ) u ( y z ) | r | x y | ( n + 2 β ) d x d y ) dz.

Then, applying the change x ˜ =xz, y ˜ =yz in the dxdy integral and using (2.3), we conclude

I ext R n ψ ( z ) ( | x ˜ y ˜ | > 1 | u ( x ˜ ) u ( y ˜ ) | r | x ˜ y ˜ | ( n + 2 β ) d x ˜ d y ˜ ) d z = | x ˜ y ˜ | > 1 | u ( x ˜ ) u ( y ˜ ) | r | x ˜ y ˜ | ( n + 2 β ) d x ˜ d y ˜ .

Using this last expression, we obtain from the assumption (1.2) that

I ext C R 2 n J(x,y)|u(x)u(y) | r dxdy.
(2.6)

Now we deal with I int . In this case, using the definition of v, we can write

I int = | x y | < 1 | R n u(z) ( ψ ( x z ) ψ ( y z ) ) dz | r | x y | ( n + 2 β ) dxdy.
(2.7)

Note that by using (2.3), we have for all x,y R n

R n u(x) ( ψ ( x z ) ψ ( y z ) ) dz=u(x) ( R n ψ ( x z ) d z R n ψ ( y z ) d z ) =0,

and then

R n u(z) ( ψ ( x z ) ψ ( y z ) ) dz= R n ( u ( z ) u ( x ) ) ( ψ ( x z ) ψ ( y z ) ) dz.

Thus, using this equality with (2.7), we get

I int = | x y | < 1 | R n ( u ( z ) u ( x ) ) ( ψ ( x z ) ψ ( y z ) ) dz | r | x y | ( n + 2 β ) dxdy.

However, note that if |xz|2 in the dz integral, since |xy|<1 necessarily |yz|>1. Then, due to the fact that ψ is supported in the unit ball, the contribution of the integrand when |xz|2 is null in the dz integral. Taking this into account, applying Hölder’s inequality into the dz-integral, we have

I int = | x y | < 1 | | x z | < 2 ( u ( z ) u ( x ) ) ( ψ ( x z ) ψ ( y z ) ) d z | r × | x y | ( n + 2 β ) d x d y | x y | < 1 ( | x z | < 2 | u ( z ) u ( x ) | r d z ) × ( | x z ˜ | < 2 | ψ ( x z ˜ ) ψ ( y z ˜ ) | r d z ˜ ) r / r | x y | ( n + 2 β ) d x d y .

By Fubini’s theorem we can write

I int = x R n ( | x z | < 2 ( u ( z ) u ( x ) ) r d z ) Ψ(x)dx,

where

Ψ(x)= | x y | < 1 ( | x z ˜ | < 2 | ψ ( x z ˜ ) ψ ( y z ˜ ) | r d z ˜ ) r / r | x y | ( n + 2 β ) dy.

Using the regularity of ψ, we have

Ψ ( x ) | x y | < 1 ( | x z ˜ | < 2 D ψ r | x y | r d z ˜ ) r / r | x y | ( n + 2 β ) d y D ψ r | B 2 | r / r | x y | < 1 | x y | r | x y | ( n + 2 β ) d y ,

and since r>2β, we conclude that the last integral is convergent, obtaining

Ψ(x) C n , β , r D ψ r | B 2 | r / r ,

which leads to the following estimate for I int :

I int C x R n | x z | < 2 |u(z)u(x) | r dzdx.

From this, it is easy to get

I int C | x z | 2 | u ( z ) u ( x ) | r ( 1 + | x z | ) n + 2 β dzdx,

which, by the use of (1.2), let us conclude that

I int C R 2 n J(x,y)|u(x)u(y) | r dxdy.

This last estimate together with (2.6) concludes the proof. □

3 Proof of Theorem 1.1

As mentioned in the introduction, existence and uniqueness of solutions to problem (1.1) follows as in [1]. In fact, the symmetry, boundedness, and integrability assumptions over J allow us to perform a fixed point argument to obtain the following result, whose proof is omitted.

Theorem 3.1Let u 0 L 1 ( R n ) L ( R n ), then there exists a unique solutionuC([0,+), L 1 ( R n ) L ( R n ))of equation (1.1). This solution satisfies u ( , t ) L 1 ( R n ) u 0 L 1 ( R n ) and u ( , t ) L ( R n ) u 0 L ( R n ) for allt0.

Now, let us introduce the main idea behind the energy methods. To clarify the exposition, let us perform these computations in the local case and next see how we can adapt them to our nonlocal problem with the help of Proposition 2.2. Let us describe briefly how the energy method can be applied to obtain decay estimates for local problems. Let us begin with the simpler case of the estimate for solutions to the p-Lapacian evolution equation in L 2 -norm. Let u be a solution to

t u= Δ p u.

If we multiply the equation by u and integrate in R n , we obtain

t R n 1 2 u 2 (x,t)dx= R n |u(x,t) | p dx.

Now we use Sobolev’s inequality

R n |u(x,t) | p dxC ( R n | u ( x , t ) | p d x ) p / p

with p =pn/(np) to obtain

t R n u 2 (x,t)dxC ( R n | u ( x , t ) | p d x ) p / p .

If we use interpolation and that u ( , t ) L 1 ( R n ) C( u 0 ) for any t>0, we have

u ( , t ) L 2 ( R n ) u ( , t ) L 1 ( R n ) α u ( , t ) L p ( R n ) 1 α C u ( , t ) L p ( R n ) 1 α

with α determined by

1 2 =α+ 1 α p ,that is, α= ( 1 2 1 p ) p ( p 1 ) .

Hence we get

t R n u 2 (x,t)dxC ( R n u 2 ( x , t ) d x ) 1 1 α

from where the decay estimate

u ( , t ) L 2 ( R n ) C t 1 2 ( n n ( p 2 ) + 2 ) ,t>0,

follows. To obtain a decay bound for u ( , t ) L q ( R n ) we can use the same idea multiplying by u q 1 at the beginning.

Now we are ready to proceed with the proof of our main result.

Proof of Theorem 1.1 The symmetry assumption on J allows us to mimic this idea and use an energy approach in order to get Theorem 1.1. Roughly speaking, this assumption allows us to ‘integrate by parts’ (1.1). For q=1 the proof is finished by Theorem 3.1. For q>1 we multiply the equation by q | u | q 2 u and integrate, obtaining the identity

t R n | u ( x ) | q d x = q 2 R 2 n J ( x , y ) | u ( y ) u ( x ) | p 2 ( u ( y ) u ( x ) ) × ( | u ( y ) | q 2 u ( y ) | u ( x ) | q 2 u ( x ) ) d y d x ,
(3.1)

where we omitted the dependence on t of the function u for simplicity.

Now we recall the following inequality (whose proof is straightforward): let q>1 and a,b0. Then there exists a constant C depending only on q, such that

(ab) ( | a | q 2 a | b | q 2 b ) C | a b | q .

Hence, using this inequality with (3.1), we conclude

t u L q ( R n ) q C R 2 n J(x,y)|u(y)u(x) | p 2 + q dydx=:CE(u).
(3.2)

Note that we find that the L q -norm of u is decreasing in t. At this point we would like to use Sobolev’s inequality, which is not available due to the lack of regularizing effect of our nonlocal operator. Instead we will use Proposition 2.2, which involves a good control of the smooth part v (but we have to take care of the rough part w).

Let us fix q2 (the case 1<q<2 will be tackled at the end of this proof using interpolation). By the definition of v and w in (2.4), we have

u L q ( R n ) q 2 q 1 ( v L q ( R n ) q + w L q ( R n ) q ) .
(3.3)

Now we note that v belongs to L p for all p. Hence, we can interpolate, obtaining

v L q ( R n ) q v L s ( R n ) q θ v L 1 ( R n ) q ( 1 θ ) ,

with

s= n ( p + q 2 ) n 2 σ ,

where θ is given by

1 q = θ s +(1θ),that is, θ= s ( q 1 ) q ( s 1 ) .

Recalling v L 1 ( R n ) u ( , t ) L 1 ( R n ) u 0 L 1 ( R n ) and the Sobolev-type inequality (2.2), we obtain

v L q ( R n ) q C [ v ] σ ˜ , p + q 2 q θ ,
(3.4)

where σ ˜ =2σ ( p + q 2 ) 1 and the constant C depends on u 0 , q, σ, and n.

Concerning w we can also interpolate and obtain

w L q ( R n ) q w L p + q 2 ( R n ) q γ w L 1 ( R n ) q ( 1 γ ) ,

with γ given by

1 q = γ p + q 2 +(1γ),that is, γ= ( p + q 2 ) ( q 1 ) q ( p + q 3 ) .

Note that we are using that p>2 here. Now we use

w L 1 ( R n ) u ( , t ) L 1 ( R n ) u 0 L 1 ( R n )

to get

w L q ( R n ) q C w L p + q 2 ( R n ) q γ ,
(3.5)

with C depending on u 0 , q, σ, and n.

From (3.3), (3.4), and (3.5) we obtain

u L q ( R n ) q C [ v ] σ ˜ , p + q 2 q θ +C w L p + q 2 ( R n ) q γ .
(3.6)

Now we use Proposition 2.2, with r=p+q2 and β=σ, to obtain

v L s ( R n ) q θ C ( E ( u ) ) q θ p + q 2 and w L p + q 2 ( R n ) q γ C ( E ( u ) ) q γ p + q 2

and we conclude that

u L q ( R n ) q C ( E ( u ) ) q θ p + q 2 +C ( E ( u ) ) q γ p + q 2 ,

that is,

H 1 ( u L q ( R n ) q ) E(u)

with H(z)=C z q θ p + q 2 +C z q γ p + q 2 . Since u L q ( R n ) q (t) u 0 L q ( R n ) q (recall that the L q -norm of the solution decreases) and q θ p + q 2 < q γ p + q 2 , we have

H 1 ( u L q ( R n ) q ) C ( u L q ( R n ) q ) p + q 2 q θ .

Then, from (3.2), we obtain

t u ( , t ) L q ( R n ) q CE(u)C H 1 ( u L q ( R n ) q ) C ( u L q ( R n ) q ) p + q 2 q θ ,

from which it follows that

u ( , t ) L q ( R n ) q C t q θ p + q 2 q θ ,

that is,

u ( , t ) L q ( R n ) q C t ( q 1 ) n q ( ( p 2 ) n + 2 σ ) ,

as we wanted to show.

Now we deal with the case 1<q<2. We can interpolate, obtaining

u ( , t ) L q ( R n ) q u ( , t ) L 2 ( R n ) q θ u ( , t ) L 1 ( R n ) q ( 1 θ ) ,

with

1 q = θ 2 +(1θ),that is, θ=2 ( 1 1 q ) .

As we have already proved that

u ( , t ) L 2 ( R n ) C t n ( p 2 ) n + 2 σ 1 2 and u ( , t ) L 1 ( R n ) C,

we conclude

u ( , t ) L q ( R n ) C t n ( p 2 ) n + 2 σ ( 1 1 q )

for 1<q<2. □

References

  1. Andreu-Vaillo F, Mazon JM, Rossi JD, Toledo-Melero JJ Mathematical Surveys and Monographs 165. In Nonlocal Diffusion Problems. Am. Math. Soc., Providence; 2010.

    Chapter  Google Scholar 

  2. Fife P: Some nonclassical trends in parabolic and parabolic-like evolutions. In Trends in Nonlinear Analysis. Springer, Berlin; 2003:153-191.

    Chapter  Google Scholar 

  3. Du, Q, Gunzburger, M, Lehoucq, R, Zhou, K: A nonlocal vector calculus, nonlocal volume-constrained problems, and nonlocal balance laws. Preprint

  4. Bates P, Chen X, Chmaj A: Heteroclinic solutions of a van der Waals model with indefinite nonlocal interactions.Calc. Var. 2005, 24: 261-281. 10.1007/s00526-005-0308-y

    Article  MathSciNet  Google Scholar 

  5. Hutson V, Martínez S, Mischaikow K, Vickers GT: The evolution of dispersal.J. Math. Biol. 2003, 47: 483-517. 10.1007/s00285-003-0210-1

    Article  MathSciNet  Google Scholar 

  6. Parks ML, Lehoucq R, Plimpton S, Silling S: Implementing peridynamics within a molecular dynamics code.Comput. Phys. Commun. 2008, 179: 777-783. 10.1016/j.cpc.2008.06.011

    Article  Google Scholar 

  7. Silling SA: Reformulation of elasticity theory for discontinuities and long-range forces.J. Mech. Phys. Solids 2000, 48: 175-209. 10.1016/S0022-5096(99)00029-0

    Article  MathSciNet  Google Scholar 

  8. Silling SA, Lehoucq R: Convergence of peridynamics to classical elasticity theory.J. Elast. 2008, 93: 13-37. 10.1007/s10659-008-9163-3

    Article  MathSciNet  Google Scholar 

  9. Ignat L, Rossi JD: Decay estimates for nonlocal problems via energy methods.J. Math. Pures Appl. (9) 2009, 92(2):163-187. 10.1016/j.matpur.2009.04.009

    Article  MathSciNet  Google Scholar 

  10. Chasseigne E, Chaves M, Rossi JD: Asymptotic behavior for nonlocal diffusion equations.J. Math. Pures Appl. (9) 2006, 86(3):271-291. 10.1016/j.matpur.2006.04.005

    Article  MathSciNet  Google Scholar 

  11. Chasseigne, E, Felmer, P, Rossi, JD, Topp, E: Fractional decay bounds for nonlocal zero order heat equations. Preprint

  12. Brändle, C, de Pablo, A: Nonlocal heat equations: decay estimates and Nash inequalities. Preprint (2014). arXiv:1312.4661v3

    Google Scholar 

  13. Ignat L, Rossi JD, San Antolin A: Lower and upper bounds for the first eigenvalue of nonlocal diffusion problems in the whole space.J. Differ. Equ. 2012, 252(12):6429-6447. 10.1016/j.jde.2012.03.011

    Article  MathSciNet  Google Scholar 

  14. Ignat L, Pinasco D, Rossi JD, San Antolin A: Decay estimates for nonlinear nonlocal diffusion problems in the whole space.J. Anal. Math. 2014, 122(1):375-401. 10.1007/s11854-014-0011-z

    Article  MathSciNet  Google Scholar 

  15. Di Neza E, Palatucci G, Valdinoci E: Hitchhiker’s guide to the fractional Sobolev spaces.Bull. Sci. Math. 2012, 136(5):521-573. 10.1016/j.bulsci.2011.12.004

    Article  MathSciNet  Google Scholar 

Download references

Acknowledgements

This work was partially supported by MEC MTM2010-18128 and MTM2011-27998 (Spain).

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Angel San Antolin.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

All authors contributed equally to the writing of this paper. All authors read and approved the final manuscript.

Rights and permissions

Open Access  This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.

The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.

To view a copy of this licence, visit https://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Esteve, C., Rossi, J.D. & San Antolin, A. Upper bounds for the decay rate in a nonlocal p-Laplacian evolution problem. Bound Value Probl 2014, 109 (2014). https://doi.org/10.1186/1687-2770-2014-109

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1687-2770-2014-109

Keywords