[論文レビュー] Relaxation Equations: Fractional Models

Ester C. F. A. Rosa, E. Capelas de Oliveira|arXiv (Cornell University)|Oct 6, 2015

Fractional Differential Equations Solutions参考文献 23被引用数 4

ひとこと要約

本稿では、Riemann-Liouville型の分数階微分（$\gamma$、$0 < \gamma \leq 1$）を用いた一般化された分数階運動方程式を導入し、Laplace変換を用いてMittag-Leffler関数で表される解を得た。主な貢献は、$\gamma = 1$ のとき、古典的Debye、Cole-Cole、Cole-Davidson、Havriliak-Negamiの緩和関数が特殊ケースとして得られることであり、異常な緩和挙動をより柔軟に記述できる分数階フレームワークによって、これらのモデルが統一されることである。

ABSTRACT

The relaxation functions introduced empirically by Debye, Cole-Cole, Cole-Davidson and Havriliak-Negami are, each of them, solutions to their respective kinetic equations. In this work, we propose a generalization of such equations by introducing a fractional differential operator written in terms of the Riemann-Liouville fractional derivative of order $γ$, $0 &lt; γ\leq 1$. In order to solve the generalized equations, the Laplace transform methodology is introduced and the corresponding solutions are then presented, in terms of Mittag-Leffler functions. In the case in which the derivative's order is $γ=1$, the traditional relaxation functions are recovered. Finally, we presente some 2D graphs of these function.

研究の動機と目的

• 運動方程式に $\gamma$（$0 < \gamma \leq 1$）の分数階微分を導入することで、古典的誘電緩和モデルを一般化すること。
• Laplace変換法を用いて、これらの分数階運動方程式の解析的解を導出すること。
• 解をMittag-Leffler関数の形で表現し、古典的緩和関数を一般化すること。
• $\gamma = 1$ のとき、古典的緩和モデル（Debye、Cole-Cole、Cole-Davidson、Havriliak-Negami）が特殊ケースとして回復されることを示すこと。
• $\gamma$、$\alpha$、$\beta$、$\sigma$ の異なるパラメータ値に対する分数階緩和関数のグラフィカルな可視化を提供すること。

提案手法

• 誘電緩和過程を記述するため、Riemann-Liouville型の分数階微分（$\gamma$）を用いた分数階運動方程式を定式化する。
• Debye、Cole-Cole、Cole-Davidson、Havriliak-Negamiの各古典的緩和モデルに対応する積分核に基づく記憶関数を構築する。
• 分数階運動方程式にLaplace変換を適用し、変換された緩和関数の代数的表現を導出する。
• 特に $\mathscr{L}[t^{\gamma-1}E_{\alpha,\gamma}^{\beta}(at^{\alpha})](s) = \frac{s^{\alpha\beta - \gamma}}{(s^{\alpha} - a)^{\beta}}$ というMittag-Leffler関数の既知のLaplace変換恒等式を用い、解を逆変換する。
• $\gamma = 1$ と特定のパラメータ値（例：$\beta = 1$、$\alpha = 1$）を設定することで、古典的モデルを回復する。
• 2次元プロットを用いて、$\gamma$、$\alpha$、$\beta$、$\sigma$ に依存する分数階緩和関数の依存関係を可視化する。

実験結果

リサーチクエスチョン

• RQ1 Riemann-Liouville微分（$\gamma \in (0,1]$）を用いた分数階微分方程式は、古典的誘電緩和モデルを一般化できるか？
• RQ2 これらの分数階方程式の解は、$\gamma = 1$ のとき、既知の古典的緩和関数（Debye、Cole-Cole、Cole-Davidson、Havriliak-Negami）に還元されるか？
• RQ3 $\gamma < 1$ のとき、Mittag-Leffler関数に基づく解はどのように振る舞い、その単調性および物理的解釈可能性はいかがなものか？
• RQ4 パラメータ $\alpha$、$\beta$、$\gamma$、$\sigma$ は、分数階フレームワークにおける緩和ダイナミクスにどのように寄与しているか？
• RQ5 分数階緩和関数から、$\tilde{\varepsilon}(s) = 1 - s^{\gamma}\tilde{\varphi}(s)$ というLaplace変換関係を用いて、複素誘電率を一貫して導出できるか？

主な発見

• Riemann-Liouville微分（$\gamma$）を用いた分数階運動方程式は、古典的緩和モデルを一般化でき、解はMittag-Leffler関数で表される。
• $\gamma = 1$ のとき、分数階Debyeモデルは古典的Debye緩和関数 $\varphi(t) = e^{-t/\sigma}$ に一致する。
• 分数階Cole-Coleモデルは、$\alpha = 0.5$ のとき、$\varphi(t)_{CCF} = \frac{t^{\gamma-1}}{\Gamma(\gamma)} - \sigma^{-\alpha}t^{\alpha + \gamma - 1}E_{\alpha,\alpha + \gamma}^{\alpha}\left(-\frac{t}{\sigma}\right)$ を与え、伸びた指数的減衰を示す。
• 分数階Cole-Davidsonモデルは $\varphi(t)_{CDF} = \frac{t^{\gamma-1}}{\Gamma(\gamma)} - \sigma^{-\beta}t^{\beta + \gamma - 1}E_{1,\beta + \gamma}^{\beta}\left(-\frac{t}{\sigma}\right)$ で与えられ、$\gamma = 1$ のとき古典的C-D関数に回復する。
• 分数階Havriliak-Negamiモデルは $\varphi(t)_{HNF} = \frac{t^{\gamma-1}}{\Gamma(\gamma)} - \sigma^{-\alpha\beta}t^{\alpha\beta + \gamma - 1}E_{\alpha,\alpha\beta + \gamma}^{\beta}\left[-\left(\frac{t}{\sigma}\right)^\alpha\right]$ で表され、$\gamma = 1$ のとき古典的H-Nモデルに簡略化される。
• グラフィカル分析により、$t > 0$ において分数階緩和関数が完全単調性を保ち、$\gamma$、$\alpha$、$\beta$ に依存する非指数的・累乗則的減衰特性を示すことが確認された。