[논문 리뷰] The probability density function for the Havriliak-Negami relaxation

K. Górska, A. Horzela외 3명|arXiv (Cornell University)|2016. 11. 19.

Nonlinear Dynamics and Pattern Formation인용 수 4

한 줄 요약

이 논문은 적분 변환을 사용하여 Havriliak-Negami 분극 모델과 관련된 시간 도메인 응답 함수 $ 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.

연구 동기 및 목표

• Havriliak-Negami 분극 모델과 관련된 시간 도메인 응답 함수와 공간 도메인 확률 밀도에 대한 정확하고 명시적인 표현을 유도하는 것.
• 기존의 Cole-Cole 모델($ \alpha < 1, \beta = 1 $)에 대한 결과를 $ \alpha = l/k < 1 $, $ \beta < k/l $ 인 더 일반적인 경우로 확장하는 것. 여기서 $ k,l $ 는 양의 정수이다.
• 모수 $ q \to 1^+ $ 로 갈 때 응답 함수의 점근적 행동을 조사하는 것. 특히 그것이 한쪽 띠 Lévy 안정 분포로 수렴하는 지를 분석하는 것.
• 확률 밀도 $ g_{\alpha, \beta}(u) $ 의 자기유사성 성질을 활용하여 적분 형태의 진화 방정식을 만족하는 이변량 밀도를 구성하는 것.

제안 방법

• 특히 Mellin 변환을 활용한 적분 변환 기법을 사용하여 시간 도메인 응답 함수 $ f_{\alpha, \beta}(t/\tau_0) $ 와 공간 도메인 확률 밀도 $ g_{\alpha, \beta}(u) $ 의 정확한 표현을 도출한다.
• 유도된 함수의 수렴성과 해석성을 보장하기 위해 $ \alpha = l/k < 1 $ 와 $ \beta < k/l $ 조건을 적용하며, $ k,l \in \mathbb{Z}^+ $ 이다.
• 응답 함수 $ f_{\alpha, \beta}(t/\tau_0) $ 의 점근적 행동을 설정하기 위해 $ \beta = (2-q)/(q-1) $ 과 $ \tau_0 = (q-1)^{1/\alpha} $ 를 설정하고, $ q \to 1^+ $ 로 갈 때의 극한을 분석한다.
• 확률 밀도 $ g_{\alpha, \beta}(u) $ 의 자기유사성 성질을 식별하여, 적분 형태의 진화 방정식을 만족하는 이변량 확률 밀도 함수를 구성할 수 있도록 한다.
• 자기유사성 구조를 활용하여 시간과 공간 변수 사이의 스케일링 대칭성을 통해 연결된 적분 형태의 진화 방정식을 유도한다.
• 특수한 경우로 $ \beta = 1 $ 이고 $ 0 < \alpha < 1 $ 일 때, 기존의 Cole-Cole 모델을 회복함으로써 일반 프레임워크의 타당성을 검증한다.

실험 결과

연구 질문

• RQ1 제약 조건 $ \alpha = l/k < 1 $, $ \beta < k/l $ 하에서 Havriliak-Negami 분극 모델의 시간 도메인 응답 함수와 공간 도메인 확률 밀도에 대해 정확하고 명시적인 표현을 유도할 수 있는가?
• 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} $ 이면, $ q \to 1^+ $ 로 갈수록 응답 함수 $ f_{\alpha, \beta}(t/\tau_0) $ 는 한쪽 띠 Lévy 안정 분포로 수렴하며, 이는 Havriliak-Negami 모델과 Lévy 통계 사이의 직접적인 연결 고리를 확립한다.
• 확률 밀도 $ g_{\alpha, \beta}(u) $ 의 자기유사성 성질을 활용하여 적분 형태의 진화 방정식을 만족하는 이변량 확률 밀도 함수를 구성한다.
• 유도된 프레임워크는 $ \alpha < 1 $ 이고 $ \beta = 1 $ 일 때 기존의 Cole-Cole 분극 모델을 특수한 경우로 포함하며, 기존 결과와의 일致성을 확인한다.
• 적분 변환 기법을 통해 시간 도메인 및 공간 도메인 함수에 대한 닫힌 형태의 표현을 성공적으로 도출하였으며, 이는 분극 과정의 심층적 분석을 가능하게 한다.
• 자기유사성 구조에서 유도된 진화 방정식은 분극 과정의 역학을 위한 새로운 적분 형태의 공식화를 제공하며, 더 넓은 스케일링 해의 클래스를 시사한다.