On a nonlinear matrix equation arising in nano research

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

*Corresponding author for this work

Research output: Contribution to journalArticlepeer-review

15 Scopus citations

Abstract

The matrix equation X + A τX -1A = Q arises in Green's function calculations in nano research, where A is a real square matrix and Q is a real symmetric matrix dependent on a parameter and is usually indefinite. In practice one is mainly interested in those values of the parameter for which the matrix equation has no stabilizing solutions. The solution of interest in this case is a special weakly stabilizing complex symmetric solution X*, which is the limit of the unique stabilizing solution X? of the perturbed equation X + A τX -1A = Q + i?I, as ? ? 0+. It has been shown that a doubling algorithm can be used to compute X? efficiently even for very small values of ?, thus providing good approximations to X *. It has been observed by nano scientists that a modified fixed-point method can sometimes be quite useful, particularly for computing X? for many different values of the parameter. We provide a rigorous analysis of this modified fixed-point method and its variant and of their generalizations. We also show that the imaginary part XI of the matrix X* is positive semidefinite and we determine the rank of XI in terms of the number of unimodular eigenvalues of the quadratic pencil λ2A τ - Q + A. Finally we present a new structure-preserving algorithm that is applied directly on the equation X + A τX -1A = Q. In doing so, we work with real arithmetic most of the time.

Original languageEnglish
Pages (from-to)235-262
Number of pages28
JournalSIAM Journal on Matrix Analysis and Applications
Volume33
Issue number1
DOIs
StatePublished - 4 Jun 2012

Keywords

  • Complex symmetric solution
  • Fixed-point iteration
  • Green's function
  • Nonlinear matrix equation
  • Structure-preserving algorithm
  • Weakly stabilizing solution

Fingerprint Dive into the research topics of 'On a nonlinear matrix equation arising in nano research'. Together they form a unique fingerprint.

Cite this