Quantum Analysis on Time Behavior of a Lengthening Pendulum
Quantum properties of a lengthening pendulum are studied under the assumption that the length of the string increases at a steady rate. Advanced analysis for various physical problems in several types of quantum states, such as propagators, Wigner distribution functions, energy eigenvalues, probability densities, and dispersions of physical quantities, is carried out using quantum wave functions of the system. In particular, the time behavior of Gaussian-type wave packets is investigated in detail. The probability density for a Gaussian wave packet displaced in the positive at oscillates back and forth from the center (). This phenomenon is very similar to the classical motion of the pendulum. As a consequence, we can confirm that there is a correspondence between its quantum and classical behaviors. When we analyze a dynamical system in view of quantum mechanics, the quantum and classical correspondence is very important in order for the associated quantum theory to be valid and viable.
If the Hamiltonian of a system is dependent on time, it is classified as a time-dependent Hamiltonian system (TDHS). Typical examples of TDHSs include a damped harmonic oscillator [1–3], a harmonic oscillator driven by a time-varying force [4, 5], and a harmonic oscillator with time-dependent parameters such as time-dependent mass [6–9] and/or frequency [10–12]. A powerful method for solving quantum solutions of a TDHS is the invariant operator method [13–17].
Among various types of TDHSs, we study quantum dynamical properties of a lengthening pendulum [18–23] in this paper. This kind of pendulum is an interesting topic as a fundamental nonstationary mechanical system and its study may enable us to acquire elementary methods for manipulating more complicated general dynamical systems. The derivation of quantum mechanical solutions of such time-varying swinging objects and the analysis of them with high precision demands exquisite theory of quantum physics and chaos based on advanced mathematical techniques.
Much attention has been paid to the classical and quantum problems of the lengthening pendulum up to now. The time variation of angular and linear amplitudes has been analyzed by Brearley for a simple lengthening pendulum under the assumption that the amplitude of the oscillation is small . The characteristics of quantum wave functions in Fock state for the lengthening pendulum have been investigated by Um et al. . The research group of McMillan et al.  studied a parametric amplification of its oscillatory motion and the physical mechanism associated with the energy transfer to the pendulum via a time-varying parameter (the length of pendulum) under the assumption that the lengthening pendulum is driven by an arbitrary force. Nevertheless, as far as we know, the exact analytical solutions have not been obtained yet due to the difficulty in mathematical treatment of the differential equation describing the motion of the pendulum, which involves a sinusoidal function. The derivation of quantum solutions of the system with higher accuracy than what we can do by hand requires the aid of computer programs such as MATLAB .
The pendulum or the lengthening thereof not only is a good mathematical model for practicing a nontrivial problem in mechanics but also can possibly be applied to many actual systems. A gravitational pendulum which has been used as a seismic sensor of translation and rotation is one of its well-known applications . Usually seismologists use a kind of sensing device (pendulum) in order to detect earthquake waves, which is more complicated than the simple pendulum. The theory of the lengthening pendulum can also be applied to analyzing a drifting spiral motion of a charged particle in mirror-confined plasma . In this case, a position-dependent magnetic field, which is a dense and weak magnetic field, is responsible for the change of the rotational radius along its spiral motion.
We will study in this work several quantum states of the lengthening pendulum beyond the Fock state. We will compare quantum motions of the pendulum with its classical behaviors. Quantized energy and its dispersion will be investigated via a rigorous development of the quantum theory associated with the pendulum.
This paper is organized as follows. The fundamental mechanics of the lengthening pendulum is surveyed in Section 2 and various quantum mechanical properties of the system are analyzed in the subsequent sections. Section 3 is devoted to the study of the time behavior of the quantized energy. Path integral formulation of quantum mechanics for the system is treated in Section 4. Wigner distribution function (WDF) of the system is investigated in Section 5. Quantum properties of the superposition of two neighboring states are discussed in Section 6. The correspondence between classical and quantum mechanics for the pendulum is demonstrated in Section 7 through Gaussian wave packet description of the system. The concluding remarks are given in Section 8 which is the last section.
2. Fundamentals of the Lengthening Pendulum
Now we are carrying out a survey of the fundamentals of the lengthening pendulum. Let us consider a pendulum of a relatively massive object hung by a vertical string from a fixed ceiling and assume that it swings from its fixed equilibrium position. Suppose that the length of the pendulum increases with a constant rate from its initial length . In case is sufficiently small, the pendulum undergoes an adiabatic change and such case will only be considered in this review.
We do not regard the quantum effects of the -component for the sake of simplicity. Then, the Hamiltonian considered up to the second order in is given by  where . A classical solution of the system can be represented as  Here, and are arbitrary constants and where and are Bessel functions of the first kind and the second kind, respectively, and . A more general solution for this system is given in .
The construction of the annihilation and creation operators is useful for developing a quantum theory of a TDHS. According to the invariant operator theory , the formulae of them are somewhat different from those of the simple harmonic oscillator. For instance, the annihilation operator is given by  Of course, the Hermitian adjoint, , of the above equation plays the role of the creation operator. These operators satisfy the boson commutation relation, such that .
If we use the property of this commutation relation, we can easily derive the wave functions in Fock state, which are represented in terms of some time functions and time-dependent phases [21, 22]:where the quantum number is given by . (All subsequent will follow this convention unless we particularly specify it.) The explicit forms of and are given by where is th-order Hermite polynomial and with The wave functions given in (6) are the basic information necessary for studying quantum features of the system. We will investigate quantum characteristics of the pendulum for more complicated cases in the next section.
3. Quantized Energy
Bohr tried to merge quantum and classical mechanics by introducing a correspondence principle between them. Even though the results of quantum and classical descriptions for a system more or less overlap under particular limits, their underlying principles are quite different from each other. There are significant differences between consequences of the two theories when the quantum number is sufficiently low, whereas those with a higher quantum number look very similar. It is well known that we can compute the probability for finding the pendulum at a particular angle from the square of the given wave function. In cases of higher energy states, the most probable angle of the pendulum will shift away from the center () and the intervals between any two adjacent peaks in the graph of probability densities in Fock state become narrow. Due to this trend associated with the probability density, the behavior of the quantum pendulum in a high energy limit may look like that of the counterpart classical one.
Another difference between them is that while the classical probability density function is confined within the two classical turning points, the quantum probabilities extend beyond the classically allowed angle. For more detailed investigations about this consequence, we consider the quantum and the classical energies of the system that are given by Let us evaluate the probability that the pendulum stays outside the classically allowed region. If we denote the amplitude of angle as , the maximum classical potential energy is given by . Then, regarding the classical turning point of the pendulum associated with the ground state wave packet, we equate the ground state energy with the maximum classical potential energy: Thus, by inserting (10) with into the left-hand side of this equation, the classical amplitude can be evaluated to be The probability density outside the classically allowed region is given by Here, factor 2 is multiplied considering the equal contribution of the probability from outside of the opposite turning angle, . In addition, we have supposed that the oscillating amplitude is sufficiently small in order to replace the integral interval in (14) with . This extension of the scope of integration will also be equally applied in further similar calculations. By substituting (6) with into (14), we have This is the probability that the quantum pendulum remains within the classically forbidden regions.
In fact, the description of quantized energy for a TDHS, such as (10), is a delicate problem. While some authors [1, 3, 11, 21, 22] including ours have investigated quantized energy levels for specific TDHSs, there is another opinion [26, 27] that such energy levels do not exist for the case of TDHSs. Hence, explicit demonstrations of the existence of quantized energy levels may be an interesting research topic for further study in the future in this field.
The path integral formulation of quantum mechanics was found by Feynman  as an alternative quantum description of dynamical systems. Historically, this achievement was partly inspired by Dirac’s idea of a quantum mechanical description . The techniques of path integrals continued to be advanced until now, leading to providing many useful tools for solving problems in quantum mechanics. One of the efficient uses of the path integral method is to evaluate a propagator. This achieved great success in both quantum mechanics and quantum field theory, because a propagator involves all the information of the quantum system. Through the analysis of the propagator, we can confirm how to propagate quantum wave packets from an initial angle and time to a final angle and time .
In terms of wave functions, the propagator is given by Let us calculate the complete analytical form of the propagator of the system using this definition. By inserting (6) with (7), we have From the use of Mehler’s formula it becomes We can use this result in estimating the evolution of quantum states under the given Hamiltonian of the pendulum and in studying the dynamics of the pendulum with varying length. More precisely, the propagator enables us to know how the initial state evolves to an arbitrary later state via the following equation The propagator plays an important role in modern physics.
5. Wigner Distribution Function
The investigation of the probability distribution is useful for understanding the novel features of quantum mechanics. The WDF  may be the most useful one among various distribution functions that are needed when we study quantum statistical mechanics. Many quantum mechanical characteristics of a dynamical system, which are absent in classical mechanics, can be derived from WDF.
The WDFs in phase space are given by We replace the range of the integration with for the same reason as that of the previous integration in (14). Using the integral formula of the form  we can evaluate (21), leading to where This can be used to study the probability distribution for and of the pendulum within a quantum mechanical view point. Figure 1 shows the WDF for several values of . Exact marginal distributions can be achievable with the use of a WDF, whereas it cannot be done by means of the Husimi distribution function . It is allowed that the WDF takes negative values as well as positive ones in phase space. The negativity of the WDF is, in fact, one of the primary nonclassical effects of quantum mechanics [34, 35]. The WDF, given in (23), enables us to study the quantum statistical properties of the system in configuration space. Some imaginary distribution functions that can be used in a wide branch of physics were also proposed in the literature by other authors [36–40]. The width of in the angle space becomes small while that in the angular momentum space becomes large as time goes by.
6. Superposition of Two Neighboring States
The general solution of the Schrödinger equation is a linear combination of separated solutions given in (6): where . To study the energy spectrum in the state, it is necessary to evaluate the eigenvalue equation of the Hamiltonian, which is
It is well known that quantum mechanics predicts only the probability for happening of particular results among a large number of mutually independent possible outcomes in physical systems. An arbitrary state function representing any dynamical state can be expanded in terms of Fock state wave functions. As a first example, let us consider the case that (25) is a superposition of two neighboring states: where phases and are arbitrary real constants. The expectation values of canonical variables in this state can be easily evaluated to be  From these equations, we can confirm that the expectation values of and oscillate with time.
Through the similar procedure, the expectation value of the energy operator can also be obtained as This is the same as the average energy of the th state and th state energies with the same proportion. This energy corresponds to that of an instantaneous eigenstate and varies with time like a classical one. However, this reduces to a constant value when , which is associated with the results of standard quantum mechanics for the simple pendulum represented in many text books in this field.
7. Gaussian Wave Packets
As a next example, we consider the case that the initial state function is given by where is some constant which is real. In view of quantum theory, the probabilistic predictions are always obtained by squaring the modulus of a probability amplitude. The probability amplitudes that the oscillator stays in th eigenstates are Then, using (6), (30), and (31) after replacing the interval of the integration with , we obtain that for . On the other hand, for . Hence, odd order amplitudes do not contribute to the state function.
One may also possibly consider other types of wave packets described as a superposition of Fock state wave functions. Now, we are going to choose a state which seems appropriate for studying whether there is a quantum and classical correspondence concerning the behavior of the lengthening pendulum. At the early era of quantum mechanics, Schrödinger tried to explain the relation between quantum and classical mechanics by introducing coherent wave packets . In order to find the possibility for the existence of a quantum solution relevant to a particular wave packet whose center oscillates with the same period as the classical motion, let us consider a state function that is expressed at as where is a real constant with the dimension of angle. This is the same as (30) except that is replaced by and the center of the initial wave packet is displaced by an amount . In this case, the probability amplitudes can also be obtained using the same method as that of the first case (the nondisplaced Gaussian case) performed with (31), such that where . Here we have supposed that the increase of the length of string is sufficiently slow. With the help of (6) and (34) the time evolution of the state function (25) becomes The corresponding probability density is obtained by squaring the above equation. Thus, a direct evaluation gives where . The wave packet equation (36) is a time-dependent Gaussian shape. We depicted this wave packet in Figure 2. The wave packet not only converges to the center () but also oscillates back and forth around the center as time goes by. The amplitude of such oscillation is determined depending on the magnitude of . This quantum behavior is very similar to the classical motion of the pendulum. The width of the wave packet is . The probabilities for the contribution of th quantum state to the state function are . Thus, using (34), we have These probabilities do not vary with time. The mean value of quantum numbers which contribute to the state function is . Using (37), we see that this can be evaluated to be Thus, for the case that is large, the contribution of high energy states is also large.
The dispersion of is given by . This can be easily calculated as Figure 3(a) reveals that increases as grows.
We can also evaluate the dispersion of energy from . Through a mathematical procedure similar to that of the calculation of , we have From Figure 3(b), we can confirm that increases as grows. However, decreases with time.
On the basis of the Schrödinger solutions that can be obtained by taking advantage of the ladder operators and , various quantum states of the lengthening pendulum have been investigated. We obtained the energy eigenvalues, propagator, WDFs, uncertainties, and probability densities in several types of quantum states. The quantum behavior of the pendulum was analyzed rigorously and compared to that of the classical one.
The propagator given in (19) entirely determines the quantum behavior of the lengthening pendulum in a gravitation. If we know the wave function at an initial angle and time, it is possible to estimate probability amplitude for another angle at a subsequent time using the propagator. Wigner distribution function of the system was derived as shown in (23). This can be a negative value in some regions of phase space as well as a positive one, which signifies a novel nonclassical feature of the quantum system. We can confirm from Figure 1 that the width of associated with decreases as time goes by while that associated with increases. This appearance is due to the fact that the uncertainty in decreases with time while that in increases. For the variation of the uncertainty for and , you can refer to .
As the length of string increases, mechanical energy cannot be conserved because some amount of work is extracted from the system . This is the reason why the amplitude decreases with time, not only for the classical pendulum but for the quantum one as well.
From (28), we see that the expectation values of canonical variables in the superposition of the two neighboring states oscillate with time. The quantum energy in this case varies with time as shown in (29). For the state that is described by (30), the odd order probability amplitudes disappear since (30) is an even function about .
To investigate the existence of the quantum and classical correspondence for the behavior of the lengthening pendulum, we have supposed that the form of the state function at is given by (33). The initial angle in this case is displaced toward a positive by an amount . The resulting probability density is represented in (36) which is a Gaussian shape. From Figure 2, we can confirm that this probability density oscillates back and forth from the center depending on the magnitude of , while its width becomes narrower over time as the quantum energy decays. So to speak, its quantum behavior is very similar to that of the classical oscillation of the pendulum. Therefore, the correspondence between quantum and classical behaviors holds.
Indeed, the quantum and classical correspondence for a certain quantum theory is a crucial requirement for the validity of the associated quantum theory. The invariant operator theory we used here is sound and enables us to obtain reasonable quantum solutions for a time-dependent dynamical system. Further, we know that it sometimes gives the exact quantum solutions for a system of which the classical solutions are completely known.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
This research was supported by the Basic Science Research Program of the year 2015 through the National Research Foundation of Korea (NRF) funded by the Ministry of Education (Grant no. NRF-2013R1A1A2062907). Some parts of this work were written and developed on the basis of Chap. 2 of doctoral thesis of Choi .
D. M. Greenberger, “A new approach to the problem of dissipation in quantum mechanics,” Journal of Mathematical Physics, vol. 20, no. 5, pp. 771–780, 1979.View at: Google Scholar
V. V. Dodonov and V. I. Man'ko, “Invariants and evolution of nonstationary quantum systems,” in Proceedings of Lebedev Physical Institute 183, M. A. Markov, Ed., Nova Science Publishers, Commack, NY, USA, 1989.View at: Google Scholar
M. L. Boas, Mathematical Methods in the Physical Sciences, John Wiley & Sons, New York, NY, USA, 2nd edition, 1983.
L. C. Andrews, Special Function of Mathematics for Engineers, SPIE Press, Oxford, UK, 2nd edition, 1998.
C. I. Um, J. R. Choi, K. H. Yeon, S. Zhang, and T. F. George, “Exact quantum theory of a lengthening pendulum,” Journal of the Korean Physical Society, vol. 41, no. 5, pp. 649–654, 2002.View at: Google Scholar
J. R. Choi, J. N. Song, and S. J. Hong, “Spectrum of quantized energy for a lengthening pendulum,” AIP Conference Proceedings, vol. 1281, no. 1, pp. 594–598, 2010.View at: Google Scholar
F. F. Chen, Introduction to Plasma Physics and Controled Fusion, vol. 1, Springer, New York, NY, USA, 2006.
R. W. Hasse, “On the quantum mechanical treatment of dissipative systems,” Journal of Mathematical Physics, vol. 16, no. 10, pp. 2005–2011, 1975.View at: Google Scholar
P. A. M. Dirac, “The Lagrangian in quantum mechanics,” Physikalische Zeitschrift der Sowjetunion, vol. 3, no. 1, pp. 64–72, 1933.View at: Google Scholar
A. Erdély, Higher Transcendental Functions, vol. II, McGraw-Hill, New York, NY, USA, 1953.
I. S. Gradshteyn and M. Ryzhik, Table of Integrals, Series and Products, Academic Press, New York, NY, USA, 1980.
K. Husimi, “Some formal properties of the density matrix,” Proceedings of the Physico-Mathematical Society of Japan, vol. 22, no. 5, pp. 264–314, 1940.View at: Google Scholar
S. Khademi and S. Nasiri, “Gauge invariant quantization of dissipative systems of charged particles in extended phase space,” Iranian Journal of Science and Technology, Transaction A: Science, vol. 26, no. 1, 2002.View at: Google Scholar
J. R. Choi, Quantum mechanical treatment of time-dependent hamiltonians [Ph.D. thesis], Korea University, Seoul, South Korea, 2001.
L. I. Schiff, Quantum Mechanics, McGraw-Hill, New York, NY, USA, 1968.