Original Research ARTICLE
Perturbation Theory for Time-Dependent Quantum Systems Involving Complex Potentials
- Department of Physics, Kyonggi University, Suwon, South Korea
We explore how to apply perturbation theory to complicated time-dependent Hamiltonian systems that involve complex potentials. To do this, we introduce a generalized time-dependent oscillator to which the complex potentials are connected through a weak coupling strength. We regard the complex potentials in the Hamiltonian as the perturbed terms. Quantum characteristics of the system, such as wave functions and expectation values of the Hamiltonian, are investigated on the basis of the perturbation theory. We apply our theory to particular systems with explicit choices of time-dependent parameters. Through such applications, the time behavior of the quantum wave packets and the spectrum of expectation values of the Hamiltonian are analyzed in detail. We confirm that the imaginary parts of expectation values of the Hamiltonian are not zero but very small, whereas the real parts deviate slightly from those of the unperturbed system.
In the case that we are unable to derive exact quantum solutions for a perturbed system, perturbation theory is a useful tool for obtaining approximate quantum solutions. When we apply perturbation theory, it should be supposed that the scale of the perturbed terms in the potential is relatively small. The perturbation theory is valid only when the quantum solutions of the system in which the perturbation potentials have been removed are exactly known or derivable. The perturbation theory was originally developed for Hermitian systems in which the potential is real. Hence, in conventional quantum mechanics, the perturbation theory has, in large, been developed for the systems in which the potentials are real Hermitian that allows only the spectrum of real expectation values for quantum observables. Whereas the eigenvalue problem and the effects of perturbations on stationary Hamiltonian systems are well known in non-relativistic quantum physics, the perturbation techniques for the time-dependent Hamiltonian systems (TDHSs) with complex potentials have been investigated much less. This is partly due to the difficulty of mathematical procedures when we apply perturbation techniques in TDHSs.
Time-independent perturbation theory is a mathematical tool for treating quantum systems whose Hamiltonian involves small static perturbing terms which do not induce transitions to other quantum states. In a case that the Hamiltonian is a function of time, transitions between quantum states may take place. Then, we should consider time-dependent perturbation theory instead of time-independent one. Provided that the perturbing Hamiltonian is differentiable with respect to time in that case, the transition probability is determined from the time derivative of the perturbing Hamiltonian . Hence, if the perturbing Hamiltonian is slowly varying, we can adopt adiabatic theorem which assumes that the quantum system remains in its instantaneous eigenstate . This allows the use of time-independent perturbation theory as an approximation in that situation. For more details of time-dependent perturbation theory and its relevance to time-independent perturbation theory, refer to Mandal and Hunt [1, 3] and Langhoff et al. .
On account of recent attention to the quantum problem of physical systems characterized by complex potentials, the necessity for the extension of perturbation theory to complex potential systems has gradually emerged. The data in some of the processes of elastic scattering such as nucleus-nucleus scattering  and electron scattering from solid crystals  fit better if one adopts complex potentials involving a non-positive imaginary part. Moreover, according to the advancement of techniques for unfolding quantum theory, the necessity of quantum manipulation with complex potentials takes place in many different areas of physics and related topics such as optics [7–11], statistical physics [12, 13], relativistic quantum field theory [14–16], nuclear physics , condensed matter physics [18–20], and biology . To analyze the stability and nonlinear modes of optical beams propagating along an asymmetric waveguide, an imaginary potential has been adopted . The structures of the energy band in elastic scattering processes governed by a periodic complex potential have also been studied .
According to such a research trend, we will investigate, in this work, how to apply perturbation theory to complex potential systems. We will adopt a general Hamiltonian that corresponds to a time-dependent harmonic oscillator to which complex potentials are coupled as perturbation terms. We will consider adiabatic evolution [23–25] of the states under slow variations of parameters while the magnitudes of perturbation terms are quite small compared to the unperturbed term at all times. Then, there will be no transition between eigenstates and, as a consequence, we can use usual time-independent perturbation theory. The wave functions corrected by perturbed potentials will be analytically evaluated. Using the wave functions involving correction terms, quantum characteristics of the system will be investigated. Furthermore, the spectrum of the expectation values of the Hamiltonian (EVH) will be rigorously evaluated. We will apply our development to some particular cases and analyze the time evolution of the corresponding wave functions and the EVH.
2. Theoretical Framework for a Complex Potential System
In this section, we outline the quantum structure for a TDHS which involves complex potentials as perturbations to the system. The equation of motion for lots of actual physical models in quantum mechanics are described not only by linear terms which allow us to have exact mathematical solutions, but also by nonlinear terms associated with interactions of the system with various subsystems. If the coupling parameters of the nonlinear parts are sufficiently small compared to others, we can treat them as perturbation terms. Regarding this, we consider the system of which Hamiltonian is represented as
where Ĥ is a time-dependent quadratic Hamiltonian and Ĥp is a perturbing Hamiltonian, while ϵ is a small coupling parameter. In this case, the perturbation theory is available only when the second term on the right hand side of Equation (1) is very small relative to the first term. We assume that Ĥ and Ĥp are given by
where f(t), g(t), h(t), kR, j(t), and kI, l(t) are time functions that are differentiable with respect to time, and f(t)≠0. m0 is the mass and ω(t) is a time-dependent angular frequency. Although this is a one-dimensional Hamiltonian, we can apply it to various different kinds of oscillators depending on the choice of the explicit formulae of time functions. As examples for such, we will apply it to two special cases at later by choosing the time functions being particular forms. In the case where the mass of the system varies over time, we can regard m0 as the initial mass. Notice that the additional term Ĥp in the Hamiltonian makes the system be an anharmonic oscillator.
The Schrödinger solutions of the system are different depending on the choice of the time functions in the Hamiltonian. If we consider a discrete spectrum of solutions, the Schrödinger equation for the overall Hamiltonian can be written as
To apply the perturbation theory at this stage, it is necessary to know complete quantum solutions associated with the unperturbed Hamiltonian Ĥ as mentioned in the introductory part. For this reason, we first derive exact solutions of the Schrödinger equation for Ĥ from
Let us take a moment to review the theory for solving this equation exactly. Because Ĥ involves time functions, it is very difficult to solve this equation relying on the conventional method. For such TDHSs, it is known based on the invariant operator theory that the Schrödinger solutions are expressed in terms of the solutions of a classical equation associated with the system [26, 27]. Hence, let us consider a time function ρ(t) which is the solution of the following nonlinear equation
Then, we can define the annihilation operator of the system in terms of ρ(t) as
where Ω is an arbitrary real positive constant, and
Of course, the Hermitian adjoint of Equation (8), â†, is the creation operator. These ladder operators obey the commutation relation [â, â†] = 1. The time evolution of â(t) can be simply represented as â(t) = â(0)e−i[η(t)−η(0)], where
We can check that the direct differentiation of Î with respect to time results in zero, which means that Î(t) is a time-invariant. If we write the eigenfunctions of Î as ϕn(q, t), they can be derived from the eigenvalue equation:
This is analogous to the time-independent Schrödinger equation in a stationary system. We can easily solve this equation from a minor evaluation; the exact results are given in Appendix A. According to the invariant operator theory, the wave functions of the unperturbed system are expressed in terms of ϕn(q, t) as [27, 28]
while is an arbitrary constant phase. Because ϕn(q, t) given in Appendix A are described by ρ(t), the wave functions are also represented in terms of it. In fact, the reason why ρ(t) has been introduced is that it is necessary in the expression of the quantum solutions of the unperturbed system of which Hamiltonian Ĥ is dependent on time.
We can easily check that the ladder operators yield
The quantum theory for a TDHS represented above will be used in the subsequent section in order to develop the perturbation theory of the system associated with the complex potentials.
3. Perturbation Theory With the Complex Potentials
Although perturbation theory does not give the spectrum of exact analytical quantum solutions for a dynamical system, it enables us to solve difficult problems in quantum mechanics from a series of routine calculations. If we regard that we cannot derive exact Schrödinger solutions in many cases of quantum mechanical problems, perturbation theory is quite useful in quantum mechanics. We will show how to manage perturbation theory for the case where the Hamiltonian involves time-dependent complex potentials on the basis of the associated theory represented in the previous section.
By expressing Equation (8) and its Hermitian adjoint â† inversely, we obtain the representation of canonical variables and in terms of â and â†. Then we can easily derive the expectation values of Ĥ using Equations (13), (15), and (16), leading to
The purpose of using perturbation techniques is to find the effects of small perturbed potentials on the whole system. In general, perturbation expansion gives a power series representation of a resultant quantity with respect to the perturbing parameter ϵ. The results of a calculation considering perturbations agree quite well with experimental data, but entails an infinite number of minor terms. However, the higher orders of such minor terms can, in general, be negligible from the perturbation corrections because their numeric scales are small enough.
Our theory is based on time-independent perturbation theory that is usually managed in terms of the eigenstates where the phases are not considered , as mentioned earlier. Hence, if we do not care about the minor phases of the wave functions contributed from the perturbing Hamiltonian, we can represent the overall wave functions of the system as
where 〈q|Φn〉 are eigenstates and θn(t) are the phases. According to the perturbation theory, the eigenstates take the form
where 〈q|ϕn〉 = ϕn(q, t), which are given in Appendix A, and
By perturbation theory, the EVH are represented as
Up until now, we have developed perturbation theory for the general time-dependent Hamiltonian system which involves complex potentials. In order to see the applicability of our theory to particular systems, let us consider the case where kR, 1(t)≠0 and kI, 3(t)≠0, while all other kR, j(t) and kI, l(t) are zero. In this case, Equation (3) can be written in the form
Now, from Equations (20)–(22), the corresponding wave functions are derived to be
where, is given by
Because in Equation (30) is a somewhat complicated form, we have represented it in Appendix B.
On the other hand, the EVH become
If large-scale perturbation expansions regarding higher order terms are required in the solutions, it is necessary to take the aid of computer algebra . The last term in Equation (35) is an imaginary part of the energy. We can elucidate the time-dependent evolution of the probability density in a quantum state based on non-Hermitian quantum theory, where imaginary terms of the energy levels take place [10, 31–34] as well as the real ones. The novel outcome of non-Hermitian analyses is excluded in the problem of conventional quantization methods relevant to Hermitian Hamiltonians. The eigenvalues of Hermitian Hamiltonians and/or Pseudo Hermitian Hamiltonians are always real from a mathematical point of view. However, some systems/subsystems can violate the requirement of Hermiticity of Hamiltonians . As can be seen from Equation (35), the expectation values of the Hamiltonian are complex. The possible existence of imaginary energy has a considerable interest from several decades ago . Zhu and Cukier  adopted an imaginary energy method in order to formulate quantum rate theories. They have obtained a rate constant that can be used to explain the WKB quantum tunneling rate at low temperature and the quantum state transition rate at high temperature. Further, they have shown that their quantum rate constant has the same structure as the Miller's theory [36, 37] of quantum transition state at high temperature limit. Meanwhile, imaginary energy gaps between the complex eigenvalues of a system have been used to analyze cooperative effects such as superconductivity and superradiance 
which can be obtained from the choice of time functions as f(t) = h−1(t) = eβt, g(t) = 0, and ω(t) = ω0, where β and ω0 are real positive constants. In addition, we choose
where a and b are real constants. Then, the solution of the nonlinear equation, Equation (6), is given by
provided that , where , and c1 and c2 are arbitrary real positive constants.
We have depicted the probability density for the system described above in Figure 1. Due to the choice of c1 = c2 = 1 in this figure, ρ(t) is a uniform function having the value 1. This figure shows that the deviation of the probability densities for the perturbed wave functions, Equation (19), at an initial time can be negligible and, as a consequence, they are almost the same as which are probabilities associated to the unperturbed wave functions. However, such deviation increases gradually as time goes by. We can also confirm from this figure that considered ϵ up to the second order are nearly the same as those considered ϵ up to the first order. This means that we can neglect the terms which are higher order in ϵ in the wave functions, as expected.
Figure 1. The probability densities (for the wave functions Equation 14) and (for the wave functions Equation 20) for the system described in Equations (37)–(39), where t is 0.5 for (A), 1.5 for (B), and 2.5 for (C). The characters 1st and 2nd in the figure legends mean that the corresponding figure has been plotted with the consideration up to the first order and up to the second order of ϵ, respectively. The parameters are chosen as ℏ = 1, β = 0.3, ω0 = 1, m0 = 1, n = 8, c1 = c2 = 1, ϵ = 0.0035, and a = b = 1.
Figure 2 shows the time evolution of the real and the imaginary parts of the expectation values of the Hamiltonian considering up to the second order of ϵ for the system described by Equations (36)–(38). From Figure 2A, we can confirm that the deviation of both the real and the imaginary parts of the EVH from those of the unperturbed Hamiltonian is, largely, not so significant, but is not zero. Figures 2B,C are enlarged graphs of the real and the imaginary parts of , respectively. If β is large, the expectation values of the real part of the Hamiltonian are also large. Real parts deviate toward the positive direction whereas imaginary parts toward the negative one and such deviations increase over time.
Figure 2. (A) The time evolutions of the real and the imaginary parts of the expectation values, Equation (36), for the system described in Equations (37)–(39), for three different values of β. The parameters are chosen as ℏ = 1, ω0 = 1, m0 = 1, n = 0, c1 = c2 = 1, ϵ = 0.0001, and a = b = 1. (B) Enlargement of the real parts in (A). The reference solid red line is the expectation value of the Hamiltonian without the perturbation term, given in Equation (18), under the choice of β = 0.5. (C) Enlargement of the imaginary parts in (A). The reference solid red line indicates zero of the vertical axis.
Let us show that the results, Equations (30)–(35), become well-known ones for a particular case. For this purpose, we put β = 0 from Equation (36) and b = 0 from Equation (37) while a ≠ 0. Then, the system becomes a simple one that is described by the Hamiltonian
where δ = ϵa. Then, under the choice c1 = c2 = 1, we have ρ = 1, W = Ω = ω0, G = 0, and . From these, we easily confirm that Equations (31) and (36) reduce in this situation to
with . The results Equations (40) and (41) are well-known in the literature [see, for example ]. If we choose c1≠1 and/or c2≠1, the quantum solutions for the unperturbed Hamiltonian Ĥ are somewhat different from the standard ones and the result given in Equation (41) becomes different; however, the quantum solutions for Ĥ still satisfy the Schrödinger equation in that case and the corresponding solutions approximated by the perturbation theory for the entire Hamiltonian are also allowed. Although we have chosen a simple reduced Hamiltonian, Equation (39), in order to show that our results become well-known ones for a particular case, the quantum solutions of the system described by this Hamiltonian can also be derived completely without using perturbation theory.
Now let us see another special case where the unperturbed Hamiltonian is given by
This corresponds to the case where the time functions are given by f(t) = h(t) = 1, g(t) = 0, and Meanwhile, we put the time functions in Equation (30) as
where c ≥ d. Then the solution of the nonlinear equation given in Equation (6) can be written in the form
provided that where ρ1, 0 and ρ2, 0 are real positive constants, Jν and Nν are Bessel functions of the first and second kind, respectively, while ν = 1/6.
In Figure 3, we have plotted the EVH for the system described in Equations (43)–(47) considering up to the second order of ϵ. The imaginary part of the EVH does not vary much over time, whereas the real part increases as time goes by. This shows that the imaginary part is small but not zero, while the real part does not significantly deviate from the unperturbed one. We have enlarged the real and the imaginary parts of the expectation value in Figures 3B,C, respectively. It is interesting that both the real and the imaginary parts of the EVH oscillate over time, provided that γ is sufficiently large. The period of such oscillation becomes short as the angular frequency γ increases. In the case of α → ∞, we have kR, 1(t) → 0, while the spectrum of the EVH is real and positive as conjectured by Bessis from numerical studies [see  and references therein].
Figure 3. (A) The time evolutions of the real and the imaginary parts of the expectation values, Equation (36), for the system described in Equations (44)–(48) for three different values of γ. The parameters are chosen as ℏ = 1, ω0 = 1, m0 = 1, n = 10, ρ1, 0 = ρ2, 0 = 1, α = 1, c = 1, d = 0.5, and ϵ = 0.001. (B) Enlargement of the real parts in (A). The reference solid red line is the expectation value of the unperturbed Hamiltonian, given in Equation (18). (C) Enlargement of the imaginary parts in (A). The reference solid red line indicates zero of the vertical axis.
The properties of the quantum states for perturbed Hamiltonian systems involving complex potentials were investigated. We have considered time-dependence of the imaginary part of the perturbing potentials as well as that of the real part, which are both coupled to a generalized harmonic oscillator through a weak coupling constant. The solutions of the Schrödinger equation of the system in the Fock state have been obtained using the perturbation theory. The effects of the perturbation on time behavior of the system have been analyzed. The perturbation corrections on the wave functions and on the EVH have been analytically investigated. Because we have considered complicated time-dependent Hamiltonian systems as a generalization, the EVH for each particular case is nonconservative. However, when we remove the time dependence of the Hamiltonian, our quantum solutions reduce to stationary states where the corresponding energy spectrum is conservative.
We have applied our theory to particular cases, such as the perturbed systems of which the unperturbed part is described by the Caldirola-Kanai Hamiltonian, and by a potential in which the angular frequency increases in proportion to the square of time. We see that the deviation of the probability densities and the EVH from those of the unperturbed system is, in large, not significant due to the weakness of the coupling. However, the deviation of the probability densities from those of the unperturbed system becomes large with the lapse of time. The imaginary part of the EVH is very small, but not zero, whereas the real part deviates slightly from that of the unperturbed system. In the case where the unperturbed part of the Hamiltonian corresponds to the Caldirola-Kanai oscillator, such deviation increases as the damping factor β grows. We have confirmed that the EVH for the system whose angular frequency is proportional to the square of t oscillate over time due to the time-variation of the chosen perturbing term in the Hamiltonian. The period of such oscillation becomes short as the angular frequency γ that appears in the perturbing term increases.
The development of the perturbation theory represented here may also be possibly applied to diverse different quantum systems which contain time-dependent complex perturbation potentials beyond those studied in the text. Some insights for characterizing complex potential systems are necessary in the light of the fact that a complex potential is receiving due attention in analyzing actual physical states, such as elastic scattering processes [5, 6], gain and loss in waveguides , band structures , and tunneling in absorptive media . In a unified description of resonance and decay phenomena in open quantum systems, imaginary energy implies decay widths of resonance states, while decay widths are inversely proportional to resonance lifetimes [31, 43]. On the other hand, in a description of quantum and/or classical waves propagating across a complex potential, the imaginary energy plays the role of gain which is the degree of amplification of the transmitted waves . Model-independent investigations of usual and unusual quantum features of complex potential systems may be a very interesting topic for subsequent research as a further generalization of conventional quantum mechanics. In addition, PT based study for a perturbation theory of time-dependent Hamiltonian systems may also be a good research topic as a next task.
Data Availability Statement
The datasets generated for this study are available on request to the corresponding author.
The author confirms being the sole contributor of this work and has approved it for publication.
This research was supported by the Basic Science Research Programs through the National Research Foundation of Korea (NRF) grants funded by the Korea government (Ministry of Education) (No. NRF-2016R1D1A1A09919503) and funded by the Korea government(MSIT) (No. NRF-2020R1F1A1075322).
Conflict of Interest
The author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
1. Mandal A, Hunt KLC. Adiabatic and nonadiabatic contributions to the energy of a system subject to a time-dependent perturbation: complete separation and physical interpretation. J Chem Phys. (2012) 137:164109. doi: 10.1063/1.4750045
3. Mandal A, Hunt KLC. Quantum transition probabilities during a perturbing pulse: Differences between the nonadiabatic results and Fermi's golden rule forms. J Chem Phys. (2018) 148:194107. doi: 10.1063/1.5019172
7. Zhang C, Bai R, Gu X, Jin XR, Zhang YQ, Lee YP. Dual-band unidirectional reflectionless phenomena in an ultracompact non-Hermitian plasmonic waveguide system based on near-field coupling. Opt Exp. (2017) 25:24281. doi: 10.1364/OE.25.024281
27. Lewis HR Jr, Riesenfeld WB. An exact quantum theory of the time-dependent harmonic oscillator and of a charged particle in a time-dependent electromagnetic field. J Math Phys. (1969) 10:1458–73. doi: 10.1063/1.1664991
31. Tellander F, Berggren KF. Spectra, current flow, and wave-function morphology in a model PT-symmetric quantum dot with external interactions. Phys Rev A. (2017) 95:042115. doi: 10.1103/PhysRevA.95.042115
32. Chávez NC, Mattiotti F, Méndez-Bermúdez JA, Borgonovi F, Celardo GL. Real and imaginary energy gaps: a comparison between single excitation superradiance and superconductivity and robustness to disorder. Eur Phys J B. (2019) 92:1–12. doi: 10.1140/epjb/e2019-100016-3
Appendix A: Solutions of Equation (12)
The eigenstates of the invariant operator for the system are obtained by directly solving Equation (12) or by using the properties of Equations (15) and (16) [28, 44, 45]. In our case, the results are given by
Appendix B: The Formula of in Equation (30)
is obtained by evaluating the summation Since Π2 given in Equation (22) is composed of two terms, let us express in the form
The second term is zero because the expectation values 〈ϕn|Ĥp|ϕn〉 which are involved in it are zero. On the other hand, a rigorous evaluation of the first term using Equation (23) gives
Keywords: perturbation theory, time-dependent Hamiltonian system, complex potential, Schrödinger equation, wave function, expectation value
Citation: Choi JR (2020) Perturbation Theory for Time-Dependent Quantum Systems Involving Complex Potentials. Front. Phys. 8:189. doi: 10.3389/fphy.2020.00189
Received: 04 February 2020; Accepted: 30 April 2020;
Published: 05 June 2020.
Edited by:Francisco Gomez, Centro Nacional de Investigación y Desarrollo Tecnológico, Mexico
Reviewed by:Silvestro Fassari, Università degli Studi Guglielmo Marconi, Italy
Fabio Rinaldi, Università degli Studi Guglielmo Marconi, Italy
Copyright © 2020 Choi. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jeong Ryeol Choi, firstname.lastname@example.org