Exact treatment of generalized modifications of finite-dimensional systems by the LRM approach

LRM (Low Rank Modification) is a mathematical method that produces eigenvalues and eigenstates of generalized eigenvalue equations. It is similar to the perturbation expansion in that it assumes the knowledge of the eigenvalues and eigenstates of some related (unperturbed) system. However, unlike perturbation expansion, LRM produces correct results however large the modification of the original system. LRM of finite-dimensional systems is here generalized to the combined (external and internal) modifications. Parent n-dimensional system An containing n eigenvalues λi and n eigenstates \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${| {\Phi_i}\rangle}$$\end{document} is described by the generalized n × n eigenvalue equation. In an external modification system An interacts with another ρ-dimensional system Bρ which is situated outside the system An. In an internal modification relatively small σ-dimensional subsystem of the parent system An is modified. Modified system Cn+ρ that contains external as well as internal modifications is described by the generalized (n + ρ) × (n + ρ) eigenvalue equation. This system has (n + ρ) eigenvalues \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varepsilon_s}$$\end{document} and (n + ρ) corresponding eigenstates \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${| {\Psi_s}\rangle}$$\end{document} . In LRM this generalized (ρ + n) × (ρ + n) eigenvalue equation is replaced with a (nonlinear) (ρ + σ) × (ρ + σ) equation which produces all eigenvalues \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varepsilon_s \notin \left\{ {\lambda_i}\right\}}$$\end{document} and all the corresponding eigenstates \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${| {\Psi_s}\rangle }$$\end{document} of Cn + ρ. Another equation produces remaining solutions (if any) that satisfy \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varepsilon_s \in \left\{ {\lambda_i}\right\}}$$\end{document} . Those two equations produce exact solution of the modified system Cn + ρ. If (ρ + σ) is small with respect to n, this approach is numerically much more efficient than a standard diagonalization of the original generalized eigenvalue equation. Unlike perturbation expansion, LRM produces exact results, however large modification of the parent system An.