Optimal probing scale for current fluctuations in a Brownian gyrator

技術要約：ブラウン・ジャイレーターにおける電流揺らぎの最適プロービング・スケール

問題設定
本論文は、駆動されたコロイド系、特にブラウン・ジャイレーターにおける循環電流の揺らぎを境界付ける問題を扱っている。熱力学的不確定性関係（TUR）やクラメール・ラオ境界から導かれるような電流揺らぎの境界は、通常、分母に散逸（コスト）項、分子に応答項をとる比の形式をとる。この境界の大きさは、摂動に対する電流の応答強度と、それに伴うオンサーガー・マッハラップ・コストに依存するが、摂動の空間的な位置（「プローブ」）は理論によって固定されていない。著者は、単位コストあたりの応答を最大化するために、摂動を適用すべき最適な半径方向の距離が存在すると仮定している。目標は、特定のクラスの摂動を受けるブラウン・ジャイレーターに対して、最適なプロービング半径（$r^*$）を決定することである。

手法
本研究では、ブラウン・ジャイレーターを表す二次元線形拡散モデルを用いた、解析的な導出と数値シミュレーションの組み合わせを用いている。この系は、トラップ剛性 $\gamma$ と拡散定数 $D$ を持つ調和トラップ内の粒子で構成され、回転率 $\Omega$ を持つ非保存的な循環力にさらされている。

1. プローブ設計: 著者は、特定の摂動場として、ガウス型エンベロープを持つ四重極せん断、$u_\lambda = \lambda g(r) (2xy, x^2-y^2)^\top$（ここで $g(r) = e^{-\alpha r^2/2\sigma_0^2}$）を導入する。このプローブは、以下の2つの重要な特性のために選択されている：

   • 閉じ込め: ガウス型エンベロープにより、正および負の両方の摂動強度に対して、摂動されたプロセスが非爆発的かつ指数関数的にエルゴード的であることが保証される。
   • 線形応答への盲目性: 角運動量の直交性により（このプローブは $m=\pm 3$ の角セクターに投影されるが、無摂動時の電流は等方的または $m=0$ である）、有限の観測窓における電流の一次線形応答は正確に消失する。その結果、主要な信号は純粋に二次的なものとなり、一次項を差し引く必要がなくなる。

2. 評価指標: 最適化の対象は、二次の応答の二乗（$\Delta \langle J \rangle^2$）と、オンサーガー・マッハラップ・コストの比である。振幅ゼロの極限において、これは二次バタチャリヤ境界 $B = \chi_2^2 / (2 E_0[Q^2])$ に帰着する（ここで $\chi_2$ は二次感受率、 $Q$ は積分されたコストである）。

3. 解析解: 問題は、以下の2つの独立したセクターに分離されることが示される：

   • 半径方向セクター: コスト項は半径座標にのみ依存し、駆動 $\Omega$ には依存しない。これは定常ガウス分布を用いて厳密に解かれる。
   • 角度方向セクター: 応答は $m=3$ の角セクターに限定される。摂動密度の支配方程式は、パラメータ $b=4$ および $a = \frac{3}{2}(1 + i\Omega/\gamma)$ を持つクンマーの微分方程式（合流型超幾何微分方程式）へと簡約される。感受率 $\chi_2$ は、エンベロープの幅パラメータ $\alpha$ の関数として、超幾何関数を用いて閉形式で導出される。

4. 最適化: 評価指標 $B(\alpha)$ をエンベロープの幅 $\alpha$ に関して最大化する。その後、最適な半径 $r^*$ は、関係式 $r^* = \sigma_0 \sqrt{2/\alpha^*}$ を通じて決定される。

主な結果
解析は、系の2つの特徴的なタイムスケール、すなわちトラップ緩和時間（$1/\gamma$）と回転時間（$1/\Omega$）の競合によって支配される閉形式の式を与える。

• 弱駆動極限 ($\Omega \ll \gamma$): 最適半径はトラップの熱的半径にロックされる。
  $$r^* \approx 1.1264 \sqrt{D/\gamma}$$
  このレジームでは、トラップ剛性 $\gamma$ がスケールを定義する。

• 強駆動極限 ($\Omega \gg \gamma$): 最適半径は、1ラジアン回転する間に拡散する距離にスケールする。極めて重要なことに、このときトラップ剛性 $\gamma$ は式の記述から完全に消失する。
  $$r^* \approx 1.5563 \sqrt{D/\Omega}$$
  係数 $1.5563$は、$\sin$関数と$\cos$関数を含む方程式の根$c_0$を用いて、$\sqrt{2/c_0}$ として解析的に導出される。

• クロスオーバー: これら2つのレジーム間の遷移は $\Omega = \gamma$ で起こる。最適半径は、2桁の駆動比にわたってわずかに変化するのみ（$1.13 \sigma_0$ から $0.92 \sigma_0$ の間）であり、単純な冪乗則をクロスオーバー全体に適合させても、いずれの極限も正確に記述することはできない。

検証
理論的予測は、オイラー・マリヤママ法およびストラトノヴィッチの中点則を用いた直接的な数値シミュレーションによって検証されている。

• シミュレーションは、弱駆動および強駆動の各極限における係数を裏付けている。
• 最適半径の駆動強度への依存性は、様々な $\Omega/\gamma$ 比におけるバタチャリヤ境界曲線のピークを測定することによって検証されている。
• 有限時間効果（観測窓 $T$）が分析されており、有限時間の測定と定常理論との間の欠損は、系の最も遅い緩和モード $|\lambda_0| = 3\sqrt{\gamma^2 + \Omega^2}$ によって制御されていることが示されている。

意義と主張
本論文は、駆動された系における電流揺らぎをプロービングするための最適な空間スケールを、初めて厳密に決定したものであると主張している。その主要な貢献は、最適なプローブの位置がトラップの幾何学的形状のみによって決まるのではなく、系の固有のタイムスケールの相互作用によって決まることを示した点にある。

特筆すべき具体的な知見は、強駆動レジームにおいて、最適なプロービング・スケールがトラップ剛性 $\gamma$ に依存しなくなることである。代わりに、それは拡散定数と回転率によってのみ決定され、粒子が1ラジアン回転する間に拡散する距離を表している。著者は、この結果の扱いやすさは、線形応答を盲目にするという特定のプローブの対称性に依存していることを強調している。これにより、問題が半径方向のコストと角度方向の応答へと明確に分離することが可能となった。本研究は、特定の定数はプローブの対称性（ここでは $m=3$）に依存するものの、半径方向のコストと角度方向の応答の分離は、堅牢な構造的特徴であることを示唆している。