Quantum Effect in a Diode Included Nonlinear Inductance-Capacitance Mesoscopic Circuit

Similar documents
Bose Description of Pauli Spin Operators and Related Coherent States

Infinite Sequence Soliton-Like Exact Solutions of (2 + 1)-Dimensional Breaking Soliton Equation

Universal Associated Legendre Polynomials and Some Useful Definite Integrals

Generalized projective synchronization between two chaotic gyros with nonlinear damping

698 Zou Yan-Li et al Vol. 14 and L 2, respectively, V 0 is the forward voltage drop across the diode, and H(u) is the Heaviside function 8 < 0 u < 0;

Synchronization and Bifurcation Analysis in Coupled Networks of Discrete-Time Systems

Chaos suppression of uncertain gyros in a given finite time

A New Kind of k-quantum Nonlinear Coherent States: Their Generation and Physical Meaning

New Integrable Decomposition of Super AKNS Equation

Multiple Image Method for the Two Conductor Spheres in a Uniform Electrostatic Field

Complete and Partial Separability Conditions of Mixed State for a Quantum Network of Any Nodes

An Improved F-Expansion Method and Its Application to Coupled Drinfel d Sokolov Wilson Equation

New Feedback Control Model in the Lattice Hydrodynamic Model Considering the Historic Optimal Velocity Difference Effect

Nonchaotic random behaviour in the second order autonomous system

POSTULATES OF QUANTUM MECHANICS

Shallow Donor Impurity Ground State in a GaAs/AlAs Spherical Quantum Dot within an Electric Field

Time Evolution, Dynamical Quantum Fluctuation and High-Order Squeezing Feature in Polariton System I

New Homoclinic and Heteroclinic Solutions for Zakharov System

Dynamical analysis and circuit simulation of a new three-dimensional chaotic system

The rotating Morse potential energy eigenvalues solved by using the analytical transfer matrix method

Time-delay feedback control in a delayed dynamical chaos system and its applications

Dynamical behaviour of a controlled vibro-impact system

A Condition for Entropy Exchange Between Atom and Field

Some exact solutions to the inhomogeneous higher-order nonlinear Schrödinger equation by a direct method

Scheme for Asymmetric and Deterministic Controlled Bidirectional Joint Remote State Preparation

Projective synchronization of a complex network with different fractional order chaos nodes

Rational Form Solitary Wave Solutions and Doubly Periodic Wave Solutions to (1+1)-Dimensional Dispersive Long Wave Equation

Superposition of two mesoscopically distinct quantum states: Coupling a Cooper-pair box to a large superconducting island

Partial factorization of wave functions for a quantum dissipative system

A Generalized Extended F -Expansion Method and Its Application in (2+1)-Dimensional Dispersive Long Wave Equation

Solving ground eigenvalue and eigenfunction of spheroidal wave equation at low frequency by supersymmetric quantum mechanics method

Bidirectional Partial Generalized Synchronization in Chaotic and Hyperchaotic Systems via a New Scheme

Quantum Theory of Light and Matter

A Realization of Yangian and Its Applications to the Bi-spin System in an External Magnetic Field

Critical entanglement and geometric phase of a two-qubit model with Dzyaloshinski Moriya anisotropic interaction

What Is Semiquantum Mechanics?

550 XU Hai-Bo, WANG Guang-Rui, and CHEN Shi-Gang Vol. 37 the denition of the domain. The map is a generalization of the standard map for which (J) = J

Formation Mechanism and Binding Energy for Icosahedral Central Structure of He + 13 Cluster

Inductance, Inductors, RL Circuits & RC Circuits, LC, and RLC Circuits

Application of Mean-Field Jordan Wigner Transformation to Antiferromagnet System

New Application of the (G /G)-Expansion Method to Excite Soliton Structures for Nonlinear Equation

Propagation of Lorentz Gaussian Beams in Strongly Nonlocal Nonlinear Media

Effects of Atomic Coherence and Injected Classical Field on Chaotic Dynamics of Non-degenerate Cascade Two-Photon Lasers

A New Integrable Couplings of Classical-Boussinesq Hierarchy with Self-Consistent Sources

A tunable corner-pumped Nd:YAG/YAG composite slab CW laser

Classical calculation of radiative lifetimes of atomic hydrogen in a homogeneous magnetic field

Absorption-Amplification Response with or Without Spontaneously Generated Coherence in a Coherent Four-Level Atomic Medium

The u(2) α and su(2) α Hahn Harmonic Oscillators

Exciton-Dependent Pre-formation Probability of Composite Particles

Radiation energy flux of Dirac field of static spherically symmetric black holes

Backstepping synchronization of uncertain chaotic systems by a single driving variable

A Piezoelectric Screw Dislocation Interacting with an Elliptical Piezoelectric Inhomogeneity Containing a Confocal Elliptical Rigid Core

Compound buried layer SOI high voltage device with a step buried oxide

Generalized projective synchronization of a class of chaotic (hyperchaotic) systems with uncertain parameters

Second Quantization Method for Bosons

Thermal quantum discord in Heisenberg models with Dzyaloshinski Moriya interaction

Elastic behaviour of an edge dislocation near a sharp crack emanating from a semi-elliptical blunt crack

Quantum secret sharing based on quantum error-correcting codes

Semi-Relativistic Reflection and Transmission Coefficients for Two Spinless Particles Separated by a Rectangular-Shaped Potential Barrier

Fano resonance and wave transmission through a chain structure with an isolated ring composed of defects

Realization of Two-Qutrit Quantum Gates with Control Pulses

Experimental observation of direct current voltage-induced phase synchronization

Anti-synchronization Between Coupled Networks with Two Active Forms

Effects of Particle Shape and Microstructure on Effective Nonlinear Response

Exact solutions of the radial Schrödinger equation for some physical potentials

A COMPACT PI-STRUCTURE DUAL BAND TRANSFORMER

Function Projective Synchronization of Fractional-Order Hyperchaotic System Based on Open-Plus-Closed-Looping

Department of Applied Mathematics, Dalian University of Technology, Dalian , China

Quantize electrical circuits

Ionization Potentials and Quantum Defects of 1s 2 np 2 P Rydberg States of Lithium Atom

Impulsive control for permanent magnet synchronous motors with uncertainties: LMI approach

Isospin and Symmetry Structure in 36 Ar

Construction of a New Fractional Chaotic System and Generalized Synchronization

The correlation between stochastic resonance and the average phase-synchronization time of a bistable system driven by colour-correlated noises

Correlation between classical Fisher information and quantum squeezing properties of Gaussian pure states

No. 5 Discrete variational principle the first integrals of the In view of the face that only the momentum integrals can be obtained by the abo

Probabilistic Teleportation of an Arbitrary Two-Qubit State via Positive Operator-Valued Measurement with Multi Parties

MP463 QUANTUM MECHANICS

Stereodynamics of the O( 3 P) with H 2 (D 2 ) (ν = 0, j = 0) reaction

Non-relativistic Limit of Dirac Equations in Gravitational Field and Quantum Effects of Gravity

Two-mode excited entangled coherent states and their entanglement properties

Lecture 5: Harmonic oscillator, Morse Oscillator, 1D Rigid Rotor

Bifurcation of Sound Waves in a Disturbed Fluid

Representation of Spin Group Spin(p, q)

Information Entropy Squeezing of a Two-Level Atom Interacting with Two-Mode Coherent Fields

The oscillator-bath model of the environment.

Isotopic effect of Cl + 2 rovibronic spectra in the A X system

arxiv: v2 [quant-ph] 6 Jun 2008

Quantum Parameter Estimation: From Experimental Design to Constructive Algorithm

NEGATIVE BINOMIAL STATES OF THE RADIATION FIELD AND THEIR EXCITATIONS ARE NONLINEAR COHERENT STATES

Deleting a marked state in quantum database in a duality computing mode

Dielectric Properties and Lattice Distortion in Rhombohedral Phase Region and Phase Coexistence Region of PZT Ceramics

Solitary Wave Solutions of KP equation, Cylindrical KP Equation and Spherical KP Equation

Effects of Different Spin-Spin Couplings and Magnetic Fields on Thermal Entanglement in Heisenberg XY Z Chain

Differential Representations of SO(4) Dynamical Group

A New Generalized Riccati Equation Rational Expansion Method to Generalized Burgers Fisher Equation with Nonlinear Terms of Any Order

Knotted Pictures of Hadamard Gate and CNOT Gate

Theoretical Analysis of Neutron Double-Differential Cross Section of n + 19 F at 14.2 MeV

Appendix A. The Particle in a Box: A Demonstration of Quantum Mechanical Principles for a Simple, One-Dimensional, One-Electron Model System

NILPOTENT QUANTUM MECHANICS AND SUSY

Transcription:

Commun. Theor. Phys. (Beijing, China) 52 (2009) pp. 534 538 c Chinese Physical Society and IOP Publishing Ltd Vol. 52, No. 3, September 15, 2009 Quantum Effect in a Diode Included Nonlinear Inductance-Capacitance Mesoscopic Circuit YAN Zhan-Yuan, ZHANG Xiao-Hong, and MA Jin-Ying School of Mathematics and Applied Physics, North China Electric Power University, Baoding 071003, China (Received August 8, 2008; Revised May 11, 2009) Abstract The mesoscopic nonlinear inductance-capacitance circuit is a typical anharmonic oscillator, due to diodes included in the circuit. In this paper, using the advanced quantum theory of mesoscopic circuits, which based on the fundamental fact that the electric charge takes discrete value, the diode included mesoscopic circuit is firstly studied. Schrödinger equation of the system is a four-order difference equation in ˆp representation. Using the extended perturbative method, the detail energy spectrum and wave functions are obtained and verified, as an application of the results, the current quantum fluctuation in the ground state is calculated. Diode is a basis component in a circuit, its quantization would popularize the quantum theory of mesoscopic circuits. The methods to solve the high order difference equation are helpful to the application of mesoscopic quantum theory. PACS numbers: 73.23.-b, 73.23.Ra, 03.65.Ca Key words: mesoscopic circuit, quantum effects 1 Introduction In recent years, research work of macroscopic quantum effects and mesoscopic physics has attracted great interest. [1 With the dramatic development in nanotechnology and nanoelectronics, the integrated circuits and components have been miniaturized toward atomic-scale dimensions in the manufacture of electronic devices. Undoubtedly, the quantum mechanical effect of circuits and devices must be taken into account. The quantum theory of mesoscopic circuits is studied as early as 1970 s, comparing with the quantization process of a harmonic oscillator, the L-C design circuit is quantized by Louisell. [2 In the middle of 90 s, much progress has been made in many aspects, such as quantum fluctuation, squeezing effect of some fundamental mesoscopic circuits. [3 6 However, they seldom consider the discreteness of the electronic charge. In Ref. [8, an improved quantum theory for mesoscopic electric circuits considering the discreteness of electric charge is proposed, which caused a new research hotpot in mesoscopic circuits, in the framework of a profound quantum theory for mesoscopic circuits, the quantum effects of RLC circuits, [9 11 inductance and capacitance coupling circuit [12,13 transmission lines [14 are studied. The circuits studied before are linear systems. If the diode is included in the mesoscopic circuit, the system acquires a nonlinear characteristic. In this paper, a mesoscopic circuit describing the anharmonic oscillator [15 studied in the framework of quantum theory developed by Li and Chen. [8 The anharmonic oscillator is quantized, and the energy spectrum, wave functions of the system are obtained. Diode is a basis component in a circuit, its quantization would popularize the quantum theory of mesoscopic circuits. 2 Quantization of Diode Included Nonlinear Inductance-Capacitance Mesoscopic Circuit Fig. 1 Diagram of mesoscopic nonlinear inductance-capacitance circuit. A diode included nonlinear inductance-capacitance mesoscopic circuit describing the anharmonic oscillator is shown in Fig. 1. Kirchhoffs voltage rule in the circuit is expressed as [16 L q + V C (q) = 0, (1) where q is the electric charge and V c is the voltage drop across both capacitors. In Eq. (1), the value for V C will vary according to the amount of charge on the capacitances C 1 and C 2. A plot of the voltage across both capacitances as a function of the charge on the capacitances is shown in Fig. 2. When the potential drop across C 2 is Supported by National Natural Science Foundation of China under Grant No. 10575028 E-mail: zhanyuanyan@tom.com

No. 3 Quantum Effect in a Diode Included Nonlinear Inductance-Capacitance Mesoscopic Circuit 535 lower than a certain voltage, the diodes do not conduct, the circuit equivalent capacitance C e is that of the two capacitances 1/C e = 1/C 1 + 1/C 2. As the charge on the capacitances increases, the potential drop across C 2 will reach a certain voltage and the diodes begin to conduct. The circuit equivalent capacitance increases and the slope of the line shown in Fig. 2 decreases to the value of 1/C 1. In this instance, the diodes act as closed switches and effectively remove the capacitance C 2 from the circuit. Fig. 2 Voltage versus charge across both capacitors. In Eq. (1), the form of V c (q) is given by the following piecewise function { 1/Ce q(t) < q b, V c (q) = q (2) 1/C 1 otherwise, where q b is the critical charge required to conduct the diodes. The three straight line segments can be approximated by a continuous of the form V c (q) = q µ νq3, (3) where µ and ν are the fitting values corresponding to experimental values. [7 As a matter of fact, the electric charges are discrete. So in the process of quantization of the mesoscopic circuits, the eigenvalues of the self-adjoint operator ˆq should take discrete values, following the quantization method in [8, ˆq q = nq e q, (4) where n Z, q stand for eigenstates of electric charge for the circuits and q e = 1.6 10 19 Coulomb, is the elementary charge. A minimum shift operators is introduced, ˆQ = exp(iq eˆp/ h), (5) the important self-adjoint momentum operator is defined as ˆP = h 2i ( q e + qe ) = h 2iq e ( ˆQ ˆQ ), (6) in which, the right and left discrete derivative operators are qe = ( ˆQ 1)/q e, qe = (1 ˆQ )/q e, (7) and the free Hamilton operator is Ĥ 0 = h2 2L ( q e qe ) = h2 2 ( ˆQ + ˆQ + 2) Quantized Hamilton of the system is = h2 2Lq e ( qe qe ). (8) Ĥ = h2 2Lq e ( qe qe ) + ˆq2 2µ νˆq4 4. (9) Now, we consider a representation in which the operator ˆp is diagonal and called ˆp presentation. using the following relations p ˆQ 2 p = 2π h3 2 q e p 2 δ(p p ), p ( qe qe ) p = 4π h q 2 e [ cos ( qe h p ) 1 δ(p p ). (10) The Schrödinger equation of nonlinear inductancecapacitance mesoscopic circuit in ˆp presentation is ν h 4 4 2 h2 p4ψ + 2µ p 2 ψ { h 2 } + qe 2L[cos(q ep/ h) 1 + E ψ = 0. (11) In above, the nonlinear inductance-capacitance mesoscopic circuit is quantized in the framework of the quantum theory for mesoscopic circuits based on discrete charge space, it is a four order difference equation. 3 Wave Function and Quantum Fluctuations of System Ordinarily, a linear inductance-capacitance-resistancepower source systems Schrödinger equations are two order difference equation, they could be solved through transforming them into a formation of Mathieu equation. [8,11,13 To solve the Schrödinger equation of the nonlinear inductance-capacitance mesoscopic circuit system, we extended the perturbative method [18 of two order difference equation to four order difference equation, this method is applicable for other high order difference equations. By introducing π q e h p = 2z, λ = 64 h2 νlq 6 e + 64E νq 4 e ξ = 32 h2 νlqe 6,, (12) the Schrödinger equation of nonlinear inductancecapacitance mesoscopic circuit is convert to: z 4ψ + 8 2 2 z2ψ + [λ 2ξ cos(2z)ψ = 0. (13)

536 YAN Zhan-Yuan, ZHANG Xiao-Hong, and MA Jin-Ying Vol. 52 For the concrete values of the parameter ν, L, Plank constant and the elementary electric charge, ξ 10 1, perturbative method is valid. ψ = ψ (0) + ξψ (1) + ξ 2 ψ (2) + ξ 3 ψ (3) +, λ = λ (0) + ξλ (1) + ξ 2 λ (2) + ξ 3 λ (3) +, (14) where ψ (0) and λ (0) = k 4 + (8/µq 2 el)k 2 are zero order approximate of wave function and eigenvalue, ψ (1,2,3,...) and λ (1,2,3,...) are high order corrections of wave function and eigenvalue. When k = 0, 1, 2,..., ψ have separate solutions. Wave function ψ (0,1,2) and eigenvalue λ (0,1,2) are determined by the following equations: ψ (0) z 4 + 8 2 ψ (0) 2 z nψ (0) = 0, z 4 + 8 2 ψ (1) 2 z nψ (1) + λ (1) ψ (0) 2 cos(2z)ψ (0) = 0, z 4 + 8 2 ψ (2) 2 z nψ (2) + λ (1) ψ (1) + λ (2) ψ (0) 2 cos(2z)ψ (1) = 0. The detail calculation process of the wave function and energy lever are in the appendix. Calculated to second order approximate, the separate periodic wave function are: 2ψ0 = 1 + 4 h2 cos2z 2 ) 2 2 cos4z 6 1 2/ 6 (2 1/ 2)(1 4/µνq2 e ) + ψ 1 + 2 h2 cos3z 2 ) 2 [ 4 cos3z = cosz + 6 5 4/ 6 (5 4/ 4 cos5z )2 (39 12/ 2)(5 4/µνq2 e ) + ψ1 2 h2 sin3z 2 ) 2 [ 4 sin 3z = sin z + 6 5 4/ 6 (5 4/) 2 4 sin5z (39 12/)(5 2 4/) [ ψ h2 2 = cos2z + 6 2/ 2 1 + 2 cos4z 2 ) 2 cos6z 15 6/ 6 (20 6/ 2)(15 6/µνq2 e ) + ψ2 2 h2 sin 4z 2 ) 2 sin 6z = sin 2z + 6 15 6/µqe 2L + 6 (20 6/µqe 2L)(15 6/µq2 e L) + (16) while the superscripts + and specify the even and odd parity solution respectively, the corresponding energy levels are E 0 = h2 Lqe 4 h4 8 1 1 4/ E 1 + = 3 h2 q2 e 8µ νq4 e 64 + h 4 L 2 νqe(5 8 4/) E1 = h2 q2 e 8µ νq4 e 64 h 4 L 2 νqe 8(5 4/µνq2 e ) + E + 2 = h2 Lq 2 e + q2 e 2µ νq4 e 4 + 2 h 4 L 2 νq 8 e (2/µνq2 e 1) + h 4 L 2 νq 8 e (15 6/µνq2 e ) + E2 = h2 Lqe q2 e 2µ νq4 e 4 + h 4 L 2 νqe(15 8 6/) In calculation, we found that the energy spectrum E + k and E k would have the same expression when k 3. This does not mean the degenerate of energy level because we only calculated to ξ 2 phase, however we can conclude that the energy level will become more and more close with the increasing of k. To verify our results, let ν = 0 and µ = C 1 in Eqs. (17) and (18), the system becomes to a LC oscillator, our results are the same as the results in Ref. [19. As an application of the results obtained above, we calculated the average values of currents P and P 2 in the ground state. P 0 = 0, (15) (17)

No. 3 Quantum Effect in a Diode Included Nonlinear Inductance-Capacitance Mesoscopic Circuit 537 P 2 h 2 0 = ( 2 2)2 (2 2 1)2 ( 2 [ 4)2 32 176µν + 346qeµ 2 2 ν 2 302qeµ 4 3 ν 3 + 257 2 q6 eµ 4 ν 4 26qeµ 8 5 ν 5 + 2qe 10 µ 6 ν 6 q 2 e + 16 h4 µ 2 92 h 4 µ 3 ν L2 qe 8 + 146 h4 µ 4 ν 2 L2 qe 6 55 h4 µ 5 ν 3 L2 qe 4 + 6 h4 µ 6 ν 4 L2 qe 2L2 qe 10 + 4 h8 µ 4 4 h 8 µ 5 ν h 8 µ 6 ν 2 + L4 L4 qe 18 qe 16 qe 14L4. (18) From Eq. (18), we can see that the average values of the currents are zero, whereas the average values of the square of the currents are not. This indicated the existence of the current quantum fluctuations, it is the shot noise in the system. In the above, using an advanced quantum theory to mesoscopic circuits developed recently by Li and Chen, which based on the fundamental fact that the electric charge takes discrete value, a diode included nonlinear inductancecapacitance mesoscopic circuit is quantized. The Schrödinger equation in ˆp representation is a four order difference equation. Using the power series and perturbative method, the detail energy spectrum and wave functions are obtained and verified, as an application of the results, the current quantum fluctuation in the ground state is calculated. This paper is an application of the profound quantum mesoscopic circuits theory to the diode included circuit, the methods of quantization and solution process would be helpful to the popularization of the mesoscopic quantum theory. Appendix A: Solution of Four Order Difference Equation In Sec. 3 Eq. (14), the Schrödinger equation in ˆp representation is a four order difference equation. To make the equation simple, let a = 8/ 2, Eq. (14) becomes 2 z4ψ + a z2ψ + [λ 2ξ cos(2z)ψ = 0. (A1) ψ(z) could be periodic function of π or 2π. We show the calculation process of ψ(z) with the periodic π, let ψ(z) = [c k sin(kz) + d k cos(kz), (A2) k=0,2,4,... insert Eq. (A2) into Eq. (A1), calculate to the order k = 8, λd 0 ξd (16c 2 4ac 2 ξc 4 )sin 2z + ( 2ξd 0 + 16d 2 4ad λd 2 ξd 4 )cos 2z + ( ξc 256c 4 16ac 4 + λc 4 ξc 6 )sin 4z + ( ξd 256d 4 16ad 4 + λd 4 ξd 6 )cos4z + ( ξc 4 + 1296c 6 36ac 6 + λc 6 ξc 8 )sin 6z + ( ξd 4 + 1296d 6 36ad 6 + λd 6 ξd 8 )cos6z + ( ξc 6 + 4096c 8 64ac 8 + λc 8 )sin 8z + ( ξd 6 + 4096d 8 64ad 8 + λd 8 )cos8z = 0. (A3) Obviously, the coefficients satisfy the following matrix equation 16 4a + λ ξ 0 0 ξ 256 16a + λ ξ 0 0 ξ 1296 36a + λ ξ 0 0 ξ 4096 64a + λ 0 0 0 ξ λ ξ 0 0 0 2ξ 16 4a + λ ξ 0 0 0 ξ 256 16a + λ ξ 0 0 0 ξ 1296 36a + λ ξ 0 0 0 ξ 4096 64a + λ 0 0 0 0 ξ c 2 c 4 c 6 c 8 = 0, d 0 d 2 d 4 d 6 d 8 (A4) = 0. (A5) The same matrix equations of 2π period ψ could be obtained with the same method. We could solve λ and c k, d k through this two equations theoretically, however, it is too difficult to calculate the roots for an equation of 5-th Order. if ξ is small, perturbative theory is valid. The zero order approximate of λ is easy obtained from the diagonal elements of coefficients matrix, λ = k 4 + ak, k = 0, 1, 2, 3, (A6)

538 YAN Zhan-Yuan, ZHANG Xiao-Hong, and MA Jin-Ying Vol. 52 According to Eqs. (15) and (16), higher order approximates of λ k and ψ k could be calculated. For example, we calculate λ + 1 and ψ+ 1. When k = 1, λ (0) 1 = a 1. The zero approximates odd solution of ψ (0) 1 could obtained from the first formula of Eq. (16) the solution is ψ (0) 1 = cosz. Insert ψ (0) 1 into the second formula of Eq. (16), ψ (0) z 4 + a 2 ψ (0) z (a 1)ψ (0) = 0, z 4 + a 2 ψ (1) z (n 1)ψ (1) + (λ (1) 1)cosz cos3z = 0, in the equation above, let λ (1) = 1, ψ (1) has a period solution, satisfy: the solution of ψ (1) is Insert ψ (0), ψ (1) into the third formula of Eq. (16), z 4 + a 2 ψ (1) z (n 1)ψ (1) = cos3z, ψ (1) = cos3z 80 8a. (A7) z 4 + a 2 ψ (2) z (a 1)ψ (2) + cos3z 80 8a + λ(2) cosz cos5z 80 8a cosz 80 8a = 0. Let λ (2) = 1/(80 8a), ψ (2) has a period solution, satisfy: the solution of ψ (2) is z 4 + a 2 ψ (2) z (a 1)ψ (2) = cos3z 80 8a + cos5z 80 8a, ψ (2) = cos3z (80 8a) cos5z (624 24a)(80 8a). (A8) So, the final results of ψ 1 + and λ+ 1 are ψ 1 + cos3z = cosz + ξ 80 8a + ξ2[ cos3z (80 8a) cos5z, (A9) (624 24a)(80 8a) λ + 1 = a 1 + ξ + 1 ξ2 80 8a. (A10) References [1 W.H. Zurek, Phys. Today 44 (1991) 36. [2 W.H. Louisell, in Quantum Statistical Properties of Radiation 1, John Wiley, New York (1973). [3 C.P. Sun and L.H. Yu, Phys. Rev. A 49 (1994) 592. [4 L.H. Yu and C.P. Sun, Phys. Rev. A 51 (1995) 1845. [5 H.Y. Fan and X.T. Liang, Chin. Phys. Lett. 17(3) (2000) 174. [6 J.S. Wang and C.Y. Sun, Int. J. Theor. Phys. 37 (1998) 1213. [7 E.L.M. Flerackers, H.J. Janssen, and L. Beerden, Am. J. Phys. 53 (1985) 574. [8 Y.Q. Li and B. Chen, Phys. Rev. B 53 (1996) 4027. [9 Y.Q. Li, in Proc. 5th Winger Symposium, eds. P. Kasperkovitz and D. Grau, World Scientific, Singpore (1998) 307. [10 T. Lu and Y.Q. Li, Mod. Phys. Lett. B 16 (2002) 975. [11 J.X. Liu and Z.Y. Yan, Commun. Theor. Phys. 45 (2006) 1126. [12 J.S. Wang, T.K. Liu, and M.S. Zhan, Phys. Lett. A 276 (2000) 155. [13 J.X. Liu and Z.Y. YAN, Commun. Theor. Phys. 44 (2005) 1091. [14 J.C. Flores, Phys. Rev. B 64 (2001) 235309. [15 R.H. Enns and G. McGuire, in Nonlinear Physics with Maple for Scientists and Engineers, Birkhauser Boston, Cambridge (1997) Chap. 2, p. 12. [16 S. Zhang, C.I. Um, and K.H. Yeon, Chin. Phys. Lett. 19 (2002) 985. [17 J.X.Liu and Z.Y. Yan, Commun. Theor. Phys. 44 (2005) 1091. [18 T.Z. Liu and L.Q. Chen, in Nonlinear Vibrations, Higher Education Press, Beijing (2001). [19 Y.Q. Li and B. Chen, Commun. Theor. Phys. 29 (1998) 137.