rspert.m

Rayleigh-Schrodinger perturbation theory to arbitrary order, Eqs 2.21-2.23 from Stefan Stoll's PhD thesis, with the typo fixed in the numerator of Eq 2.21.

Syntax

    [Ep,Vp]=rspert(E0,H1,order)

Parameters

    E0    - eigenvalues of H0, a column vector of real 
            numbers

    H1    - perturbation, written in the basis that di-
            agonalises H0

    order - order of perturbation theory to be used, 6
            is the sensible maximum

Outputs

    E     - eigenvalues of H0+H1 to the specified order
            in perturbation theory, a vector of reals 

    V     - normalised eigenvectors of H0+H1 to the spe-
            cified order in perturbation theory, a squa-
            re unitary matrix with eigenvectors in cols
            in the same order as the eigenvalues in E

Examples

Below is the output of examples/fundamentals/perturb_theory.m example file.

Notes

1. There must be no degeneracies in H0.
2. H1 must be Hermitian.
3. The theory only converges when H1 << H0 in 2-norm.
4. Numerical artefacts appear beyond sixth order.
5. Complexity is linear in the order and cubic in matrix dimension.

See also

vvpert.m, acomm.m, arnoldi.m, atranspose.m, aux_mat.m, binpack.m, cheap_norm.m, cheb_coeff.m, clean_up.m, dirdiff.m, eigenfields.m, expdrop.m, expmint.m, expmint2.m, fftdiff.m, fourdif.m, fourlap.m, frob_chop.m, gaussfun.m, hdot.m, herm_spline.m, jacobianest.m, keep_rank.m, krondelta.m, kronm_new.m, logfactorial.m, lorentzcon.m, lorentzfun.m, md5_hash.m, mprealloc.m, remncomm.m, remtrace.m, rspt_eig.m, snormpdf.m, svd_shrink.m, tikhoind.m, tikhonov.m, trapdiff.m, unit_oper.m, unit_state.m, Kernel_utilities

Version 2.6, authors: Ilya Kuprov