[논문 리뷰] A second-order numerical method for Landau-Lifshitz-Gilbert equation with large damping parameters

Yongyong Cai, Jingrun Chen외 2명|arXiv (Cornell University)|2021. 05. 08.

Magnetic properties of thin films참고 문헌 28인용 수 11

한 줄 요약

이 논문은 큰 감쇠 계수를 가진 랑두-리프시츠-기르베르트(LLG) 방정식을 위한 이阶 정확도를 가지며 조건부 안정인 수치적 방법을 제안한다. 감쇠 항을 조화 사상 흐름으로 재구성하고, 비선형 항을 완전히 명시적으로 처리한 두 번째 차수의 BDF 시간 이산화를 사용함으로써, 각 시간 단계에서 일정 계수를 가진 대칭적이고 양의 정부호 선형 시스템만을 풀게 된다—이를 통해 효율적인 FFT 기반 솔버를 적용할 수 있다. 핵심 기여는 큰 감쇠 조건에서도 안정성과 물리적 정확성을 유지하는 계산적으로 효율적인 두 번째 차수 정확도를 가진 방법을 제공한다는 점이며, 1D 및 3D 시뮬레이션(예: 도메인 벽 동역학 포함)을 통해 검증되었다.

ABSTRACT

A second order accurate numerical scheme is proposed and implemented for the Landau-Lifshitz-Gilbert equation, which models magnetization dynamics in ferromagnetic materials, with large damping parameters. The main advantages of this method are associated with the following features: (1) It only solves linear systems of equations with constant coefficients where fast solvers are available, so that the numerical efficiency has been greatly improved, in comparison with the existing Gauss-Seidel project method. (2) The second-order accuracy in time is achieved, and it is unconditionally stable for large damping parameters. Moreover, both the second-order accuracy and the great efficiency improvement will be verified by several numerical examples in the 1D and 3D simulations. In the presence of large damping parameters, it is observed that this method is unconditionally stable and finds physically reasonable structures while many existing methods have failed. For the domain wall dynamics, the linear dependence of wall velocity with respect to the damping parameter and the external magnetic field will be obtained through the reported simulations.

연구 동기 및 목표

• 현대 자성 물질에서 흔한 큰 감쇠 계수를 가진 랑두-리프시츠-기르베르트(LLG) 방정식을 수치적으로 효율적이고 안정적으로 풀 수 있는 방법을 개발하는 것.
• 기존의 반음성적 방법이 각 시간 단계에서 비대칭 선형 시스템을 풀어야 하여 GMRES와 같은 반복 솔버에 의존함으로써 효율성이 제한되는 계산적 병목 현상을 해결하는 것.
• 1차 정확도를 가지는 가우스-세이델 투영 방법(GSPM)을 개선하여 동일한 계산 복잡도를 유지하면서도 시간에 대해 두 번째 차수 정확도를 확보하고 빠른 솔버를 가능하게 하는 것.
• 기존 방법이 실패하거나 불안정해지는 큰 감쇠 조건에서도 도메인 벽 이동 및 에너지 소산과 같은 물리적으로 의미 있는 역학을 정확히 포착할 수 있는 견고한 수치적 프레임워크를 제공하는 것.

제안 방법

• 감쇠 항을 조화 사상 흐름으로 표현함으로써, 확산 유사 행동을 더 체계적으로 다룰 수 있도록 LLG 방정식을 재구성한다.
• 시간 도함수에 대해 두 번째 차수의 후방 미분 공식(BDF2)을 적용하여 시간에 대해 두 번째 차수 정확도를 확보한다.
• 일정 계수를 가진 라플라스 연산자(확산) 항은 BDF2 스킴을 사용해 암시적으로 처리하여 대칭적이고 양의 정부호 선형 시스템을 형성한다. 이는 빠른 솔버에 적합하다.
• 모든 비선형 항—특히 기립모멘트 항과 조화 사상 흐름의 비선형 부분—은 두 번째 차수 보간 공식을 사용해 명시적으로 처리함으로써 안정성과 효율성을 유지한다.
• 각 시간 단계에서 유도된 선형 시스템은 대칭적이고 양의 정부호이므로 FFT 기반 빠른 솔버를 적용할 수 있으며, 이는 계산 효율성을 크게 향상시킨다.
• 조화 사상 흐름에 의한 안정화 소산 효과로 인해, 큰 감쇠 계수 조건에서도 이 방법은 조건부 안정성을 확보한다.

실험 결과

연구 질문

• RQ1 큰 감쇠 계수를 가진 LLG 방정식에 대해 두 번째 차수 정확도를 가지며 조건부 안정인 수치적 스킴을 개발할 수 있는가? 기존 방법은 종종 실패하거나 불안정해지는 상황이다.
• RQ2 벡터성과 비선형성을 가진 LLG 방정식임에도 불구하고, FFT와 같은 빠른 솔버의 사용을 가능하게 함으로써 이 방법이 높은 계산 효율성을 유지할 수 있는가?
• RQ3 도메인 벽 속도는 감쇠 계수 α와 외부 자기장 he에 어떻게 의존하는가? 그리고 새로운 스킴은 이 의존성을 정확히 포착할 수 있는가?
• RQ4 특히 큰 감쇠 조건에서 3D 시뮬레이션을 수행할 때, 제안된 방법은 GSPM 및 반음성적 투영 방법(SIPM)보다 더 효율적이고 정확한가?

주요 결과

• 제안된 방법은 1D 및 3D 시뮬레이션 모두에서 두 번째 차수 시간 정확도를 달성하였으며, 테스트된 모든 오차 노름에 대해 수렴 차수 약 2.0을 보였다.
• GSPM 및 SIPM이 동일한 조건에서 불안정성 또는 물리적으로 비합리적인 행동을 보이는 반면, 제안된 방법은 큰 감쇠 계수(최대 α = 40)에서도 조건부 안정성을 유지한다.
• 도메인 벽 속도는 감쇠 계수 α와 외부 자기장 he에 대해 선형 의존성을 보이며, 최소 squares 피팅 결과 기준 기울기는 α당 약 0.91~1.02 m/s, he당 약 0.91~0.95 m/s/mT 범위였다.
• 제안된 방법은 GSPM 및 SIPM보다 상당히 더 효율적이다: 1D에서는 SIPM보다 CPU 시간이 적게 소요되고 GSPM보다 훨씬 적게 소요되며, 3D에서는 GSPM과 유사한 성능이지만 동일한 정확도에서 SIPM보다 훨씬 빠르다.
• 에너지 변화 곡선은 α = 2, 5, 8일 경우 세 방법 모두 일관된 소산 패턴을 보였으며, α = 10일 경우 SIPM에서 略로 다른 패턴을 보였는데, 이는 제안된 방법이 극한의 감쇠 조건에서도 뛰어난 견고성을 확보하고 있음을 시사한다.
• 기존 방법이 실패하는 α = 40 조건에서도 안정된 도메인 벽 및 균일한 자화 상태와 같은 물리적으로 합리적인 자화 구조를 성공적으로 포착하였다.