[논문 리뷰] Transport of Charged Particles: Entropy Production and Maximum Dissipation Principle

Chia-Yu Hsieh, YunKyong Hyon외 3명|arXiv (Cornell University)|2014. 07. 31.

Spectroscopy and Quantum Chemical Studies참고 문헌 18인용 수 4

한 줄 요약

이 논문은 이온 종 간 상대 속도 차이로 인한 추가적인 에너지 손실을 포함하는 수정된 포아송-너른스트-플랭크(Poisson-Nernst-Planck, PNP) 시스템을 제안한다. 에너지적 변분 접근법을 사용하여 유도된 이 시스템은 원래 PNP와 동일한 평형 상태를 유지하지만, 추가적인 엔트로피 생성으로 인해 다를 바 있는 역학을 보이며, 갈레르킨 방법과 슈아우더 고정점 정리에 기반한 엄밀한 국소 존재 정리가 증명되었다. 유한요소 시뮬레이션을 통해 동일한 평형 상태를 보이지만 임시적 행동이 다름을 확인하였다.

ABSTRACT

In order to describe the dynamics of crowded ions (charged particles), we use an energetic variation approach to derive a modified Poisson-Nernst-Planck (PNP) system which includes an extra dissipation due to the effective velocity differences between ion species. Such a system is more complicated than the original PNP system but with the same equilibrium states. Using Schauder's fixed-point theorem, we develop a local existence theorem of classical solutions for the modified PNP system. Different dynamics (but same equilibrium states) between the original and modified PNP systems can be represented by numerical simulations using finite element method techniques.

연구 동기 및 목표

• 표준 PNP 이론이 배제 체적 효과와 높은 이온 농도로 인해 실패하는 혼잡하고 비이상적인 환경에서의 이온 이동을 모델링하기 위해.
• 열역학적으로 일관된 변분 원리에 기반해 이종 간 상대 운동으로 인한 추가적인 손실을 PNP 프레임워크에 통합하기 위해.
• 강화된 비선형성을 지닌 수정된 PNP 시스템에 대한 엄밀한 국소 존재 이론을 수립하기 위해.
• 수정된 PNP 시스템과 원래 PNP 시스템의 역학을 수치적으로 비교하여, 동일한 평형 상태임에도 불구하고 임시적 행동의 차이를 부각하기 위해.

제안 방법

• 보존력에 대해 최소 작용 원리(Least Action Principle, LAP)와 소산력에 대해 최대 소산 원리(Maximum Dissipation Principle, MDP)를 조합한 에너지적 변분 접근법을 사용한다.
• 양이온과 음이온 간 상대 속도의 제곱에 비례하는 항을 소산 함수에 추가하여 소산 함수를 수정한다: $ \frac{k_B T}{D_{n,p}} c_n c_p |\vec{u}_n^* - \vec{u}_p^*|^2 $, 이는 종 간 마찰을 나타낸다.
• 수정된 PNP 시스템을 이온 농도와 전기적 위치 에너지에 대한 비선형이고 탈기성 있는 편미분방정식(PDE)의 쌍으로 유도한다.
• 적절한 함수 공간 프레임워크 내에서 갈레르킨 방법과 슈아우더 고정점 정리를 사용하여 수정된 PNP 시스템의 고전적 해에 대한 국소 존재성을 증명한다.
• 안정성과 에너지 법칙 근사화를 보장하기 위해 암시적 시간 이산화를 사용하는 유한요소법(Finite Element Method, FEM) 스킴을 구현하여 시스템을 수치적으로 해석한다.
• 경계 조건과 초기 자료를 동일하게 설정하여 수정된 PNP 시스템과 원래 PNP 시스템의 수치적 해를 비교함으로써 추가 소산의 영향을 분리한다.

실험 결과

연구 질문

• RQ1 상대 속도에 의존하는 소산 항의 포함이 혼잡한 이온계에서의 이온 이동 역학에 어떤 영향을 미치는가?
• RQ2 추가적인 이종 간 마찰 항을 포함한 수정된 PNP 시스템이 원래 PNP 시스템과 동일한 평형 상태를 유지할 수 있는가?
• RQ3 수정된 PNP 시스템의 수학적 구조는 어떠한가? 비록 비선형성과 잠재적 탈기성으로 인해 복잡해졌지만, 국소 고전적 해가 존재하는가?
• RQ4 동일한 평형 상태와 경계 조건을 가진 수정된 PNP 시스템과 원래 PNP 시스템의 임시적 역학은 어떻게 다른가?
• RQ5 추가 소산 항은 이온 이동 시뮬레이션에서 엔트로피 생성률과 에너지 소산 프로파일에 어느 정도의 영향을 미치는가?

주요 결과

• 수정된 PNP 시스템은 양이온과 음이온 간 상대 속도에 의존하는 새로운 소산 항을 도입하여 원래 PNP 모델에 존재하지 않는 종 간 마찰 효과를 캡처한다.
• 추가적인 복잡성과 비선형성에도 불구하고, 수정된 PNP 시스템은 원래 PNP 시스템과 동일한 평형 상태를 유지한다. 이는 양자 모두 동일한 자유 에너지 함수에 의해 결정되기 때문이다.
• 갈레르킨 방법과 슈아우더 고정점 정리를 사용하여 수정된 PNP 시스템의 고전적 해에 대한 국소 존재 정리를 엄밀히 증명하였다.
• 유한요소법을 사용한 수치 시뮬레이션을 통해 동일한 경계 조건 하에서 수정된 PNP 시스템과 원래 PNP 시스템이 동일한 평형 상태에 도달하는 것으로 확인되었다.
• 두 시스템은 임시적 역학에서 크게 다른 행동을 보이며, 특히 시간에 따라 변화하는 소산 프로파일에서 수정된 시스템이 추가 소산 항으로 인해 더 높은 엔트로피 생성률을 보인다.
• 그림 2에서 수정된 시스템의 소산률 $\triangle^*$ 과 원래 시스템의 소산률 $\triangle$ 을 비교한 결과, 추가 소산 항이 평형 상태에 도달하는 데 필요한 회복 경로를 변화시키지만 최종 상태에는 영향을 주지 않는 것으로 나타났다.