[논문 리뷰] Solving the stochastic Landau-Lifshitz-Gilbert-Slonczewski equation for monodomain nanomagnets : A survey and analysis of numerical techniques

Sebastian Ament, Nikhil Rangarajan외 2명|arXiv (Cornell University)|2016. 07. 15.

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

한 줄 요약

이 논문은 단일도메인 나노자석에서 스토케스틱 랑다우-리프시츠-길버트-슬론체프스키(s-LLGS) 방정식을 해결하기 위한 수치적 방법을 평가하며, 스트라토노비치 미적분과 수렴 정확도에 중점을 둔다. 소음이 작은 SDE에 적합한 RK4-Heun 방법을 제안하고, 암시적 중간점, 헤운, 오일러-헤운 방법을 비교 평가하며, 표준 적분기법을 사용하는 SPICE 기반 시뮬레이터가 곱셈형 소음으로 인해 실패할 수 있음을 밝히며, 삼각법적 방법은 제한된 시간 간격 정확도를 갖지만 스트라토노비치 해를 정확히 도출함을 보여준다.

ABSTRACT

The stochastic Landau-Lifshitz-Gilbert-Slonczewski (s-LLGS) equation is widely used to study the temporal evolution of the macrospin subject to spin torque and thermal noise. The numerical simulation of the s-LLGS equation requires an appropriate choice of stochastic calculus and numerical integration scheme. In this paper, we comprehensively evaluate the accuracy and complexity of various numerical techniques to solve the s-LLGS equation. We focus on implicit midpoint, Heun, and Euler-Heun methods that converge to the Stratonovich solution of the s-LLGS equation. By performing numerical tests for both strong (path-wise) and weak (statistical) convergence, we quantify the accuracy of various numerical schemes used to solve the s-LLGS equation. We demonstrate a new method intended to solve Stochastic Differential Equations (SDEs) with small noise (RK4-Heun), and test its capability to handle the s-LLGS equation. We also discuss the circuit implementation of nanomagnets for large-scale SPICE-based simulations. We evaluate the efficacy of SPICE in handling the stochastic dynamics of the multiplicative noise in the s-LLGS equation. Numerical schemes such as Euler and Gear, typically used by SPICE-based circuit simulators do not yield the expected outcome when solving the Stratonovich s-LLGS equation. While the trapezoidal method in SPICE does solve for the Stratonovich solution, its accuracy is limited by the minimum time step of integration in SPICE. We implement the s-LLGS equation in both its cartesian and spherical coordinates form in SPICE and compare the stability and accuracy of the two implementations. The results in this paper will serve as guidelines for researchers to understand the tradeoffs between accuracy and complexity of various numerical methods and the choice of appropriate calculus to solve the s-LLGS equation.

연구 동기 및 목표

• 단일도메인 나노자석에서 s-LLGS 방정식을 해결하기 위한 수치적 방법의 정확도와 계산 복잡도를 체계적으로 평가하는 것.
• s-LLGS 방정식의 스트라토노비치 해에 수렴하는 가장 적합한 수치적 방법을 규명하는 것.
• s-LLGS 방정식에 내재된 곱셈형 소음 처리에 있어 SPICE 기반 회로 시뮬레이터의 한계를 조사하는 것.
• 나노자기 동역학에 관련된 소음이 작은 SDE에서 더 높은 정확도를 확보하기 위해 새로운 하이브리드 RK4-Heun 방법을 제안하고 시험하는 것.
• SPICE 내에서 카르테시안 좌표계와 구면좌표계 간의 s-LLGS 구현에서의 안정성과 정확도를 비교하는 것.

제안 방법

• s-LLGS 방정식의 물리적 일관성을 확보하기 위해 스트라토노비치 스토케스틱 미적분을 적용하여 자화 동역학에서의 소음 상호작용을 정확히 모델링한다.
• 정확도와 안정성을 평가하기 위해 강한(경로 기반) 및 약한(통계적) 수렴을 기반으로 세 가지 수치적 방법—암시적 중간점, 헤운, 오일러-헤운—을 시험한다.
• 소음이 작은 SDE를 해결하기 위한 새로운 RK4-Heun 방법을 제안하며, RK4의 4차 정확도와 헤운의 안정성의 장점을 융합하여 나노자기 시스템에서의 수렴성을 향상시킨다.
• SPICE 내에서 카르테시안 좌표계와 구면좌표계 양쪽 모두에 s-LLGS 방정식을 구현하여 동일한 스토케스틱 역학 하에서 수치적 행동과 안정성을 비교한다.
• 표준 SPICE 적분기법(Euler, Gear)과 삼각법적 방법이 스트라토노비치 해를 어떻게 처리하는지 분석하며, 시간 간격 제약으로 인한 한계를 규명한다.
• 크기 없는 s-LLGS 방정식을 자화, 등效장, 스핀 전류의 정규화를 통해 유도하여 다양한 나노자기 기하학 및 재료 매개변수에 대해 일관되고 확장 가능한 수치적 테스트를 가능하게 한다.

실험 결과

연구 질문

• RQ1 스트라토노비치 해에 수렴하는 s-LLGS 방정식을 해결하기 위해 어떤 수치적 통합 방법이 정확도와 계산 복잡도 사이의 최적의 균형을 이루는가?
• RQ2 소음이 작은 나노자기 시스템에서 제안된 RK4-Heun 방법은 전통적 방법과 비교해 약한 수렴과 강한 수렴 측면에서 어떻게 성능을 발휘하는가?
• RQ3 왜 표준 SPICE 적분기법(Euler, Gear)은 회로 수준 시뮬레이션에서 널리 사용되지만 스트라토노비치 s-LLGS 방정식을 정확히 시뮬레이션하지 못하는가?
• RQ4 SPICE 기반 시뮬레이션에서 카르테시안 좌표계와 구면좌표계로 표현된 s-LLGS 방정식 간의 안정성과 정확도의 상호 교환은 어떻게 이루어지는가?
• RQ5 SPICE 내 삼각법적 방법이 s-LLGS 방정식을 정확히 해결할 수 있는 정도는 어느 정도이며, 그 정확도에 제약을 가하는 요소는 무엇인가?

주요 결과

• 암시적 중간점과 헤운 방법은 오일러 기반 방법보다 우수한 약한 수렴과 강한 수렴을 보이며, 특히 고감쇠 영역에서 헤운 방법이 더 뛰어난 안정성을 확보한다.
• 제안된 RK4-Heun 방법은 고차 정확도의 약한 수렴을 달성하며, 통계적 모멘트 측면에서 표준 오일러 및 헤운 방법보다 더 높은 정확도를 보인다.
• 표준 SPICE 적분기법(Euler 및 Gear)은 곱셈형 소음을 잘못 처리하여 정확한 스트라토노비치 해에 수렴하지 못하며, 물리적으로 비합리적인 동역학을 유도한다.
• SPICE 내 삼각법적 방법은 스트라토노비치 s-LLGS 방정식을 정확히 해결하지만, 최소 통합 시간 간격의 제약으로 인해 고속 시뮬레이션에서 정확도가 제한된다.
• SPICE 내에서 s-LLGS 방정식의 구면좌표계 표현은 카르테시안 표현보다 더 뛰어난 수치적 안정성을 보이며, 특히 카르테시안 좌표계에서 특이점이 발생하는 극 부근에서 두드러진다.
• 논문에서 유도한 차원 없는 s-LLGS 방정식은 다양한 나노자기 기하학과 재료 매개변수 간에 일관되고 확장 가능한 벤치마킹을 가능하게 한다.