[논문 리뷰] Mean exit time and escape probability for dynamical systems driven by Levy noise

Ting Gao, Jinqiao Duan외 2명|arXiv (Cornell University)|2012. 01. 29.

Ecosystem dynamics and resilience참고 문헌 13인용 수 4

한 줄 요약

이 논문은 $\alpha$-stable Lévy 소음에 의해 구동되는 확률적 동역계 시스템의 평균 도착 시간과 탈출 확률을 계산하기 위한 정확한 수치적 방법을 개발한다. 이는 국소적이지 않은 적분미분방정식으로 특징지어진다. 주요 발견은 평균 도착 시간이 도메인 크기와 $\alpha$에 크게 의존하며, 더 작은 $\alpha$는 큰 점프를 유도하고 큰 도메인에서의 탈출을 가속화하는 반면, 탈출 확률은 비선형성에 따라 크게 달라지며 $\alpha \in (0,1)$일 경우 경계에서 불연속성이 나타난다는 것이다.

ABSTRACT

The mean first exit time and escape probability are utilized to quantify dynamical behaviors of stochastic differential equations with non-Gaussian alpha-stable type Levy motions. Both deterministic quantities are characterized by differential-integral equations(i.e.,differential equations with non local terms) but with different exterior conditions. The non-Gaussianity of noises manifests as nonlocality at the level of mean exit time and escape probability. An objective of this paper is to make mean exit time and escape probability as efficient computational tools, to the applied probability community, for quantifying stochastic dynamics. An accurate numerical scheme is developed and validated for computing the mean exit time and escape probability. Asymptotic solution for the mean exit time is given when the pure jump measure in the Levy motion is small. From both the analytical and numerical results, it is observed that the mean exit time depends strongly on the domain size and the value of alpha in the alpha-stable Levy jump measure. The mean exit time can measure which of the two competing factors in alpha-stable Levy motion, i.e. the jump frequency or the jump size, is dominant in helping a process exit a bounded domain. The escape probability is shown to vary with the underlying vector field(i.e.,drift). The mean exit time and escape probability could become discontinuous at the boundary of the domain, when the process is subject to certain deterministic potential and the value of alpha is in (0,1).

연구 동기 및 목표

• 비정규 확률적 동역계 시스템에서 평균 도착 시간과 탈출 확률을 위한 효율적인 수치 계산 프레임워크를 개발하기 위해.
• 특히 점프 빈도와 크기의 영향을 분석하기 위해 비정규 Lévy 소음이 유한 도메인에서의 탈출 행동에 미치는 영향을 연구하기 위해.
• 특히 작은 점프 조건에서의 점근적 및 해석적 해와의 비교를 통해 수치적 방법의 타당성을 검증하기 위해.
• 벡터장(비선형성)과 $\alpha$-stable Lévy 파ameter가 탈출 가능성과 도착 시간에 미치는 영향을 조사하기 위해.

제안 방법

• 비선형성 $f(x)$와 점프 측도 $\nu$를 갖는 $\alpha$-stable Lévy 운동에 의해 구동되는 확률적 미분 방정식을 사용하여 시스템을 모델링하기 위해.
• Lévy 생성자에서 유래하는 비국소적 항을 포함한 미분-적분 방정식으로 평균 도착 시간 $u(x)$와 탈출 확률 $P_E(x)$를 수식화하기 위해.
• 원점 근처에서 Lévy 측도의 특이성을 처리하기 위해 수치적 방법에 보정 항을 통합하기 위해.
• 비국소적 적분 항에 대해 적응형 사다리꼴 적분법을 사용한 유한 차분 방법을 적용하였으며, 점근적 및 해석적 해와의 비교를 통해 검증하였다.
• 이중 웰 포텐셜 $f(x) = x - x^3$에 대한 수치 시뮬레이션을 수행하여 안정 상태 간의 전이 확률을 연구하기 위해.
• 도메인 크기 $D$, $\alpha$, 소음 강도 $\varepsilon$의 변화가 도메인 탈출 동역학에 미치는 영향을 분석하기 위해.

실험 결과

연구 질문

• RQ1 α-안정 Lévy 과정에서 $\alpha$의 값이 유한 도메인에서의 평균 도착 시간에 어떤 영향을 미치는가?
• RQ2 α-안정 Lévy 소음에서 탈출 동역학을 결정하는 데 있어 점프 빈도와 점프 크기 중 어느 것이 더 중요한가?
• RQ3 기본적인 비선형성 벡터장 $f(x)$는 도메인 보완부분의 목표 부분집합 $E$로의 탈출 확률에 어떤 영향을 미치는가?
• RQ4 평균 도착 시간 또는 탈출 확률이 도메인 경계에서 불연속성이 되는 조건은 무엇인가?
• RQ5 Gaussian 소음($d > 0$)을 추가할 경우 순수 $\alpha$-안정 Lévy 운동과 비교해 평균 도착 시간과 탈출 확률의 행동은 어떻게 변화하는가?

주요 결과

• α ∈ (0,1)일 경우, 시스템이 결정론적 포텐셜에 의해 영향을 받을 때 평균 도착 시간은 도메인 경계에서 불연속성이 나타난다.
• 더 작은 $\alpha$ 값(예: $\alpha = 0.5$)은 경계에서 멀어질수록 평균 도착 시간 프로파일이 더 평탄하게 만들며, 이는 도메인 크기보다 큰 점프가 지배적임을 시사한다.
• α가 2에 가까울 경우, 높은 빈도의 작은 점프는 작은 도메인에서의 탈출 가능성을 높이며, α가 0에 가까워질수록 큰 도메인에서 큰 점프가 유리하다.
• 탈출 확률 $P_E(x)$는 Gaussian 소음($\alpha = 2$)일 때 가장 크며, $\alpha$가 감소함에 따라 감소한다. 특히 도메인 중심 근처의 시작점일 경우 두드러진 감소가 관찰된다.
• α < 1일 경우, 이중 웰 포텐셜 설정에서 도메인의 왼쪽 경계에서 탈출 확률이 불연속성을 보인다.
• 작은 양의 Gaussian 소음($d = 0.1$)을 추가하면 $\alpha$ 값에 따른 탈출 확률의 차이가 감소하여 전이 행동이 매끄럽게 개선된다.