Evaluated Excited-State Time-Independent Correlation Function and Eigenfunction of the Harmonics Oscillator Cosine Asymmetric Potential via Numerical Shooting Method

We aimed to evaluate the ground-state and excite-state energy eigenvalue (En),wave function, and the time-independent correlation function of the atomic density fluctuation of a particle under the harmonics oscillator Cosine asymmetric potential (Saad et al. 2013). Instead of using the 6-point kernel of 4 Green’s function (Cherroret and Skipetrov, 2008), averaged over disorder, we use the numerical shooting method (NSM) to solve the Schrödinger equation of quantum mechanics system with Cosine asymmetric potential. Since our approach does not use complicated formulas, it requires much less computational effort when compared to the Green functions techniques (Cherroret and Skipetrov, 2008).We show that the idea of the program of evaluating time-independent correlation function of atomic density is underdamped motion for the Cosine asymmetric potential from the numerical shooting method of this problem. Comparison of the time-independent correlation function obtained from numerical shooting method by Boonchui andHutem (2012) and correlation function experiment by Kasprzak et al. (2008).We show the intensity of atomic density fluctuation (δn(x) = n(x) − ?̃?(x)) in harmonics oscillator Cosine asymmetric potential by numerical shooting method.


Introduction
Most problems encountered in quantum mechanics cannot be solved exactly. Exact solutions of the Schrödinger equation exist only for a few idealized systems. To solve general problems, one must resort to approximation methods. A variety of such methods have been developed, and each has its own area of applicability. There exist several means to study them, for example, Wentzel-Kramers-Brillouin [1], perturbation [2], the quasilinearization method [3], the variational method [4], function analysis [5,6], the eigenvalue moment method [7], the analytical transfer matrix method [8][9][10], and numerical shooting method [11,12]. Grobe et al. 1994 [13] proposed a criterion to determine the numerical degree of global correlation function of the multiparities quantum system. They applied this method to several situations, including electron-atom scattering and strong-field photoionization. Lye et al. 2005 [14] discussed the effect of a weak random potential, indicated by stripes in the expanded density profile of the Bose-Einstein Condensate and damped dipole oscillations. Henseler and Shapiro 2008 [15] defined the disorder-induced intensity-intensity correlation function, (r, r ) = | * (r) (r )| 2 , for the Bose-Einstein Condensate for Fermi gas. Cherroret and Skipetrov 2008 [16] showed decay of the average atomic density ( (r, ) = | (r, )| 2 ) as a function of time.  [19] numerically studied the dynamics regimes of classical transport of cold atoms gases in a twodimensional anisotropic disorder potential. In this paper, we consider approximation methods that deal with stationary states corresponding to time-independent Hamiltonian.
To study problem of stationary states, we focus on one approximation method: numerical shooting method useful to evaluate wave function and time-independent correlation function of a particle around attraction by the harmonics oscillator with Cosine asymmetric potential. The scheme of the paper is as follows. In Section 2, we write the basic time-independent Schrödinger equation in terms of finite difference and the harmonics oscillator Cosine asymmetric potential in terms of the new variable is given by where 2 +2 cos ( ) is the Cosine asymmetric potential. In Section 3, we show the idea of writing a program for evaluating energy eigenvalue wave function and correlation function of atomic density for the Cosine asymmetric potential via the numerical shooting method (Asaithambi, Ledoux and van Daele, Boonchui and Hutem [12,20,21]). Section 4 contains our conclusions.

Time-Independent Schrödinger Equation in Finite Difference Formula for Harmonics Oscillator Cosine Asymmetric Potential
We consider a particle of mass moving on the -axis in a time-independent potential V( ). The time-independent Schrödinger equation corresponding to this one-dimensional motion is where E is the total energy eigenvalues of the particle. The solution of this equation yields the allowed energy eigenvalues E and the corresponding wave function ( ). To solve this partial differential equation, we need to specify the potential V( ) as well as the boundary condition; the boundary condition can be obtained from the physical requirement of the system. Suppose a particle is bound state to around of attraction by the harmonics oscillator Cosine asymmetric potential (see Figure 1): with + cos ( ) that is called the Cosine asymmetric potential [22], where , , and are positive constants. Substituting the harmonics oscillator Cosine asymmetric potential from (3) into (2) leads to the following equation: For our approach (4) can be solved in the numerical shooting method. It is convenient to simplify the arithmetic involved in the shooting solution. We define some new dimensionless variables. Then the position variable is replaced with the dimensionless variable : With this definition, the second-derivative term can be written as Substituting in for in terms of and setting = 2 /ℏ and setting ℏ = = = 1 into (4), we can rewrite the timeindependent Schrödinger equation completely in terms of as follows: Also, the time-independent potential in terms of the new variable is given by For the finite difference method, (7) is rewritten in the form of many small segment Δ in the domain length. The secondderivative of the first term in (7) can be approximated in finite difference form as follows:  We can obtain the form of the time-independent Schrödinger equation in terms of finite difference by substituting (9) into (7), and we get where +1 = Δ + . The special potential given by harmonics oscillator Cosine asymmetric potential has been used in calculating (10) in the mathematica program (see Section 3).

Numerical Shooting Method and Results
We construe the new variable to be used in calculating the ground-state energy eigenvalue, wave function, and the timeindependent correlation function of the harmonics oscillator Cosine asymmetric potential.  (1) min is the start position in the analysis range.
(2) max is the ultimate position in the analysis range.
(3) is any position in the analysis range.
(4) is a number of very small bars in the analysis range.
(5) Δ is the length of very small bars so that The logic of the numerical shooting method evaluation of energy eigenvalue, eigenfunction, and time-independent correlation function for the harmonics oscillator Cosine asymmetric potential is as follows.
(i) Input values min and max in mathematica program for the harmonics oscillator Cosine asymmetric potential. Find the initial value for calculation. Input the initial condition by setting 1 = 0 for the position imprisons and set / = 1 from the slope of positions 1 and 2, so that By inputting 1 and 2 as two initial values for calculation, we can find 3 from (10). In the same way, we can find 4 by substituting 2 and 3 in the equation. As we keep doing this, we can find (see Figure 2 in [12]).
(i) The next task is to calculate wave function in (10) ( +1 ) so that it approaches zero as closely as desired. Normally, we assign a small value as the standard to make sure that wave function in (10)    (viii) Plot the wave function that is normalized for the harmonics oscillator potential by the graph related to .

Conclusion
In conclusion, we then represented the method by obtaining numerical solution of the one-dimensional harmonic oscillator, perturbed from a set of the Cosine asymmetric potentials.
Physics Research International  In our calculation, we can obtain the time-independent correlation function corresponding with the Green functions techniques [16]. Although the numerical shooting method does not use complicated formulas, it requires much less computational effort when compared to the Green functions techniques. Generally, regarded as one of the most efficient methods, the numerical shooting method [12] gives very accurate results because it integrates the Schrdinger equation directly, though in the numerical sense.
In this case, the wave function of the harmonics oscillator Cosine asymmetric potential is different from that in case of a typical harmonics oscillator potential (V( ) = (1/2) 2 2 ) (see . In this case, the time-independent correlation function ( ( )) of the harmonics oscillator Cosine asymmetric potential via numerical shooting method and the intensity correlation experiment by reference [23] have the same appearance (see