matlab 磷烯Hubbard自洽平均场计算程序完善版本

最新推荐文章于 2023-03-21 12:49:49 发布

wwxy261

最新推荐文章于 2023-03-21 12:49:49 发布

阅读量949

点赞数 3

分类专栏：算法

本文链接：https://blog.csdn.net/wwxy1995/article/details/82846797

版权

设置四种磁性构型

function [Nupavg,Ndownavg]=Initialization_1D_FM(m,Ny)
% m 为原子总数
% Ny 为y方向原子个数
Nupavg=0.5*ones(m,1);   % 将所有的<Nupi>取为1作初始值,自旋向上
Ndownavg=0.5*ones(m,1); % 将所有的<Ndowni>取为1作初始值，自旋向下
Nupavg(1:Ny:m,1)=   0.8; %上边界上自旋电子密度
Ndownavg(1:Ny:m,1)= 0.2; %上边界下自旋电子密度
Nupavg(Ny:Ny:m,1)=  0.8; %下边界上自旋电子密度
Ndownavg(Ny:Ny:m,1)=0.2; %下边界下自旋电子密度
Nupavg=diag(Nupavg);
Ndownavg=diag(Ndownavg);

AFM1

function [Nupavg,Ndownavg]=Initialization_1D_AM(m,Ny)
% m 为原子总数
% Ny 为y方向原子个数
Nupavg=0.5*ones(m,1);   % 将所有的<Nupi>取为1作初始值,自旋向上
Ndownavg=0.5*ones(m,1); % 将所有的<Ndowni>取为1作初始值，自旋向下
Nupavg(1:Ny:m,1)=   0.8; %上边界上自旋电子密度
Ndownavg(1:Ny:m,1)= 0.2; %上边界下自旋电子密度
Nupavg(Ny:Ny:m,1)=  0.2; %下边界上自旋电子密度
Ndownavg(Ny:Ny:m,1)=0.8; %下边界下自旋电子密度
Nupavg=diag(Nupavg);
Ndownavg=diag(Ndownavg);

AFM2

function [Nupavg,Ndownavg]=Initialization_1D_AM2(m,Ny)
% m 为原子总数
% Ny 为y方向原子个数
Nupavg=0.5*ones(m,1);   % 将所有的<Nupi>取为1作初始值,自旋向上
Ndownavg=0.5*ones(m,1); % 将所有的<Ndowni>取为1作初始值，自旋向下
Ndownavg(1:2*Ny:m,1)=0.2; %上边界前半段下自旋电子密度
Nupavg(Ny+1:2*Ny:m,1)=0.2; %上边界后半段上自旋电子密度
Ndownavg(Ny+1:2*Ny:m,1)=0.8; %上边界后半段下自旋电子密度
Nupavg(Ny:2*Ny:m,1)=0