標題: Numerical solution of nonlinear matrix equations arising from Green's function calculations in nano research 作者: Guo, Chun-HuaKuo, Yueh-ChengLin, Wen-Wei應用數學系Department of Applied Mathematics 關鍵字: Nonlinear matrix equation;Weakly stabilizing solution;Structure-preserving algorithm;Green's function 公開日期: 1-Nov-2012 摘要: 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 eta D*)X-1(C + i eta D) = R+i eta P. where R and P are Hermitian, P + lambda D* + lambda D-1 is positive definite for all lambda on the unit circle, and eta -> 0(+). For each fixed eta > 0, we show that the required solution is the unique stabilizing solution X-eta. Then X-center dot = lim(eta -> 0+) X-eta 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 g for which the equation X + C*X-1C = R has no stabilizing solutions or, equivalently, the quadratic matrix polynomial P(lambda) = lambda C-2* - lambda R + C has eigenvalues on the unit circle. We point out that a doubling algorithm can be used to compute X-eta efficiently even for very small values eta, 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(lambda) 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 X-eta. (C) 2012 Elsevier B.V. All rights reserved. URI: http://dx.doi.org/10.1016/j.cam.2012.05.012http://hdl.handle.net/11536/16801 ISSN: 0377-0427 DOI: 10.1016/j.cam.2012.05.012 期刊: JOURNAL OF COMPUTATIONAL AND APPLIED MATHEMATICS Volume: 236 Issue: 17 起始頁: 4166 結束頁: 4180 Appears in Collections: Articles

Files in This Item:

If it is a zip file, please download the file and unzip it, then open index.html in a browser to view the full text content.