This thread is locked.Only browsing is available.
Top Page > Browsing
Projected band structure
Date: 2016/08/17 20:46
Name: Riemann   <riemann.derakhshan@gmail.com>

Dear Prof.Ozaki

I have a question about the projection of band structure with OpenMX.

My question is : The projection and extraction of the band structure of a 2D layer supported

on the substrate (same as -----> Y. Cai, C. P. Chuu, C. M. Wei, and M. Y. Chou, “Stability

and electronic properties of two-dimensional silicene and germanene on graphene,” Phys. Rev.

B - Condens. Matter Mater. Phys., vol. 88, no. 24, pp. 2&#8211;7, 2013.) can be done with Openmx

package?

I would be grateful if You give me any guidance about it.


Sincerely Yours
Riemann

メンテ
Page: [1]

Re: Projected band structure ( No.1 )
Date: 2016/08/18 08:56
Name: T. Ozaki

Hi,

Yes, you can perform the same analysis using OpenMX.
How to do that is found in the manual.

Regards,

TO
メンテ
Re: Projected band structure ( No.2 )
Date: 2016/08/18 13:48
Name: Riemann  <riemann.derakhshan@gmail.com>

Dear Prof.Ozaki

Many thanks for your prompt reply.

I've searched in the manual and within section Interface with Wannier90 I've found the

phrase " OpenMX interfaces with Wannier90 [99] which constructs maximally localized Wannier

functions, and calculates physical properties such as Wannier projected DOS and

bandstructure,".

Did You mean that? should I do the band structure projection with Wannier function?


Sincerely Yours
Riemann
メンテ
Re: Projected band structure ( No.3 )
Date: 2016/08/18 22:00
Name: T. Ozaki

Hi,

Since the bands can be decomposed to pseudo-atomic orbital contributions as shown in

http://www.openmx-square.org/openmx_man3.8/node154.html
http://www.openmx-square.org/openmx_man3.8/node155.html

you only have to collect the contributions associated with pseudo-atomic orbitals in the
two-dimensional layer. Then, it is easy to see how the original band structure of the
free-standing two-dimensional layer is perturbed by the interaction with the substrate.
The analysis can be performed in both the unfolding and non-unfolding cases, where the
non-unfolding means that the unit cell used for the SCF calculation is used in the analysis
as the 'reference' cell.

Regards,

TO
メンテ
Re: Projected band structure ( No.4 )
Date: 2016/08/19 04:18
Name: Riemann  <riemann.derakhshan@gmail.com>

Dear Prof.Ozaki

Many thanks for Your kind and prompt reply.

I've done the band structure analysis for Bilayer silicene (as toy example) to find out

which pseudo-atomic orbitals are corresponding to which atom in structure. Herewith I've

attached the input file for Bilayer silicene. As can be seen from input file I've divided

the Silicon atom to two types. First, the (Si1) atoms belong to the lower layer and second,

(Si2) to the upper layer. After finishing the run, when I've seen the "silicene.unfold_orb"

file, it contains 784 rows and 54 columns. Now the question is this: Which row and columns

is belong to which pseudo-atomic orbitals and corresponds to which atoms in the structure?




#
# File Name
#

System.CurrrentDirectory ./ # default=./
System.Name silicene
level.of.stdout 1 # default=1 (1-3)
level.of.fileout 1 # default=1 (0-2)
DATA.PATH /home/riemann/DFT/openmx3.7/DFT_DATA13
#
# Definition of Atomic Species
#

Species.Number 2
<Definition.of.Atomic.Species
Si1 Si7.0-s2p2d1 Si_PBE13
Si2 Si7.0-s2p2d1 Si_PBE13
Definition.of.Atomic.Species>

#
# Atoms
#

Atoms.Number 4
Atoms.SpeciesAndCoordinates.Unit FRAC # Ang|AU
<Atoms.SpeciesAndCoordinates # Unit=Ang.
1 Si1 0.3333333333300033 0.6666666666669983 0.5000000000000000 2 2
2 Si1 0.6666666666699967 0.3333333333330017 0.4668390626705445 2 2
3 Si2 0.0000000000000000 0.0000000000000000 0.6592873730830462 2 2
4 Si2 0.3333333333300033 0.6666666666669983 0.6261262380352548 2 2
Atoms.SpeciesAndCoordinates>
Atoms.UnitVectors.Unit Ang # Ang|AU
<Atoms.UnitVectors # unit=Ang.
3.342858059 1.930000000 0.000000000
-3.342858059 1.930000000 0.000000000
0.000000000 0.000000000 20.000000000
Atoms.UnitVectors>

#
# SCF or Electronic System
#

scf.XcType GGA-PBE # LDA|LSDA-CA|LSDA-PW|GGA-PBE
scf.SpinPolarization nc # On|Off|NC
scf.SpinOrbit.Coupling on
scf.ElectronicTemperature 300.0 # default=300 (K)
scf.energycutoff 200.0 # default=150 (Ry)
scf.maxIter 200 # default=40
scf.EigenvalueSolver Band # DC|GDC|Cluster|Band
scf.Kgrid 12 12 3 # means 4x4x4
scf.Mixing.Type Rmm-Diisk # Simple|Rmm-Diis|Gr-Pulay|Kerker|Rmm-Diisk
scf.Init.Mixing.Weight 0.30 # default=0.30
scf.Min.Mixing.Weight 0.001 # default=0.001
scf.Max.Mixing.Weight 0.400 # default=0.40
scf.Mixing.History 15 # default=5
scf.Mixing.StartPulay 8 # default=6
scf.criterion 1.0e-6 # default=1.0e-6 (Hartree)
scf.lapack.dste dstevx # dstegr|dstedc|dstevx, default=dstevx


#
# Band structure
#

Band.dispersion on # on|off, default=off
## if <Band.KPath.UnitCell does not exist,
## the reciprical lattice vector is employed.
Band.Nkpath 3
<Band.kpath
200 0.33333333333 0.33333333333 0.0000000000 0.00000000000 0.00000000000 0.0000000000 K G
200 0.00000000000 0.00000000000 0.0000000000 0.50000000000 0.00000000000 0.0000000000 G M
200 0.50000000000 0.00000000000 0.0000000000 0.33333333333 0.33333333333 0.0000000000 M K
Band.kpath>

#
# DOS and PDOS
#

Dos.fileout on # on|off, default=off
Dos.Erange -20.0 20.0 # default = -20 (eV) 20 (eV)
Dos.Kgrid 30 30 1 # default = Kgrid1 Kgrid2 Kgrid3



#
# Band unfolding
#

Unfolding.Electronic.Band on # on|off, default=off
Unfolding.LowerBound -10.0 # default=-10 eV
Unfolding.UpperBound 6.0 # default= 10 eV

Unfolding.desired_totalnkpt 30
Unfolding.Nkpoint 4

<Unfolding.kpoint
K 0.33333333333 0.33333333333 0.0000000000
G 0.00000000000 0.00000000000 0.0000000000
M 0.50000000000 0.00000000000 0.0000000000
K 0.33333333333 0.33333333333 0.0000000000
Unfolding.kpoint>

メンテ
Re: Projected band structure ( No.5 )
Date: 2016/08/19 15:15
Name: T. Ozaki

Hi,

The index is the same as for orbitally decomposed Mulliken populations and
LCAO coefficients found at

http://www.openmx-square.org/openmx_man3.8/node19.html
http://www.openmx-square.org/openmx_man3.8/node93.html


Also, in Fig. 50 at
http://www.openmx-square.org/openmx_man3.8/node154.html

you will find a similar analysis.

Hope this helps you.

Regards,

TO
メンテ
Re: Projected band structure ( No.6 )
Date: 2016/08/19 15:45
Name: Riemann  <riemann.derakhshan@gmail.com>

Dear Prof.Ozaki

Many thanks for Your reply.

It helped me a lot and I've found out which pseudo-atomic orbitals are corresponding to

which atom in the structure in 'System.Name.out' file.

Sincerely Yours
Riemann
メンテ

Page: [1]