Equivalence between free and harmonically trapped quantum particles
It is shown that general solutions of the free-particle Schrödinger equation can be mapped onto solutions of the Schrödinger equation for the harmonic oscillator. This is done in such a way that the time evolution of a free particle subjected to a sudden transition to a harmonic potential can be described by a simple coordinate transformation applied at the transition time. This procedure is computationally more efficient than either state-projection or propagator techniques. A concatenation of the map and its inverse allows us to map from one harmonic oscillator to another with a different spring constant.
pacs:03.65.Ca, 03.65.Db, 03.65.Fd, 03.65.Ge
It is known that laser beam modes in the paraxial approximation are equivalent to the eigenmodes of the harmonic oscillator Nienhuis and Allen (1993); Steuernagel (2005). The corresponding (one-dimensional, scalar) wave equation
arises as an approximation to Maxwell’s equations Haus (2000). In this optical context the parameter is proportional to the (stationary, monochromatic beam’s) wave number or inversely proportional to the light’s wavelength .
The Schrödinger equation for a free particle in one spatial dimension is given by
where the parameter contains the particle mass and Planck’s constant , namely . It is obviously equivalent to eq. (1) above.
Since the wave functions for free particles and those subjected to harmonic potentials factorize with respect to their spatial coordinates we will only discuss the one-dimensional case.
It has recently been realized that a simple coordinate transformation Steuernagel (2005) maps paraxial beams onto two-dimensional harmonic oscillator wave functions. Clearly, the same coordinate transformation can be used for the mapping of a free onto a harmonically trapped particle; in a more general setting this was noted before Takagi (1990); Bluman and Shtelen (1996). This case has added meaning since a non-adiabatic physical transition constitutes an experimentally implementable sequence of changing environmental conditions. It turns out that modelling such a transition using the standard eigenstate-projection or quantum propagator techniques is more cumbersome. State projection typically leads to infinite sums over eigenfunctions which have to be truncated and are difficult to simplify, propagator techniques involve non-trivial integrals. The results reported here may therefore not only be of fundamental but also of technical interest.
It is noteworthy that the discontinuities of the mapping of the time-coordinate in eq. (3) maps from a Schrödinger equation with a continuous to another with a discrete energy spectrum.
ii.1 Map from free to trapped case
The coordinate transformations Takagi (1990)
applied to the wave function mapping
yield the solution for the Schrödinger equation
of a harmonic oscillator with mass , spring constant and resonance frequency .
Although this transformation maps onto , compare with the Gouy-phase of optics Steuernagel (2005), the periodicity arising through the use of trigonometric functions meaningfully represents the oscillator’s motion for all times .
In beam optics the confocal parameter Pampaloni and Enderlein (2004) parameterizes the strength of the beam’s focussing and the curvature of its hyperbolic flow lines Steuernagel (2005). Here, it serves as a rescaling parameter of the transverse coordinate transformation and thus allows us to stretch or squeeze the width of the wave packet we want to map.
In order to preserve the wave function normalization we have to determine the spatial coordinate stretching at the mapping time . This yields the normalization factor which the wave function has to be multiplied with. The specific choice
yields the natural mapping and which we will use from now on.
ii.2 Inverse map (trapped to free)
ii.3 Concatenated maps (trapped to trapped)
The concatenated coordinate transformations from an initial harmonic trapping potential with spring constant and wave function , via the free particle-case , to a final harmonic potential with spring constant and wave function is given by
Here , and solves Schrödinger eq. (6) with substituted by .
Iii An Example
A well-known textbook example is the freely evolving Gaussian wave-packet with initial position spread
Here parameterizes spatial, and momentum displacement of the wave functions. If either of these two quantities is non-zero the mapping onto a harmonically trapped state results in a state with oscillating center-of-mass. In general, although a wave function of Gaussian shape, this freely evolving wave packet will also not ‘fit’ the width of the harmonic potential and therefore be squeezed. In short, the state of eq. (13) trapped in a harmonic potential becomes a squeezed coherent state Schleich (2001). For a superposition of such states see Fig. 1 and for this superposition state being trapped consult Fig. 2.
It is shown that general solutions of the free-particle Schrödinger equation can be mapped onto solutions of the Schrödinger equation for the harmonic oscillator using a simple coordinate transformation in conjunction with a multiplication of the wave function by a suitable phase factor. This map is invertible and a concatenation of two such maps allows us to map from one harmonic oscillator to another with a different spring constant. The simplicity of the approach described here makes it a tool of choice for the description of the wave function of a particle experiencing instantaneous transitions from a free to a harmonically trapped state, the sudden release from a harmonic trap, or the sudden change of the strength of its harmonic trapping potential.
The mapping introduced here is computationally more efficient than state-projection or propagator techniques and conceptually simpler than mapping techniques such as those used for supersymmetric potentials Gendenshteïn (1983); Cooper et al. (1989). It may even help with numerical calculations because it allows for the determination of the behaviour of wave functions trapped in a harmonic potential using the lower computational overheads of wave functions in free space.
- Nienhuis and Allen (1993) G. Nienhuis and L. Allen, Phys. Rev. A 48, 656 (1993).
- Steuernagel (2005) O. Steuernagel, Am. J. Phys. 73, 625 (2005), eprint arXiv:physics/0312116v2.
- Haus (2000) H. A. Haus, Electromagnetic Noise and Quantum Optical Measurements (Springer, Heidelberg, 2000).
- Takagi (1990) S. Takagi, Prog. Theor. Phys. 84, 1019 (1990).
- Bluman and Shtelen (1996) G. Bluman and V. Shtelen, J. Phys. A: Math. Gen. 29, 4473 (1996).
- Pampaloni and Enderlein (2004) F. Pampaloni and J. Enderlein (2004), eprint arXiv:physics/0410021.
- Schleich (2001) W. P. Schleich, Quantum Optics in Phase Space (Wiley-VCH, Berlin, 2001).
- Gendenshteïn (1983) L. É. Gendenshteïn, Sov. J. Exp. Theor. Phys. Lett. 38, 356 (1983).
- Cooper et al. (1989) F. Cooper, J. N. Ginocchio, and A. Wipf, J. Phys. A: Math. Gen. 22, 3707 (1989).