[논문 리뷰] The Spectrum of the Force-Based Quasicontinuum Operator for a Homogeneous Periodic Chain

Matthew Dobson, Christoph Ortner외 1명|Warwick Research Archive Portal (University of Warwick)|2010. 04. 20.

Numerical methods in engineering참고 문헌 16인용 수 7

한 줄 요약

이 논문은 동일한 주기적 사슬에서 두 번째 이웃 원자 상호작용의 경우 선형화된 力기반 준연속체(QCF) 연산자의 스펙트럼이 양의 값이며, 에너지 기반 준국소(QNL) 연산자와 정확히 일치함을 엄밀하게 증명한다. 또한 시스템 크기와 원자적 영역 크기에 관계없이 고유벡터 행렬의 조건수에 대해 균일한 상한을 확립한다. 이러한 결과는 이전의 수치적 추측을 검증하며, GMRES 솔버의 엄밀한 수렴 추정과 QCF 방법에 대한 새로운 안정성 한계를 가능하게 한다.

ABSTRACT

We show under general conditions that the linearized force-based quasicontinuum (QCF) operator has a positive spectrum, which is identical to the spectrum of the quasinonlocal quasicontinuum (QNL) operator in the case of second-neighbour interactions. Moreover, we establish a bound on the condition number of a matrix of eigenvectors that is uniform in the number of atoms and the size of the atomistic region. These results establish the validity of and improve upon recent conjectures ([arXiv:0907.3861, Conjecture 2] and [arXiv:0910.2013, Conjecture 8]) which were based on numerical experiments. As immediate consequences of our results we obtain rigorous estimates for convergence rates of (preconditioned) GMRES algorithms, as well as a new stability estimate for the QCF method.

연구 동기 및 목표

• 선형화된 力기반 준연속체(QCF) 연산자의 스펙트럼과 고유벡터 조건수에 대한 이전의 수치적 추측을 엄밀히 검증하기 위해.
• 시스템 크기와 원자적 영역 크기에 관계없이 QCF 방법에 대한 균일한 안정성 한계를 설정하기 위해.
• QCF 평형 방정식에 적용된 반복적 해법(GMRES 등)의 수렴에 대한 이론적 기반을 제공하기 위해.
• 다음 이웃 원자 상호작용에서 일반적인 유한 범위 상호작용으로의 스펙트럼 분석을 확장하기 위해.
• 시스템 크기 $N$과 원자적 영역 $τ$에 관계없이 $\upsilon^{2,2}$ 노름에서 QCF 방법에 대한 새로운 안정성 추정을 유도하기 위해.

제안 방법

• 주기적 사슬에서 총 $N$개의 원자를 가진 균일한 탈형 상태를 기준으로 QCF 연산자를 선형화한다.
• 스펙트럼 이론과 $\times^N$ 내의 행렬 표현을 사용하여 QCF 연산자의 스펙트럼을 에너지 기반 준국소(QNL) 연산자와 비교한다.
• 고유벡터 행렬을 명시적으로 구성하고, 스펙트럼 분해와 연산자 노름 추정을 통해 $N$과 원자적 영역 크기에 대해 균일하게 조건수를 상한화한다.
• 유한 범위 상호작용의 경우, 잠재 에너지의 이阶도수에 대한 추정을 통해 QCF 연산자의 스펙트럼이 원자적 연산자와 연속체 연산자의 스펙트럼 사이에 있음을 보인다.
• 균일한 조건수 상한과 조건화된 연산자의 스펙트럼 성질을 이용하여, 조건화된 GMRES 수렴 속도를 유도한다.
• 평균이 0인 이산 소볼레프 공간에서 QCF 연산자의 역행렬을 분석함으로써 $\mathcal{U}^{2,2}$ 노름에서 새로운 안정성 추정을 도출한다.

실험 결과

연구 질문

• RQ1 두 번째 이웃 원자 상호작용의 경우 선형화된 QCF 연산자의 스펙트럼이 QNL 연산자와 일치하는가?
• RQ2 QCF 연산자의 고유벡터 행렬 조건수는 원자 수와 원자적 영역 크기에 대해 균일하게 유계인가?
• RQ3 QCF 연산자의 스펙트럼 성질은 다음 이웃 원자 상호작용에서 일반적인 유한 범위 상호작용으로 확장될 수 있는가?
• RQ4 QCF 평형 방정식을 풀기 위한 조건화된 GMRES 방법의 수렴 속도는 무엇인가?
• RQ5 QCF 연산자는 시스템 크기와 원자적 영역 크기에 관계없이 이산 $\upsilon^{2,2}$ 노름에서 안정적인가?

주요 결과

• 두 번째 이웃 원자 상호작용의 경우 선형화된 QCF 연산자의 스펙트럼은 양이며, QNL 연산자와 정확히 일치한다.
• QCF 연산자의 고유벡터 행렬의 조건수는 $N$과 원자적 영역 크기에 대해 균일하게 유계이며, 이는 이전의 추측을 향상시킨다.
• 유한 범위 상호작용의 경우 QCF 연산자의 스펙트럼은 원자적 연산자와 연속체 연산자의 스펙트럼 사이에 있으며, 이에 대한 상한은 $N$과 $\tau$에 관계없이 일정하다.
• QCF 시스템을 풀기 위한 조건화된 GMRES 방법은 잔차 상한 $\|L_{1}^{-1}r^{(m)}\| \leq C N^3 q^m \|L_{1}^{-1}r^{(0)}\|$ 및 $\|L_{1}^{-1/2}r^{(m)}\| \leq C N^2 q^m \|L_{1}^{-1/2}r^{(0)}\|$을 가진다. 여기서 $C>0$, $q\in(0,1)$은 $N$과 $\mathcal{A}$에 관계없이 일정하다.
• QCF 연산자는 $\mathcal{U}^{2,2}$ 노름에서 안정적이다: $\|(L^{\rm qcf}_0)^{-1}\|_{L(\mathcal{U}^{0,2},\mathcal{U}^{2,2})}$는 $N$이나 $\mathcal{A}$에 관계없이 잠재 에너지 매개수에만 의존하는 상수로 유계이다.