This thread is locked.Only browsing is available.
Top Page > Browsing
Bi2Te3 topological insulator
Date: 2014/06/19 19:29
Name: Artem   <artem.pulkin@epfl.ch>

Dear Taisuke,

Exploring the possibilities of the localized basis sets in openmx, I tried to reproduce the band structure of (topological insulator) Bi2Te3. I face the 2 problems now.

1) I cannot reproduce zero gap at gamma point. I tried to increase integration cutoff, k-point mesh but that did not help. I will probably go to a larger supercell but I am not sure if it helps. Any suggestions?

2) Here are the 4 band structures obtained with 4 similar input files. The only difference is the basis function set used for calculation.

Bi8.0 Te7.0 http://postimg.org/image/85qw5lib5/
Bi8.0 Te9.0 http://postimg.org/image/9g1nrvl0n/
Bi10.0 Te7.0 http://postimg.org/image/orbgshicn/
Bi10.0 Te9.0 http://postimg.org/image/dmznd4dmv/

First 3 ones are very similar. The last one is totally different (and wrong). I would expect that the choice of the basis set would not influence calculation that much but it does!

For your convenience, I post here the input file for the wrong band structure. As I said, all of them differ only by the choice of the basis set. All calculations converged to 1e-10 in ~40 iterations.

Yours faithfully,

Artem Pulkin
EPFL

___________________________________________________________

System.CurrrentDirectory ./
System.Name bite_3ql
level.of.stdout 1
level.of.fileout 0
Species.Number 2
<Definition.of.Atomic.Species
Bi Bi10.0-s4p4d3f2 Bi_PBE13
Te Te9.0-s4p3d3f2 Te_PBE13
Definition.of.Atomic.Species>
Atoms.UnitVectors.Unit Ang
<Atoms.UnitVectors
4.386 0.000 0.000
2.193 3.798387421 0.00
0.00000000 0.00000000 50.8283333968
Atoms.UnitVectors>
Atoms.Number 15
Atoms.SpeciesAndCoordinates.Unit Frac
<Atoms.SpeciesAndCoordinates
1 Te 0.000000 0.000000 0.000000 8.0 8.0
2 Bi 0.333333 0.333333 0.034180 7.5 7.5
3 Te 0.666667 0.666667 0.074180 8.0 8.0
4 Bi 0.000000 0.000000 0.114180 7.5 7.5
5 Te 0.333333 0.333333 0.148360 8.0 8.0
6 Te 0.666667 0.666667 0.200000 8.0 8.0
7 Bi 0.000000 0.000000 0.234180 7.5 7.5
8 Te 0.333333 0.333333 0.274180 8.0 8.0
9 Bi 0.666667 0.666667 0.314180 7.5 7.5
10 Te 0.000000 0.000000 0.348360 8.0 8.0
11 Te 0.333333 0.333333 0.400000 8.0 8.0
12 Bi 0.666667 0.666667 0.434180 7.5 7.5
13 Te 0.000000 0.000000 0.474180 8.0 8.0
14 Bi 0.333333 0.333333 0.514180 7.5 7.5
15 Te 0.666667 0.666667 0.548360 8.0 8.0
Atoms.SpeciesAndCoordinates>
scf.XcType GGA-PBE # LDA|LSDA-CA|LSDA-PW|GGA-PBE
scf.SpinPolarization nc # On|Off|NC
scf.ElectronicTemperature 300.0 # default=300 (K)
scf.energycutoff 600.0 # default=150 (Ry)
scf.maxIter 1000 # default=40
scf.EigenvalueSolver band # DC|GDC|Cluster|Band
scf.Kgrid 5 5 1 # means n1 x n2 x n3
scf.Mixing.Type rmm-diisk # Simple|Rmm-Diis|Gr-Pulay|Kerker|Rmm-Diisk
scf.Mixing.History 20 # default=5
scf.criterion 1.0e-10 # default=1.0e-6 (Hartree)
scf.spinorbit.coupling on
MD.Type nomd # Nomd|Opt|NVE|NVT_VS|NVT_NH
Band.dispersion on
Band.Nkpath 3
<Band.kpath
100 0.3333333333 0.6666666667 0.0000000000 0.0 0.0 0.0 M G
100 0.0 0.0 0.0 0.5 0.5 0.0 G K
50 0.5 0.5 0.0 0.3333333333 0.6666666667 0.0 K M
Band.kpath>
メンテ
Page: [1]

Re: Bi2Te3 topological insulator ( No.1 )
Date: 2014/06/19 22:33
Name: T. Ozaki

Hi,

Your first three results may be correct, while the result with "Bi10.0 Te9.0" is
due to overcompleteness.

Similar calculations by OpenMX have been reported in
http://www.pnas.org/content/108/1/24.full
http://journals.aps.org/prl/abstract/10.1103/PhysRevLett.110.016403

They, who published the first paper, used 40 quintuple-layers to avoid
the surface and surface interaction.

Your system seems to be too thin compared to their system, which may be
the reason why you didn't get what you expect.

Regards,

TO
メンテ
Re: Bi2Te3 topological insulator ( No.2 )
Date: 2014/06/20 18:03
Name: Artem  <artem.pulkin@epfl.ch>

Dear Taisuke,

Thank you for the references. I understand the problem with interacting surface states. In Quantum Espresso, however, it is enough to have the unit cell as large as I demonstrated: the gap is of order of 1 meV. Instead, with the localized basis set I obtain ~100 meV. I wonder if I can obtain desired result in the same unit cell by introducing ghost atoms instead of increasing the cell size?

As for overcompleteness, I do not get it: the basis size is the same, isn't it? We just change matrix elements and obtain something very different. Or you meant that the initial basis set is too large? Then why does this problem reveal itself only for this combination of basis functions?

Thank you in advance,

Yours faithfully,

Artem
メンテ
Re: Bi2Te3 topological insulator ( No.3 )
Date: 2014/06/21 10:11
Name: T. Ozaki

Hi,

> In Quantum Espresso, however, it is enough to have the unit cell as large as
> I demonstrated: the gap is of order of 1 meV.

A VASP calculation shows a band gap of 148 and 36 meV at Gamma point for 2 and 3 QLs, respectively,
published in https://journals.aps.org/prb/abstract/10.1103/PhysRevB.81.041307

There seems to be a large difference even among PW codes.

> Instead, with the localized basis set I obtain ~100 meV. I wonder if I can obtain
> desired result in the same unit cell by introducing ghost atoms instead of increasing the
> cell size?

It is easy to see how the band gap at Gamma point is affected by introducing ghost atoms.
Why don't you try it?

> As for overcompleteness, I do not get it: the basis size is the same, isn't it? We just
> change matrix elements and obtain something very different. Or you meant that the initial
> basis set is too large? Then why does this problem reveal itself only for this combination
> of basis functions?

The overcompleteness means a condition that some of the eigenvalues of k-dependent overlap matrix
become negative. You changed cutoff radii of the basis functions, thus leading to change of the
spectrum of overlap matrix. The setting you made for basis functions must be close to the condition,
since your basis functions look very rich.

Regards,

TO
メンテ

Page: [1]