Interface crossing behavior of prolate microswimmers: thermo and hydrodynamics

기술 요약: 편평한 미세 유영자(Prolate Microswimmers)의 계면 교차 거동

문제 정의
수성 환경 내 미생물의 운동성은 주변 환경으로부터 오는 피드백 메커니즘, 특히 액체-액체 계면에 의해 크게 영향을 받는다. 미세 유영자와 계면 사이의 유체 역학적 상호작용은 잘 문서화되어 있으나, 깨끗한 액체-액체 계면을 가로지르는 활성 입자의 역학에서 열역학적 트래핑 힘(특히 피커링 효과)의 역할은 상대적으로 덜 탐구되었다. 최근 Cheon 등(2024)의 실험에서는 막대 모양의 Bacillus subtilis 박테리아가 속도와 초기 투영 각도에 따라 이방성-액정(isotropic-nematic) 액체-액체 계면을 통과하거나 그곳에 갇힐 수 있음을 관찰했다. 본 연구는 깨끗한 액체-액체 계면에서 편평한(prolate, 막대형) 미세 유영자의 트래핑 역학을 조사하고, 계면의 열역학적 힘과 활성 유체 역학적 힘 사이의 경쟁을 규명하는 것을 목표로 한다.

방법론
저자들은 편평한 미세 유영자와 깨끗한 액체-액체 계면 사이의 상호작용을 모델링하기 위해 대규모 3D 격자 볼츠만 시뮬레이션(Lattice Boltzmann simulations)을 채택하였다.

• 유영자 모델: 미세 유영자는 스토크스 흐름(Stokes flow) 하에서 구형 스커머(spheroidal squirmer) 모델을 사용하여 고체 편평 입자(종횡비 $AR = 3$)로 모델링되었다. 유영자는 Bacillus subtilis를 모사하기 위해 정해진 경계 슬립 속도를 갖는 "푸셔(pusher, $\beta = -1$)"로 정의되었다.
• 유체 계면: 이성분 유체 혼합물은 Cahn-Hilliard 방정식을 통해 상 진화(phase evolution)를 수행하고 나비에-스토크스 방정식(Navier-Stokes equation)을 통해 유체 속도를 계산하는 Ginzburg-Landau 자유 에너지 범함수를 사용하여 모델링되었다. 계면은 중립적 젖음성(접촉각 $\approx 90^\circ$)을 갖도록 처리되었다.
• 시뮬레이션 파라미터: 시뮬레이션은 낮은 레이놀즈 수($Re \sim O(10^{-2})$)에서 수행되었다. 모세관 수($Ca = \eta u_0 / \sigma_0$)는 계면 장력($\sigma_0$)을 변경하거나 유영 속도($u_0$)를 변경함으로써 변화되었다. 계면에 대한 유영자의 초기 투영 각도($\phi_0$) 또한 체계적으로 변화되었다.
• 분석: 연구에서는 입자 궤적, 방향 역학 및 속도 변동을 분석하였다. 계면의 열역학적 에너지(계면 면적 감소로부터 유도됨)와 계면을 가로지르는 데 필요한 활성 에너지 사이의 에너지 균형을 기반으로 한 이론적 프레임워크가 개발되었다.

주요 기여 및 결과

1. 교차 vs 트래핑의 상도표: 시뮬레이션은 모세관 수($Ca$)와 초기 정렬 각도($\phi_0$) 공간에서의 상도표를 구축하였다. 결과는 중립적 젖음성을 가진 유영자가 활성 유영 힘과 계면 장력 사이의 경쟁에 따라 계면을 통과하거나 갇힐 수 있음을 보여준다.

   • 트래핑(Trapping): 계면(열역학적) 힘이 지배할 때 발생하며, 일반적으로 낮은 $Ca$ 또는 얕은 초기 각도에서 나타난다. 유영자는 계면에 평행하게 정렬된다.
   • 교차(Crossing): 활성 힘이 지배할 때 발생하며(높은 $Ca$), 유영자가 계면을 횡단할 수 있게 한다.
   • 검증: 수치적 상경계는 두 연구에서 제안된 트래핑 힘의 물리적 메커니즘이 다름에도 불구하고, Bacillus subtilis가 이방성-액정 계면을 통과하는 Cheon 등(2024)의 실험적 관찰 결과와 놀라울 정도로 잘 일치한다.

2. 열역학적 트래핑 메커니즘: 본 연구는 트래핑이 피커링 효과에 의해 구동된다고 제안한다. 편평한 유영자의 경우, 계면에서 제거되는 데 드는 에너지 비용은 정렬 각도 $\phi$에 따라 변하는 접촉 면적에 따라 달라진다. 저자들은 유영자에 의한 활성 작업량과 열역학적 트래핑 에너지($\Delta E_{th} \propto \sigma_0 \cdot \text{Area}(\phi)$)를 균형 맞춤으로써 임계 모세관 수($Ca^*$)를 도출하였다. 이 해석적 예측($Ca^*$)은 시뮬레이션에서 관찰된 트래핑 상태와 교차 상태 사이의 경계를 성공적으로 정의한다.

3. 재정렬 역학(Reorientation Dynamics): 계면에 평행하게 정렬되는 유영자의 재정렬은 유체 역학적 성분과 열역학적 성분을 모두 포함하는 결합된 토크에 의해 구동된다는 점을 밝혀냈다.

   • 유체 역학적 토크: 푸셔 유영자의 힘 쌍극자($B_2$)에서 발생하며, 모세관 수에 비례한다.
   • 열역학적 토크: 계면 자유 에너지의 기울기($\tau_{thermo} \sim -\partial \Delta E_p / \partial \phi$)에서 발생하며, 정렬 각도에 의존한다.
   • 관찰: 얕은 각도로 초기 정렬된 유영자는 (원거리 유체 역학에 의해) 계면 법선 방향으로 초기 재정렬을 보인 후, 계면에 갇혀 평행하게 재정렬되는 양상을 보일 수 있다. 재정렬 속도는 $Ca$에 선형적으로 증가한다.

4. 퍼텐셜 에너지 지형: 속도 변동을 통합하고 자유 에너지 차이를 계산함으로써, 저자들은 트래핑 퍼텐셜을 매핑하였다. 퍼텐셜은 계면 근처에서 포물선 형태를 띠며, 고정된 계면 장력에 대해 그 깊이가 비교적 일정하게 유지됨을 확인하였는데, 이는 낮은 $Ca$에서 열역학적 트래핑의 지배력을 확인시켜 준다.

의의 및 주장
본 논문은 막대 모양 미세 유영자의 계면 교차 거동에 대한 통합적인 이론적 및 수치적 설명을 제공한다고 주장한다. 주요 의의는 다음과 같다:

• 메커니즘 식별: 중립적 젖음성을 가진 편평한 유영자의 트래핑이 (계면 변형에 의한 스프링 같은 힘에 기인한다고 해석했던 이전의 실험적 해석과 달리) 피커링 효과(계면 에너지의 열역학적 감소)에 의해 지배됨을 입증하였다.
• 예측 능력: 에너지 균형 방정식을 기반으로 도출된 임서 임계 모세관 수는 트래핑 상태와 교차 상태 사이의 전이를 정확하게 예측하며, 저자의 시뮬레이션 및 독립적인 박테리아 실험과 매우 잘 일치한다.
• 이중 토크 기원: 계면에서의 미세 유영자 재정렬이 (구형 스커머에서 흔히 가정되는 것처럼) 단지 유체 역학적 현상만이 아니라, 비구형 입자의 경우 상당한 열역학적 성분을 포함한다는 것을 확립하였다.

저자들은 자신들의 모델이 트래핑 힘의 구체적인 물리적 기원에 있어서 Cheon 등(2024)의 실험과는 다르지만, 액체-액체 계면에서의 미세 유영자의 거시적 역학 및 상태 전이를 성공적으로 포착한다고 결론짓는다.