Variation in Magnetic Moments depending on Supercell Size
Posted: Mon Feb 21, 2022 1:13 pm
Dear all,
I am looking at running some bulk calculations for some antiperovskite materials (A3BX). I have only just started with some very preliminary calculations, so I am running some tests to determine suitable parameters. I got a rather odd result - starting from the crystallographic unit cell (i.e. .cif file), the wavefunctions were converged, but oddly I got different local magnetic moments for each A atom in my unit cell, despite the fact that they ought to be symmetry equivalent (the experimentally determined unit cell is undistorted), although they were all ferromagnetically aligned. The difference between the largest and smallest local moment was 0.75 BM. I haven't attempted any optimisation yet, so the atoms should be precisely fixed on their crystallographic lattice sites.
I repeated the calculation but multiplied the unit cell to a 2x2x2 supercell, and adjusted the MP k-point mesh accordingly. This time, the A atom local magnetic moments were almost the same (identical to 1 decimal point, but varying with a range of 0.02 BM). The B and X atom local moments were unchanged from the 1x1x1 cell.
I am puzzled as to why this is the case - it would make more sense if there were a complex antiferromagnetic ordering that could not be accommodated within the symmetry imposed by the smaller cell. However, in both cases, the calculation converged on a ferromagnetic solution. Moreover, the total magnetic moment per formula unit was more or less unchanged, i.e. the average of the local magnetic moments for the A atoms in the 1x1x1 roughly equal that of each atom in the 2x2x2 cell.
I am probably missing something obvious here. Any comments or suggestions would be most appreciated.
I am looking at running some bulk calculations for some antiperovskite materials (A3BX). I have only just started with some very preliminary calculations, so I am running some tests to determine suitable parameters. I got a rather odd result - starting from the crystallographic unit cell (i.e. .cif file), the wavefunctions were converged, but oddly I got different local magnetic moments for each A atom in my unit cell, despite the fact that they ought to be symmetry equivalent (the experimentally determined unit cell is undistorted), although they were all ferromagnetically aligned. The difference between the largest and smallest local moment was 0.75 BM. I haven't attempted any optimisation yet, so the atoms should be precisely fixed on their crystallographic lattice sites.
I repeated the calculation but multiplied the unit cell to a 2x2x2 supercell, and adjusted the MP k-point mesh accordingly. This time, the A atom local magnetic moments were almost the same (identical to 1 decimal point, but varying with a range of 0.02 BM). The B and X atom local moments were unchanged from the 1x1x1 cell.
I am puzzled as to why this is the case - it would make more sense if there were a complex antiferromagnetic ordering that could not be accommodated within the symmetry imposed by the smaller cell. However, in both cases, the calculation converged on a ferromagnetic solution. Moreover, the total magnetic moment per formula unit was more or less unchanged, i.e. the average of the local magnetic moments for the A atoms in the 1x1x1 roughly equal that of each atom in the 2x2x2 cell.
I am probably missing something obvious here. Any comments or suggestions would be most appreciated.