Thinning-by-spinning: shear rheology of dense chiral fluids

技術要約：スピンによる希薄化：高密度カイラル流体の剪断レオロジー

問題提起
自己回転粒子からなるカイラル流体は、回転が奇（横方向）の力を発生させ、相挙動やダイナミクスを再形成する非平衡状態の物質の多用途なプラットフォームを代表するものである。自己回転が六角形秩序を破壊し欠陥の増殖を促進することで、2次元固体を流動化させることは確立されているが、自己回転と局所的な剪断との等価性、およびこのようなカイラル流体が外部から課される剪断にどのように応答するかについては、依然として根本的な疑問が残っている。特に、高密度系（固体および流体領域の両方にわたる）における固有のカイラリティと外部駆動の相互作用には、統一されたレオロジーの枠組みが必要である。

手法
著者らは、2次元における粒子ベースのシミュレーションを用いて、モデルカイラル流体の線形および非線形レオロジーを調査している。系は、$N=4096$個の球状粒子（ディスク）で構成され、非保存的な横方向の力が加わったレナード＝ジョーンズ・ポテンシャルによって相互作用する。ダイナミクスは、以下の要素を含む過減衰ランジュバン方程式によって支配される：

1. 保存的相互作用： レナード＝ジョーンズ・ポテンシャル $U(r_{ij})$ から導かれる。
2. 横方向の力： $f_{ij} = -\Gamma u_{ij}$ と定義される。ここで、$u_{ij}$ は分離ベクトルと平面に垂直な単位ベクトルの外積に比例し、カイラリティパラメータ $\Omega$ に依存するトルクを導入する。
3. 外部剪断： 剪断速度 $\dot{\gamma}$ を伴うリースの・エドワーズ境界条件を介して実装される。

シミュレーションは、広範な剪断速度（$10^{-5} \leq \dot{\gamma} \leq 1$）およびカイラリティの強さ（$\Omega$）を探索し、正（剪断渦と一致）および負（剪断渦と逆）の両方の値の範囲をカバーしている。系は、固体領域（高密度 $\rho=0.7$, 低温 $k_BT=0.35$）と流体領域（$\rho=0.6$, $k_BT=0.47$）の両方で分析される。

主要な観測量には以下が含まれる：

• 応力テンソル： アーヴィング・カークウッド応力テンソル $\sigma_{\alpha\beta}$ から、対称化された剪断応力 $\sigma_s$ および粘度 $\eta$ が導出される。
• 構造的秩序： 局所的および全域的な六角形秩序パラメータ（$\psi_6$ および $\Psi_6$）を用い、相転移と欠陥の増殖を監視する。
• 有効温度 ($T_{eff}$)： ゆらぎ散逸定理（FDT）の破れ、具体的にはストークス・アインシュタインの関係式 $k_B T_{eff} = D/\mu$ を用いて算出される。ここで $D$ は拡散係数、$\mu$ は微小な外部力に対する線形応答から得られる移動度である。

主要な貢献および結果

1. 固体の流動化と剪断希薄化：
   固体領域において、系は降伏応力挙動を示す（平衡状態 $\Omega=0$ では低剪断速度において流動しない）。しかし、カイラリティパラメータ $|\Omega|$ を増加させると、系は流動化し、ゼロ剪断極限において有限の粘度を持つようになる。剪断速度およびカイラリティの増加の両方が粘度を減少させ、系を剪断希薄化領域へと駆動する。粘度のスケーリングは $\eta \sim 1/\dot{\gamma}$（モード結合理論と一致）および $\eta \sim 1/|\Omega|$ に従い、カイラリティが固有の剪断応力の源として機能することを示している。

2. スピンによる希薄化メカニズム：
   著者らは、「スピンによる希薄化」メカニズムを特定した。これは、固定された剪断速度においても、粘度がカイラリティパラメータ $|\Omega|$ の大きさに伴って単調に減少する現象である。この効果はマクロなものであるが、純粋に微視的な対カイラル相互作用によって駆動されている。

3. 有効温度による統一的枠組み：
   線形応答領域において、著者らはグリーン・クボ関係の一般化を提案している。熱力学的温度 $T$ をカイラリティ依存の有効温度 $T_{eff}(\Omega)$ に置き換えることで、平衡状態のグリーン・クボ積分が非平衡剪断粘度を定量的に予測できる。これは、カイラリティが高次元のスケールにおいて、上昇した有効温度として捉えられる形で、系にエネルギーを注入していることを示唆している。

4. フローカーブの崩壊（コラップス）：
   線形応答を超えた領域では、粘度を $T_{eff}$ でスケールし、剪断速度を $|\Omega|$ でスケールすることで、異なる剪断速度およびスピン速度に対するフローカーブが単一のマスタカーブへと崩壊（コラップス）する。この崩壊は、自己回転が外部変形と直接競合する固有の応力源として機能し、剪断とスピンによる流動化の経路を統一していることを実証している。

5. ハンドネスの非対称性と層形成：
   大きな強制力が加わる際、スピンと剪断の対応関係は崩れ、顕著なハンドネスの非対称性が現れる。

• 逆向きのハンドネス ($\Omega < 0$)： 横方向の力が課された剪断に反対する場合、系は臨界剪断速度 $\dot{\gamma}^*$ において構造転移を起こす。流体は、紐状の準1次元的な秩序ドメイン（層形成）へと自己組織化する。この構造は流動を促進し、応力の著しい緩和と剪断希薄化の強化をもたらす。
• 一致するハンドネス ($\Omega > 0$)： 系は剪断方向に適度に伸長した六角形パッチを示すが、非常に高い剪断速度に達するまでは、上述のような明確な紐状のチャネルを形成しない。
• この層形成領域への転移は、平均六角形秩序パラメータ $|\psi_6|$ の極小値と、構造因子のピークの急激な増加によって特徴付けられる。

意義
本論文は、高密度粒子系における流動化の一般的なメカニズムとしてカイラリティを確立している。著者らは、カイラリティと外部駆動が互いに補強し合い、固体および流体の両相にわたって剪断希薄挙動を生み出す、統一されたレオロジーの枠組みを提供している。これらの知見は、線形応答が有効温度によって記述できる一方で、非線形領域は外部変形と固有のカイラリティとの間の競合によって支配され、相対的なハンドネスに依存した特異な構造状態（紐状の流動チャネルなど）を生成することを実証している。これらの結果は、横方向の相互作用が材料特性を根本的に変容させる、高密度カイラル懸濁液、アクティブマター、および奇（オッド）機械的メタマテリアルのレオロジーを理解するための参照となるものである。