[论文解读] The probability density function for the Havriliak-Negami relaxation

K. Górska, A. Horzela+3 more|arXiv (Cornell University)|Nov 19, 2016

Nonlinear Dynamics and Pattern Formation被引用 4

一句话总结

本文利用积分变换推导出与哈夫里利阿科-内加米弛豫模型相关的时域响应函数 $ f_{\alpha, \beta}(t/\tau_0) $ 和空域概率密度 $ g_{\alpha, \beta}(u) $ 的精确、显式表达式。关键贡献在于证明：当 $ \beta = (2-q)/(q-1) $ 且 $ \tau_0 = (q-1)^{1/\alpha} $ 时，响应函数在 $ q \to 1^+ $ 时收敛于单边 Lévy 稳定分布，并通过自相似性引入了一种满足积分演化方程的双变量密度函数。

ABSTRACT

We study the functions related to the Havriliak-Negami frequency relaxation $\sim [1+(i\omega au_{0})^{\alpha}]^{-\beta}$ with $ au_{0}$ characteristic time, measured in many experiments. We furnish exact and explicit expression for the response function $f_{\alpha, \beta}(t/ au_{0})$ in time domain and a probability density $g_{\alpha, \beta}(u)$ in space domain for $\alpha = l/k < 1$ and $\beta < k/l$, with $k$ and $l$ positive integers. For $0 < \alpha < 1$ and $\beta=1$ we reproduce the functions related to the Cole-Cole relaxation. We use the method of integral transforms. We show that $f_{\alpha, \beta}(t/ au_{0})$ with $\beta = (2-q)/(q-1)$ and $ au_{0} = (q-1)^{1/\alpha}$, $1 < q < 2$, goes over to the one-sided Levy stable distribution when $q$ tends to one. Moreover, applying the self-similar property of $g_{\alpha, \beta}(u)$ we introduce two-variable density which satisfies the integral form of evolution equation.

研究动机与目标

• 推导与哈夫里利阿科-内加米弛豫模型相关的时域响应函数和空域概率密度的精确、显式表达式。
• 将已知的 Cole-Cole 模型结果（$ \alpha < 1, \beta = 1 $）推广至更一般的情形，其中 $ \alpha = l/k < 1 $，$ \beta < k/l $，且 $ k,l $ 为正整数。
• 研究当参数 $ q \to 1^+ $ 时响应函数的极限行为，特别是其收敛于单边 Lévy 稳定分布的性质。
• 利用概率密度 $ g_{\alpha, \beta}(u) $ 的自相似性，构造满足积分形式演化方程的双变量密度函数。

提出的方法

• 利用积分变换，特别是梅林变换，推导时域响应函数 $ f_{\alpha, \beta}(t/\tau_0) $ 和空域概率密度 $ g_{\alpha, \beta}(u) $ 的精确表达式。
• 应用条件 $ \alpha = l/k < 1 $ 和 $ \beta < k/l $，其中 $ k,l \in \mathbb{Z}^+ $，以确保所导函数的收敛性与解析性。
• 通过设定 $ \beta = (2-q)/(q-1) $ 和 $ \tau_0 = (q-1)^{1/\alpha} $，建立 $ f_{\alpha, \beta}(t/\tau_0) $ 的极限行为，并分析 $ q \to 1^+ $ 时的极限。
• 识别出 $ g_{\alpha, \beta}(u) $ 的自相似性，从而实现满足积分形式演化方程的双变量概率密度函数的构造。
• 利用自相似结构推导出以时间与空间变量通过标度对称性关联的积分形式演化方程。
• 通过当 $ \beta = 1 $ 且 $ 0 < \alpha < 1 $ 时恢复 Cole-Cole 模型作为特例，验证了通用框架的有效性。

实验结果

研究问题

• RQ1 在约束 $ \alpha = l/k < 1 $，$ \beta < k/l $ 下，能否为哈夫里利阿科-内加米弛豫模型的时域响应函数与空域概率密度推导出精确、显式的表达式？
• RQ2 当 $ \beta = (2-q)/(q-1) $ 且 $ \tau_0 = (q-1)^{1/\alpha} $ 时，响应函数在 $ q \to 1^+ $ 时的行为如何？
• RQ3 概率密度 $ g_{\alpha, \beta}(u) $ 的自相似性质是否可用于构造满足积分形式演化方程的双变量密度函数？
• RQ4 当 $ \beta = 1 $ 时，所推导的响应函数是否退化为已知的 Cole-Cole 弛豫函数？
• RQ5 当 $ q \to 1^+ $ 时，响应函数的极限形式是否等价于一个单边 Lévy 稳定分布？

主要发现

• 本文在 $ \alpha = l/k < 1 $，$ \beta < k/l $，且 $ k,l $ 为正整数的条件下，推导出时域响应函数 $ f_{\alpha, \beta}(t/\tau_0) $ 和空域概率密度 $ g_{\alpha, \beta}(u) $ 的精确、显式表达式。
• 当 $ \beta = (2-q)/(q-1) $ 且 $ \tau_0 = (q-1)^{1/\alpha} $ 时，响应函数 $ f_{\alpha, \beta}(t/\tau_0) $ 在 $ q \to 1^+ $ 时收敛于单边 Lévy 稳定分布，建立了哈夫里利阿科-内加米模型与 Lévy 统计之间的直接联系。
• 利用 $ g_{\alpha, \beta}(u) $ 的自相似性质，构造出满足积分形式演化方程的双变量概率密度函数。
• 所推导的框架在 $ \alpha < 1 $ 且 $ \beta = 1 $ 时可还原为已知的 Cole-Cole 弛豫模型，验证了与已有结果的一致性。
• 积分变换方法成功获得了时域与空域函数的闭式表达式，为弛豫过程的深入分析研究提供了可能。
• 由自相似结构导出的演化方程提供了一类新的积分形式，用于描述弛豫过程的动力学，暗示存在更广泛的标度解类。