应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » vasp 论坛自旋轨道耦合效应（SOC ）问题集锦,有问题的时候大家可以看看这个
181/2返回列表12下一页
查看: 5057  |  回复: 17
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

wuli8

荣誉版主 (知名作家)

…………

优秀版主优秀版主

[交流] vasp 论坛自旋轨道耦合效应（SOC ）问题集锦,有问题的时候大家可以看看这个 已有6人参与

Starting NonCollinear NEB calculations

Postby Neutrino » Thu Sep 15, 2011 1:39 am

Hi ,
I posted this question in the "physics questions" forum , but then realized that it is more appropriate to be posted here.

As the manual recommended, I started with non-spin calculations (ISPIN=1) to generate WAVECAR and CHGCAR for each image of my NEB calculations. Then I tried to restart the job and turnning on the non-collinear calculations. I set
ISTART=2 ; ICHARG = 1
LNONCOLLINEAR=.TRUE.
I also set NBANDS, NGX,NGY,NGZ,NGXF, NGYF,NGZF manually so that both calculations match each other. In spite of that the job stops with message:
ERROR: while reading WAVECAR, plane wave coefficients changed.
Any hints are appreciated!

As another question, does it hurt if I start the noncollinear calculations from scratch by setting:
ISTART=0 ; ICHARG = 2

Last edited by Neutrino on Thu Sep 15, 2011 1:39 am, edited 1 time in total.
Neutrino Jr. Member
Jr. Member Posts: 59Joined: Wed Mar 16, 2011 3:53 am

Top



--------------------------------------------------------------------------------






Starting NonCollinear NEB calculations

Postby admin » Fri Sep 16, 2011 4:50 pm

as in the non-collinear run all symmetry is reduced, it may well be that some of the PW coefficients have changed. In this case it would be much safer to start with ISTART=0 in order to avoid that coefficients are read from WAVECAR that do not match the basis functions generated for the non-collinear run

Last edited by admin on Fri Sep 16, 2011 4:50 pm, edited 1 time in total.
admin Administrator
Administrator Posts: 2626Joined: Tue Aug 03, 2004 8:18 am

http://cms.mpi.univie.ac.at/vasp ... 150&p=10150

[ Last edited by wuli8 on 2015-4-20 at 14:28 ]
回复此楼

» 收录本帖的淘帖专辑推荐

第一性原理相关文档 第一性原理计算-tg 自旋轨道耦合 VASP软件从入门到高深
第一性原理和电化学 第一性原理 第一性 VASP学习

» 本帖已获得的红花（最新10朵）

hy1987007 仲_______ linjiahe

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
…………
1楼 2015-04-17 21:45:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

wuli8

荣誉版主 (知名作家)

…………

优秀版主优秀版主
引用回帖:
14楼: Originally posted by linjiahe at 2015-05-26 11:36:54
楼主你好，我想请教加自旋轨道耦合，如何控制自旋轨道耦合的大小，比如百分三十，是否是要重新编译，需要的话，该怎么改怎么编译，万分感谢！！！

请参考这里http://muchong.com/bbs/taotie.php?action=view&ttid=66284
赞一下(1人) 回复此楼
…………
15楼2015-05-26 15:55:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

ljw4010

荣誉版主 (职业作家)

小木虫从头派教主
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
大神一出手就知道有没有，赞个
赞一下 回复此楼
不要说话
2楼2015-04-17 22:26:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wuli8

荣誉版主 (知名作家)

…………

优秀版主优秀版主
DOS in noncollinear calculation?

Postby zlujpdavid » Mon Oct 17, 2005 4:17 am

When I tried to look at the density of state after noncollinear caculation, the DOSCAR just gave half of the density of state information(only upspin). And it was said that my calculation is non- spinpolarized when I ran split_dos. It is not true. I did set ISPIN=2 and OUTCAR gave magnetic moment for all the atoms. What is the reason for DOSCAR only giving part of the DOS information? How to get the complete DOS information in noncollinear calculation?

Last edited by zlujpdavid on Mon Oct 17, 2005 4:17 am, edited 1 time in total.
zlujpdavid Newbie
Newbie Posts: 5Joined: Sat Oct 15, 2005 11:11 pm





Hm, which LORBIT did you set? Maybe try LORBIT = 11 or 12 ...

Last edited by Veronika on Mon Oct 24, 2005 7:52 am, edited 1 time in total.
Veronika Jr. Member
Jr. Member Posts: 70Joined: Wed Feb 02, 2005 4:14 pm

Top





DOS in noncollinear calculation?

Postby zlujpdavid » Tue Oct 25, 2005 10:32 pm

Usually, I set LORBIT =.TRUE. if I want to incude spin-obital coupling. For many times, I just did noncollinear calculation by setting LNONCOLLINEAR=.TRUE. with setting LORBIT. It is not necessary to set LORBIT, isn't it?

Last edited by zlujpdavid on Tue Oct 25, 2005 10:32 pm, edited 1 time in total.
zlujpdavid Newbie
Newbie Posts: 5Joined: Sat Oct 15, 2005 11:11 pm

Top


DOS in noncollinear calculation?

Postby Veronika » Fri Nov 04, 2005 7:01 am

If you run with vasp 4.6 or higher you should use the integer numbers for LORBIT as described in the manual (http://cms.mpi.univie.ac.at/vasp/vasp/node113.html) ... The exact interplay between LORBIT and LNONCOLLINEAR I do not know but I suggest to set LORBIT if you want to get more information on the DOS ...

Last edited by Veronika on Fri Nov 04, 2005 7:01 am, edited 1 time in total.
Veronika Jr. Member
Jr. Member Posts: 70Joined: Wed Feb 02, 2005 4:14 pm


DOS in noncollinear calculation?

Postby windiks » Fri Nov 04, 2005 10:43 am

Dear all,

the terms "spin up" and "spin down" do not exist in the concept of noncolinear magnetism or spin-orbit calculations. You obtain a magnetisation density only as defined in the paper of Hobbs, Kresse and Hafner PRB 62(17), 1156 (2000). The paper explaines how the concept of noncolinear magnetism within the local spin density theory is implemented in VASP.
As far as I know, the DOS is calculated from that magnetization density and, hence, you do not get a DOS for spin-up and spin-down states. Furthermore, the tag ISPIN has nothing to do with the tags LORBIT and LNONCOLLINEAR!
But you can calculated noncolinear magnetic moments without taking into account spin orbit coupling (i.e. LNONCOLLINEAR=.TRUE. and LORBIT=.FALSE. is possible).

Again, do not mix the terms "spin-polarized calculations" with "calculations of noncolinear magnetism". Both calculations are based on completely different concepts!

Regards,
René

Last edited by windiks on Fri Nov 04, 2005 10:43 am, edited 1 time in total.
windiks Newbie
Newbie Posts: 7Joined: Tue Feb 01, 2005 9:07 pmLocation: Switzerland







DOS in noncollinear calculation?

Postby admin » Wed Nov 16, 2005 1:13 pm

for non-collinear magnetism, the second block of data in DOSCAR
contains
E, and 4 columns for each, s,p,d, giving:
rho, m_x, m_y, m_z
with m....magnetisation

Last edited by admin on Wed Nov 16, 2005 1:13 pm, edited 1 time in total.
admin Administrator
Administrator Posts: 2626Joined: Tue Aug 03, 2004 8:18 am

Top





DOS in noncollinear calculation?

Postby lahaye » Tue Dec 26, 2006 12:20 am

I did a NON-Collinear / LSORBIT=.TRUE. / ISPIN=2 calculation,
and the second block in DOSCAR has for each atom NEDOS
rows with 37 columns. The first column is clearly the DOS energy,
but what are the other 36 columns?

Thanks,
Rob.

Last edited by lahaye on Tue Dec 26, 2006 12:20 am, edited 1 time in total.
lahaye Jr. Member
Jr. Member Posts: 98Joined: Fri Apr 14, 2006 5:08 amLocation: Suwon - Korea




DOS in noncollinear calculation?

Postby admin » Wed Dec 27, 2006 10:57 am

my last posting answers your question.
Energy , 9x4 entries for the decomposed local DOS of each atom.

As already pointed out by Rene (2 postings above), it makes absolutely NO SENSE to set ISPIN=2 (up and down) for non-clollinear runs, therfore this tag is ignored when it s read from INCAR.

Last edited by admin on Wed Dec 27, 2006 10:57 am, edited 1 time in total.
admin Administrator
Administrator Posts: 2626Joined: Tue Aug 03, 2004 8:18 am


DOS in noncollinear calculation?

Postby wellsdm » Mon Feb 11, 2008 9:28 pm

I have been doing f-electron LS-coupling runs can you explain the make-up of the DOS in that case? i have what seems to be 37 columns (E+rho,mx,my,mz for each s,p,d) then a second line of another 28 before i get another line starting with an understandable energy value. Can you please explain this. if the second line is the f-electron components, why is the first number formatted more like energy (0.000 instead of 0.0000E+00 like the others)

Last edited by wellsdm on Mon Feb 11, 2008 9:28 pm, edited 1 time in total.
wellsdm Newbie
Newbie Posts: 3Joined: Thu Nov 01, 2007 9:09 pm



DOS in noncollinear calculation?

Postby jermysmart » Wed May 25, 2011 11:27 am

whether I can get a spin-up-DOS and spin-down-DOS invovling SOC by projecting DOS to up and down

Last edited by jermysmart on Wed May 25, 2011 11:27 am, edited 1 time in total.
jermysmart Newbie
Newbie Posts: 1Joined: Tue Jun 16, 2009 1:45 pm





DOS in noncollinear calculation?

Postby chelman » Thu Oct 27, 2011 5:44 pm

Dear head admin: Following this thread, I would like to know if the file vasprun.xml contain the same information that DOSCAR. if it's true, then p4v can read it? how?
Thanks!

Last edited by chelman on Thu Oct 27, 2011 5:44 pm, edited 1 time in total.
chelman Newbie
Newbie Posts: 10Joined: Fri Mar 11, 2011 4:58 pmLocation: Buenos Aires, Argentina





DOS in noncollinear calculation?

Postby admin » Fri Oct 28, 2011 1:04 pm

Yes. p4v/control/Electronic/LocalDOS+bands control

Last edited by admin on Fri Oct 28, 2011 1:04 pm, edited 1 time in total.
admin Administrator
Administrator Posts: 2626Joined: Tue Aug 03, 2004 8:18 am




DOS in noncollinear calculation?

Postby yywang » Mon Apr 01, 2013 6:40 am

[quote="jermysmart"]whether I can get a spin-up-DOS and spin-down-DOS invovling SOC by projecting DOS to up and down[/quote]

I also want to know the answer, and how to project DOS to up and down. Has your problem been solved?

Last edited by yywang on Mon Apr 01, 2013 6:40 am, edited 1 time in total.

http://cms.mpi.univie.ac.at/vasp ... =466&p=9414
赞一下 回复此楼
…………
3楼2015-04-19 21:51:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wuli8

荣誉版主 (知名作家)

…………

优秀版主优秀版主
Spin orbit coupling energy

Postby ido » Tue Aug 27, 2013 9:59 am

Hello,

Sorry for my following naive question...

1) If I do a spin-orbit coupling calculation, does the line:

energy without entropy= -64.80983412 energy(sigma->0) = -64.80998344

inculdes the spin-orbit coupling energy?

2) How can I get only the spin-orbit coupling energy of the system?

3) While doing a calculation with ISPIN=2 from scratch and then use ICHARG = 11 with LSORBIT = .TRUE. is equivalent to a self-consistent calculation with ISPIN = 1 and LSORBIT = .TRUE. from scratch?

THANKS FROM ADVANCED!!!

Last edited by ido on Tue Aug 27, 2013 9:59 am, edited 1 time in total.
ido Newbie
Newbie Posts: 7Joined: Mon Oct 29, 2012 8:44 am


Spin orbit coupling energy

Postby admin » Tue Sep 03, 2013 9:56 am

1) no, this includes the extrapolation of the electronic entropy term to T=0
2) I suppose by doing the -otherwise- identical calc with/without SOC.
3) if LSORBIT=.True. non-collinear magnetism is used, i.e. there are
4 spinors and ISPIN=2 is ignored.
http://cms.mpi.univie.ac.at/vasp ... 661&p=13720
赞一下 回复此楼
…………
4楼2015-04-20 09:21:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wuli8

荣誉版主 (知名作家)

…………

优秀版主优秀版主
Orbital moment of a single isolated Terbium (Tb) atom: LDA+U+SOC

Postby saurabh_2802 » Thu Sep 16, 2010 8:13 pm

Dear All,

I am doing a spin orbit calculation of a single isolated Tb atom (at the centre of a 16x16x16 A3 cell and the calculation has been done at Gama point)

I follow the usual procedure i.e. first I have done the spin polarized calculation and kept the self consistence CHGCAR and WAVECAR file, and then I start spin orbit calculation form these two files.

INCAR file is given as:  

SYSTEM = Tb atom
ENCUT = 400.0
PREC = Normal
ISPIN = 2
MAGMOM = 0 0 6
ISTART = 1
ICHARG = 1
NELM = 80
NELMIN = 4
#NELMDL = -5
EDIFF = 0.0001
EDIFFG = -0.005
##############################
ISYM = 0
NSW = 0
#NBLOCK = 1
#KBLOCK = 0
IBRION = 1
POTIM = 0.1
ISIF = 2
IALGO = 48
LDIAG = .TRUE.
NSIM = 4
LPLANE= .TRUE.
#=TIME = 0.1
###################################
LDAU = .TRUE.
LDAUTYPE = 2
LDAUL = 3
LDAUU = 4.0
LDAUJ = 1.0
#LDAUPRINT=
LMAXMIX = 6
###################
LSORBIT = .TRUE.
SAXIS = 0 0 1
LORBMOM = .TRUE.
NBANDS = 32
################################
#RWIGS =
LORBIT = 11
#EMIN =
#EMAX =
ISMEAR = 0
SIGMA = 0.1
LREAL= Auto
#ROPT =
#NEDOS =
#ADDGRID = .TRUE.
#########################
#LWAVE = .FALSE.
#LCHARG = .FALSE.
#LVTOT = .FALSE.
#LELF = .FALSE.
#NPAR= 8
################################

POTCAR, KPOINTS and POSACR file is okay.. The calculation converged and the orbital momentum that I got is zero. But the value has to be 3 meuB. Although I am getting the correct spin moment, as given below.....
-------------------------------------------------------------------------
orbital moment (x)

# of ion p d f tot
------------------------------------------------
1 0.000 0.000 0.000 0.000
------------------------------------------------
0.000 0.000 0.000 0.000



orbital moment (y)

# of ion p d f tot
------------------------------------------------
1 0.000 0.000 0.000 0.000
------------------------------------------------
0.000 0.000 0.000 0.000



orbital moment (z)

# of ion p d f tot
------------------------------------------------
1 0.000 0.000 0.000 0.000
------------------------------------------------
0.000 0.000 0.000 0.000



total charge

# of ion s p d f tot
------------------------------------------------
1 2.908 6.032 0.000 9.428 18.368



magnetization (x)

# of ion s p d f tot
------------------------------------------------
1 0.000 0.000 0.000 0.000 0.000



magnetization (y)

# of ion s p d f tot
------------------------------------------------
1 0.000 0.000 0.000 0.000 0.000



magnetization (z)

# of ion s p d f tot
------------------------------------------------
1 0.072 -0.021 0.000 5.017 5.069

---------------------------------------------------------------
My question: Is my calculation right?? or I am doing something wrong?

Looking forward for your kind reply.

Thanking you and with regards,

Saurabh Ghosh

Last edited by saurabh_2802 on Thu Sep 16, 2010 8:13 pm, edited 1 time in total.
saurabh_2802 Newbie
Newbie Posts: 2Joined: Sat Feb 23, 2008 11:36 amLocation: Koltaka, India


Orbital moment of a single isolated Terbium (Tb) atom: LDA+U+SOC

Postby boris » Fri Sep 17, 2010 8:19 am

Hi,

Starting from the CHGCAR and WAVECAR files without SOC is not an appropriate method. I know that it is the method described in the vasp manual, but unfortunately it does not work.

Your calculation results show that the orbital moment is zero. The reason for this is that the nondiagonal spin components, that give rise to the orbital moment, are probably empty. You can check this by setting LDAUPRINT = 1. This will allow vasp to print the occupation matrix of the Tb f shells at the end of the OUTCAR file. There should be 8 occupation matrices (4 spin components, each with a real and imaginary part) with only the 1st and 4th spin components filled. The 2nd and 3rd spin components should be empty, hence no orbital moment.

I guess the right way of doing this is not to use the CHGCAR and WAVECAR files and to start directly the calculation with SOC (ISTART=0, ICHARG=2).

Also, you should use a non cubic supercell, something like 15x16x17 A3. This will break all symmetries in the system and will allow the calculation to reach the ground state more easily.

Finally, beware of the metastable states created by the LDA+U formalism. With SOC, it will be even worse.

Regards
http://cms.mpi.univie.ac.at/vasp ... 521&p=10521
赞一下 回复此楼
…………
5楼2015-04-20 11:54:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wuli8

荣誉版主 (知名作家)

…………

优秀版主优秀版主
How to accelerate SOC + HSE06 calculation?


How to accelerate SOC + HSE06 calculation?

Postby jzzhao » Fri Jun 04, 2010 7:47 am

Hi, all vasp users and developers,

I'm now tring to calculate some semiconductors materials with Spin-Orbital Coupling (SOC) and
Hybride Functional (HSE06) method in Vasp 5.2. I'm in trouble about this type calculation,
and ask for your help.


My calculation process is:

First, general static DFT calculation.
Second, standart SOC calculation (LSORBIT=T), and obtain the WAVECAR file.
Third, based on this WAVECAR file, add LHFCALC = .TRUE. and other related HF parameters to perform HF type calculation.


Here are my SOC + HF input files:

======== INCAR ====================================================
ENCUT = 350
LMAXMIX=6

LSORBIT=T

ISTART = 1
LHFCALC = .TRUE.
HFSCREEN = 0.2
ALGO = Damped
TIME = 0.4
#ENCUTFOCK = 0
#NKRED = 2

======== KPOINTS ==================================================
Auto
0
Gamma
8 8 8
0 0 0

======= POSCAR ====================================================
Cd As Ge
1.00000000000000
-2.9596082483964499 2.9596082483964552 6.0150253167398420
2.9596082483964552 -2.9596082483964525 6.0150253167398402
2.9596082483964539 2.9596082483964539 -6.0150253167398411
2 4 2
Direct
0.0000000000000000 0.0000000000000000 0.0000000000000000
0.7500000000000000 0.2500000000000000 0.5000000000000000
0.3750000000000000 0.3566726093521447 0.4816726093521447
0.8750000000000000 0.8933273906478553 0.5183273906478553
0.6433273906478553 0.1250000000000000 0.0183273906478553
0.1066726093521447 0.6250000000000000 0.9816726093521447
0.5000000000000000 0.5000000000000000 0.0000000000000000
0.2500000000000000 0.7500000000000000 0.5000000000000000

========= POTCAR ==================================================
potpaw_PBE for every element.




My calculating enviroment is a blade-based system:

Each node contains four AMD Opteron Quad-Core 64-bit processors (16 cores in all) on a
single board, as an SMP unit. The core frequency is 2.3 GHz and supports 4 floating-point
operations per clock period with a peak performance of 9.2 GFLOPS/core or 128 GFLOPS/node.
Each node contains 32 GB of memory. The memory subsystem has a 1.0 GHz HyperTransport
system Bus, and 2 channels with 667 MHz DDR2 DIMMS. Each socket possesses an independent
memory controller connected directly to an L3 cache.



I used 4 nodes(total 64 CPUs) to perform this calculation, but one electronic step is
NOT done with more than 20 hours. In order to accelerate the calculation, I also try
to uncomment the #ENCUTFOCK = 0 and #NKRED = 2, but there is an error about memory:

----------------------------------------------------------------------------------
forrtl: severe (41): insufficient virtual memory
Image PC Routine Line Source
libintlc.so.5 00002B102C3703A6 Unknown Unknown Unknown
libintlc.so.5 00002B102C36F338 Unknown Unknown Unknown
libifcore.so.5 00002B102BD2829E Unknown Unknown Unknown
libifcore.so.5 00002B102BCB2DA4 Unknown Unknown Unknown
libifcore.so.5 00002B102BCEF6EF Unknown Unknown Unknown
libifcore.so.5 00002B102BCEF5DF Unknown Unknown Unknown
vasp5.2 00000000005CC44A Unknown Unknown Unknown
vasp5.2 00000000005CBC26 Unknown Unknown Unknown
vasp5.2 00000000004293D0 Unknown Unknown Unknown
vasp5.2 000000000040A3E2 Unknown Unknown Unknown
libc.so.6 0000003D58A1C3FB Unknown Unknown Unknown
vasp5.2 000000000040A32A Unknown Unknown Unknown
MPI process terminated unexpectedly
----------------------------------------------------------------------------------



So, I have two questions:

1. Are there some parameters I don't sensibly used? In order to accelerate the calculation
and reduce the demand for computing resources, which parameters should be adjusted?

2. In Vasp 5.2, is the calculation with a combination of SOC and HF stable and reliable?
If yes, how much the computing resources are needed for this type calculation?


I'm looking forward to your reply, Thank you very much!

Last edited by jzzhao on Fri Jun 04, 2010 7:47 am, edited 1 time in total.
jzzhao Newbie
Newbie Posts: 12Joined: Mon Dec 11, 2006 8:04 am


How to accelerate SOC + HSE06 calculation?

Postby admin » Thu Jul 01, 2010 9:02 am

sorry, NKRED and ENCUTFOCK cannot lead to the errors you show if they are set properly (it seems that an error message cannot be written due to problems with the libraries )
1) the most efficient way to speed up HF calculations in fact is to set NKRED
2) yes, this depends on the size of your system. pelase note that SOC calculations do not use any symmetry, which increases the memory demand significanlty, especially if you compare to to runs with highly symmetric cells.

http://cms.mpi.univie.ac.at/vasp ... ?f=4&t=7516
赞一下 回复此楼
…………
6楼2015-04-20 13:28:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wuli8

荣誉版主 (知名作家)

…………

优秀版主优秀版主
LDA+U, Spin-Orbit Coupling and Orbital Moment
Moderators: Moderator, Global Moderator
Post a reply
LDA+U, Spin-Orbit Coupling and Orbital Moment

Postby boris » Fri Apr 16, 2010 12:55 pm

Dear all,

I'm using vasp 5.2 and I want to calculate the orbital moment carried by uranium atoms in uranium dioxide UO2.

I'm using the LDA+U formalism and I included the spin-orbit coupling in the calculations. I'm printing the orbital moment by setting LORBMOM=.TRUE.

The result is, however, inconsistent with the value given in the literature. I'm supposed to find something like -3.5 MuB for the orbital moment and vasp yields 0.0 MuB.

Is there any known issue with the spin-orbit coupling implementation in vasp together with the LDA+U formalism?

Thank you for your answers.

Regards

Last edited by boris on Fri Apr 16, 2010 12:55 pm, edited 1 time in total.
boris Full Member
Full Member Posts: 106Joined: Wed Jan 16, 2008 1:27 pmLocation: France

Top

LDA+U, Spin-Orbit Coupling and Orbital Moment

Postby Danny » Fri Apr 16, 2010 2:13 pm

Maybe you need to initialise the orbital momenta yourself (similar to magmom) on a value close to what it should be?

I don't seem to be able to find any reference to LORBMOM in the online manual...

cheers
Danny

Last edited by Danny on Fri Apr 16, 2010 2:13 pm, edited 1 time in total.

HIVE-STM program for VASP
Danny Full Member
Full Member Posts: 201Joined: Thu Nov 02, 2006 4:35 pmLocation: Ghent, Belgium


LDA+U, Spin-Orbit Coupling and Orbital Moment

Postby boris » Fri Apr 16, 2010 2:50 pm

There is indeed no reference to LORBMOM in the manual. I managed to find this keyword thanks to a question asked on this forum.

Unfortunately, there was nothing more about it.

Last edited by boris on Fri Apr 16, 2010 2:50 pm, edited 1 time in total.
boris Full Member
Full Member Posts: 106Joined: Wed Jan 16, 2008 1:27 pmLocation: France


LDA+U, Spin-Orbit Coupling and Orbital Moment

Postby admin » Mon Apr 19, 2010 11:17 am

just a silly question: have you given all other parameters to do SOC as well (LSORBIT, LNONCOLLINEAR, SAXIS, ISYM=-1)?

Last edited by admin on Mon Apr 19, 2010 11:17 am, edited 1 time in total.
admin Administrator
Administrator Posts: 2626Joined: Tue Aug 03, 2004 8:18 am


LDA+U, Spin-Orbit Coupling and Orbital Moment

Postby boris » Mon Apr 19, 2010 11:52 am

Yes, I have set LSORBIT = .TRUE., LNONCOLLINEAR is automatically set to TRUE as well, SAXIS = 0 0 1 but I set ISYM to 0, not to -1. Does this make any difference?

Thank you for your answer

Last edited by boris on Mon Apr 19, 2010 11:52 am, edited 1 time in total.
boris Full Member
Full Member Posts: 106Joined: Wed Jan 16, 2008 1:27 pmLocation: France
LDA+U, Spin-Orbit Coupling and Orbital Moment

Postby boris » Mon Apr 19, 2010 1:01 pm

I have performed another calculation using ISYM=-1. There is no difference compared to the ISYM=0 case.

Here is the result:



orbital moment (z)

 # of ion     p       d       f       tot
 ------------------------------------------------
  1        -0.002  -0.002  -0.368  -0.372
  2        -0.002  -0.002  -0.368  -0.372
  3         0.002   0.002   0.368   0.372
  4         0.002   0.002   0.368   0.372
  5         0.000   0.000   0.000   0.000
  6         0.000   0.000   0.000   0.000
  7         0.000   0.000   0.000   0.000
  8         0.000   0.000   0.000   0.000
  9         0.000   0.000   0.000   0.000
 10         0.000   0.000   0.000   0.000
 11         0.000   0.000   0.000   0.000
 12         0.000   0.000   0.000   0.000
 ------------------------------------------------
            0.000   0.000   0.000   0.000



There is no orbital moment in the other directions. The spin moment is given by :



magnetization (z)

# of ion     s       p       d       f       tot
------------------------------------------------
  1        0.013  -0.001   0.030   1.967   2.010
  2        0.013  -0.001   0.030   1.967   2.010
  3       -0.013   0.001  -0.030  -1.967  -2.010
  4       -0.013   0.001  -0.030  -1.967  -2.010
  5        0.000   0.000   0.000   0.000   0.000
  6        0.000   0.000   0.000   0.000   0.000
  7        0.000   0.000   0.000   0.000   0.000
  8        0.000   0.000   0.000   0.000   0.000
  9        0.000   0.000   0.000   0.000   0.000
 10        0.000   0.000   0.000   0.000   0.000
 11        0.000   0.000   0.000   0.000   0.000
 12        0.000   0.000   0.000   0.000   0.000
------------------------------------------------
tot        0.000   0.000   0.000   0.000   0.000



Therefore, it makes a total magnetic moment of 1.64 MuB, which I think is in very good agreement with the experimental value (1.74 MuB).

However, other people have calculated it using first principles DFT+U calculations and found the orbital moment to be -3.5.

I can't figure out what is wrong.

Thank you




<span class='smallblacktext'>[ Edited Mon Apr 19 2010, 03:08PM ]</span>
Last edited by boris on Mon Apr 19, 2010 1:01 pm, edited 1 time in total.
boris Full Member
Full Member Posts: 106Joined: Wed Jan 16, 2008 1:27 pmLocation: France

LDA+U, Spin-Orbit Coupling and Orbital Moment

Postby boris &raquo; Tue Apr 20, 2010 5:38 pm

Ok, I found what was going wrong.

It is written in the vasp manual that one should first calculate a collinear WAVECAR and CHGCAR file, and then add the spin-orbit coupling while restarting with the two files.

Actually,it appears that starting from a charge density without any spin-orbit coupling terms is not a good method.

I started from scratch directly with the spin-orbit coupling and it gave an orbital moment of -3 MuB.

Last edited by boris on Tue Apr 20, 2010 5:38 pm, edited 1 time in total.

http://cms.mpi.univie.ac.at/vasp ... 7104&p=7116
赞一下 回复此楼
…………
7楼2015-04-20 14:22:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

hy1987007
送红花一朵
谢谢分享！
8楼2015-04-21 14:31:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

仲_______

金虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
送红花一朵
楼主您好，感谢您的分享。
我对加自旋轨道耦合后的DOSCAR很感兴趣，想知道其中的37列数据分别是什么。但是您文中的回复是“my last posting answers your question”。我在VASP论坛上没有搜到这个帖子，请问您能解释一下这37列数据或者给个链接吗？
非常感谢！

引用：
DOS in noncollinear calculation?

Postby admin &raquo; Wed Dec 27, 2006 10:57 am

my last posting answers your question.
Energy , 9x4 entries for the decomposed local DOS of each atom.

As already pointed out by Rene (2 postings above), it makes absolutely NO SENSE to set ISPIN=2 (up and down) for non-clollinear runs, therfore this tag is ignored when it s read from INCAR.

Last edited by admin on Wed Dec 27, 2006 10:57 am, edited 1 time in total.
admin Administrator
Administrator Posts: 2626Joined: Tue Aug 03, 2004 8:18 am
赞一下 回复此楼
9楼2015-04-24 09:48:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wuli8

荣誉版主 (知名作家)

…………

优秀版主优秀版主
引用回帖:
9楼: Originally posted by 仲_______ at 2015-04-24 09:48:33
楼主您好，感谢您的分享。
我对加自旋轨道耦合后的DOSCAR很感兴趣，想知道其中的37列数据分别是什么。但是您文中的回复是“my last posting answers your question”。我在VASP论坛上没有搜到这个帖子，请问您能解 ...

http://cms.mpi.univie.ac.at/vasp ... =466&p=9414
赞一下 回复此楼
…………
10楼2015-04-24 11:06:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 wuli8 的主题更新
181/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭