[Pw_forum] SOC didn't make difference in GaAs band structure

"Masakatsu ITO(伊藤)" m-ito at aist.go.jp
Thu Mar 13 06:35:06 CET 2014


Dear Alexandra,

Thank you for your advice. I changed my input files
to use fully relativistic pseudopotentials,
Ga.rel-pbe-dn-kjpaw_psl.0.2.UPF
As.rel-pbe-n-kjpaw_psl.0.2.UPF  .

But the resultant band plotting with spin orbit coupling
is not different from that without SOC.

Probably I made other mistakes. I'd be really grateful
if you or anybody further helps me correct my input files.

with my best regards,
Masa

Masakatsu ITO

Green Nanoelectronics Collaborative Research Center
National Institute of Advanced Industrial Science and Technology
(AIST)


(3/13/14, 12:58 PM), "Alexandra Carvalho (雅丽)" wrote:
> Dear Masakatsu,
>
> My understanding is that you should use a fully relativistic 
> pseudopotential.
> Ga.pbe-dn-kjpaw_psl.0.2.UPF is scalar relativistic.
>
> Alexandra
> On 13/03/2014 10:43, "Masakatsu ITO(伊藤)" wrote:
>> Dear all,
>>
>> I'm trying to calculate GaAs band structure with
>> spin orbit coupling (SOC), but am sorry that I'm
>> relatively a novice about density functional theory.
>> I still don't get the right band structure of GaAs.
>> Could you please help me correct my mistakes?
>>
>> The problem is that SOC didn't make difference
>> in the band structure plot of GaAs. As you see
>> in the attached ps files, "gaas_bands.ps" with SOC is
>> quite the same to the band structure plot without SOC.
>>
>>  I plotted those bands along L - Gamma - X path and
>> guess that the Fermi level is around 7.0 eV.
>> I did this band calculation using the attached input files
>> like the following.
>>
>> $ pw.x < gaas.scf.david.in > gaas.scf.david.out
>> $ pw.x < gaas.bands.david.in > gaas.bands.david.out
>> $ bands.x < gaas.post.bands.in > gaas.post.bands.out
>> $ plotband.x < gaas.plotband.in > gaas.plotband.out
>>
>> Probably I miss some important parameters to
>> deal with SOC, but I cannot figure it out from
>> the documents in the network.  So I'd be grateful
>> if you tell me about my possible mistakes or
>> give me any advice.
>>
>> with my best regards,
>> Masa
>>
>> Masakatsu ITO
>>
>> Green Nanoelectronics Collaborative Research Center
>> National Institute of Advanced Industrial Science and Technology
>> (AIST)
>>
>>
>>
>>
>>
>> _______________________________________________
>> Pw_forum mailing list
>> Pw_forum at pwscf.org
>> http://pwscf.org/mailman/listinfo/pw_forum
>
>
> -- 
> Alexandra Carvalho (雅丽)
> Research Fellow
> NUS Graphene Research Centre, Office S16-06-07
> Singapore
>
>
> _______________________________________________
> Pw_forum mailing list
> Pw_forum at pwscf.org
> http://pwscf.org/mailman/listinfo/pw_forum

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.quantum-espresso.org/pipermail/users/attachments/20140313/871976cb/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: gaas_bands.ps
Type: application/postscript
Size: 81864 bytes
Desc: not available
URL: <http://lists.quantum-espresso.org/pipermail/users/attachments/20140313/871976cb/attachment.ps>
-------------- next part --------------
 &control
    calculation = 'scf'
    restart_mode='from_scratch',
    prefix='gaas',
    pseudo_dir = '/home/mito/espresso-5.0.2/pseudo',
    outdir='./temp'
 /
 &system
    ibrav=  2, celldm(1) =10.6827, nat=  2, ntyp= 2,
    ecutwfc =80.0,
    noncolin = .true.,
    lspinorb=.true.,
    starting_magnetization(1)=0.0,
    input_dft = 'PBE'
 /
 &electrons
    diagonalization='david'
    mixing_mode = 'plain'
    mixing_beta = 0.7
    conv_thr =  1.0d-8
 /
ATOMIC_SPECIES
 Ga  69.72  Ga.rel-pbe-dn-kjpaw_psl.0.2.UPF
 As  74.92  As.rel-pbe-n-kjpaw_psl.0.2.UPF 
ATOMIC_POSITIONS
 Ga 0.00 0.00 0.00
 As 0.25 0.25 0.25
K_POINTS {automatic}
6 6 6 1 1 1
-------------- next part --------------
 &control
    calculation='bands'
    pseudo_dir = '/home/mito/espresso-5.0.2/pseudo',
    outdir='./temp',
    prefix='gaas'
    verbosity='high'
 /
 &system
    ibrav=  2, celldm(1) =10.6827, nat=  2, ntyp= 2,
    ecutwfc =80.0,
    nbnd = 30, 
    noncolin = .true.,
    lspinorb=.true.,
    starting_magnetization(1)=0.0,
    input_dft = 'PBE'
 /
 &electrons
    diagonalization='david'
 /
ATOMIC_SPECIES
 Ga  69.72  Ga.rel-pbe-dn-kjpaw_psl.0.2.UPF
 As  74.92  As.rel-pbe-n-kjpaw_psl.0.2.UPF 
ATOMIC_POSITIONS
 Ga 0.00 0.00 0.00
 As 0.25 0.25 0.25
K_POINTS
51
  0.500000  0.500000  0.500000 1.0
  0.480000  0.480000  0.480000 1.0
  0.460000  0.460000  0.460000 1.0
  0.440000  0.440000  0.440000 1.0
  0.420000  0.420000  0.420000 1.0
  0.400000  0.400000  0.400000 1.0
  0.380000  0.380000  0.380000 1.0
  0.360000  0.360000  0.360000 1.0
  0.340000  0.340000  0.340000 1.0
  0.320000  0.320000  0.320000 1.0
  0.300000  0.300000  0.300000 1.0
  0.280000  0.280000  0.280000 1.0
  0.260000  0.260000  0.260000 1.0
  0.240000  0.240000  0.240000 1.0
  0.220000  0.220000  0.220000 1.0
  0.200000  0.200000  0.200000 1.0
  0.180000  0.180000  0.180000 1.0
  0.160000  0.160000  0.160000 1.0
  0.140000  0.140000  0.140000 1.0
  0.120000  0.120000  0.120000 1.0
  0.100000  0.100000  0.100000 1.0
  0.080000  0.080000  0.080000 1.0
  0.060000  0.060000  0.060000 1.0
  0.040000  0.040000  0.040000 1.0
  0.020000  0.020000  0.020000 1.0
  0.000000  0.000000  0.000000 1.0
  0.020000  0.000000  0.020000 1.0
  0.040000  0.000000  0.040000 1.0
  0.060000  0.000000  0.060000 1.0
  0.080000  0.000000  0.080000 1.0
  0.100000  0.000000  0.100000 1.0
  0.120000  0.000000  0.120000 1.0
  0.140000  0.000000  0.140000 1.0
  0.160000  0.000000  0.160000 1.0
  0.180000  0.000000  0.180000 1.0
  0.200000  0.000000  0.200000 1.0
  0.220000  0.000000  0.220000 1.0
  0.240000  0.000000  0.240000 1.0
  0.260000  0.000000  0.260000 1.0
  0.280000  0.000000  0.280000 1.0
  0.300000  0.000000  0.300000 1.0
  0.320000  0.000000  0.320000 1.0
  0.340000  0.000000  0.340000 1.0
  0.360000  0.000000  0.360000 1.0
  0.380000  0.000000  0.380000 1.0
  0.400000  0.000000  0.400000 1.0
  0.420000  0.000000  0.420000 1.0
  0.440000  0.000000  0.440000 1.0
  0.460000  0.000000  0.460000 1.0
  0.480000  0.000000  0.480000 1.0
  0.500000  0.000000  0.500000 1.0


More information about the users mailing list