[논문 리뷰] Analysis of an Energy-based Atomistic/Continuum Coupling Approximation of a Vacancy in the 2D Triangular Lattice

Christoph Ortner, Alexander V. Shapeev|arXiv (Cornell University)|2011. 04. 02.

Graphene research and applications참고 문헌 23인용 수 27

한 줄 요약

이 논문은 유한 범위의 두체 잠재력으로 이루어진 2D 삼각 격자에서 흐린 결함을 모델링하기 위한 에너지 기반 원자역학/연속체 결합(a/c) 방법에 대한 엄밀한 사전 오차 분석을 제시한다. 이 방법은 결함 근처에서는 원자역학 해상도를 유지하고 원거리 영역에서는 연속체 모델을 사용하며, 새로운 일致성 및 안정성 분석을 통해 메쉬 크기와 해의 부드러움에 따라 H¹-노름과 에너지에서 수렴 속도를 도출한다. 실용적인 메쉬 및 원자역학 영역 설계를 위한 힌트도 제시한다.

ABSTRACT

We present a comprehensive a priori error analysis of a practical energy based atomistic/continuum coupling method (Shapeev, arXiv:1010.0512) in two dimensions, for finite-range pair-potential interactions, in the presence of vacancy defects. The majority of the work is devoted to the analysis of consistency and stability of the method. These yield a priori error estimates in the H1-norm and the energy, which depend on the mesh size and the "smoothness" of the atomistic solution in the continuum region. Based on these error estimates, we present heuristics for an optimal choice of the atomistic region and the finite element mesh, which yields convergence rates in terms of the number of degrees of freedom. The analytical predictions are supported by extensive numerical tests.

연구 동기 및 목표

• 2D 격자에서 흐린 결함가 있는 실용적인 에너지 기반 원자역학/연속체 결합 방법에 대한 종합적인 사전 오차 분석을 제공하는 것.
• 원자역학-연속체 경계에서 발생하는 고스트 힘 문제를 엄밀한 일치성 및 안정성 분석을 통해 다루는 것.
• 메쉬 크기와 연속체 영역 내 원자역학 해의 부드러움에 따라 H¹-노름과 에너지에서 오차 추정치를 도출하는 것.
• 최적의 수렴 속도를 달성하기 위해 연속체 영역의 유한요소 메쉬와 원자역학 영역 크기의 최적 선택을 위한 힌트를 개발하는 것.

제안 방법

• 방법은 Shapeev(2010)의 ECC(에너지 기반 결합) 형식을 2D 삼각 격자에 적용하여 주기적 경계 조건과 이체 잠재력 상호작용을 고려한다.
• 연속체 영역 내 원자역학 해의 정규성을 측정하기 위해 새로운 이산 부드러움 척도를 도입하며, 이는 오차 추정에 핵심적이다.
• 잠재력 도함수의 국소 리프시츠 추정과 결합 밀도 레미마를 사용하여 이산 W⁻¹ᵖ-노름에서 일치 오차를 분석한다.
• 결함이 있는 격자의 안정성 분석을 균일 격자로의 환원을 가능하게 하는 '흐린 결함 안정성 인덱스'를 통해 안정성을 확립한다.
• 일치성 및 안정성 추정치를 통합하여 H¹-노름과 에너지에서 사전 오차 한계를 도출한다.
• 오차 한계를 자유도 수와 연관지어 유한요소 메쉬 및 원자역학 영역 설계의 최적화를 안내하는 수치적 힌트를 유도한다.

실험 결과

연구 질문

• RQ1 2D 에너지 기반 원자역학/연속체 결합 방법에서 흐린 결함에 대해 일치성 및 안정성을 어떻게 엄밀하게 분석할 수 있는가?
• RQ2 사전 오차 한계가 메쉬 크기와 연속체 영역 내 원자역학 해의 부드러움에 따라 어떻게 의존하는가?
• RQ3 흐린 결함 안정성 인덱스는 어떻게 정의되고, 이를 통해 결함이 있는 격자의 안정성 분석을 균일 격자로 환원할 수 있는가?
• RQ4 최적의 수렴 속도를 달성하기 위해 연속체 영역의 원자역학 영역 크기와 유한요소 메쉬를 최적 선택하기 위한 힌트는 무엇인가?
• RQ5 분석적 오차 추정치가 실질적 시뮬레이션에서의 수치 수렴 행동과 얼마나 잘 일치하는가?

주요 결과

• 논문은 메쉬 크기와 연속체 영역 내 원자역학 해의 부드러움에 따라 H¹-노름과 에너지에서 사전 오차 추정치를 확립한다.
• 새로운 이산 부드러움 척도를 도입하여 이전 연구에 비해 더 날카운 일치 오차 추정치를 도출할 수 있으며, 더 적은 기술적 가정을 필요로 한다.
• 흐린 결함 안정성 인덱스를 통해 결함이 있는 격자의 안정성을 균일 격자의 안정성 분석으로 환원할 수 있다.
• 국소 리프시츠 추정과 결합 밀도 레미마를 사용하여 이산 W⁻¹ᵖ-노름에서 일치 오차를 유계로 제시한다.
• 일치성 및 안정성 분석을 통합하여 자유도 수에 따라 H¹-노름과 에너지에서 수렴 속도를 도출하며, 최적의 메쉬 및 원자역학 영역 설계를 위한 힌트를 제공한다.
• 광범위한 수치 실험을 통해 분석 예측이 잘 일치하며, 이론적 수렴 속도와 양호한 일치를 보인다.