Computing the Enclosures Eigenvalues Using the Quadratic Method

In this article, We compute the enclosures eigenvalues (upper and lower bounds) using the quadratic method. The Schrodinger operator (A) (harmonic and anharmonic oscillator model) has used as an example. We study a new technique to get more accurate bounds. We compare our results with Boulton and Strauss method.


Introduction
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 [1] [2], and Boulton & Strauss method in [3] [4].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 [5].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 [6] [7].It was then suggested by Shargrodsky and subsequently by Levitin and Shargorodsky (2000) in [8] that the second order relative spectra can also be employed for the pollu-S.M. Abusamra DOI: 10.4236/am.2019.104016213 Applied Mathematics tion-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 [9] and then by Bolton & Strauss (2007Strauss ( , 2011) in [3] [4].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 ( 1 2 3 4 5 , , , , λ λ λ λ λ ) (i.e.; around λ 1 we will choose 1 1 In this paper we studied two models (harmonic and anharmonic oscillator) which are: 1):

Notation
• Below  denotes a generic separable Hilbert space with inner product .and norm . .

Let the operator ( ( )
 ) be self-adjoint.We will write

( )
Spec  to denote the spectrum of (  ).
•  is real inner product space means that  is real vector space on which there is an inner product , x y associating a real number to each pair of elements x, y.The norm is real function such that • For any pair of elements x, y of  satisfies the following properties: • 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: - .,

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 [6] [7], it was then suggested by Shargrodsky and subsequently by Levitin and Shargorodsky (2000) in [8], Properties of second order relative spectra have been studied recently by Boulton & Levitin in [9] and then by Boluton & Strauss (2007Strauss ( , 2011) ) in [3] [4].For the benefit of the reader, we include here some definitions.
Let  be self-adjoint operator on Hilbert space  .Let ( ) number z ∈  belongs to second order spectrum ( ) For all u is nonzero.Means that the Second Order Spectrum usually contain complex numbers, but it turns out that if , , n e e e =  .
Be a basis of  [  is finite dimensional subspace of ( ) Dom  ], 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 z ∈  such that The standard way of finding Spec(Q) is to reach to the linear pencil eigenvalue problem: 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 2 1 0 2 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: (2.5) Indeed, the assertion that Q(z) is singular is equivalent to the existence of 0 u ≠ 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: be given: ( ) ( ) Then ( ) z  is an upper bound for the distance from z to the spectrum of  , that means: (2.12) Assume that z not in ( ) Therefore ( ) x  can be small only when x ∈  is close to ( ) Spec  , so we can make a connection between ( ) Spec  and ( ) Suppose [a, b] be an interval of the spectrum of  which lies in it, and , , , , , Thus the zeros of G(z) appear in conjugate pair.( ) Then either β = 0, u ∈  , u = v we get:

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 [2].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 [3] [4].The spectrum of  ,

( )
Spec  , may be expressed as the union of the discrete spectrum consisting of all isolated eigenvalues of finite multiplicity, ( ) disc Spec  , and the essential spectrum, where.

( ) ( ) ( )
In most standard situations the essential spectrum can be found analytically, but points in ( ) disc Spec  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 Strauss method depends on sub-interval around the eigenvalues in ( ) Spec  to find more accurate approximation enclosure eigenvalues than the Quadratic method around λ. Let: By corresponding limit sets which has been stimulated by the following prop-

New Technique to Get More Accurate Results
we described a new technique to get more accurate bounds for the eigenvalues in ( ) Spec  .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 know the accuracy and efficiency for this method in [5].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 Our improvement depends on domain expansion around the first five eigenvalues in

( )
Spec  , 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 1 1 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.

Schrodinger Operator
Consider the trial subspace  constructed via the finite element method on fi- nite segment.Let 0 >  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: ( ) , har har Spec L  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.6) in all cases as n increases (Figure 3).2, Figure 5).

Conclusions
We compute ( )  , anh anh Spec L  for value n = 200, clearly in both the second-order relative spectra is not the same.
In Table 1 and Table 2 we  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.
is called triangular inequality.

( ) 2 ,
Spec   are closed to the real line, so the conjugate pairs will give small intervals (a, b)s enclosing points in ( ) 2 Spec  .After that Boulton in [9] [10] [11], was modified this studying and gave a value for each a, b within a specific domain for the first five eigenvalues in ( )

Figure 1 and 1 : 6 har
Figure 1 and Figure 2 show the conjugate pair for each eigenvalue with upper and lower bounds of eigenvalues in

Figure 1 .
Figure 1.Second order spectra relative to L har .The horizontal axis is the real part of the points in

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 6

Figure 4 .
Figure 4. Second order spectra relative to 6 anh L .The horizontal axis is the real part of the

Figure 5 .
Figure 5. Approximating enclosures for the first five eigenvalues of 6 anh  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.

1  and 2 
described before and calculate the eigenvalues enclosure.All the coefficient of the matrices 0  , were found analytically.

Figure 1
Figure 1 and Figure 4 show show the approximation of the first five eigenvan = 200 by three different methods (Quadratic method, Boulton & Strauss method, and Our improvement)

Table 1 .
Approximating enclosures for the first five eigenvalues of 6 low is the lower bound of the enclosing λ and λ upp is the upper bound.

Table 2 .
Approximating enclosures for the first five eigenvalues of 6 low is the lower bound of the enclosing λ and λ upp is the upper bound.