Brownian motion in terms of Lifson and Jackson (LJ) formula has been widely explored in periodic systems and it has been believed for a long time that the LJ formula only applies to periodic potentials. Recently we show that for the following Brownian motion $γ\dot{x} = -U'(x) + ξ$, where $U(x)$ is the quasi-periodic potential, the effective diffusion constant can still be described by the LJ formula $D^* = D/(\langle \exp(βU)\rangle \langle \exp(-βU)\rangle)$, where the average is redefined as $\langle \exp(βU)\rangle = \lim_{L\rightarrow \infty} L^{-1} \int_0^L \exp(βU(x))dx$. In this manuscript we prove this result exactly using the mean first passage time $Ï(x)$, with boundary conditions $Ï(\pm L) = 0$, and show that the effective diffusion constant can be determined using $D^* =\lim_{L \rightarrow \infty} L^2/(2Ï(0))$, where $\pm L$ is the two positions of the absorbing boundary. We exactly solve the equation of motion of $Ï(x)$ and obtain the above result with the aid of Jacobi-Anger expansion method. Our result can be generalized to the other potentials and even higher dimensions, which can greatly broaden our understanding of Brownian motion in more general circumstances. The requirement for a well-defined effective diffusion constant $D^*$ in more general potentials is also discussed.
論文ID : 2510.10435タイトル : Proof of the exact diffusion constant via first passage time in quasi-periodic potentials著者 : Ming Gong(中国科学技術大学量子情報重点実験室)分類 : cond-mat.stat-mech(統計力学)発表日 : 2025年10月14日論文リンク : https://arxiv.org/abs/2510.10435 本論文は、Lifson-Jackson (LJ)公式が準周期ポテンシャルにおいて厳密に適用可能であることを証明している。従来、LJ公式は周期ポテンシャルにのみ適用可能と考えられていたが、著者は平均第一通過時間(MFPT)の厳密な数学的証明を通じて、準周期ポテンシャルU ( x ) U(x) U ( x ) において、有効拡散定数が再定義されたLJ公式で記述できることを示した:D ∗ = D / ( ⟨ exp ( β U ) ⟩ ⟨ exp ( − β U ) ⟩ ) D^* = D/(\langle \exp(\beta U)\rangle \langle \exp(-\beta U)\rangle) D ∗ = D / (⟨ exp ( β U )⟩ ⟨ exp ( − β U )⟩) 。ここで平均値は⟨ exp ( β U ) ⟩ = lim L → ∞ L − 1 ∫ 0 L exp ( β U ( x ) ) d x \langle \exp(\beta U)\rangle = \lim_{L\rightarrow \infty} L^{-1} \int_0^L \exp(\beta U(x))dx ⟨ exp ( β U )⟩ = lim L → ∞ L − 1 ∫ 0 L exp ( β U ( x )) d x と定義される。Jacobi-Anger展開法を用いて、著者はMFPT方程式を厳密に解き、D ∗ = lim L → ∞ L 2 / ( 2 τ ( 0 ) ) D^* =\lim_{L \rightarrow \infty} L^2/(2\tau(0)) D ∗ = lim L → ∞ L 2 / ( 2 τ ( 0 )) を証明した。
古典理論の限界 :1962年に提唱されたLifson-Jackson公式は、厳密に周期的なポテンシャル場におけるブラウン運動にのみ適用可能と考えられてきた。この制限は、より一般的な物理系への公式の応用を大きく制約している。実際的応用の必要性 :実際の物理系では、完全に周期的なポテンシャル場は実現困難であり、準周期的または非周期的ポテンシャル場がより一般的である。したがって、LJ公式の適用範囲を拡張することは、理論的および実践的に重要な意義を持つ。理論的争点 :Zwanzigは、系が2つの異なる周期a a a とb b b を有する場合、LJ公式はa ≫ b a \gg b a ≫ b の条件下でのみ近似的に成立すると主張し、その適用範囲をさらに制限している。著者は先行研究(参考文献13 )において、数値的方法によってLJ公式の準周期ポテンシャルにおける有効性を示したが、査読者から厳密な数学的証明の欠如を指摘された。本論文は、この指摘に応えるため、平均第一通過時間に基づいた厳密な数学的証明を提供するものである。
理論的突破 :LJ公式が準周期ポテンシャルにおいて厳密に適用可能であることを初めて証明し、従来の周期性に関する理論的制限を打破した。数学的方法の革新 :Jacobi-Anger展開に基づいた解析的方法を開発し、準周期ポテンシャルにおけるMFPT方程式を厳密に解くことが可能になった。公式の一般化 :古典的なLJ公式をより一般的なポテンシャル場に拡張し、平均値の計算方法を再定義するだけで対応可能にした。普遍的結果 :本方法が多重準周期ポテンシャルおよび高次元系に推広可能であることを証明し、ブラウン運動理論の適用範囲を大幅に拡大した。準周期ポテンシャルU ( x ) U(x) U ( x ) における過減衰ブラウン運動を研究する:
γ x ˙ = − U ′ ( x ) + ξ \gamma\dot{x} = -U'(x) + \xi γ x ˙ = − U ′ ( x ) + ξ
ここでγ \gamma γ は粘性係数、ξ \xi ξ は⟨ ξ ( t ) ⟩ = 0 \langle\xi(t)\rangle = 0 ⟨ ξ ( t )⟩ = 0 および⟨ ξ ( t ) ξ ( t ′ ) ⟩ = 2 γ k B T δ ( t − t ′ ) \langle\xi(t)\xi(t')\rangle = 2\gamma k_BT\delta(t-t') ⟨ ξ ( t ) ξ ( t ′ )⟩ = 2 γ k B T δ ( t − t ′ ) を満たすランダム力である。
目標は、有効拡散定数が以下のように表現できることを証明することである:
D ∗ = D ⟨ exp ( Φ ( x ) ) ⟩ ⟨ exp ( − Φ ( x ) ) ⟩ D^* = \frac{D}{\langle\exp(\Phi(x))\rangle\langle\exp(-\Phi(x))\rangle} D ∗ = ⟨ e x p ( Φ ( x ))⟩ ⟨ e x p ( − Φ ( x ))⟩ D
ここでΦ ( x ) = β U ( x ) \Phi(x) = \beta U(x) Φ ( x ) = β U ( x ) 、β = 1 / k B T \beta = 1/k_BT β = 1/ k B T である。
境界条件τ ( ± L ) = 0 \tau(\pm L) = 0 τ ( ± L ) = 0 を有するMFPT問題に対する支配方程式は:
D e − Φ ( x ) ∂ ∂ x e Φ ( x ) ∂ ∂ x τ ( x ) = − 1 De^{-\Phi(x)} \frac{\partial}{\partial x}e^{\Phi(x)} \frac{\partial}{\partial x}\tau(x) = -1 D e − Φ ( x ) ∂ x ∂ e Φ ( x ) ∂ x ∂ τ ( x ) = − 1
周期ポテンシャルΦ ( x ) = U cos ( x ) \Phi(x) = U\cos(x) Φ ( x ) = U cos ( x ) に対して、展開式を利用する:
exp ( U cos ( x ) ) = ∑ n I n ( U ) exp ( i n x ) \exp(U\cos(x)) = \sum_n I_n(U) \exp(inx) exp ( U cos ( x )) = ∑ n I n ( U ) exp ( in x )
ここでI n ( U ) I_n(U) I n ( U ) は修正ベッセル関数である。
L → ∞ L \to \infty L → ∞ 極限における各項の振る舞いを分析することで、主導項と次主導項を分離する:
主導項:∝ L 2 \propto L^2 ∝ L 2 次主導項:∝ L \propto L ∝ L 高次補正項:o ( L ) o(L) o ( L ) 従来の方法は周期的境界条件に依存していたが、本論文はL → ∞ L \to \infty L → ∞ 極限により境界効果を無視できるようにした。これが重要な技術的突破である。
準周期ポテンシャルΦ ( x ) = U a cos ( 2 π x / a ) + U b cos ( 2 π x / b ) \Phi(x) = U_a\cos(2\pi x/a) + U_b\cos(2\pi x/b) Φ ( x ) = U a cos ( 2 π x / a ) + U b cos ( 2 π x / b ) (ここでa / b a/b a / b は無理数)に対して、二重Jacobi-Anger展開を使用する:
exp ( Φ ( x ) ) = ∑ n , m I n ( U a ) I m ( U b ) e i ( n / a + m / b ) x \exp(\Phi(x)) = \sum_{n,m} I_n(U_a)I_m(U_b)e^{i(n/a+m/b)x} exp ( Φ ( x )) = ∑ n , m I n ( U a ) I m ( U b ) e i ( n / a + m / b ) x
振動積分の界を厳密に評価することで:
∣ ∫ 0 L e i k ⋅ r x g k ( x ) d x ∣ ≤ 2 L ∣ k ⋅ r ∣ \left|\int_0^L e^{ik \cdot rx}g_k(x)dx\right| \leq \frac{2L}{|k \cdot r|} ∫ 0 L e ik ⋅ r x g k ( x ) d x ≤ ∣ k ⋅ r ∣ 2 L L → ∞ L \to \infty L → ∞ における振動項の次主導的地位を証明した。
Φ ( x ) = U cos ( x ) \Phi(x) = U\cos(x) Φ ( x ) = U cos ( x ) に対して、厳密に以下を証明した:
D ∗ = D I 0 2 ( U ) D^* = \frac{D}{I_0^2(U)} D ∗ = I 0 2 ( U ) D
これは古典的なLJ公式と完全に一致する。
Φ ( x ) = U a cos ( 2 π x / a ) + U b cos ( 2 π x / b ) \Phi(x) = U_a\cos(2\pi x/a) + U_b\cos(2\pi x/b) Φ ( x ) = U a cos ( 2 π x / a ) + U b cos ( 2 π x / b ) に対して、以下を証明した:
D ∗ = D I 0 2 ( U a ) I 0 2 ( U b ) D^* = \frac{D}{I_0^2(U_a)I_0^2(U_b)} D ∗ = I 0 2 ( U a ) I 0 2 ( U b ) D
N N N 個の通約不可能な周期を有する一般的な場合:
U ( x ) = ∑ i = 1 N U i cos ( 2 π x a i + θ i ) U(x) = \sum_{i=1}^N U_i\cos\left(\frac{2\pi x}{a_i} + \theta_i\right) U ( x ) = ∑ i = 1 N U i cos ( a i 2 π x + θ i )
に対して、以下を得た:
D ∗ = D ∏ i = 1 N I 0 2 ( U i ) D^* = \frac{D}{\prod_{i=1}^N I_0^2(U_i)} D ∗ = ∏ i = 1 N I 0 2 ( U i ) D
まずU ( x ) = 0 U(x) = 0 U ( x ) = 0 の場合、τ ( x ) = ( L 2 − x 2 ) / ( 2 D ) \tau(x) = (L^2-x^2)/(2D) τ ( x ) = ( L 2 − x 2 ) / ( 2 D ) を検証し、D ∗ = D D^* = D D ∗ = D を確認する。
Jacobi-Anger展開と項別分析を通じて、以下を証明する:
− D τ ( 0 ) = − I 0 ( U ) I 0 ( − U ) L 2 2 + C L + ⋯ -D\tau(0) = -I_0(U)I_0(-U)\frac{L^2}{2} + CL + \cdots − D τ ( 0 ) = − I 0 ( U ) I 0 ( − U ) 2 L 2 + C L + ⋯
ここでC C C は有界定数であり、省略記号はo ( L ) o(L) o ( L ) 項を表す。
同一の分析枠組みを使用して、交差項と振動積分を処理し、最終的に同一の漸近挙動を得る。
概念的革新 :LJ公式が周期ポテンシャルにのみ適用可能という従来の認識を打破した方法的革新 :準周期ポテンシャルを処理するための体系的な数学的枠組みを確立した普遍性 :より一般的なポテンシャル場における拡散問題に理論的基礎を提供した実験的実現 :実験系が厳密な周期性に対する要求を低減した複雑系 :無秩序系、ガラス状態などの複雑な凝縮系の研究に道具を提供した多次元推広 :方法は高次元系に直接推広可能である平均値⟨ exp ( ± Φ ( x ) ) ⟩ = lim L → ∞ L − 1 ∫ 0 L exp ( ± Φ ( x ) ) d x \langle\exp(\pm\Phi(x))\rangle = \lim_{L\to\infty} L^{-1}\int_0^L \exp(\pm\Phi(x))dx ⟨ exp ( ± Φ ( x ))⟩ = lim L → ∞ L − 1 ∫ 0 L exp ( ± Φ ( x )) d x が存在し有限であることが必要である。
調和ポテンシャルΦ ( x ) = k x 2 \Phi(x) = kx^2 Φ ( x ) = k x 2 に対しては、上記の極限が存在しないため、D ∗ = 0 D^* = 0 D ∗ = 0 となり、本方法は束縛ポテンシャルには適用不可能である。
高次元への推広と、より複雑なポテンシャル場の処理には、さらなる研究が必要である。
本論文は厳密な数学的証明を通じて、LJ公式が準周期ポテンシャルにおいて厳密に適用可能であることを確立した。この結果は重要な理論的価値を有するだけでなく、実験研究に新たな可能性をもたらす。本研究は、より一般的な条件下でのブラウン運動の理解に堅実な理論的基礎を提供し、統計物理学、凝縮系物理学などの分野に重要な影響を与えることが予想される。
今後の研究方向には以下が含まれる:
高次元系への推広 異常拡散への応用 傾斜ポテンシャル場における巨大拡散現象 実験的検証と具体的物理系への応用 数学的厳密性 :完全な解析的証明を提供し、理論的空白を埋めた方法的革新性 :Jacobi-Anger展開法の応用は独創的である結果の重要性 :古典理論の限界を突破し、里程碑的意義を有する適用性の広さ :方法は多種の複雑系に推広可能である適用条件の制限 :特定の数学的条件(平均値の存在性)が必要である実験的検証の欠如 :理論結果はまだ実験的検証が必要である次主導項の物理的意義 :次主導項に関する物理的意義の議論が不十分であるこれは統計物理学分野における重要な理論的突破であり、当該分野のさらなる研究熱潮を引き起こすことが予想され、ブラウン運動理論の発展に深遠な影響を与えるであろう。