This thread is locked.Only browsing is available.
Top Page > Browsing
Total energy in NEGF
Date: 2015/01/14 20:30
Name: Artem Pulkin   <artem.pulkin@epfl.ch>

Is the total energy calculated by NEGF OpenMX formally correct? Is it possible to use it to compare formation energies of a scattering region?

I ask this question because caluclation of a total energy of a scattering region obviously involves some changes compared to the periodic boundary conditions case. However it seems like the Total_Energy.c does it in a uniform way (regardless of a 'Solver' variable).
メンテ
Page: [1]

Re: Total energy in NEGF ( No.1 )
Date: 2015/01/19 17:44
Name: T. Ozaki

Hi,

The total energy appearing in the NEGF calculations in OpenMX is calculated based on the
formulation used for the conventional cluster or band methods. Though the theoretical basis
is obscure, one might be able to use the energy to discuss stability of the scattering region
of system from a practical sense. Even the structural optimization might be possible by fixing
position of atoms around the boundaries, as several groups do such calculations, while I have
never checked the issue carefully, however.

Regards,

TO
メンテ
Re: Total energy in NEGF ( No.2 )
Date: 2015/01/26 18:27
Name: Artem Pulkin  <artem.pulkin@epfl.ch>

Dear Taisuke,

I am working with polar MoS2 and I want to compare the boundary/edge formation energies. Conventionally it is done using 0D structures, see for example here:

http://arxiv.org/pdf/1101.4457.pdf

This approach has a number of drawbacks, mainly due to an unbalanced number of cations and anions in a cluster. I was interested in calculating this kind of energy using NEGF. For example, I create a following symmetric NEGF geometry with identical edges:

| + - | + - | + - | vacuum | - + | - + | - + |

or

| - + | - + | - + | vacuum | + - | + - | + - |

The question is: does the total energy calculated by NEGF includes contribution from a spurious defect coming from connected left and right leads in periodic geometry?

I would like to emphasize that this is not a question about comparing forces or stability. This is a question about comparing 2 different(!) NEGF calculations with right and left leads swapped.

Regards,

Artem
メンテ
Re: Total energy in NEGF ( No.3 )
Date: 2015/01/28 12:02
Name: T. Ozaki

Hi,

> The question is: does the total energy calculated by NEGF includes contribution from
> a spurious defect coming from connected left and right leads in periodic geometry?

Do you mean the Coulombic interaction between defects in ......| + - | vacuum | - + |.......
by "contribution from a spurious defect"?

If so, the interaction is actually included in the total energy calculated by NEGF through
the Poisson equation.

Regards,

TO
メンテ
Re: Total energy in NEGF ( No.4 )
Date: 2015/01/28 18:11
Name: Artem Pulkin  <artem.pulkin@epfl.ch>

Hi,

I was not clear enough: I mean that there may be some component total energy which is calculated in Fourier space where the scattering region is periodic:

| + - | + - | + - | vacuum | - + | - + | - + | <- spurious interaction -> | + - | + - | + - | vacuum | - + | - + | - + |

I.e. the spurious defect in this case is

| - + | + - |

Of course, I am not a specialist here that's why I ask you if it is the case in OpenMX. If not - good, I will continue to study whether this approach is feasible for polar materials.

Regards,

Artem
メンテ
Re: Total energy in NEGF ( No.5 )
Date: 2015/01/28 19:20
Name: T. Ozaki

Hi,

> | + - | + - | + - | vacuum | - + | - + | - + | <- spurious interaction -> | + - | + - | + - | vacuum | - + | - + | - + |
>
> I.e. the spurious defect in this case is
>
> | - + | + - |


The system treated by NEGF is not periodic along the transport direction, but semi-infinite.
The semi-infinity in evaluating Coulombic interaction is taken into account by solving
the Poisson equation with boundary condition that values of the Hartree potential at
the boundaries between scattering region and leads are fixed at those obtained by
the calculation of the first step. Thus, the interaction "| - + | + - |" is not included
in the NEGF calculation.

Also, the ESM method can be used to eliminate such a spurious interaction
in case of the slab model. The details can be found at
http://www.openmx-square.org/openmx_man3.7/node118.html
and an application can be found in Phys. Rev. B 90, 075422 (2014).

I hope it helps you.

Regards,

TO
メンテ
Re: Total energy in NEGF ( No.6 )
Date: 2015/01/29 18:56
Name: Artem Pulkin  <artem.pulkin@epfl.ch>

Thank you, it is a very relevant information.
メンテ

Page: [1]