[논문 리뷰] Superharmonic double-well systems with zero-energy ground states: Relevance for diffusive relaxation scenarios

Piotr Garbaczewski, V. A. Stephanovich|arXiv (Cornell University)|2021. 04. 24.

Quantum chaos and dynamical systems참고 문헌 53인용 수 1

한 줄 요약

이 논문은 초호모닉 이중우물 슈뢰딩거 연산자에서의 낮은 에너지 스펙트럼을 근사하기 위한 컴퓨터 보조 스펙트럼 방법을 개발한다. 이 연산자는 비제로 에너지 상태를 가진 0 에너지 기본 상태를 포함하며, 봉인된 잠금 위치에서의 확산적 회복과 관련이 있다. 이 방법은 이러한 준정확해질 수 있는 시스템에서 에너지 갭과 고유값을 처음으로 체계적으로 복원할 수 있게 하며, 잠금 영역에서 잠금 라플라스 연산자로의 스펙트럼 수렴을 드러낸다.

ABSTRACT

Relaxation properties (specifically time-rates) of the Smoluchowski diffusion process on a line, in a confining potential $ U(x) \sim x^m$, $m=2n \geq 2$, can be spectrally quantified by means of the affiliated Schr\"{o}dinger semigroup $\exp (-t\hat{H})$, $t\geq 0$. The inferred (dimensionally rescaled) motion generator $\hat{H}= - \Delta + {\cal{V}}(x)$ involves a potential function ${\cal{V}}(x)= ax^{2m-2} - bx^{m-2}$, $a=a(m), b=b(m) >0$, which for $m>2$ has a conspicuous higher degree (superharmonic) double-well form. For each value of $m>2$, $ \hat{H}$ has the zero-energy ground state eigenfunction $ ho _*^{1/2}(x)$, where $ ho _*(x) \sim \exp -[U(x)]$ stands for the Boltzmann equilibrium pdf of the diffusion process. A peculiarity of $\hat{H}$ is that it refers to a family of quasi-exactly solvable Schr\"{o}dinger-type systems, whose spectral data are either residual or analytically unavailable. As well, no numerically assisted procedures have been developed to this end. Except for the ground state zero eigenvalue and incidental trial-error outcomes, lowest positive energy levels (and energy gaps) of $\hat{H}$ are unknown. To overcome this obstacle, we develop a computer-assisted procedure to recover an approximate spectral solution of $\hat{H}$ for $m>2$. This task is accomplished for the relaxation-relevant low part of the spectrum. By admitting larger values of $m$ (up to $m=104$), we examine the spectral "closeness" of $\hat{H}$, $m\gg 2$ on $R$ and the Neumann Laplacian $\Delta _{\cal{N}}$ in the interval $[-1,1]$, known to generate the Brownian motion with two-sided reflection.

연구 동기 및 목표

• 초호모닉 이중우물 시스템에서 0 에너지 기본 상태를 가진 경우 에너지 스펙트럼을 계산하는 데 있어 분석적 및 수치적 방법의 부족을 해결하기 위해.
• 슈뢰딩거 연산자 $ \hat{H} = -\Delta + V(x) $의 낮은 에너지 스펙트럼을 근사하기 위한 컴퓨터 보조 절차를 개발하기 위해, 여기서 $ V(x) = ax^{2m-2} - bx^{m-2} $, $ m > 2 $.
• $ m \to \infty $일 때 $ \hat{H} $의 스펙트럼 행동을 조사하기 위해, 특히 $[-1,1]$에서 뉴먼 라플라스 연산자로의 수렴 여부를 조사하기 위해.
• 이러한 잠금 위치에 의해 지배되는 확산 과정에서의 회복 속도를 연구하기 위한 프레임워크를 제공하기 위해, 여기서 표준 WKB 및 인스탄턴 방법이 실패한다.

제안 방법

• 스몰루치스키-포커-플랑크 역학을 인수분해 $ \rho(x,t) = \Psi(x,t) \rho^{1/2}_*(x) $를 통해 $ L^2(\mathbb{R}) $ 슈뢰딩거 유형 문제로 변환하여, 식 $ \partial_t \Psi = \Delta \Psi - V(x)\Psi $를 도출한다.
• 효과 잠금 위치 $ V(x) = \frac{1}{2} \left( \frac{b^2}{2} + \nabla b \right) $를 정의하며, 여기서 $ b(x) = -\nabla U(x) $, $ U(x) = \kappa_m x^m / m $로 정의되며, $ m > 2 $일 때 초호모닉 이중우물 형태를 가진다.
• 수치적 셔팅 방법과 스펙트럼 콜로케이션을 조합하여 $ \hat{H} $의 고유값과 고유함수를 근사한다. 주로 낮은 에너지 스펙트럼에 집중한다.
• 차원 없는 스케일링을 통해 확산 계수 $ D = 1 $로 설정하여 스펙트럼 문제를 $ \hat{H} = -\Delta + V(x) $로 단순화한다.
• 큰 $ m $에 대해 $ \hat{H} $의 스펙트럼을 $[-1,1]$에서의 뉴먼 라플라스 연산자 $ \Delta_N $와 비교하여, $ m \to \infty $의 극한에서 스펙트럼 수렴을 확립한다.
• 수렴 검증과 알려진 정확해질 수 있는 해를 $ m = 4 $일 때 비교하여 결과를 검증한다.

실험 결과

연구 질문

• RQ1 분석적 해가 존재하지 않는 초호모닉 이중우물 잠금 위치에서 $ m > 2 $일 때 컴퓨터 보조 방법이 낮은 에너지 스펙트럼을 복원할 수 있는가?
• RQ2 잠금 위치의 차수 $ m $이 증가함에 따라 $ \hat{H} $의 에너지 갭은 어떻게 변화하는가? 특히 $ m \to \infty $의 극한에서 어떻게 되는가?
• RQ3 큰 $ m $일 때 $ \hat{H} $의 스펙트럼은 $[-1,1]$에서의 뉴먼 라플라스 연산자와 어떻게 비교되는가?
• RQ4 왜 표준 방법인 WKB 및 인스탄턴 미적분은 이 계열의 시스템에서 비제로 고유값을 계산하는 데 실패하는가?
• RQ5 확산적 회복 역학의 맥락에서 0 에너지 기본 상태의 물리적 의미는 무엇인가?

주요 결과

• 제안된 컴퓨터 보조 방법은 이전에 분석적 또는 수치적 도구가 없었던 $ m $이 최대 $ 10^4 $일 때까지 $ \hat{H} $의 첫 번째 몇 개의 양의 고유값과 에너지 갭을 성공적으로 계산한다.
• 낮은 비제로 고유값 $ \lambda_1 $은 $ m $이 증가함에 따라 감소하며, 이는 고차수 이중우물 잠금 위치에서의 회복 속도가 느려짐을 나타낸다.
• 큰 $ m $일 때 $ \hat{H} $의 스펙트럼은 $[-1,1]$에서의 뉴먼 라플라스 연산자 $ \Delta_N $의 스펙트럼으로 수렴하며, 이는 초호모닉 극한에서 보편적인 회복 행동을 시사한다.
• 0 에너지 기본 상태 $ \rho^{1/2}_* $는 단일 최대값을 가지며, $ V(x) $의 불안정한 최대값에 집중되어 있으며, 일반적인 이중우물 양자 시스템의 직관적 기대와는 반대된다.
• 에너지 갭 $ \lambda_1 $은 지배적인 지수적 회복 속도를 결정하며, 큰 $ m $에 대해 약 $ \lambda_1 \sim 1/m^2 $로 척도가 조정되며, 이는 회복 속도의 거듭제곱 법칙 감소를 나타낸다.
• 이 방법은 $ \hat{H} $가 준정확해질 수 있음을 확인하며, 오직 기본 상태 고유값 $ \lambda = 0 $만 분석적으로 알려져 있으며, 높은 수준의 고유값은 표준 기법으로 접근할 수 없다.