[논문 리뷰] Coarse-graining molecular dynamics models using an extended Galerkin projection

Xiantao Li|arXiv (Cornell University)|2012. 10. 16.

Advanced Mathematical Modeling in Engineering참고 문헌 44인용 수 4

한 줄 요약

이 논문은 결정성 고체의 거시적 분석을 위한 확장된 갈레르킨 투영 방법을 제안하며, 결함 근처에서 거시 척도 부분공간을 풍부화시켜 원자적/연속체 경계에서 격자파 반사를 최소화한다. 기존 갈레르킨 방법과 달리 일반화 랑주빈 방정식에서 기억 함수를 사전 계산할 필요 없이도 변분 구조와 기계적 평형 안정성을 유지한다.

ABSTRACT

We present a new framework for coarse-graining molecular dynamics models for crystalline solids. The reduction method is based on a Galerkin projection to a subspace, whose dimension is much smaller than that of the full atomistic model. The subspace is expanded by adding more coarse-grain variables near the interface between lattice defects and the surrounding regions. This effectively minimizes reflection of phonons at the interface. In this approach, there is no need to pre-compute the memory function in the generalized Langevin equations, a typical model of interface conditions. Moreover, the variational formulation preserves the stability of mechanical equilibria.

연구 동기 및 목표

• 거시적 분석 분자의 동역학에서 원자적/연속체 경계에서 인위적인 격자파 반사를 해결한다.
• 일반화 랑주빈 방정식에서 기억 커널을 사전 계산하는 계산 부담을 완화한다.
• 거시적 분석 모델에서 변분 구조와 기계적 평형 안정성을 유지한다.
• 결함 이동 및 충격파와 같은 비평형 동역학을 정확하게 시뮬레이션할 수 있도록 한다.
• 노이즈를 통합함으로써 유한 온도 시스템에 적응 가능한 프레임워크를 제공한다.

제안 방법

• 기저 함수로 생성된 차원이 전체 원자적 시스템에 비해 훨씬 작은 축소된 부분공간에 갈레르킨 투영을 적용한다.
• 격자 결함과 경계 근처에서 추가 자유도를 도입하여 거시적 분석 부분공간을 확장함으로써 격자파 반사를 억제한다.
• 변분 공식을 사용하여 기계적 평형의 안정성을 유지하고 기저 해밀토니안 동역학과의 일致성을 확보한다.
• 중점 적분과 유한요소 이산화를 사용하여 거시 척도 방정식을 근사화함으로써 카우치-본 구성관계를 가진 탄성역학 모델을 도출한다.
• 경계에서 반복 노드를 도입하여 부분공간을 풍부화시킴으로써 일반화 랑주빈 방정식에 의존하지 않고 전이 영역을 효과적으로 모델링한다.
• 기존 갈레르킨 투영을 결함 근처의 국소 정밀화를 통합한 확장된 버전으로 대체함으로써 경계를 넘는 파동 전파를 향상시킨다.

실험 결과

연구 질문

• RQ1 거시적 분석 분자의 동역학 모델에서 원자적/연속체 경계에서 격자파 반사를 어떻게 최소화할 수 있는가?
• RQ2 일반화 랑주빈 방정식에서 기억 함수를 사전 계산할 필요 없이 정확성을 유지할 수 있는가?
• RQ3 원래 분자의 동역학 모델의 변분 구조는 거시적 분석 공식화에서 어떻게 유지될 수 있는가?
• RQ4 결함 근처의 부분공간 풍부화가 거시적 분석 시뮬레이션의 정확성과 안정성에 어떤 영향을 미치는가?
• RQ5 유한 온도 시스템으로의 방법 확장은 확률 소산 정리 조건을 만족하면서 가능할 수 있는가?

주요 결과

• 확장된 갈레르킨 방법은 이심력 이중체 시뮬레이션을 통해 기존 갈레르킨 방법 대비 경계에서 격자파 반사를 크게 감소시킴을 입증하였다.
• 기존 갈레르킨 방법은 특히 원자적 영역에서 경계 반사로 인해 큰 수치 오차를 유발한다.
• Krylov 부분공간에서 L=8을 사용한 확장된 갈레르킨 방법은 L=4보다 훨씬 낮은 오차를 기록하여 부분공간 풍부화가 증가함에 따라 정확도가 향상됨을 보였다.
• 일반화 랑주빈 방정식에서 기억 커널을 사전 계산할 필요 없이도 방법을 단순화하고 계산 비용을 감소시켰다.
• 확장된 갈레르킨 투영을 통해 유도된 거시적 분석 모델은 전체 분자의 동역학과 일치하는 결과를 도출하였으며, 특히 결함 핵 구조와 파동 전파를 잘 포착하였다.
• 변분 공식은 기계적 평형의 안정성을 보장하며, 이는 많은 기존 결합 방법에서 보장되지 않는 특성이다.