Numerical solution of nonlinear matrix equations arising from Green's function calculations in nano research

Chun Hua Guo*, Yueh Cheng Kuo, Wen-Wei Lin

*Corresponding author for this work

Research output: Contribution to journalArticlepeer-review

11 Scopus citations


The Green's function approach for treating quantum transport in nano devices requires the solution of nonlinear matrix equations of the form X+( C*+iη D*)X- 1(C+iηD)= R+iηP, where R and P are Hermitian, P+λ D*+λ- 1D is positive definite for all λ on the unit circle, and η→0 +. For each fixed η>0, we show that the required solution is the unique stabilizing solution . Then X*=limη→0 + Xη is a particular weakly stabilizing solution of the matrix equation X+ C*X- 1C=R. In nano applications, the matrices R and C are dependent on a parameter, which is the system energy E. In practice one is mainly interested in those values of E for which the equation X+ C*X- 1C=R has no stabilizing solutions or, equivalently, the quadratic matrix polynomial P(λ)= λ2C*-λR+C has eigenvalues on the unit circle. We point out that a doubling algorithm can be used to compute efficiently even for very small values of η, thus providing good approximations to X*. We also explain how the solution X* can be computed directly using subspace methods such as the QZ algorithm by determining which unimodular eigenvalues of P(λ) should be included in the computation. In some applications the matrices C,D,R,P have very special sparsity structures. We show how these special structures can be exploited to drastically reduce the complexity of the doubling algorithm for computing .

Original languageEnglish
Pages (from-to)4166-4180
Number of pages15
JournalJournal of Computational and Applied Mathematics
Issue number17
StatePublished - 1 Nov 2012


  • Green's function
  • Nonlinear matrix equation
  • Structure-preserving algorithm
  • Weakly stabilizing solution

Fingerprint Dive into the research topics of 'Numerical solution of nonlinear matrix equations arising from Green's function calculations in nano research'. Together they form a unique fingerprint.

Cite this