24小时热门版块排行榜    

查看: 5702  |  回复: 1
【悬赏金币】回答本帖问题,作者@邓邓同学将赠送您 10 个金币

@邓邓同学

新虫 (初入文坛)

[求助] vasp计算 已有1人参与

磁各向异性能如何计算,求助各位大佬
发自小木虫手机客户端
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tfang

木虫之王 (文学泰斗)

木虫008之大懒虫

【答案】应助回帖

感谢参与,应助指数 +1
磁各向异性能(Magnetic Anisotropy Energy, MAE)的计算没有一个单一的公式,它取决于各向异性的类型和采用的理论方法。但其核心思想是一致的:磁各向异性能是沿不同方向磁化时体系总能量的差值。

基本概念
磁各向异性能,是指磁性材料沿不同方向磁化时,其自由能会有所不同。其中,能量最低的方向叫易磁化方向,能量最高的则叫难磁化方向。磁各向异性能的大小,就定义为这两个方向上自由能的差值。它的单位是 J/m³(焦耳每立方米)。

唯象模型:用各向异性常数描述,这是最常用的方法,通过引入磁各向异性常数(如 K)来构建能量的解析表达式。

单轴各向异性:这是最常见的情况,公式为:
E = K₁·sin²θ + K₂·sin⁴θ + ...
E:各向异性能。
K₁, K₂:各向异性常数(J/m³)。
θ:磁化强度方向与易磁化轴方向的夹角。
对于只取一阶近似的简单单轴各向异性,公式可简化为 E = K·sin²θ。此时,易磁化轴(θ=0°)能量为0,难磁化轴(θ=90°)能量为K。
立方各向异性:对于立方晶体,其表达式更复杂,通常包含方向余弦的高次项。

第一性原理计算:从电子结构出发
这种方法基于密度泛函理论(DFT),通过计算不同磁化方向的体系总能量差来得到MAE。由于MAE通常只有总能量的百万分之一量级,计算精度要求极高。
主流的第一性原理计算软件(如VASP)操作思路如下:
设置自旋方向:在VASP中通过 SAXIS 标签来设定磁化方向。
计算总能:分别计算沿易磁化方向(如θ=0°)和难磁化方向(如θ=90°)的体系总能 E(θ=0°) 和 E(θ=90°)。
计算能量差:MAE即为这两个能量的差值,MAE = E(θ=90°) - E(θ=0°)。
注意:计算时必须开启自旋-轨道耦合(SOC) 效应,因为它是磁晶各向异性的主要物理起源。

其他计算方法
除了上述两种主流方法,还有一些特定场景下的计算方法:
扭矩法(Torque Method):通过计算磁矩偏离平衡位置时受到的扭矩来间接获得MAE。这种方法在某些情况下收敛更快。
微磁学模拟:在宏观尺度下,将磁各向异性能作为总能量的一项(包含形状各向异性、磁晶各向异性等)进行模拟。
计算磁各向异性能的核心是求能量差。具体采用唯象模型还是第一性原理计算,取决于你的研究尺度和精度要求。
如果是宏观、唯象的研究,使用各向异性常数的解析公式(如 E=K·sin²θ)会更方便。
如果是微观、机理层面的探索,则需要借助第一性原理计算(如VASP)来获得精确的MAE值。
同志们:别炒股,风险太大了,还是做豆腐最安全!做硬了是豆腐干,做稀了是豆腐脑,做薄了是豆腐皮,做没了是豆浆,放臭了是臭豆腐!稳赚不亏呀!何乐而不为!
2楼2026-07-05 23:06:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 @邓邓同学 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 zSbvVX9JFl9J 2026-11-09 6/300 2026-11-11 07:06 by Z7qfCsdIiD1T
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 FChKX8meu1ds 2026-11-09 9/450 2026-11-11 04:54 by Z7qfCsdIiD1T
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 0TEMgxdT1eIw 2026-11-09 6/300 2026-11-11 04:07 by Z7qfCsdIiD1T
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 zSbvVX9JFl9J 2026-11-09 6/300 2026-11-11 03:56 by Z7qfCsdIiD1T
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 FChKX8meu1ds 2026-11-09 6/300 2026-11-11 01:54 by Z7qfCsdIiD1T
[考研] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 cqv1zKnzd7hr 2026-11-09 5/250 2026-11-11 01:24 by Z7qfCsdIiD1T
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 cqv1zKnzd7hr 2026-11-09 6/300 2026-11-11 01:21 by Z7qfCsdIiD1T
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 zSbvVX9JFl9J 2026-11-09 5/250 2026-11-11 00:59 by Z7qfCsdIiD1T
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 FChKX8meu1ds 2026-11-09 4/200 2026-11-10 20:51 by Z7qfCsdIiD1T
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 cqv1zKnzd7hr 2026-11-09 5/250 2026-11-10 20:45 by Z7qfCsdIiD1T
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 kvewwGtKktqB 2026-11-04 12/600 2026-11-09 13:26 by Equinoxhua
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +5 5sJFDKp2kwuI 2026-11-06 15/750 2026-11-09 11:41 by GwrQVVR6Oenx
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 5sJFDKp2kwuI 2026-11-06 10/500 2026-11-09 08:00 by IhnEQauGmT3K
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 5sJFDKp2kwuI 2026-11-06 12/600 2026-11-09 06:32 by l48hFd8Nwae2
[公派出国] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 5sJFDKp2kwuI 2026-11-06 8/400 2026-11-09 00:14 by l48hFd8Nwae2
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +6 ttH3f1oA6ace 2026-11-04 12/600 2026-11-06 16:10 by 5sJFDKp2kwuI
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 ttH3f1oA6ace 2026-11-04 11/550 2026-11-06 16:10 by 5sJFDKp2kwuI
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 kvewwGtKktqB 2026-11-04 14/700 2026-11-06 05:10 by qL1cU46YIdK0
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 7/350 2026-11-06 02:13 by qL1cU46YIdK0
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 ttH3f1oA6ace 2026-11-04 10/500 2026-11-06 02:04 by qL1cU46YIdK0
信息提示
请填处理意见