이동력법

Shifted force method

입자 집합 내에 포함된 색인 을(를) 사용하여 충전된 입자에 작용하는 순 정전기력은 다음과 같이 주어진다.

여기서 는) 공간 좌표, j (는) 입자 i (와) 사이의 분리 거리 입자 로부터의 단위 벡터다.입자 대한 의 힘 크기, i 전하 정전력이 - 에 비례하는 상황에서 개별 입자-입자 간 상호작용은 자연적으로 장기적이며, 미립자 시스템 시뮬레이션에서 어려운 연산 문제를 제시한다 입자에 작용하는 순력을 결정하기 위해 일반적으로 에발트 또는 레크너 합계 방법을 사용한다. 하나의 대안이며 일반적으로 장거리에서의 상호작용(예: 1nm >)은 특정 시스템에서 작용하는 순힘에 미미하다는 개념에 기초하여 계산적으로 더 빠른 기법이 구형 절단 방법이다.[1] 기본 절단 방정식은 다음과 같다.

여기서 컷오프 거리다. 단순히 이 컷오프 방법을 적용하면 r 에서 힘이 불연속성을 도입하여 다른 입자가 각각의 상호작용 구 경계를 통과할 때 갑작스런 자극을 경험하게 된다. 정전기력의 특별한 경우, 힘의 크기가 경계에서 크기 때문에 이 비물리적 특성은 시뮬레이션 정확도를 손상시킬 수 있다. 이 문제를 해결하는 방법은 r 에서 힘을 0으로 전환하여 불연속성을 제거하는 것이다[2] 이는 다양한 기능으로 달성할 수 있지만, 가장 단순하고 컴퓨터적으로 효율적인 접근방식은 다음과 같이 컷오프 거리에서 정전기력 크기 값을 단순히 빼는 것이다.

앞에서 언급한 바와 같이, 이동력(SF) 방법은 일반적으로 장기적 성격인 순정전기 상호작용이 없는 시스템에 적합하다. 전기장 선별 효과를 나타내는 응축 시스템의 경우가 이에 해당한다. 최근에 인터페이스에 대한 SF 방법의 채택이 제안되었지만,[3] 비등방성 시스템(예: 인터페이스)은 SF 방법으로 정확하게 시뮬레이션되지 않을 수 있다는 점에 유의한다.[4] 또한 특정 시스템 속성(예: 에너지 의존적 관측 가능)은 다른 방법보다 SF 방법의 사용에 더 큰 영향을 받을 것이라는 점에 유의한다. SF 방법을 사용하여 특정 시스템의 특정 속성을 정확하게 판단할 수 있다고 가정하는 것은 타당성 있는 주장 없이 안전하지 않다. 만약 SF 방법의 정확도를 시험할 필요가 있다면, 이것은 정합성 시험(즉, 시뮬레이션 결과가 컷오프를 증가해도 크게 변하지 않는다는 것을 보여준다)이나, 성능이 좋다고 알려진 다른 전기 공학 기법(예: Ewald)을 통해 얻은 결과와 비교함으로써 이루어질 수 있다.[5] 대략적인 경험으로 볼 때, SF 방법으로 얻은 결과는 컷오프가 근린상호 상호작용 거리보다 최소 5배 이상 클 때 충분히 정확한 경향이 있다.

SF 방법으로, 힘의 파생상품에는 불연속성이 여전히 존재하며, 이온성 액체가 그 불연속성을 제거하기 위해 힘 방정식을 더 변경하는 것이 바람직할 수 있다.[6]

참조

  1. ^ Fennell, C. J.; Gezelter, J. D. (2006). "Is the Ewald summation still necessary? Pairwise alternatives to the accepted standard for long-range electrostatics". The Journal of Chemical Physics. 124 (23): 234104. Bibcode:2006JChPh.124w4104F. doi:10.1063/1.2206581. PMID 16821904.
  2. ^ Toxvaerd, S. R.; Dyre, J. C. (2011). "Communication: Shifted forces in molecular dynamics". The Journal of Chemical Physics. 134 (8): 081102. arXiv:1012.1116. Bibcode:2011JChPh.134h1102T. doi:10.1063/1.3558787. PMID 21361519.
  3. ^ Spohr, E. (1997). "Effect of electrostatic boundary conditions and system size on the interfacial properties of water and aqueous solutions". The Journal of Chemical Physics. 107 (16): 6342–6348. Bibcode:1997JChPh.107.6342S. doi:10.1063/1.474295.
  4. ^ Welch, D. A.; Mehdi, B. L.; Hatchell, H. J.; Faller, R.; Evans, J. E.; Browning, N. D. (2015). "Using molecular dynamics to quantify the electrical double layer and examine the potential for its direct observation in the in-situ TEM". Advanced Structural and Chemical Imaging. 1. doi:10.1186/s40679-014-0002-2.
  5. ^ Hansen, J. S.; Schrøder, T. B.; Dyre, J. C. (2012). "Simplistic Coulomb Forces in Molecular Dynamics: Comparing the Wolf and Shifted-Force Approximations". The Journal of Physical Chemistry B. 116 (19): 5738. arXiv:1108.5267. doi:10.1021/jp300750g.
  6. ^ Kale, S.; Herzfeld, J. (2011). "Pairwise Long-Range Compensation for Strongly Ionic Systems". Journal of Chemical Theory and Computation. 7 (11): 3620–3624. doi:10.1021/ct200392u. PMC 3254088. PMID 22247701.