Galerkin method is one of the best methods for determining upper bounds for the eigenvalues of semi-definite operators, unfortunately this method cannot find enclosures eigenvalue. This paper shows how to compute enclosures of the eigenvalues of self-adjoint operators by the Quadratic method. At first, we study the second-order relative spectrum (The Quadratic method) in   , and Boulton & Strauss method in   . These methods have used for computing eigenvalue enclosures (upper and lower bounds) of the eigenvalues of self-adjoin operators. The quadratic method, which relies on calculation of the second-order which is providing, certified a priori intervals of spectral enclosure. Then we study our new technique which gives more accurate results, we also follow the results that have been published by Boulton & Hobiny in  . The method will be examined by harmonic and anharmonic oscillator models.
Second-order relative spectra were first considered by Davies (1998) in the context of resonances for general self-adjoint operators in   . It was then suggested by Shargrodsky and subsequently by Levitin and Shargorodsky (2000) in  that the second order relative spectra can also be employed for the pollution-free computation of eigenvalues in gaps of the essential spectrum. Various implementations, including on models from elasticity, solid state, physics, relativistic quantum mechanics and magneto hydrodynamics confirm that the Quadratic method is a reliable tool for eigenvalue approximation in the spectral pollution regime. Properties of second order relative spectra have been studied recently by Bolton & Leviton in  and then by Bolton & Strauss (2007, 2011) in   . Boulton and Strauss extended this method to normal operators and optimal convergence rates for eigenvalues and estimated that by an order of magnitude for the harmonic & anharmonic oscillator models, by cut the interval into sub-interval around the eigenvalues in the Spectrum, and the approximation enclosure eigenvalues results of these models are more accurate than the Quadratic method around λ.
Our improvement depends on domain expansion around the first five eigenvalues in the spectrum, we will take a value for a less than λ1 and a value for b greater than λ5 and calculate the conjugate pairs of eigenvalues ( ) (i.e.; around λ1 we will choose ). In this paper we studied two models (harmonic and anharmonic oscillator) which are:
・ Below denotes a generic separable Hilbert space with inner product and norm .
Let the operator ( ) be self-adjoint. We will write to denote the spectrum of ( ).
・ is real inner product space means that is real vector space on which there is an inner product associating a real number to each pair of elements x, y. The norm is real function such that , and the distance is: . And: ,which is called triangular inequality.
・ For any pair of elements x, y of satisfies the following properties:
・ Given a subspaces of dimension n such that:
we will write
・ The discrete spectrum: is the set of eigenvalues of finite multiplicity.
・ The essential spectrum of : is the remaining part of .
・ Theorem: Let be a self-adjoint operator:
- ,there exists a sequence such that:
- as [This sequence is called Welys sequence].
- is not Fredholm.
・ The min-max theorem tells us if is bounded above and below by the mini-mum and maximum eigenvalues respectively ,and we have: ,where denoted by Rayleigh-Ritz quotient on a self-adjoint operator as:
2. The Quadratic Method
We begin by describing the basic framework of the Quadratic method associated to a self-adjoint operator. First considered by Davis (1998) in the context of resonances for general self-adjoint operator in   , it was then suggested by Shargrodsky and subsequently by Levitin and Shargorodsky (2000) in  , Properties of second order relative spectra have been studied recently by Boulton & Levitin in  and then by Boluton & Strauss (2007, 2011) in   . For the benefit of the reader, we include here some definitions.
Let be self-adjoint operator on Hilbert space . Let ,a number belongs to second order spectrum if there exists such that;
For all u is nonzero. Means that the Second Order Spectrum usually contain complex numbers, but it turns out that if then: .
Be a basis of [ is finite dimensional subspace of ], then the Quadratic matrix polynomial is:
where the mass, stiffness and bending matrices are;
We define the spectrum of Q(z) as the set of such that
for some . (2.3)
The standard way of finding Spec(Q) is to reach to the linear pencil eigenvalue problem:
for some . (2.4)
In this method we need to construct companion matrices which depend on the quadratic matrix polynomial to find the enclosure eigenvalues.
We have two possible companion matrices form which are given by
The eigenvalues of the matrix polynomial can be determined from one of this companion matrices form
LEMMA 2. 1: Let Q(z) is defined and be singular then:
Indeed, the assertion that Q(z) is singular is equivalent to the existence of such that
Denoting v = zu, this can be rewritten as
In turn, the latter is equivalent to:
As needed for the verification of (2.3).
LEMMA 2.2: Let and is the orthogonal projection onto ,then:
Typically contains non-real points. From the lemma it is easy to see that if and only if .
We now discuss a strategy suggested by Davis and Plum for computing the spectrum of self-adjoint operator.
Let and for to ,consider the function be given:
Then is an upper bound for the distance from z to the spectrum of ,that means:
So . (2.12)
Assume that z not in ,since then:
Therefore can be small only when is close to ,so we can make a connection between and .
Suppose [a, b] be an interval of the spectrum of which lies in it, and
Clearly, . For ,let Q(z) be as in (2.1),
Thus the zeros of G(z) appear in conjugate pair.
THEOREM 2.1: (Shargorodsky):
For ,and ,let
Let for ,then:
Then either β = 0, ,u = v we get:
So: and (2.16)
3. The Boulton and Strauss Method
The method of second order relative spectra has been shown to reliably approximate the discrete spectrum for a self-adjoint operator in  . Boulton and Strauss extended this method to normal operators and find optimal convergence rates for eigenvalues and estimated that by an order of magnitude in   . The spectrum of , ,may be expressed as the union of the discrete spectrum consisting of all isolated eigenvalues of finite multiplicity, ,and the essential spectrum, where.
In most standard situations the essential spectrum can be found analytically, but points in are usually estimated by numerical procedures. The standard numerical techniques, such as the Quadratic method, aim at solving Galerkin approximate problems posed in weak form: find and such that ( is finite dimensional). Boulton and Strauss method depends on sub-interval around the eigenvalues in to find more accurate approximation enclosure eigenvalues than the Quadratic method around λ. Let:
for any sequence ( ).
By corresponding limit sets which has been stimulated by the following property: if ,then:
Thus inclusions of points in the spectrum of are achieved from Re(z) with a two-sided explicit residual given by |Im(z)|. Indeed, if and with ,then:
4. New Technique to Get More Accurate Results
we described a new technique to get more accurate bounds for the eigenvalues in . Hobiny was discussed the Quadratic method to estimate enclosure eigenvalues for self-adjoint operator and illustrated the results on the harmonic and anharmonic oscillator models to know the accuracy and efficiency for this method in  . We discussed also the bounds of the size of the enclosure eigenvalues and studied them in the context of one dimensional Schrodinger operators, and illustrated that the conjugate pairs of the eigenvalues in are closed to the real line, so the conjugate pairs will give small intervals (a, b)s enclosing points in .
After that Boulton in    , was modified this studying and gave a value for each a, b within a specific domain for the first five eigenvalues in ,(i.e.: , , , , ).
Our improvement depends on domain expansion around the first five eigenvalues in ,we will take a value for a less than λ1 and a value for b greater than λ5 and calculate the conjugate pairs of eigenvalues (λ1, λ2, λ3, λ4, λ5) (i.e.; around λ1 we will choose where:
And the same technique for other eigenvalues. Then we will apply a, b in our programs to find the first five enclosure eigenvalues and produce our comparison by numerical experiments on the harmonic and anharmonic oscillator models.
5. Schrodinger Operator
Consider the trial subspace constructed via the finite element method on finite segment.
Let and set then, the general Schrodinger equation is:
where V(x) is called the potential function, and it must be bounded below ; which is called the Dirichlet Boundary Conditions.
If we take u(x) as a common factor then:
so the operator is a differential equation.
then the equation is related to harmonic oscillator model.
then the equation is related to anharmonic oscillator model.
We compare between the Quadratic method, Boulton & Strauss and our development on these models to calculate the enclosure eigenvalues to know which one is the best in this field.
6. Harmonic Oscillator Model
The Harmonic Oscillator is one of the most important models of quantum theory.
Let ,for ,then the exact eigenvalue is ; .
By Schrodinger equation we have:
This equation can be solved explicitly and we can find the approximation eigenvalues using Matlab program by matrices M, N, R where,
Now we compute as described before and calculate the eigenvalues enclosure. All the coefficients of the matrices were found analytically.
NOTE: In our improvement method we choose (n = 200, L = 6) to compare our results with Boulton and Strauss method, which approximate the first five eigenvalues of harmonic and anharmonic models with n = 200, L = 6.
Trial 1: In this trial we use the Quadratic method, Boulton & Strauss method, and our improvement with: n = 200, L = 6, to find the first five approximation
We calculate the Error between the exact and the lower bound for the first five approximation enclosure eigenvalues of . . The slope of the graphs is close to the value (6) in all cases as n increases (Figure 3).
Table 1. Approximating enclosures for the first five eigenvalues of with n = 200, λlow is the lower bound of the enclosing λ and λupp is the upper bound.
Figure 1. Second order spectra relative to Lhar. The horizontal axis is the real part of the points in and the vertical axis is imaginary part.
Figure 2. and illustration on the end pints of the segment given in Theorem 2.1, and .
Figure 3. Approximating enclosures for the first five eigenvalues of by the: Quadratic, Boulton & Strauss, and Our improvement methods with n = 200, λlow is the lower bound of the enclosing λ and λupp is the upper bound. (a) This sub-figure shows the error between the exact and lower bound for the first five approximation eigenvalues of using the quadratic method; (b) This sub-figure shows the error between the exact and lower bound for the first five approximation eigenvalues of using the Boulton & Strauss method; (c) This sub-figure shows the error between the exact and lower bound for the first five approximation eigenvalues of using our improvement method.
7. Anharmonic Oscillator Model
The anharmonic oscillator is another model of Schrodinger equation in one dimension. This model is one of the most important problems of quantum mechanics.
We choose this model to show that the technique for both methods can be applied to operator where the exact spectrum is not known.
Let ,for then; by Schrodinger equation we have
The exact eigenvalue is unknown.
This equation can be solved explicitly and we can find the approximation eigenvalues using Matlab program by matrices M, N, R where:
We compute ,as in Figure 4.
Trial 2: In this trial we use the Quadratic method, Boulton & Strauss method, and Our improvement technique with: n = 200, L = 6, to find the first five approximation enclosure eigenvalues (upper and lower eigenvalues) of ,
Table 2. Approximating enclosures for the first five eigenvalues of with n = 200, λlow is the lower bound of the enclosing λ and λupp is the upper bound.
Figure 4. Second order spectra relative to . The horizontal axis is the real part of the points in and the vertical axis is imaginary part.
Figure 5. Approximating enclosures for the first five eigenvalues of by the Quadratic, Boulton & Strauss, and Our improvement methods with n = 200, λlow is the lower bound of the enclosing λ and λupp is the upper bound.
We compute and ,as we described before and calculate the eigenvalues enclosure. All the coefficient of the matrices , and were found analytically.
In Table 1 and Table 2 we show the approximation of the first five eigenvalues enclosures of , respectively, with n = 200 by three different methods (Quadratic method, Boulton & Strauss method, and Our improvement) to compare between these methods and identify which one is more accurate for computing eigenvalues.
For the Harmonic oscillator model, the error between the approximation results and the first five exact eigenvalues (1, 3, 5, 7, 9) in our technique is less than the error between the exact and the approximation eigenvalues by the Boulton and Strauss method, and the same thing between of the Quadratic method results and the exact eigenvalues, so the approximation enclosure eigenvalues by our technique is more accurate and effective than the Boulton & Strauss method and the Quadratic method, also Figure 3 shows that clearly.
For the anharmonic oscillator model, the exact eigenvalues are unknown but by results and Figure 5, we can confirm the previous result.
 Boulton, L. and Strauss, M. (2011) On the Convergence of second order spectra and Multiplicity. Proceedings of the Royal Society A Mathematical Physical and Engineering Sciences, 467, 264-284.
 Levitin, M. and Shargorodsky, E. (2004) Spectral Pollution and second order relative Spectra for Self-Adjoint Operators. IMA Journal of Numerical Analysis, 24, 393-416.
 Boulton, L. and Levitin, M. (2007) On Approximation of the Eigenvalue of Perturbed Periodic Schrodinger Operators. Journal of Physics A Mathematical and Theoretical, 40, 9319.
 Boulton, L., Garcia del Moral, M.P. and Restuccia, A. (2012) Spectral Properties in Super Symmetric Matrix Models. Nuclear Physics B, 856, 716-747.