Spectral Analysis of Brownian Motion with its Rheological Analogues

技術要約：ブラウン運動とそのレオロジー的類似体に関するスペクトル解析

問題提起
本論文は、線形かつ等方的な粘弾性材料中に浸漬されたプローブ微粒子のブラウン運動のパワースペクトル（パワースペクトル密度、PSD）を特徴付ける際の課題に取り組んでいる。記憶のないニュートン流体におけるブラウン運動の長期的な拡散挙動は確立されており（Einstein, 1905）、そのような流体における速度自己相関関数も既知であるが（Ornstein, 1917; Wang and Uhlenbeck, 1945）、周囲の媒体が粘弾性、流体力学的メモリ、または劣拡散性を示す場合、スペクトル解析は著しく複雑になる。従来のアプローチでは、速度相関を導出するために、時間領域で複雑な積分方程式や一般化ランジュバン方程式を直接解いてからフーリエ変換を行う必要があることが多い。本論文は、マックスウェル流体、ジェフリーズ流体、劣拡散性材料、および流体力学的メモリが顕著な高密度粘性流体を含む、多様なレオロジー環境に対して、PSDを算出するためのより合理的な手法を求めている。

手法
核心となる手法は、ブラウン運動に対する粘性–粘弾性対応原理を採用している。著者は、物理系を表現するために特定のレオロジー的類似体を合成している：

1. レオロジーネットワークの構築： 本論文は、質量 $m$ および半径 $R$ の粒子の粘弾性材料中でのブラウン運動を、以下の2つの要素の並列接続としてモデル化できると仮定している：
   • 線形粘弾性材料自体（複素動的弾性率 $G_{ve}(\omega)$ によって特徴付けられる）。
   • イナーター（inerter）（力が相対加速度に比例する機械的要素）であり、分散慣性 $m_R = m / (6\pi R)$ を持つ。
2. 複素動的流動性： 解析は、複素動的流動性（または複素移動度）に焦点を当てている。これは、合成された並列ネットワーク（粘弾性材料 + イナーター）の複素動的弾性率 $G(\omega)$ を用いて $\phi(\omega) = i\omega / G(\omega)$ と定義される。
3. スペクトルの導出： ブラウン運動のパワースペクトル密度 $S(\omega)$ は、この複素動的流動性の実部と直接比例することが示されている：
   $$S(\omega) = \frac{N k_B T}{3\pi R} \Re\{\phi(\omega)\}$$
   ここで、$N$ は空間次元数、$k_B$ はボルツマン定数、$T$ は温度である。
4. 特定モデルへの適用： この枠組みは、以下の4つのケースにおけるPSDを導出するために適用される：
   • 記憶のない粘性流体： 古典的なローレンツ型スペクトルに対する手法の妥当性を検証する。
   • 調和トラップ（ケルビン固体）： 調和ポテンシャル内にある粒子をモデル化する。
   • マックスウェル流体： 直列に接続されたスプリングとダッシュポット。
   • ジェフリーズ流体： ダッシュポットと並列に接続されたマックスウェル要素。
   • 劣拡散性材料： べき乗則の平均二乗変位をモデル化するために、分数微積分と「スプリングポット（springpot）」要素（スコット・ブレア流体）を利用する。
   • 流体力学的メモリ： 変位した流体が長距離相関を生み出す高密度流体をモデル化する。これは、ランジュバン方程式における分数階微分項（次数1/2）によって表され、特定の分数要素を含むレオロジーネットワークに対応する。

主な貢献と結果

• 統一された枠組み： 本論文は、任意の線形等方粘弾性材料におけるブラウン運動のPSDは、粘弾性材料とイナーターが並列に接続されたレオロジーネットワークの複素動的流動性の実部に比例することを確立している。これにより、様々な複雑な流体の扱いを単一の対応原理の下に統一している。
• 解析解： 本論文は、いくつかの複雑なシナリオに対して明示的な閉形式の式を提供している：
  • マックスウェル流体： 直列のスプリングの剛性が増加するにつれて、スペクトルが記憶のないニュートン限界に収束することを示す。
  • ジェフリーズ流体： 無次元の粘性比 $\xi = \eta_\infty / \eta$ の関数としてスペクトルを導出し、高周波挙動が並列ダッシュポットによって支配される様子を示す。
  • 劣拡散性材料： PSDは分数指数 $\alpha$ と特性時間スケール $\lambda$ を用いて表現され、$\alpha = 1$ のときにニュートン流体の結果を回収する。
  • 流体力学的メモリ： バセット（Basset）履歴力を表す分数階微分項をレオロジー的類似体に組み込むことで、高密度粘性流体のPSDを導出する。得られるスペクトルは、無次元密度比パラメータ $\gamma$ に依存する。
• 可視化： 結果は、無次元周波数に対してプロットされた正規化パワースペクトルとして提示され、材料パラメータ（緩和時間、剛性、粘性比、密度比など）によってスペクトルの形状がどのように変化するかを示している。

意義と主張
本論文は、この特定のレオロジー的類似体（粘弾性材料 + イナーターの並列接続）の合成が、ブラウン運動のパワースペクトルの計算を大幅に簡略化すると主張している。時間領域で速度自己相関関数の複雑な積分方程式を解いてからフーリエ変換を行う代わりに、この手法を用いることで、複素動的流動性を通じてPSDを直接計算することが可能となる。

著者は、この対応原理が、ブラウン粒子が媒体と力を交換する具体的なメカニズム（粘性、慣性、流体力学的メモリ、または粘弾性）に関わらず、「包括的な妥当性」を持つことを強調している。本論文は新しい実験装置や将来の応用を提案するものではなく、マックスウェル、ジェフリーズ、劣拡散、および流体力学的メモリを持つ流体における既知の結果（ニュートン流体や調和トラップなど）を回収しつつ、複雑な流体中の微粒子追跡のスペクトル解析を効率化するための理論的ツールを提供するものである。