[논문 리뷰] Coupling Lattice Boltzmann with Atomistic Dynamics for the multiscale simulation of nano-biological flows

Maria Fyta, Simone Melchionna외 2명|ArXiv.org|2007. 10. 05.

Lattice Boltzmann Simulation Studies참고 문헌 24인용 수 6

한 줄 요약

이 논문은 메조스코프적 용매 유체역학을 위한 격자볼츠만(LB)과 장수 생고분자에 대한 제약 조건이 부여된 분자역학(MD)을 동시적으로 결합하는 다스케일 방법을 제시한다. 이는 나노생물학적 유량의 효율적 시뮬레이션을 가능하게 하며, 국소적이고 명시적인 상호작용을 통해 유체역학적 상관관계를 포착한다. 이는 전이 스케일 법칙을 정확히 재현하고, 계산 비용이 용매 부피와 용질 크기와 선형으로 증가하는 방식으로 협동 동역학을 드러낸다.

ABSTRACT

We describe a recent multiscale approach based on the concurrent coupling of constrained molecular dynamics for long biomolecules with a mesoscopic lattice Boltzmann treatment of solvent hydrodynamics. The multiscale approach is based on a simple scheme of exchange of space-time information between the atomistic and mesoscopic scales and is capable of describing self-consistent hydrodynamic effects on molecular motion at a computational cost which scales linearly with both solute size and solvent volume. For an application of our multiscale method, we consider the much studied problem of biopolymer translocation through nanopores: we find that the method reproduces with remarkable accuracy the statistical scaling behavior of the translocation process and provides valuable insight into the cooperative aspects of biopolymer and hydrodynamic motion.

연구 동기 및 목표

• 원자적 및 메조스코프적 동역학이 결합된 나노생물학적 시스템을 효율적으로 시뮬레이션하기 위한 동시 다스케일 방법을 개발하는 것.
• 비국소적 근사에 의존하지 않고 나노포를 통과하는 생고분자에서의 유체역학적 상호작용을 시뮬레이션하는 데 도전하는 것.
• 분자 크기와 용매 부피에 따라 계산 비용이 선형적으로 증가하도록 보장하여 장수 생고분자의 대규모 시뮬레이션을 가능하게 하는 것.
• 용매와 용질 간의 국소적이고 시공간적으로 연결된 상호작용을 통해 유체역학적 상관관계를 명시적으로 포함하는 프레임워크를 제공하는 것.
• 기존 통계적 스케일링 행동을 재현함으로써 표준 문제인 생고분자 전이 시뮬레이션을 통해 방법을 검증하는 것.

제안 방법

• 방법은 동시적이고 반복이 없는 방식으로 메조스코프적 유체역학을 위한 격자볼츠만(LB)과 장수 생고분자를 위한 제약 조건이 부여된 분자역학(MD)을 결합한다.
• 유체역학적 힘은 고분자 입자에 대한 LB 분포 함수의 국소적 보간을 통해 계산되며, 운동량 교환을 통해 힘이 LB 격자로 피드백된다.
• 시간에 따른 하향분할(subcycling)을 사용하여 MD 단계가 LB 단계보다 더 높은 빈도로 수행되며, 이는 시간적 일致성을 확보한다.
• 비국소적 유체역학 표현(예: 오세인 텐서)을 피하기 위해, 이론은 고체-고분자 상호작용을 이산 격자 동역학을 통해 직접 모델링한다.
• LB 해법기는 이산 격자를 기반으로 하며, 단일 입자 분포 함수를 발전시켜 연속체 근사에서 라우-나비에-스토크스 유사 행동을 회복한다.
• 이 방법은 LB 알고리즘의 본질적인 국소성과 효율성을 활용하여 고도로 병렬화 가능하도록 설계되어 있다.

실험 결과

연구 질문

• RQ1 LB와 MD의 동시적이고 반복이 없는 결합이 나노포를 통과하는 생고분자의 통계적 스케일링 행동을 정확히 재현할 수 있는가?
• RQ2 유체역학적 상호작용은 전이 과정 중 고분자와 용매의 협동 동역학에 어떻게 영향을 미치는가?
• RQ3 LB-MD 결합의 계산 스케일링은 고분자 길이와 용매 부피에 대해 어떻게 되는가?
• RQ4 비국소적 근사 없이도 협동적이고 유체역학적으로 향상된 전이 과정의 발생을 방법이 포착할 수 있는가?
• RQ5 이 방법은 원자적 척도와 메조스코프적 척도 간의 운동량과 에너지 교환에서 물리적 일致성을 어느 정도 유지하는가?

주요 결과

• 방법은 생고분자 전이의 알려진 통계적 스케일링 지수를 정확히 재현하여 물리적 정확성을 검증한다.
• 유체역학적 상호작용이 전이 속도를 크게 향상시킴을 보여주며, 이는 용매가 협동적 환경으로 작용함을 시사한다.
• 전이 좌표 $ S_E(t) $ 의 분포는 작은 음의 尾를 보이며, 이는 전기장 반대 방향으로 움직이는 입자들이 드물게 존재함을 나타내며, 복잡한 동역학을 드러낸다.
• 용매 부피를 일정하게 유지할 때 계산 비용은 고분자 입자 수에 따라 선형적으로 증가하며, 각 단계당 실행 시간은 각각 50, 100, 400개의 입자를 가진 경우 0.433, 0.489, 0.882 초이다.
• 500개의 입자를 가진 전이 사건이 30,000개의 LB 단계 동안 지속될 경우 총 시뮬레이션 비용은 약 6시간이며, 관측된 런타임(~7시간)과 매우 유사하다.
• 이 방법은 중간 정도의 계산 비용으로 장수 생고분자(N > 1000)의 시뮬레이션을 가능하게 하여, 향후 화학적으로 특정된 고분자를 포함한 시뮬레이션의 길을 열어 놓는다.