- Research
- Open access
- Published:
Large time behavior for the fractional Ginzburg-Landau equations near the BCS-BEC crossover regime of Fermi gases
Boundary Value Problems volume 2017, Article number: 8 (2017)
Abstract
In this paper, we consider the fractional Ginzburg-Landau equations near the Bardeen-Cooper-Schrieffer-Bose-Einstein-condensate (BCS-BEC) crossover of atomic Fermi gases. This fractional Ginzburg-Landau equations can be viewed as a generalization of the integral differential equations proposed by Machida and Koyama (Phys. Rev. A 74:033603, 2006). By using the Galerkin method and a priori estimates, together with the properties of Sobolev spaces, we first establish the existence and uniqueness of weak solutions to these equations and then we prove the existence of global attractors.
1 Introduction
In this paper, we consider the fractional Ginzburg-Landau equations for atomic Fermi gases near the BCS-BEC crossover as follows:
where \(u(x,t)\) is the fermion-pair field and \(\varphi(x,t)\) is the condensed boson field, μ is the chemical potential, 2ν is the threshold energy of the Feshbach resonance, g is the coupling constant describes the process arising from the Feshbach resonance, in which a boson is created from two fermionic atoms and vice versa. d is generally complex, in the BCS limit d can be considered to be pure imaginary, while in the BEC region, the imaginary part of d vanished. \(U>0\), \(a, b,c, \beta\) are real coefficients, the square root of the Laplacian, \(\Lambda=(-\Delta)^{\frac{1}{2}}\), is the so-called Zygmund operator. In this paper, we restrict ourselves to the fractional order \(\alpha\in(\frac{1}{2},1]\) and the space dimension \(n=1\).
The BCS-BEC crossover phenomenon has been experimentally realized by using ultracold gases of \(^{6}Li\) and 40 K atoms and it has opened a new era for the study of some longstanding theoretical proposals in many fermion systems. Since the conventional perturbation theory is no longer valid, the equation of state and dynamic properties of the BCS-BEC crossover have become a big challenge for quantum theory. In recent years, a number of groups have studied the atomic Fermi gases from various views which help us to deeply understand the physics of the pseudogap and Berezinskii-Kosterlitz-Thouless transitions in fermionic systems.
In 2006, on the basis of the functional integral formalism, Machida and Koyama [1] constructed a time-dependent Ginzburg-Landau theory for the superfluid atomic Fermi gases near the Feshbach resonance from the fermion-boson model. The time-dependent Ginzburg-Landau theory just can be described as (1)-(2). In the case of \(g=0,\alpha=1\), (1)-(2) are decoupled and reduced to the conventional time-dependent Ginzburg-Landau (TDGL) equation and the linearized Gross-Pitaevskii (GP) equation, which intensively were studied in the past decades. In the case of \(g\ne0,\alpha=1\), Chen and Guo [2–4] obtained the global existence and uniqueness of weak solutions to (1)-(2) with periodic boundary conditions, and then by the properties of Besov and Sobolev spaces and matrix theory, together with the energy method, they established the global existence of classical solutions. The existence of global attractors to (1)-(4) was proved by Fang, Jin and Guo [5]. In [6], Guo et al. considered the numerical solution of Ginzburg-Landau equations near the BCS-BEC crossover through finite difference method and analyzed its convergence and stability. In [7], the Hopf bifurcation of the above equations was studied and numerical simulations were also given.
In this paper, motivated by [8] and [9], we extend the result to fractional Ginzburg-Landau equations for atomic Fermi gases near the BCS-BEC crossover. When α is not an integer, the fractional dissipation operator \(\Lambda^{2\alpha}\) is nonlocal and can be regarded as the infinitesimal generators of Lévy stable diffusion processes. More and more researchers have found that fractional differential equations play an important role in mathematical physics and can be used to describe some physical phenomenon more exactly than integral differential equations. There have been extensive study of fractional differential equations including the fractional Ginzburg-Landau equation [8], the fractional Schrödinger equation [9], the fractional Landau-Lifshitz-Gilbert equation [10], the fractional Landau-Lifshitz equation [11] etc. For more details, see [12–15].
The rest of this paper is organized as follows. In the next section, we give some notations. In Section 3, we give a priori estimates. In Section 4, we prove the global existence of weak solutions to (1)-(4). Finally, the existence of global attractors is obtained.
2 Notations
Let \(\Omega=[0,2\pi]\), \(d=d_{r}+id_{i}\), \(m=m_{r}+im_{i}\), and \(|d|^{2}=d_{r}^{2}+d_{i}^{2}\), \(|m|^{2}=m_{r}^{2}+m_{i}^{2}\). Denote by \(L^{p}(\Omega)\) the usual Sobolev space of the pth-power integrable functions normed by
When no confusion arises, we set \(L^{p}:=L^{p}(\Omega)\) for \(1\le p\le \infty\). Similarly, we define the space \(L^{p}(0,T;X)\) with the norm
If u is a periodic function, then it can be identified with the Fourier series
Moreover, \(\Lambda^{\alpha}u\) can be defined by
Define
and \(H^{\alpha}\) denotes a complete space of E with the induced norm
Then \(H^{\alpha}\) is a Banach space. It is easy to show that \(H^{\alpha }\) is a Hilbert space with the inner product
where \(\Lambda^{\alpha} v=\sum_{j\in\mathbb{Z}}|j|^{\alpha }b_{j}e^{ij\cdot x}\). For more details, see [9].
We denote positive constants by C; they may change from one line to the next line.
3 A priori estimates
Lemma 1
Assume that \(u_{0}\in L^{2}(\Omega), \varphi_{0}\in L^{2}(\Omega), f(x)\in L^{2}(\Omega), h(x)\in L^{2}(\Omega)\). Let \((u,\varphi)\) be the solution to (1)-(4). Then, for any \(0< T<\infty, U>0, c>0, m_{i}d_{r}+m_{r}d_{i}>0, b>0\), \(m_{i}>0\), \(d_{i}>0\), and \(\beta>0\), we have
where the constant C depending on the initial data and T.
Proof
Multiplying (5) by \(\overline{(u+g\varphi)}\), integrating over Ω, and taking the real part, we have
Multiplying (2) by φ̅, integrating over Ω, and taking the imaginary part, we have
Combining the above two equations, we obtain
where \(C_{1}=\max\{|\frac{g}{dU}|+\frac{|g|}{U}+|\beta g|+1,|\frac {g}{dU}|+\frac{|g|}{U}+|\beta g|+\frac{2ad_{i}}{|d|^{2}}-\frac {2d_{i}}{|d|^{2}U}+1\}\).
Noticing that \(c>0, m_{i}d_{r}+m_{r}d_{i}>0, b>0\), \(m_{i}>0\), \(d_{i}>0\), and \(\beta >0\), we have
By Gronwall’s inequality, we have
From (8), we can also deduce that
Thus we complete the proof. □
Lemma 2
Let \(u_{0}\in L^{2}(\Omega), \varphi_{0}\in L^{2}(\Omega), f\in L^{2}(\Omega), h\in L^{2}(\Omega)\). Assume that \(U>0\), \(c>0\), \(m_{i}d_{r}+m_{r}d_{i}>0\), \(b>0\), \(m_{i}>0\), \(d_{i}>0\) and \(\beta>0\), we have
Proof
It follows from Lemma 1 that
For any \(\psi\in H^{\alpha}(\Omega)\), we have
By Hölder’s inequality, we have
Applying the Sobolev embedding
we have
Therefore
Thus the proof of Lemma 2 is completed. □
Let \(I_{\phi}(t)=((u+g\varphi,\phi_{1}),(\varphi,\phi_{2}))\), \(\phi=(\phi _{1},\phi_{2})\).
Lemma 3
Under the conditions as in Lemma 1, for any \(\phi_{1}\in L^{2}(\Omega)\), \(\phi_{2}\in L^{2}(\Omega)\), \(I_{\phi}(t)\) is a continuous function with respect to t.
Proof
First, let \(\phi=(\phi_{1},\phi_{2})\in C^{\infty}(\Omega)\times C^{\infty}(\Omega)\). We have
For \(0\le t_{1}, t_{2}\le T\) and \(|t_{2}-t_{1}|<1\), we have
Then
Therefore, the continuity of \(I_{\phi}(t)\) follows.
Next we use a density argument to extend the result for \(\phi\in L^{2}(\Omega)\times L^{2}(\Omega)\). Let \(\epsilon>0\) be an arbitrary positive number, for \(\phi\in L^{2}(\Omega)\times L^{2}(\Omega)\), we may select some \(\phi^{\epsilon}\in C^{\infty}(\Omega)\times C^{\infty}(\Omega)\) such that \(\|\phi^{\epsilon}-\phi\|_{L^{2}(\Omega)}\le\epsilon\). By the triangle inequality and Hölder’s inequality, we have
Since \(I_{\phi^{\epsilon}}\) is continuous in t and ϵ is arbitrary, the continuity of \(I_{\phi}(t)\) follows for \(\phi\in L^{2}\times L^{2}\). □
4 Global existence of weak solutions
In this section, we will show the global existence of the weak solutions. We first give the definition of weak solutions.
Definition
Let \(u_{0}\in L^{2}(\Omega) ,\varphi_{0}\in L^{2}(\Omega ), f(x)\in L^{2}(\Omega), h(x)\in L^{2}(\Omega)\), we say that \((u,\varphi)\) is a weak solution of (1)-(4) if
-
(i)
for all \(T>0\), \(u\in L^{\infty}(0,T;L^{2}(\Omega))\cap L^{2}(0,T;H^{\alpha}(\Omega) ),\varphi\in L^{\infty}(0,T;L^{2}(\Omega))\cap L^{2}(0,T; H^{\alpha}(\Omega) )\);
-
(ii)
for all \(\psi\in C^{\infty}(Q_{T})\), we have
$$\begin{aligned}& \begin{aligned}[b] &{-}id(u,\psi)+id(u_{0},\psi) \\ &\quad= \biggl(- \frac{dg^{2}+1}{U}+a\biggr) \int_{0}^{t}(u,\psi )\,dt+g\bigl[a+d(2\nu-2\mu)\bigr] \int_{0}^{t}(\varphi,\psi)\,dt\\ &\qquad{} -\frac{c}{4m} \int_{0}^{t}\bigl(\Lambda^{\alpha}u, \Lambda^{\alpha}\psi\bigr)\,dt -\frac{g}{4m}(c-d) \int_{0}^{t}\bigl(\Lambda^{\alpha}\varphi, \Lambda^{\alpha}\psi \bigr)\,dt\\ &\qquad{} -b \int_{0}^{t}\bigl(|u+g\varphi|^{2}(u+g \varphi),\psi\bigr)\,dt+id \int _{0}^{t}\bigl(f(x),\psi\bigr)\,dt, \end{aligned} \\& \begin{aligned}[b] i(\varphi,\psi)-i(\varphi_{0},\psi) ={}& {-}i \beta \int_{0}^{t}(\varphi,\psi )\,dt-\frac{g}{U} \int_{0}^{t}(u,\psi)\,dt+(2\nu-2\mu) \int_{0}^{t}(\varphi,\psi )\,dt \\ &{} +\frac{1}{4m} \int_{0}^{t}\bigl(\Lambda^{\alpha}\varphi, \Lambda^{\alpha }\psi\bigr)\,dt+i \int_{0}^{t}\bigl(h(x),\psi\bigr)\,dt, \end{aligned} \end{aligned}$$
with initial conditions
where \(Q_{T}=(0,T)\times\Omega\).
Next, we recall the following lemmas which will be used later.
Lemma 4
Let \(B_{0}, B, B_{1}\) be three Banach spaces such that
where the injections are continuous and \(B_{0}, B_{1}\) are reflexive, and the injection \(B_{0} \to B\) is compact. Denote
where T is finite and \(1< p_{0}, p_{1}<\infty\). Then W equipped with the norm
is a Banach space and the embedding \(W\hookrightarrow L^{p_{0}}(0,T;B)\) is compact.
Lemma 5
Assume that D is a bounded domain in \(\mathbb{R}_{x}^{n}\times\mathbb {R}_{t}\), functions \(f_{l},f \in L^{q}(D)\) (\(1< q<\infty\)) and
Then \(f_{l} \to f\) weakly in \(L^{q}(D)\).
Lemma 6
X is a Banach space, suppose that \(g\in L^{p}(0,T;X),\frac{\partial g}{\partial t}\in L^{p}(0,T;X)\) (\(1\le p\le\infty\)). Then \(g\in C([0,T];X)\) (after possibly being redefined on a set of measure zero).
Now, we state our main result as follows.
Theorem 1
Let \(\frac{1}{2}< \alpha\le1\), \(u_{0}\in L^{2}(\Omega), \varphi_{0}\in L^{2}(\Omega), f(x)\in L^{2}(\Omega), h(x)\in L^{2}(\Omega)\), \(U>0, c>0, m_{i}d_{r}+m_{r}d_{i}>0, b>0\), \(m_{i}>0\), \(d_{i}>0\), and \(\beta>0\), then there exists a unique weak solution \((u,\varphi)\) to (1)-(4) such that
and
where we say that \(u+g\varphi\in C([0,T];w-L^{2}(\Omega))\) and \(\varphi \in C([0,T];w-L^{2}(\Omega))\), if \((u+g\varphi,\phi_{1})\in C([0,T])\) and \((\varphi,\phi_{2})\in C([0,T])\) for any \(\phi=(\phi_{1},\phi_{2})\in L^{2}\times L^{2}\).
Proof
By Galerkin’s method, we look for approximate solutions \((u_{N}(x,t),\varphi_{N}(x,t))\) for equations (1)-(4) in the form
where \(\omega_{j}(x)=e^{ij\cdot x}\) and \(\alpha_{jN}(t),\beta_{jN}(t)\) satisfy the following system of ordinary differential equations:
The local existence theory for nonlinear ordinary differential equations ensures that the initial value problem (9)-(11) has at least one solution on \([0, t_{m}]\). By a priori estimates, we know that there exists a global solution for the initial value problem of the nonlinear ordinary differential system (9)-(11) on \([0,T]\).
Similar to the proof of Lemma 1 and Lemma 2, we have
Applying the compactness lemma (Lemma 4), there exist two subsequences, still denoted by \(\{u_{N}+g\varphi_{N}\},\{\varphi_{N}\}\), such that
Since \(u_{N}+g\varphi_{N}\in L^{4}(0,T;L^{4}(\Omega))\), we can deduce that
From Lemma 5, we have
For any \(\phi\in L^{2}\times L^{2}\), we see that \(\{(u+g\varphi,\phi_{1})\} _{N}\) and \(\{(\varphi,\phi_{2})\}_{N}\) are equicontinuous in \(C([0,T])\) from Lemma 3. On the other hand, it follows from Lemma 1 that \(\{(u+g\varphi ,\phi_{1})\}_{N}\) and \(\{(\varphi,\phi_{2})\}_{N}\) are uniformly bounded in \(C([0,T])\). Therefore, by Arzela-Ascoli theorem, we have \(\{(u+g\varphi ,\phi_{1})\}_{N}\) and \(\{(\varphi,\phi_{2})\}_{N}\) are compact in \(C([0,T])\).
Letting \(N\to\infty\), we immediately obtain
Therefore, we have for any \(\psi\in C^{\infty}(Q_{T})\)
Next we show the initial conditions hold.
By \(u_{N}+g\varphi_{N}\in L^{2}(0,T;H^{\alpha}(\Omega)),u_{Nt}+g\varphi_{Nt}\in L^{2}(0,T;H^{-\alpha}(\Omega))\), and Lemma 6, we have
Similarly, \(\varphi_{N}\in C([0,T],H^{-\alpha}(\Omega))\). Then
But from (12), we have
Therefore
Assume that \((u,\varphi)\) and \((u^{\star},\varphi^{\star})\) are the two different solutions of (1)-(4). For \(w=u-u^{\star},v=\varphi-\varphi^{\star}\), we have
Multiplying (12)-(13) by \(\overline{(w+gv)}\) and v̅, respectively, then integrating over Ω and taking the real part, we have
Combining the above two equations yields
where we have used the Gagliardo-Nirenberg inequality, the Young inequality, and the elementary inequality
Note that \(w(x,0)=0\) and \(v(x,0)=0\). By Gronwall’s inequality, we have
Thus we complete the proof. □
5 Global attractors
Finally, we consider the large time behavior of the solution to (1)-(4) in \(L^{2}\times L^{2}\). Before doing so, we give the following lemma which will be used later.
Lemma 7
[16]
Let \(g, h\) and y be three nonnegative locally integrable functions on \((t_{0},\infty)\) such that \(y'\) is locally integrable on \((t_{0},\infty)\), and which satisfy
and
where \(r,a_{1},a_{2},a_{3} \) are positive constants. Then
Theorem 2
Assume that \(\alpha\in(\frac{1}{2},1]\), \(U>0, c>0, m_{i}d_{r}+m_{r}d_{i}>0, b>0\), \(m_{i}>0\), \(d_{i}>0\), \(0< g<2\) and \(\beta>\frac{\frac{g}{|d|U}+\frac {g}{U}+1}{2-g}\), then the solution operator \(S(t):S(t)(u_{0}+g\varphi _{0},\varphi_{0})=(u+g\varphi,\varphi)\), for all \(t>0\), well defines a semigroup in the space \(L^{2}\times L^{2}\) and the following statements hold:
-
(i)
For any \(t>0\), \(S(t)\) is continuous in \(L^{2}\times L^{2}\).
-
(ii)
For any \((u_{0},\varphi_{0})\in L^{2}\times L^{2}\), S(t) is continuous from \([0,T]\) to \(L^{2}\times L^{2}\).
-
(iii)
For any \(t>0\), \(S(t)\) is compact in \(L^{2}\times L^{2}\).
-
(iv)
The semigroup \(\{S(t)\}_{t\ge0}\) possesses a global attractor \(\mathcal{A}\) in \(L^{2}\times L^{2}\).
Proof
It is easy to show that the solution \((u,\varphi)\) to (1)-(4) well defines a semigroup \(S(t)\) on \(L^{2}\times L^{2}\).
First, we consider the absorbing set in \(L^{2}\times L^{2}\). Similarly to the proof of Lemma 1, we have
Since
Thus we can choose \(\epsilon_{0}\) small enough such that
Define
we have
Applying Gronwall’s inequality yields
Therefore
where \(\rho_{0}^{2}=\frac{2bd_{i}}{|d|^{2}\epsilon_{0}^{2}}|\Omega |+\|f\|_{L^{2}}^{2}+(1+g^{2})\|h\|_{L^{2}}^{2}\).
It follows that the existence of an absorbing ball in \(L^{2}\times L^{2}\). Indeed, for any \(\rho>\rho_{0}\), denote by \(\mathcal{B}_{0}\) the ball \(B(0,\rho)\), for any bounded set \(\mathcal{B}\), there exists a positive \(t_{0}=\frac{1}{R_{1}}\log\frac{\rho^{2}}{\rho^{2}-\rho_{0}^{2}}\) such that \(S(t)\mathcal{B}\subset\mathcal{B}_{0}\) for all \(t\ge t_{0}\).
Next we consider the absorbing set in \(H^{\alpha}\times H^{\alpha}\).
Integrating (15) from t to \(t+1\), and using the definition of \(R_{1}\), we have
So for all \(t\ge t_{0}\), we see that
are uniformly bounded.
Multiplying (5) by \(\overline{\Lambda^{2\alpha}(u+g\varphi)}\), then integrating over Ω and taking the real part, we have
Since
we can choose δ small enough such that
Therefore
Multiplying (2) by \(\overline{\Lambda^{2\alpha}\varphi}\), integrating over Ω, and taking the imaginary part, we have
Combining (16) and (17) yields
where \(R_{2}=\max\{\frac{g}{|d|U}+\frac{2ad_{i}}{|d|^{2}}-\frac {2d_{i}}{|d|^{2}U}+\beta g+\frac{g^{2}}{\beta U^{2}},\frac{g}{|d|U}+\beta g-\beta\}\).
Since
applying the uniform Gronwall’s inequality (Lemma 7), we have
The above inequality can ensure us the existence of an absorbing ball in \(H^{\alpha}\times H^{\alpha}\). In fact, let \(\mathcal{B}\) be a bounded set in \(H^{\alpha}\times H^{\alpha}\). Obviously, it is also a bounded set in \(L^{2}\times L^{2}\) and \(S(t)\mathcal{B}\subset\mathcal{B}_{0}\) for \(t\ge t_{0}\). From (18), it follows that \(S(t)\mathcal{B}\subset\mathcal {B}_{1}\), where \(\mathcal{B}_{1}=B(0,\rho_{1})\) is a ball with radius \(\rho _{1}^{2}=\rho_{0}^{2}+(a_{1}+a_{2})e^{R_{2}}\) in \(H^{\alpha}\times H^{\alpha}\). Since the embedding \(H^{\alpha}\times H^{\alpha}\hookrightarrow L^{2}\times L^{2}\) is compact, we obtain
Thus the assertion of item (iii) has been proved.
Notice that items (i) and (ii) can be verified similar to the proof of Lemmas 1-4 and item (iv) is a direct corollary of items (i)-(iii). Thus the proof is completed. □
References
Machida, M, Koyama, T: Time-dependent Ginzburg-Landau theory for atomic Fermi gases near the BCS-BEC crossover. Phys. Rev. A 74, 033603 (2006)
Chen, S, Guo, B: On the Cauchy problem of the Ginzburg-Landau equations for atomic Fermi gases near the BCS-BEC crossover. J. Partial Differ. Equ. 22(3), 218-233 (2009)
Chen, S, Guo, B: Existence of the weak solution of coupled time-dependent Ginzburg-Landau equations. J. Math. Phys. 51(3), 033507 (2010)
Chen, S, Guo, B: Classical solutions of time-dependent Ginzburg-Landau theory for atomic Fermi gases near the BCS-BEC crossover. J. Differ. Equ. 251, 1415-1427 (2011)
Fang, S, Jin, L, Guo, B: Global attractor for the initial-boundary value problems for Ginzburg-Landau equations for atomic Fermi gases near the BCS-BEC crossover. Nonlinear Anal. 72, 4063-4070 (2010)
Guo, C, Chen, R, Fang, S, Wang, X: Finite difference method for Ginzburg-Landau equations for atomic Fermi gases near the BCS-BEC crossover. In: World Congress on Engineering and Technology, pp. 526-529 (2011)
Guan, J, Fang, S, Wang, X, Guo, C: Hopf bifurcations of traveling wave solutions for time-dependent Ginzburg-Landau equation for atomic Fermi gases near the BCS-BEC crossover. Commun. Nonlinear Sci. Numer. Simul. 18, 124-135 (2013)
Pu, X, Guo, B: Well-posedness and dynamics for the fractional Ginzburg-Landau equation. Appl. Anal. 1, 1-17 (2011)
Guo, B, Han, Y, Xin, J: Existence of the global smooth solution to the period boundary value problem of fractional nonlinear Schrödinger equation. Appl. Math. Comput. 204, 468-477 (2008)
Pu, X, Guo, B: The fractional Landau-Lifshitz-Gilbert equation and the heat flow of harmonic maps. Calc. Var. 42, 1-19 (2011)
Pu, X, Guo, B: Well-posedness for the fractional Landau-Lifshitz equation without Gilbert damping. Calc. Var. 46, 441-460 (2013)
Tarasov, V: Fractional Dynamics: Applications of Fractional Calculus to Dynamics of Particles, Fields and Media. Springer, Berlin (2010)
Kilbas, AA, Srivastava, HM, Trujillo, JJ: Theory and Applications of Fractional Differential Equations. Elsevier, Amsterdam (2006)
Podlubny, I: Fractional Differential Equations. Academic Press, San Diego (1999)
Miller, KS, Ross, B: An Introduction to the Fractional Calculus and Fractional Differential Equations. Wiley, New York (1993)
Temam, R: Infinite-Dimensional Dynamical Systems in Mechanics and Physics. Springer, New York (1998)
Acknowledgements
This work was supported by the National Natural Science Foundation of China (Grant No. 11271141; No. 11101160; No. 11426069; No. 61375006).
Author information
Authors and Affiliations
Corresponding author
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
The authors contributed equally to the writing of this paper. All authors read and approved the final manuscript.
Lang Li and Lingyu Jin contributed equally to this work.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
About this article
Cite this article
Li, L., Jin, L. & Fang, S. Large time behavior for the fractional Ginzburg-Landau equations near the BCS-BEC crossover regime of Fermi gases. Bound Value Probl 2017, 8 (2017). https://doi.org/10.1186/s13661-016-0738-9
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s13661-016-0738-9