[논문 리뷰] Multiscale modeling of a rectifying bipolar nanopore: explicit-water versus implicit-water simulations

Zoltán Ható, Mónika Valiskó외 3명|Repository of the Academy's Library (Library of the Hungarian Academy of Sciences)|2017. 11. 30.

Nanopore and Nanochannel Transport Studies참고 문헌 116인용 수 13

한 줄 요약

이 연구는 정류 기능을 가진 이중극성 나노포어를 대상으로 명시적 물 분자 동역학(MD)과 암시적 물 분자 Nernst-Planck/국소 평형 몽테카를로(NP+LEMC) 시뮬레이션을 비교한다. 분자의 스케일에서의 상당한 차이에도 불구하고, 두 모델은 정량적으로 유사한 전류-전압 응답과 축방향 이온 농도 프로파일을 생성하며, 이는 분자의 현실성을 기여하지 않더라도 감소된 모델이 핵심적인 장치 수준의 물리적 거동을 잘 포착할 수 있음을 보여준다.

ABSTRACT

In a multiscale modeling approach, we present computer simulation results for a rectifying bipolar nanopore on two modeling levels. In an all-atom model, we use explicit water to simulate ion transport directly with the molecular dynamics technique. In a reduced model, we use implicit water and apply the Local Equilibrium Monte Carlo method together with the Nernst-Planck transport equation. This hybrid method makes the fast calculation of ion transport possible at the price of lost details. We show that the implicit-water model is an appropriate representation of the explicit-water model when we look at the system at the device (i.e., input vs.\ output) level. The two models produce qualitatively similar behavior of the electrical current for different voltages and model parameters. Looking at details of concentration and potential profiles, we find profound differences between the two models. These differences, however, do not influence the basic behavior of the model as a device because they do not influence the $z$-dependence of the concentration profiles which are the main determinants of current. These results then address an old paradox: how do reduced models, whose assumptions should break down in a nanoscale device, predict experimental data? Our simulations show that reduced models can still capture the overall device physics correctly, even though they get some important aspects of the molecular-scale physics quite wrong; reduced models work because they include the physics that is necessary from the point of view of device function. Therefore, reduced models can suffice for general device understanding and device design, but more detailed models might be needed for molecular level understanding.

연구 동기 및 목표

• 정류 기능을 가진 이중극성 나노포어를 통과하는 이온 이동을 모의하는 데 있어 암시적 물 원자 모델과 모든 원자 수준의 명시적 물 분자 동역학(MD) 모델 간의 비교를 수행하고자 한다.
• 단순화된 모델이 분자의 가정에 결함이 있음에도 불구하고 실험적 전류-전압 특성을 잘 재현하는 데 오랫동안 남아있는 역설을 해결하고자 한다.
• 분자의 스케일 물리적 정확성이 떨어지더라도 장치 수준의 예측에 대해 감소된 모델이 어떤 조건에서 유효한지를 규명하고자 한다.
• 나노스케일 이온 이동 장치의 다스케일 모델링에서 계산 효율성과 분자의 세부 정보 사이의 상충 관계를 평가하고자 한다.

제안 방법

• GROMACS를 사용한 명시적 물 분자를 포함한 모든 원자 수준의 분자 동역학(MD) 시뮬레이션으로, CHARMM27 힘장과 주기적 경계 조건을 적용하였다.
• 암시적 물 모델과 Nernst-Planck(NP) 운반 방정식을 국소 평형 몽테카를로(LEMC)와 결합하여 이온 확산 및 전기화학적 포텐셜을 계산하는 감소된 모델.
• NP+LEMC 모델은 이온을 유전율(ε = 78.5)의 연속체 내에 있는 점입자로 간주하며, 확산 계수를 기초 상태의 값에서 유도하였다.
• MD에서는 이온 수를 세어 전류를 계산하고, NP+LEMC에서는 확산 유량 통합을 통해 전류를 계산하며, 축방향 농도 및 전위 프로파일을 두 모델 간 비교하였다.
• 경계 조건은 다를 수 있다: MD는 주기적 경계 조건과 외부 전기장을 사용하는 반면, NP+LEMC는 실린더 형태의 구획에 고정된 기초 농도 및 전위를 사용한다.
• 전위 프로파일은 두 가지 방법으로 검증되었으며, 시험 전하 샘플링과 포아송 방정식 해법을 사용하였고, NP+LEMC에서는 뛰어난 일치를, MD에서는 합리적인 일치를 보였다.

실험 결과

연구 질문

• RQ1 정류 기능을 가진 이중극성 나노포어의 전류-전압 응답을 예측하는 데 있어 명시적 물 분자 MD와 암시적 물 분자 NP+LEMC 시뮬레이션 간의 비교는 어떻게 이루어지나?
• RQ2 두 모델 간의 이온 및 물의 구조에 대한 분자의 스케일에서의 차이가 장치 수준의 전기적 특성에 어느 정도 영향을 미치는가?
• RQ3 점입자, 암시적 용매 등의 단순화된 가정을 가진 감소된 모델이 알려진 물리적 근사에도 불구하고 실험적 정류 행동을 여전히 재현하는 이유는 무엇인가?
• RQ4 암시적 물 모델은 명시적 물 모델 대비 이온 농도 및 전위 프로파일의 어떤 요소를 유지하거나 왜곡하는가?
• RQ5 암시적 물 모델이 세부적인 분자 현상을 포착하지 못하더라도, 장치 설계에 신뢰할 수 있는가?

주요 결과

• 두 모델은 다양한 전압, 표면 전하 및 이온 농도 조건에서 정량적으로 동일한 전류-전압 응답을 보였다.
• 명시적 물과 암시적 물 모델 간의 축방향(z-의존성) 농도 프로파일은 국소 구조 및 용매화의 차이에도 불구하고 놀랍도록 유사하였다.
• 암시적 물 모델은 전류의 핵심 결정 요소인 z-의존성 이온 농도 프로파일을 정확히 포착하고 있어, 이 모델의 예측 성공을 설명한다.
• 특히 좁은 포어 영역에서 국소 전위 및 이온 분포 프로파일에 상당한 차이가 있으며, 이는 감소된 모델이 분자의 스케일 물리적 거동을 잘못 묘사하고 있음을 시사한다.
• NP+LEMC 모델의 계산 비용은 매우 낮아 몇 시간 내에 수행되며, MD의 수일에서 수주에 이르는 계산 대비 효율적인 파arameter 스캔이 가능하다.
• NP+LEMC 모델은 두 가지 전위 계산 방법(시험 전하 샘플링 및 포아송 방정식 해법) 간에 뛰어난 일치를 보였으며, 이는 수치적 안정성을 검증한다.