如何确定蒙卡模拟中vortex位置？ 第1页

一个格子上的所有点下标顺序记为1,2,..m（一般逆时针）

如vortex在格子中，会有 大于0（是2pi的整数倍）。antivortex这个值应该是2pi的负整数倍。

然后vortex的位置就大致认为是在格点中心了...虽然这样不准，但如果要算的是vortex对间距的话也还可以接受。毕竟距离小于晶格常数，根本观测不到的vortex对都被忽视掉了嘛...

(这个具体位置确实应该可以进一步优化，不过一时没有想到更简洁的办法。)

Vortex对之间距离的计算...

那我就引自己文章打广告了，vortex 对之间距离统计使用的优化算法也写了，简单说就是最小化系统每对间的能量（~ln(d)）。

这个方法不敢称最佳，但应该已经比较简便。是考虑了如果一对很大vortex中间如果夹了一对很紧密的vortex，物理上是倾向将其分为一大一小而不是两个都很大的（考虑到它们对周围自旋的影响，小的那一对影响范围不会大过外面的）。

文章用的是立方格子，推广到三角格子，六角格子也都可以，不过就得考虑有的格子上的vortex的charge可能不是1。

另外一提，Hungarian Method 复杂度为O(n^3)，临界点（及以上，虽然这时vortex对的能量定义似乎有些不够严格）vortex多的时候配对还是挺烦人的。