Publication Cover
Numerical Heat Transfer, Part B: Fundamentals
An International Journal of Computation and Methodology
Volume 70, 2016 - Issue 3
116
Views
3
CrossRef citations to date
0
Altmetric
Original Articles

An atomistic Green’s function method hybrid with a substructure technique for coherent phonon analysis

, &
Pages 200-214 | Received 07 Jan 2016, Accepted 01 Apr 2016, Published online: 17 Aug 2016
 

ABSTRACT

By virtue of substructure technique, the atomistic Green’s function method can be applied to efficiently evaluate the coherent phonon transport in a large device. Comprised of several substructures without inner atoms, the device can be described by a small system of algebraic equations, which can give a small but crucial submatrix of the retarded Green’s function of the open device without loss in accuracy. While avoiding solving the entire device problem, the method can basically give all the information to describe the quantum transport within the ballistic limit, such as the transmission function and local density of states. The validity and efficiency of the proposed method is demonstrated by the one- and three-dimensional cases. The multilevel substructure technique is also feasible.

Nomenclature

A=

spectral function matrix

a=

lattice constant

D=

effective stiffness matrix

G=

Green’s function matrix

GD1=

first column of the device Green’s function block matrix, defined in Eq. (10)

GD3=

third column of the device Green’s function block matrix, defined in Eq. (10)

g=

uncoupled Green’s function of contacts or the substructure’s inner part

I=

identity matrix

i=

unit of imaginary number

K=

harmonic matrix, defined in Eq. (1)

K(j, j′)=

coupling matrix between substructure j and j

K(LC,LD)=

coupling matrix between left contact and left device region

K(RD,RC)=

coupling matrix between right contact and right device region

k=

element of the harmonic matrix

m=

mass

N=

number of substructures

n=

number of the degrees of freedom

q=

wave number

qmax=

maximum wave number

U=

bonding energy

u=

displacement

V=

coupling harmonic matrix

0+=

broadening constant, which is a positive infinitesimal number

β=

spring constant

Γ=

broadening function, defined in Eq. (7)

Ξ=

transmission function

ρ=

local density of states

Σ=

self-energy matrix

τ=

nonzero submatrix of self-energy matrix, defined in Eq. (5)

ω=

angular frequency of phonons

=

complex angular frequency of phonons

Subscripts=
C=

contact region

D=

device region

ID=

intermediate device region not directly bonding with either contact

I, E=

the inner part and external part

i, j=

index of a degree of freedom

LC=

left contact region

LC0=

0th layer in left contact region

LD=

left device region bonding with the left contact

RC=

right contact region

RC0=

0th layer in right contact region

RD=

right device region bonding with the right contact

Superscripts=
T=

matrix transpose

R=

reduced matrix

(j)=

substructure number

*=

conjugate transpose of a matrix

Nomenclature

A=

spectral function matrix

a=

lattice constant

D=

effective stiffness matrix

G=

Green’s function matrix

GD1=

first column of the device Green’s function block matrix, defined in Eq. (10)

GD3=

third column of the device Green’s function block matrix, defined in Eq. (10)

g=

uncoupled Green’s function of contacts or the substructure’s inner part

I=

identity matrix

i=

unit of imaginary number

K=

harmonic matrix, defined in Eq. (1)

K(j, j′)=

coupling matrix between substructure j and j

K(LC,LD)=

coupling matrix between left contact and left device region

K(RD,RC)=

coupling matrix between right contact and right device region

k=

element of the harmonic matrix

m=

mass

N=

number of substructures

n=

number of the degrees of freedom

q=

wave number

qmax=

maximum wave number

U=

bonding energy

u=

displacement

V=

coupling harmonic matrix

0+=

broadening constant, which is a positive infinitesimal number

β=

spring constant

Γ=

broadening function, defined in Eq. (7)

Ξ=

transmission function

ρ=

local density of states

Σ=

self-energy matrix

τ=

nonzero submatrix of self-energy matrix, defined in Eq. (5)

ω=

angular frequency of phonons

=

complex angular frequency of phonons

Subscripts=
C=

contact region

D=

device region

ID=

intermediate device region not directly bonding with either contact

I, E=

the inner part and external part

i, j=

index of a degree of freedom

LC=

left contact region

LC0=

0th layer in left contact region

LD=

left device region bonding with the left contact

RC=

right contact region

RC0=

0th layer in right contact region

RD=

right device region bonding with the right contact

Superscripts=
T=

matrix transpose

R=

reduced matrix

(j)=

substructure number

*=

conjugate transpose of a matrix

Reprints and Corporate Permissions

Please note: Selecting permissions does not provide access to the full text of the article, please see our help page How do I view content?

To request a reprint or corporate permissions for this article, please click on the relevant link below:

Academic Permissions

Please note: Selecting permissions does not provide access to the full text of the article, please see our help page How do I view content?

Obtain permissions instantly via Rightslink by clicking on the button below:

If you are unable to obtain permissions via Rightslink, please complete and submit this Permissions form. For more information, please visit our Permissions help page.