16

Two-band model for graphene. To simplify the discussion of Dirac points, it is sufficient to consider a nearest-neighbor tight-binding 2-band model for graphene. This is a spinless model because spin-orbit is negligible in graphene. The Bloch Hamiltonian is: $$ \hat{H}(\mathbf{k})= \begin{pmatrix} 0 & h(\mathbf{k}) \\ h^{\dagger}(\mathbf{k}) & 0 \end{...


10

One way of determining this is using the projected density of states (P-DOS) This resolves the DOS into specific orbitals thereby allowing you to discretize each orbitals weight for a specific energy. $ \mathrm{PDOS}_\nu(E) = \sum_i \psi^*_{i,\nu} [\mathbf S | \psi_{i}\rangle]_\nu D(E-\epsilon_i) $ Note here that $|\psi_i\rangle$ is the $i$th eigenvector and ...


10

Spherical harmonics are not themselves full atomic orbitals. Consider the Hydrogen wave function, which separates into a radial part and an angular part. The latter is a spherical harmonic, but the former is some other function (in the case of Hydrogen it's a Laguerre polynomial). In general, we can approximate the angular part for other atoms with the same ...


9

Oh! I found the answer already. Since I'm doing the calculation along the High symmetry point. Im doing the calculation along G -> K -> M -> G, that's why the right hand sides show such phenomena. In tight binding I plot the band along G -> K -> G. I try to reproduce the band along this direction, it matched with the tight binding solution. Thanks guys. It's ...


8

Yes! It is definitely possible and it is useful for calculating other things like electronic transport and first-principles values of Hubbard U (i.e. the ACBN0 method, which I have used a bit). Some good papers to read are from Prof. Marco Buongiorno Nardelli's group at UNT: L. A. Agapito, A. Ferretti, A. Calzolari, S. Curtarolo, and M. Buongiorno Nardelli, ...


7

Not really sure about benchmarks for solid state TS geometries. A couple of scattered piece of references are: (i) This paper by Ceder's group has several related TS barriers: aip.scitation.org/doi/pdf/10.1063/1.4960790 (ii) Then there is Henkelman's paper for bulk unit cell TS calculations: aip.scitation.org/doi/abs/10.1063/1.3684549. Maybe you are already ...


7

The following code is a python function to read the wannier90_hr.dat, from which you should figure out its data structure. def read_hamiltonian(path): """ Read hopping matrix element from the wannier90 output file:wannier90_hr.dat wan_num: number of wannier functions wsc_num: number of wigner-sitez cells wsc_count: ...


6

You can read the wanniertools code. In wanniertools, to calculate surface state, they write a surfstate subroutine in surfstate.f90 file. The slab Hamiltonian is restructured from a bulk Hamiltonian in ham_qlayer2qlayer.f90 file. enter code here ! This is a fortran code. ! H00 Hamiltonian between nearest neighbour-quintuple-layers ! the factor 2 is ...


5

The answer by Anyon does not address how to calculate the form of the Slater-Koster two-center integrals (2CI) but rather how to determine the direction cosines that appear in the 2CI's (the $l$,$m$,$n$'s in things like $\sqrt{3}lm*(sd\sigma)$, but not the $\sqrt{3}$). I am currently trying to explicitly perform the same calculation to include in my comps so ...


5

Consider the energy eigenstate $|\psi_{n\mathbf{k}}\rangle$, where $n$ is the band you are interested in and $\mathbf{k}$ is wave vector of the TRIM point. Then, if the system has inversion symmetry, these energy eigenstates are also eigenstates of the parity operator $\hat{\pi}$. This means that to determine the parity of that state, you can calculate the ...


4

Wannier90 might not be good at preserving the symmetry. But they probably include a few new methods to enforce symmetry in Wannier90.v.3.1.0. Maybe you can check this. http://www.wannier.org/features/ Also, WannierTools can symmetrize the hr.dat, but from my personal experience it sometimes gives you worse results than the original hr.dat. http://www....


Only top voted, non community-wiki answers of a minimum length are eligible