Dynamics of Particle in Confined-Harmonic Potential in External Static Electric Field and Strong Laser Field

Dynamics of a particle in confined-harmonic potential, subjected to external static electric and time-dependent laser fields is studied. The energy levels and wave functions of unperturbed harmonic oscillator are evaluated using B-polynomial Galerkin method. Matrix formulation is used throughout the procedure. This procedure is very simple and efficient in comparison with other methods. Modifications of wave functions and energy levels due to static electric field are also calculated. Finally, absorption spectra of such a driven oscillator are studied and explained.


Introduction
The systems for which exact quantum mechanical solutions for Schrödinger equation can be found are few in number, for example, the harmonic oscillator potential and nonrelativistic hydrogen atom.The harmonic oscillator potential is a model of great practical importance, as it approximates any arbitrary potential close to equilibrium.In nanotechnology, potentials of simple shape such as quantum dots are often well approximated by such parabolic potentials.In fact, almost all exactly solvable problems in Quantum Mechanics are harmonic oscillator problems in disguise.
The confined-harmonic oscillator potential plays an important role in many applications of Quantum Mechanics.Such a potential is extensively used to describe the bound states of nonrelativistic systems.It also plays a basic role in chemical and molecular physics.In quantum chemistry, simple harmonic potential is used as a simplified model to describe vibrational motion of two atoms, where, more precise model is the Morse potential.In nonrelativistic quantum mechanics, the Schrödinger equation for this potential has been studied for systems ranging from 1-Dimensional to D-Dimensional Space [1][2][3][4].Such a system has been widely studied as it can be exactly solved and is a very relevant system [5].
The perturbation of quantum harmonic oscillators with external fields has recently attracted a renewed interest due to different aspects of the problem, catalysed by recent developments as follows: 1) quantum dynamics of ion in a Paul trap [6], 2) confining potentials for various quantum heterostructures, which leads to modifications of various physical properties of the media they are composed of [7,8] 3) dynamics of a harmonic oscillator with time-dependent force constant and perturbed by weak quartic anharmonicity [9], 4) need for exact propagators for the anisotropic two-dimensional charged harmonic oscillator in presence of external fields [10].
The effects of external fields on systems under the effect of other types of potentials like pseudo-harmonic oscillator potential have also been explored in literature.For example, the effect on energy levels of a 2D Klein Gordon particle under pseudo-harmonic oscillator interaction has been studied [11].The Schrödinger equation has been solved for a particle in the general 1D time-dependent linear potential [12].The quantum motion of an electron driven by a strong time-dependent linear potential in a 1D quantum wire has been investigated and interesting physical properties studied [13].The possibility of exactly manipulating the quantum motional states of a single particle held in a double cosine potential by using laser beams has been explored [14].
Time-dependent perturbations of such systems have also been studied extensively [15,16].Explicit wave functions and geometric phases of time-dependent harmonic oscillator in external time-dependent magnetic and electric field have been derived [17].The exact wave functions and eigenvalues of a 2D time-dependent harmonic oscillator under the influence of a static magnetic field have been calculated [18].The time evolution of a 2D harmonic oscillator, with time-dependent mass and frequency, in a static magnetic field has also been studied analytically [19].
An electron in confined-harmonic oscillator potential exposed to an external electric field is equivalent to a charged harmonic oscillator in a uniform electric field or a harmonic oscillator in an external dipole field.Such a system has an important role in quantum chemical applications [20].Recently, O. Kidun and D. Bauer [21] have studied two interacting electrons in harmonic potential driven by a strong laser field.They have studied population dynamics of the system.They have further shown the conditions of complete survival and complete depletion of the ground state of "harmonium".C. Liang et al. have studied the properties of Hooke's atom (two electrons interacting with Coulomb potential in an external harmonic oscillator potential) in an arbitrary time-dependent electric field [22].The dynamics of a perturbed quantum Hooke's atom exposed to intense ultrashort laser pulses has been studied by Torres and Vicario [23].
The traditional techniques of studying such quantum mechanical systems have lately been supplemented by finite basis set methods like B-spline [3,[24][25][26] and Bernstein-polynomial (B-polynomial) methods [27][28][29].Recently, Heidari et al. [30] have investigated the case of Hydrogen atom in spherical cavity using B-spline basis functions.The energy spectra of one-and two-electron atoms centered in an impenetrable spherical box have been calculated by Shi Ting Yun et al. by applying Bspline method [31].The B-spline basis set is highly flexible and localized which leads to very accurate results.The B-spline basis functions of degree are piecewise polynomials defined on a knot sequence.When the number of B-splines is taken as , the basis set becomes a set of continuous B-polynomials over the range under consideration [32].These B-polynomials are independent of the grid defined by knots and are simple algebraic polynomials.Each of these polynomials is positive and their sum is unity.
Polynomials are incredibly useful mathematical tools as they can be calculated very easily and accurately on computer systems.Their evaluation is also fast.They are capable of representing a tremendous variety of functions, can be differentiated and integrated quite easily, and can be pieced together to form spline curves that can approximate any function to any desired accuracy.The B-polynomial method is, therefore, much simpler and efficient.Recently, J. Liu et al. have proposed a new numerical method based on B-polynomials expansion for solving one dimensional elliptic interface problems [33].B-polynomial basis has also been used for numerically solving differential equations [34][35][36].
In this paper, the dynamics of an electron in a confined-harmonic potential in static electric and strong laser fields is studied.We have used B-polynomial Galerkin method to solve static field modified harmonic oscillator system.The populations of states modified by static electric field are calculated.The eigenenergies, eigenfunctions and dipole matrix elements of the system are also calculated.The interaction of static field modified confined-harmonic oscillator system with the laser field is taken into account by non-perturbative quasi-energy technique [37][38][39][40].The sequence of the paper is as follows.In Section 2, necessary description of B-polynomials is given.In Section 3, the model under consideration is defined and methods adopted for solving the time-independent as well as time-dependent Schrödinger equation are given.Section 4 deals with interpretation of results and finally, in Section 5, concluding remarks are made.

Bernstein-Polynomial Basis
The B-polynomials [41] of degree over an interval [a, b] are defined as [27,32] There are 1 n  , n-th degree B-polynomials.For mathematical convenience, we usually set More details of these polynomials are available in literature [24,28,29,32,35,42].
The B-polynomial Galerkin method is employed to solve Schrödinger equation for the present case.In the area of numerical analysis, Galerkin methods are a class of methods for converting a problem such as a differential equation to a linear system of equations.A few of the related formulas used are mentioned here for reference.
, , (7) where the l  are expressed as

Problem Formulation and Method of Solution
Consider an electron under the effect of a confined-harmonic oscillator potential subjected to an external static electric field 0 E , where 0 is the strength of the electric field.The units used throughout are the atomic units, i.e., . The confining potential is given by where is a positive constant representing the strength of the potential called the force constant.If the electric dipole moment of the electron is denoted by , the potential energy of the electron due to electric field is given by .The electric dipole moment of the electron is given by , representing the position vector of the electron with respect to the origin and , the charge.

r q
Assume that the electric field is along x direction, therefore the potential energy term becomes 0 xE  .The Hamiltonian for the system can be written as Therefore, the Schrödinger equation for the system becomes A fixed interval   , a b is chosen to study the system.The desired solution may be expanded in terms of a set of continuous polynomials over the closed interval and is given by where s are the coefficients of expansion and are B-polynomials of degree as defined in Section 2. Substituting Equation (12) into Equation (11), taking scalar product with n   , j n B x on both sides and using Equation (7), Equation (11) becomes where the matrix elements , i j a , , i j , , i j b f and , i j assume closed forms by applying the formulas in Section 2 [29].Equation (13) where the column matrix can be determined by solving this symmetric generalized eigenvalue problem.
and the number of B-polynomials is taken to be 26.The accuracy and efficiency of the method depend on the number of B-polynomials chosen to construct the approximate solutions.In the present case, the number of B-polynomials is taken to be 26 as there is not much gain in accuracy beyond this value.A , , B F and in Equation ( 14) are 26 × 26 matrices.The standard Fortran EISPACK library is used to solve the generalized eigenvalue problem and find the eigenvalues and eigenvectors.The eigenvalues D  give the energy levels of the system.The initial eigenvalues for 0 0 E  have been found to be correct to five places of decimal.The eigenvectors are used to calculate the corresponding wave functions using Equation (12).These wave functions are the dressed states of the system and are denoted by where H is given by Equation (10).The time-dependent Schrödinger equation for the system is now written as The solution of Equation ( 16) in quasi-energy formalism can be written as [43]   where i  are defined as quasi-energies and are time-independent eigenvectors to be determined. and applying the exact rotating wave approximation [44], these equations assume the following form The , i j s are the dipole matrix elements and i V  s are the energies of first six levels.The s are defined as and can be easily evaluated using the calculated wave functions  .The set of Equations ( 18) can be solved to determine the quasi-energies i  and the corresponding eigenvectors .These eigenvectors can be used to determine the new dressed state wave functions i i m a   given by Equation (17).In order to solve the set of Equations (18), it is written in matrix form and the corresponding matrix, called the quasi-energy matrix, is diagonalized using standard Fortran subroutines.The calculated eigenvectors are used to determine the transition probabilities to study the absorption spectra.The transition probability from ground state 0 to final state j can be computed from the eigenvectors of the quasi-energy matrix as [45,46] The photoionization probability, ion , i.e., the probability of electron to come out of bound states, is given by where bound states is the sum of the probabilities of the system being in various bound states.Using Equation ( 21) the phenomenon of photoionization is also studied.

Results and Discussion
A single electron in a confined-harmonic oscillator potential is considered to be under the effect of a static electric field.The B-polynomial Galerkin method is used to calculate the dressed states of the confined electron as discussed in Section 3. The variation of eigenvalues for the first six energy states with the static electric field 0 and force constant has been studied for this perturbed system.The values have been plotted in Figure 1 relative to those for 0 E k 0 E  a.u.so that the changes are evident.It is observed that with the increase in the strength of electric field, the energy values are deviated more from the corresponding values for 0 a.u.For a higher force constant, the change in energy values is less.

E 
According to the standard result from perturbation theory for a charged harmonic oscillator in electric field, the energy levels are always lowered by an amount 2 0 2 E k (in atomic units) due to the field.The "dressed" potential [47] in this case is written as which is just a shift of the harmonic potential.From Figure 1 it can be observed that the first two energy levels follow this pattern for low strengths of applied electric field but with the increase in field value the perturbation theory result is not exactly valid and the deviation is found to increase.For the third level it is observed that with increase in k, the relative value is first positive and gradually it becomes negative.The pattern followed by the energy levels is due to the change in wave functions for the system.As a check on the calculations it has been verified that for 1 k  , the energy values for the first few levels, for the range of electric field considered in Figure 1, are in accordance with Equation (22).This is due to the fact that in this case perturbation is small.The increase in energy values with is clearly seen in Figure 2 for The effect of and on the dipole matrix elements 12 , 23 , 34 , 45 and 56 V can be seen from Figure 3.The plots with respect to 0 are for different values of as mentioned in the respective graphs.The values plotted are relative to the corresponding ones for It may be mentioned that the dipole matrix elements for the harmonic oscillator potential are given as Since with the introduction of electric field the system is perturbed, this relation would not be valid.With increase in the value of 0 , the dipole elements diverge from the corresponding values for 0 .For example, increases marginally for some values but  decreases for all .The effect of electric field is much less for higher values.The pattern followed by these values is again related to the change in wave functions.

k k
The dipole matrix elements have been plotted with respect to in Figure 4 for 0 a.u. 12 and 23 V decrease with but , and increase with .The reason for this difference is the fact that in the system only three levels are bound.
The system is now exposed to laser field cos The response of the perturbed system is now investigated by varying different control parameters like force constant , static electric field 0 , laser field strength l and laser frequency   The phenomenon of photoionization also shows up for some 0 values.The peaks for the fourth state, i.e., the first level in the continuum, represent photoionization probability.It is observed that the blue shift for this case is much more than that for the second excited state.It is evident that the first excited state peaks show exact resonance as the probability reaches 0.5.For the peaks corresponding to the second excited state and the fourth level, there is variation in peak strength.This is because the particular frequencies do not represent the condition of exact resonance, i.e., they are slightly off-resonant.

E
By keeping 0 as a.u. and l as a.u., the variation of transition probabilities with respect to a.u. and a.u. and also show blue shift.It may be inferred that for these particular frequencies of the laser field, photoionization probability is more than transition probability.

0.06
The probability for photoionization can be seen more clearly if total probability of bound states and continuum are represented separately.

Summary and Conclusion
The dynamics of an electron in confined-harmonic oscillator potential under the effect of static electric field and strong laser field is studied.The method based on Bpolynomial basis set is employed to solve the Schrödinger equation for the charged confined-harmonic oscillator.The static electric field modifies the wave functions and energies of such confined oscillator and the response of the oscillator to external applied laser field gets affected.Photoionization probabilities show strong dependence on the applied static as well as laser field parameters.
polynomials of degree n form a complete basis over the interval   , a b .The B-polynomials can be generated by a recursive relation[33] system is now exposed to a time-dependent laser field cos l E t  polarized along x-axis, where l is the strength and E  is the frequency of the laser field.The corresponding Hamiltonian becomes

1
lowest energy level of the system under the effect of static electric field and is the number of levels considered.The i N   are the dressed states of the system in presence of laser field.The first six energy levels are taken into account and the range of and is cho-k 0 E sen such that three of them are bound.Substituting the above form of the solution into Equation (16), multiplying by * m  and integrating over dx for 1 m  to results in a set of six homogeneous coupled equations in .Using the orthogonality property of wave func tions 6 i m a

Figure 1 .
Figure 1.Variation of eigenvalues relative to those for 0 0 E  a.u. with respect to static electric field 0 E for various values of force constant .k

Figure 2 .
Figure 2. Variation of eigenvalues with respect to force constant for a.u.k 0 0 E 


. The variation of transition probabilities for first four energy states with respect to  has been depicted in Figure5for force constant 0.06 k  a.u. and laser field a.u.The plots have been made for different values of static electric field 0 .The values chosen to represent variation of transition probabilities have no special significance.These are just some typical values to show relevant effects.It is observed that the resonant frequency for the first excited state shows red shift with increase in electric field but that for the second excited state shows blue shift.The resonance for the first excited state occurs exactly for 0corresponding to the difference in the first two energy levels.The resonance for the second excited state is a two-photon process and occurs at exactly half the energy difference between the ground and the second excited state.

Figure 3 .
Figure 3. Variation of dipole matrix elements relative to those for 0 0 E  a.u. with respect to 0 E for various values of force constant .k

Figure 4 .
Figure 4. Variation of dipole matrix elements with respect to force constant for a.u.k 0 0 E 

Figure 6 .
has been shown for different values of force constant in The figure shows blue shift in resonant frequency for the first as well as the second excited state with increase in value of .The blue shift for the second excited state is less as compared to that for the first excited state.The peaks for the fourth state, representing the probability of photoionization, are very prominently

Figure 5 .Figure 6 .Figure 7 .ωFigure 8 .
Figure 5. Variation of transition probabilities for the first four energy states with respect to the laser frequency  for force constant a.u. and laser field a.u.for different values of static electric field 0.06 k  0.0004 l E 