[논문 리뷰] Diffusion on a lattice: transition rates, interactions and memory effects

M. A. Di Muro, M. Hoyuelos|arXiv (Cornell University)|2022. 04. 03.

Advanced Thermodynamics and Statistical Mechanics참고 문헌 35인용 수 1

한 줄 요약

이 논문은 초과 화학적 포텐셜에서 유도된 전이 비율을 통해 기억 효과를 고려하여 격자 확산에서 상관 계수 f에 대한 일반적인 해석적 표현을 제안한다. 부드러운 코어 및 확장된 하드코어 제약 조건을 통해 입자 상호작용을 모델링함으로써 농도에 따라 변하는 f를 유도하여 평균장 트레이서 확산계수 예측을 향상시켰으며, 다양한 입자 배제 크기와 밀도에서 몬테카를로 시뮬레이션과 강한 일치를 보였다.

ABSTRACT

We analyze diffusion of particles on a two dimensional square lattice. Each lattice site contains an arbitrary number of particles. Interactions affect particles only in the same site, and are macroscopically represented by the excess chemical potential. In a recent work, a general expression for transition rates between neighboring cells as functions of the excess chemical potential was derived. With transition rates, the mean field tracer diffusivity, $D^ ext{MF}$, is immediately obtained. The tracer diffusivity, $D = D^ ext{MF} f$, contains the correlation factor $f$, representing memory effects. An analysis of the joint probability of having given numbers of particles at different sites when a force is applied to a tagged particle allows an approximate expression for $f$ to be derived. The expression is applied to soft core interaction (different values for the maximum number of particles in a site are considered) and extended hard core.

연구 동기 및 목표

• 메모리 효과를 고려한 2차원 격자에서의 트레이서 확산에 대한 일반적인 해석적 프레임워크를 개발하는 것.
• 입자 상호작용으로 인한 공간적 상관관계를 포함시켜 평균장 확산계수 이론을 확장하는 것.
• 초과 화학적 포텐셜과 열역학적 요인을 사용하여 부드러운 코어 및 확장된 하드코어 상호작용을 매크로스코픽적으로 모델링하는 것.
• 다양한 입자 배제 영역과 밀도에서 이론적 상관 계수를 몬테카를로 시뮬레이션과 비교 검증하는 것.

제안 방법

• 초과 화학적 포텐셜로 파arameter화된 상세 평형 조건과 유효한 위도 삽입 공식을 사용하여 격자 위치 간 전이 비율을 유도한다.
• 열역학적 요인 Γ를 입자 농도와 화학적 포텐셜에서 유도한 평균장 트레이서 확산계수 DMF = D0 / Γ 를 도입한다.
• 입자 수의 공동 확률 분포를 사용하여 소외 외부 힘 하에서 전이 비율 비대칭성에 기반한 상관 계수 f를 전이 비율의 함수로 기술한다.
• Ω 입자당 격자 위치를 갖는 부드러운 코어 및 b 기반 배제 영역을 갖는 확장된 하드코어 상호작용에 이 방법을 적용하고, 유효한 포장 분율 ξ를 정의한다.
• 이차 근사 f = 1 − c₁ξ + c₂ξ² 를 사용하고, b 값에 따라 다른 몬테카를로 데이터에 맞추어 c₁, c₂ 계수를 피팅한다.
• 전체 트레이서 확산계수를 D/D₀ = f / Γ = (1 − c₁ξ + c₂ξ²)(1 − ξ/ξ_max) 로 유도하며, 여기서 ξ_max = ρ_max b / 4 이다.

실험 결과

연구 질문

• RQ1 격자 확산에서의 기억 효과는 어떻게 이론적 표현의 상관 계수 f에 체계적으로 통합될 수 있는가?
• RQ2 2차원 정사각형 격자에서 입자 농도와 상호작용 유형(부드러운 코어 대 비해 확장된 하드코어)에 따라 상관 계수 f는 어떻게 달라지는가?
• RQ3 초기 열역학적 양만을 사용하여 다양한 상호작용 모델에 적용 가능한 일반적인 f 표현식을 유도할 수 있는가?
• RQ4 조정 가능한 매개변수를 가진 이론적 f는 다양한 입자 배제 영역(b)과 농도에서 시뮬레이션 결과와 얼마나 잘 일치하는가?

주요 결과

• 조정된 c₁ 및 c₂를 사용하여 몬테카를로 데이터에 맞춘 이론적 표현 D/D₀ = (1 − c₁ξ + c₂ξ²)(1 − ξ/ξ_max) 는 다양한 b 값(배제 영역 크기)에서 트레이서 확산계수를 정확하게 예측한다.
• 부드러운 코어 상호작용의 경우 Ω이 증가할수록 상관 계수 f는 1에 수렴하며, 이는 입자 이동 시 위치에 대한 영향이 최소화되어 기억 효과가 감소함을 나타낸다.
• 확장된 하드코어 상호작용의 경우 최대 포장 분율 ξ_max 는 b 증가에 따라 감소하며, b=1일 땐 1에서 b=37일 땐 1/12로 감소하여 더 큰 배제 영역을 반영한다.
• 상관 계수 f는 포장 분율 ξ 증가에 따라 감소하며, 이는 고밀도에서 더 강한 기억 효과를 나타내며, 특히 작은 b(예: b=1,5)에서 두드러진다.
• 몬테카를로 시뮬레이션은 이론 모델과 뛰어난 일치를 보이며, 모든 테스트된 상호작용 유형과 밀도에서 접근 방식의 타당성을 검증한다.
• 이 방법은 f의 비단조화적 행동을 성공적으로 포착하며, 예를 들어 b=21일 때 c₂=-0.41과 같이 음수인 c₂ 값을 포함함으로써 고밀도에서 f(ξ)의 곡률 반전을 나타낸다.