百科问答小站 logo
百科问答小站 font logo



Python 有哪些能够应用在物理(尤其是凝聚态当中)的 package?主要功能是什么? 第1页

  

user avatar   alex-34-8 网友的相关建议: 
      

做微磁学模拟可能会用到oommf这样一个开源软件,它有一个python版本joommf = jupyter + oommf, 最近推出了船新版本

官网提供了一个快速上手的例子,计算一个vortex在磁场下的演化,最后计算了vortex的winding number。有兴趣可以试着模拟一下skyrmion。

首先定义了一个涡旋态:

       # initial magnetisation state def m_init(point):     x, y, z = point     c = 1e9  # (1/m)     return (-c*y, c*x, 0.1)     

然后求解能量最小值来弛豫系统的磁矩构型:

       import oommfc as mc  # Micromagnetic Calculator  md = mc.MinDriver() md.drive(system)  system.m.plane('z').mpl()     

winding number的计算公式是

       import math  m = system.m.orientation.plane('z') q = m @ (m.derivative('x') & m.derivative('y')) S = 1/(4*math.pi) * df.integral(q * df.dx*df.dy) S     

得到0.40036,与理论值0.5比较接近。对于skyrmion,winding number为1。

求解器的功能与oommf相似,但是设置起来要简单一些,由于笔者已经很久没用过oommf,故不再赘述,有兴趣的读者可自行探索。


user avatar   yongle-li-86 网友的相关建议: 
      

PySCF

Psi4




  

相关话题

  为什么自然科学里很多公式符合平方反比定律,为什么指数是 2 不是 3? 
  刚被LHC发现的Xicc++这样的重子还有多少种?为什么重子物质在自然条件下观测不到? 
  阿尔法解决蛋白质折叠的突破,是否预示人类需要放弃寻找精确的物理公式改用AI方法预测自然规律? 
  量子化学和计算凝聚态物理有什么共同点和不同点? 
  为什么宇宙当中的温度会有上下限? 
  光子怎么产生的? 
  一维谐振子的解中,为什么量子数 n 不能为小数? 
  对物理学研究而言,发一篇 Science 或 Nature 是什么水平? 
  普朗克长度是否在瓦解文明的意义? 
  E²=p²+m² 这个方程是如何推导出来的? 

前一个讨论
有哪些能轻松穿透人类真皮层的小生物?
下一个讨论
蝗虫等小型植食动物是如何消化食物的?





© 2024-06-01 - tinynew.org. All Rights Reserved.
© 2024-06-01 - tinynew.org. 保留所有权利