2025-11-14T06:52:10.379364

Flocking transition in phoretically interacting active particles with pinning disorder

Adhikary, Subramaniam, Singh

Recent studies in the collective behavior of active colloids have shown that a global polar order may emerge due to long-ranged chemo-repulsive interactions between them. Here, we report the role of pinning disorder in the flocking transition for such a system. To this end, we study the problem of chemically interacting active colloids with some fraction of the colloids randomly pinned over space such that they can only rotate while phoretically interacting with other particles. Using this model, we investigate the sustenance of global polar order in the presence of quenched disorder. We quantify the flocking transition by studying the global polarization, and the role of finite-size effects. We find that in the crystalline flocking phase, even a small fraction of pinning can destroy spatial crystalline order, although polar order in the form of a liquid phase is maintained. It is observed that polar order is sustained in a system with a higher pinning fraction if the long-ranged repulsive force is subsequently increased. However, in absence of chemo-repulsive forces between particles, polar order drastically decreases even with a smaller pinning fraction. Overall, this work suggests a novel route of solid-to-liquid transition that can be induced via "translationally inert" obstacles, that rotate but do not translate whilst interacting with the bulk.

academic

化学相互作用する活性粒子における固定無秩序を伴うフロッキング転移

基本情報

論文ID: 2511.10386
タイトル: Flocking transition in phoretically interacting active particles with pinning disorder
著者: Sagarika Adhikary, Arvin Gopal Subramaniam, Rajesh Singh (インド工科大学マドラス校)
分類: cond-mat.soft, cond-mat.stat-mech
投稿日時: 2025年11月13日
論文リンク: https://arxiv.org/abs/2511.10386

要旨

本研究は、化学相互作用する活性コロイド系における固定無秩序(pinning disorder)がフロッキング転移(flocking transition)に及ぼす影響を検討している。結晶フロッキング相において、わずかな比率の固定であっても空間結晶秩序を破壊できるが、液相形態で極性秩序を維持することが明らかになった。長距離反発力が増強されると、システムはより高い固定比率でも極性秩序を維持できる。本研究は、「並進惰性」障害物によって誘起される固-液転移の新しい経路を提案している。これらの障害物は回転のみが可能で平行移動はできず、同時に全体と相互作用する。

研究背景と動機

中核的問題

主要研究課題：淬火無秩序(quenched disorder)が存在する場合、化学反発相互作用を持つ活性コロイド系はいかにして全体的な極性秩序とフロッキング行動を維持するのか？
問題の重要性：
- 実際の活性系は通常、理想的な均一環境ではなく、複雑性と異質性を示す
- 幾何学的無秩序などの複雑な条件下で集団運動がいかに出現し持続するかを理解することは基礎科学の問題である
- 制御可能な活性材料と相転移材料の設計に応用価値がある
既存研究の限界：
- ほとんどの自駆動粒子研究は均一環境で実施されている
- 最近発見された化学反発結晶フロッキング(CCF)の無秩序環境での堅牢性は未知である
- 異なるタイプの無秩序が活性物質の集団行動にいかに影響するかについて、体系的な理解が欠けている
研究動機：
- 新しい形態の無秩序を探索する：ランダムに固定された粒子で、回転のみが可能で平行移動はできない
- この「並進惰性」障害物がフロッキング転移に及ぼす影響メカニズムを研究する
- 相互作用、無秩序、涌現秩序間の相互作用を理解するための完全な相図を構築する

中核的貢献

新しい相転移メカニズムの発見：固定無秩序によって誘起される固-液転移の新しい経路を提案し、わずかな固定比率(np ~ 0.01)でも結晶フロッキングを液体フロッキングに変換できることを示した
体系的な相図の構築：(np, Λt)および(Λr, Λt)パラメータ空間における完全な相図を構築し、極性秩序、六角秩序、密度分散、磁化率などの複数の秩序パラメータを定量化した
相互作用-無秩序競争の解明：長距離反発力(Λt)が固定無秩序の破壊作用を補償でき、より高いΛt値がより高い固定比率でも極性秩序を維持できることを発見した
異なる相互作用効果の区別：長距離反発力が欠ける場合(Λt=0)、システムが固定無秩序に対してより敏感であり、固定が密度帯の形成を抑制することを証明した
理論的メカニズムの説明：衝突滑動メカニズムの理解に基づき、固定粒子がいかに増強反発を減少させることで結晶局在化を破壊するかを説明した

方法論の詳細

タスク定義

二次元空間における化学相互作用する活性コロイドN個の集団動力学を研究し、そのうちNp個の粒子がランダムに固定され(回転のみ可能)、Nm個の粒子が自由に運動できる(Nm + Np = N)。目標は、固定比率np = Np/Nがフロッキング転移と空間構造に及ぼす影響を定量化することである。

モデルアーキテクチャ

1. 運動方程式

第i番目の粒子の位置riと方向eiの進化は以下に従う：

運動粒子(i = 1,...,Nm)の場合： $\dot{r}_i = v_s e_i + \zeta_t J_i + \mu F_i$

すべての粒子の方向進化： $\dot{e}_i = [\zeta_r (e_i \times J_i) + \eta^r_i] \times e_i$

固定粒子(i = Nm+1,...,N)の場合： $\dot{r}_i = 0$

2. 相互作用メカニズム

化学泳動： $J_i = -D_c[\nabla c(r,t)]_{r=r_i}$ 、ここで濃度場は以下を満たす： $D_c\nabla^2 c(r,t) + \sum_{i=1}^N \lambda_0 \delta(r-r_i) = 0$
解析解は以下の通り： $J_i = \frac{\lambda_0}{4\pi} \sum_{j\neq i} \frac{r_i - r_j}{|r_i - r_j|^3}$
反発トルク：ζr > 0の場合、粒子は化学相互作用を通じて相互に遠ざかる
反発力：ζt > 0の場合、位置動力学に反発を提供する
短距離反発： $F_i = -\nabla_i U$ 、ここで $U = \sum_{i<j} U_e(r_i, r_j)$ $U_e = \kappa(r_{ij} - 2b)^2, \quad \text{if } |r_i - r_j| < 2b$

3. 主要な無次元パラメータ

固定比率： $n_p = N_p/N$
回転反発強度： $\Lambda_r = \zeta_r/(b^3 v_s)$
並進反発強度： $\Lambda_t = \zeta_t/(b^4 v_s)$

技術的革新点

新しい無秩序モデル：従来の静的障害物と異なり、固定粒子は回転能力を保持し、化学場を通じてシステムと相互作用し、実験的実現可能性がより高い(光ピンセット技術など)
多スケール秩序パラメータ分析：
- 全体的極性秩序： $m = \langle M \rangle_{ss}$ 、 $M = |\frac{1}{N}\sum_i e_i|$
- 六角秩序： $\psi_6 = \frac{1}{N}\sum_i \psi_i$ 、 $\psi_i = \frac{1}{N^n_i}\sum_j e^{i6\theta_{ij}}$
- 密度分散： $\sigma = \langle \rho_l^2 \rangle - \langle \rho_l \rangle^2$
- 磁化率： $\chi = L^2[\langle M^2 \rangle - \langle M \rangle^2]$
極性対相関関数： $g(r,\phi)$ は位置と方向情報の両方を含み、 $\cos(\phi_{ij}) = e_i \cdot \hat{r}_{ij}$ 、結晶フロッキングと液体フロッキングを区別できる

実験設定

シミュレーションパラメータ

パラメータ	値	説明
粒子数N	117-7514	異なる実験要件に応じて
システムサイズL	32-256	有限サイズ効果研究
面積分率ϕ	0.36 (主要)	$\phi = N\pi b^2/L^2$
自駆動速度vs	50	固定
粒子半径b	1	固定単位
剛性係数κ	175	短距離反発強度
回転拡散Dr	10⁻⁴	ノイズ強度
時間ステップdt	0.01	Euler-Maruyama積分器

数値方法

初期条件：粒子位置と方向はランダムに分布し、Np個の粒子がランダムに選択されて固定される
境界条件：周期的境界条件
定常状態平均：瞬時的状態を除外した後、2×10⁵時間ステップで平均化
ペクレ数：Pe = vs/(bDr) ~ 10⁵、決定論的効果が支配的であることを確認

評価指標

相転移の識別：磁化率χのピーク位置を通じて臨界点を決定
構造特性：
- 動径分布関数g(r)のピーク高さと位置
- 六角秩序パラメータψ₆
- 極性対相関g(r,φ)の強度分布
有限サイズスケーリング：m(np,L)のシステムサイズに対する収束挙動を研究

実験結果

主要な結果

1. 3つの主要な相状態の識別(図2)

(np, Λt)相図において3つの顕著な相状態を識別：

結晶フロッキング(Crystalline Flock, np=0)：
- 極性秩序m ≈ 0.95
- 六角秩序ψ₆ ≈ 0.8
- g(r)は複数の鋭いピークを示す
- 粒子は六角格子構造を形成
極性液体(Polar Liquid, np=0.15)：
- 極性秩序m ≈ 0.75
- 六角秩序ψ₆ ≈ 0.1(結晶秩序が破壊される)
- g(r)ピークは幅広化し、強度が低下
- 集団運動を保持するが空間周期性を失う
無秩序相(Random Phase, np=0.36)：
- 極性秩序m ≈ 0.2
- 磁化率χが最大値に達する(臨界変動)
- 粒子運動がランダム化

2. 固-液転移の感度(図3)

np = 0.01(わずか1%の固定)でも結晶秩序を破壊：
- ψ₆が0.8から0.2に低下
- g(r,φ)が局在化した鋭いピークから分散分布へ転移
- 衝突に対する極性相関が結晶パターンから液体パターンへ変化
動径分布関数の変化：
- np = 0：複数の鋭いピーク、長距離結晶秩序を反映
- np = 0.01：ピーク幅広化、高さ低下
- np = 0.15：最初のピークが小さいrへ移動、粒子間距離が減少

3. 相互作用-無秩序競争(図4c-e)

(Λr, Λt)相図においてnp=0.15を固定：

フロッキングを維持するにはΛtとΛrの協調作用が必要
臨界線の比較：
- np = 0(青線)：低いΛtでもフロッキング可能
- np = 0.15：臨界線が右へシフト
- np = 0.5(赤線)：より高いΛtとΛrが必要
密度分散σ：
- 秩序相：σ < 0.05
- 転移領域：σが急速に増加
- 無秩序相：σ > 0.15(特に低Λr領域)

有限サイズ効果分析(図2g-h)

システムサイズ依存性(L = 32-256、ϕ=0.36固定)：

小さいシステム(L=32)：転移が滑らか、臨界領域が広い
大きいシステム(L=256)：転移がより鋭い、臨界点nc ≈ 0.35
有限サイズスケーリングは明確な熱力学相転移の存在を示唆

密度依存性(ϕ = 0.12-0.48、L=128固定)：

相転移位置は基本的に変わらない
転移特性は密度に対して不敏感
固定効果が支配的であり、密度効果ではないことを示す

長距離力がない場合(Λt=0)の分析(図5)

主要な発見：

臨界固定比率の大幅な低下：
- Λtがある場合：nc ≈ 0.35
- Λtがない場合：nc ≈ 0.16(約50%低下)
密度帯の抑制：
- np = 0：運動方向に垂直な密度帯を形成
- np = 0.03：密度帯が消失、液体フロッキングへ転移
- Vicsekモデルにおける障害物が帯状構造を抑制するという発見と一致
磁化率分布：
- 転移領域でχが最大
- 無秩序相側へより広く拡張
- 臨界変動の非対称性を反映

アブレーション実験

相互作用成分の役割：

条件	臨界np	主要相状態	構造特性
完全モデル(Λt>0, Λr>0)	~0.35	結晶→液体→無秩序	六角格子→液体フロッキング
長距離力なし(Λt=0)	~0.16	密度帯→液体→無秩序	帯状→均一液体
固定のみ(np↑, 固定Λ)	連続的に低下	秩序パラメータの単調低下	段階的無秩序化

物理的メカニズムの説明：

Λtが提供する増強反発が固定による局在化損失を補償
固定粒子は衝突時の増強反発に寄与しない(固定されているため)
滑動長 ∝ b/Λr、固定が有効滑動時間を増加させ、結晶を破壊

結論と議論

主要な結論

新しい相転移経路：固定無秩序誘起の結晶-液体転移、従来の熱力学相転移または密度駆動相転移と異なる
極めて高い感度：わずか1%の固定でも空間結晶秩序を破壊できるが、極性秩序はより堅牢(約35%の固定が必要)
相互作用補償メカニズム：長距離反発力Λtが固定の破壊作用を部分的に補償でき、制御手段を提供
構造依存性：
- 長距離力あり：結晶→液体→無秩序
- 長距離力なし：密度帯→液体→無秩序(固定に対してより敏感)
微視的メカニズム：衝突滑動長の理解に基づき、固定が増強反発を減少させることで結晶局在化を破壊

限界

二次元制限：モデルは二次元システムのみを考慮し、三次元の場合は異なる挙動の可能性
特定の相互作用形式：化学泳1/r²相互作用、他の長距離相互作用(1/rなど)の挙動は未知
固定モデルの簡略化：
- 固定粒子は完全に固定(実際には微小変位の可能性)
- 固定強度の分布を考慮していない
- 固定位置は完全にランダム(実際には相関の可能性)
有限サイズ効果：
- 最大システムN~7500、有限サイズ効果を完全に排除するには不十分の可能性
- 臨界指数の体系的分析が未実施
- 普遍性クラスの帰属が不明確
動力学分析の不足：
- 主に定常状態特性に焦点
- 転移動力学と緩和時間の詳細研究が不足
- 瞬時的挙動が重要な情報を含む可能性
実験検証の欠落：純粋な理論/シミュレーション研究、実験との定量的比較が欠ける

将来の方向

実験的実現：
- マイクロフルイディクス技術による障害物製造
- 光ピンセット技術による制御可能な固定実現
- 固定-活性粒子相互作用の調整が必要
理論の深化：
- 有限サイズスケーリング分析による普遍性クラスの決定
- 臨界指数の精密測定
- 連続場理論による記述
研究の拡張：
- 三次元システム
- 他の長距離相互作用形式
- 固定強度と空間分布の影響
- 動的固定(時間依存)
応用探索：
- 活性「相転移材料」の設計
- 調整可能な固-液転移
- マイクロフルイディクスにおける流動制御

深度評価

利点

1. 科学的革新性

新規な無秩序モデル：回転可能だが平行移動できない固定粒子を提案し、従来の静的障害物と異なり、物理的実現可能性がより高い
重要な発見：わずか1%の固定でも結晶秩序を破壊するが極性秩序を保持、2つの秩序の異なる堅牢性を明らかにした
メカニズム洞察：衝突滑動メカニズムに基づいて相転移を説明し、微視的動力学と巨視的相挙動を結びつけた

2. 研究の体系性

完全な相図：多次元パラメータ空間(np, Λt, Λr)の体系的走査
複数の秩序パラメータ表現：極性秩序、六角秩序、密度分散、磁化率など多角的分析
対比研究：長距離力の有無による対比が異なる相互作用の役割を明らかにした

3. 技術的厳密性

有限サイズ効果：L=32-256のシステムサイズ研究
統計的充分性：2×10⁵ステップの定常状態平均、1000個の独立構成
パラメータの透明性：詳細なパラメータ表(表I)により再現性を確保

4. 物理的図像の明確性

優れた可視化(粒子構成、相図、相関関数)
時間進化が相転移動力学を示す
補足動画材料が理解を深める

不足点

1. 理論的深さ

解析理論の欠落：主に数値シミュレーションに依存し、連続場理論または有効理論が欠ける
臨界現象分析の不足：相転移の種類(一次/二次)と普遍性クラスが未決定
スケーリング則の欠落：臨界指数β, γ, νなどが抽出されていない

2. 実験との関連

実験検証の欠落：純粋な理論研究、実験との定量的比較が欠ける
パラメータマッピングの曖昧性：シミュレーションパラメータと実際の実験系の対応が不明確
実現可能性の議論不足：光ピンセットやマイクロフルイディクスでの実現方法の詳細な議論が不足

3. 方法の限界

二次元制限：三次元効果が大きく異なる可能性(位相欠陥など)
特定の相互作用：1/r²化学泳相互作用、普遍性の検討が必要
ノイズ強度：Dr=10⁻⁴は非常に小さく、より大きなノイズ下の挙動は未知

4. 分析の詳細

動力学情報の不足：緩和時間、特性時間スケールの体系的研究が不足
変動スペクトルの欠落：密度または秩序パラメータの時空相関分析が欠ける
位相欠陥の未検討：相転移における欠陥の役割が議論されていない

5. 執筆面

一部の章の構成が最適化可能(III.H節のタイトルが長い)
一部の記号定義がより早く導入できる
文献34との関係がより明確に述べられるべき

影響力評価

1. 学術的貢献

分野の重要性：★★★★☆
- 複雑な環境における活性物質の挙動に新しい視点を提供
- 軟物質物理、非平衡統計、相転移理論を結びつける
方法の革新性：★★★★☆
- 固定-活性相互作用モデルが新規
- 複数秩序パラメータ分析フレームワークが完善

2. 実用価値

応用可能性：★★★☆☆
- 活性相転移材料設計の理論基礎
- マイクロフルイディクスにおける活性流体制御
- さらなる実験検証が必要

3. 再現性：★★★★★

詳細なパラメータ表とシミュレーション設定
方程式とアルゴリズムの記述が明確
補足動画材料を提供

4. 後続研究価値

高い：複数の明確な将来方向
- 臨界現象の深い研究
- 三次元への拡張
- 実験的実現
- 他の相互作用形式

適用場面

1. 直接的適用

化学活性コロイド系の理論予測
マイクロフルイディクス装置における粒子操作設計
生物系における障害物が集団運動に及ぼす影響(細胞遊走など)

2. 潜在的拡張

多孔質媒体における活性物質の輸送
障害環境における群ロボットの協調
交通流における固定障害物の影響

3. 方法論の参考

複数秩序パラメータ相図構築方法
有限サイズスケーリング分析フレームワーク
活性-無秩序相互作用のモデリング思考

推奨読者

強く推奨：軟物質物理、活性物質、非平衡統計物理研究者
推奨：複雑系、集団行動、相転移理論研究者
参考価値：実験軟物質物理学者、マイクロフルイディクス技術者

参考文献(精選)

7 J. Toner, The Physics of Flocking: Birth, Death, and Flight in Active Matter (Cambridge, 2024) - フロッキング物理の総説

30 S. Das et al., Flocking by turning away, Phys. Rev. X 14, 031008 (2024) - 反発トルクフロッキングメカニズム

34 A. G. Subramaniam et al., Minimal mechanism for fluidic flocks, Soft Matter (2025) - 本研究の基礎モデル

36 A. Morin et al., Distortion and destruction of colloidal flocks in disordered environments, Nat. Phys. 13, 63 (2017) - 無秩序環境におけるフロッキング

57 J. Codina et al., Small obstacle in a large polar flock, Phys. Rev. Lett. 128, 218001 (2022) - 障害物がフロッキングに及ぼす影響

総合評価：これは活性物質分野における高品質な理論/計算物理論文である。体系的な数値シミュレーションを通じて、固定無秩序誘起の新しい固-液転移を明らかにし、物理的図像が明確で分析が包括的である。主な不足は解析理論の支援と実験検証の欠落である。Soft MatterまたはPhysical Review E等の学術誌への掲載を推奨する。推奨指数：★★★★☆