Trilobites revived: fragile Rydberg molecules, Coulomb Green’s function, Lambert’s theorem

The trilobite state
The trilobite Rydberg molecule can be modeled by the Coulomb Green’s function, which represents the quantized version of Lambert’s orbit determination problem.

The recent experimental realization observation of giant Rydberg molecules by Bendkowsky, Butscher, Nipper, Shaffer, Löw, Pfau [theoretically studied by Greene and coworkers, see for example Phys. Rev. Lett. 85, 2458 (2000)] shows Coulombic forces at work at large atomic distances to form a fragile molecule. The simplest approach to Rydberg molecules employs the Fermi contact potential (also called zero range potential), where the Coulomb Green’s function plays a central role. The quantum mechanical expression for the Coulomb Green’s function was derived in position space by Hostler and in momentum space by Schwinger. The quantum mechanical expression does not provide immediate insights into the peculiar nodal structure shown on the left side and thus it is time again to look for a semiclassical interpretation, which requires to translate an astronomical theorem into the Schrödinger world, one of my favorite topics.

Johann Heinrich Lambert was a true “Universalgelehrter”, exchanging letters with Kant about philosophy, devising a new color pyramid, proving that π is an irrational number, and doing physics. His career did not proceed without difficulties since he had to educate himself after working hours in his father’s tailor shop. After a long journey Lambert ended up in Berlin at the academy (and Euler choose to “escape” to St. Petersburg).

Lambert followed Kepler’s footsteps and tackled one of the most challenging problems of the time: the determination of celestial orbits from observations. In 1761 Lambert did solve the problem of orbit determination from two positions measurements. Lambert’s Theorem is a cornerstone of astronavigation (see for example the determination of Sputnik’s orbit using radar range measurements and Lambert’s theorem). Orbit determination from angular information alone (without known distances) is another problem and requires more observations.

Lambert poses the following question [Insigniores orbitae cometarum proprietates (Augsburg, 1761), p. 120, Lemma XXV, Problema XL]: Data longitudine axis maioris & situ foci F nec non situ punctorum N, M, construere ellipsin [Given the length of the semi-major axis, the location of one focal point, the points N,M, construct the two possible elliptical orbits connecting both points.]

Lambert's construction of two ellipses.
Lambert’s construction to find all possible trajectories from N to M and to map them to a ficticious 1D motion from n to m.

Lambert finds the two elliptic orbits [Fig. XXI] with an ingenious construction: he maps the rather complicated two-dimensional problem to the fictitious motion along a degenerate linear ellipse. Some physicists may know how to relate the three-dimensional Kepler problem to a four-dimensional oscillator via the Kustaanheimo–Stiefel transformation [see for example The harmonic oscillator in modern physics by Moshinsky and Smirnov]. But Lambert’s quite different procedure has its advantages for constructing the semiclassical Coulomb Green’s function, as we will see in a moment.

Shown are two ellipses with the same lengths of the semimajor axes 1/2 A1B1=1/2 A2 B2 and a common focus located at F. The centers of the two ellipses are denoted by C1 and C2. Lambert’s lemma allows to relate the motion from N to M on both ellipses to a common collinear motion on the degenerate linear ellipse Fb, where the points n and m are chosen such that the time of flight (TOF) along nm equals the TOF
along the elliptical arc NM on the first ellipse. On the second ellipse the TOF along the arc NB2M equals the TOF along nbm. The points n and m are found by marking the point G halfway between N and M. Then the major axis Fb=A1 B1=A2 B2 of the linear ellipse is drawn starting at F and running through G. On this line the point g is placed at the distance Fg=1/2(FN+FM). Finally n and m are given by the intersection points of a circle around g with radius GN=GM. This construction shows that the sum of the lengths of the shaded triangle α±=FN + FM ± NM is equal to α±=fn+ fm ± nm. The travel time depends only on the distances entering α±, and all calculations of the travel times etc. are given by one-dimensional integrations along the ficticious linear ellipse.

Lambert did find all the four possible trajectories from N to M which have the same energy (=semimajor axis a), regardless of their eccentricity (=angular momentum). The elimination of the angular momentum from Kepler’s equation is a tremendous achievement and the expression for the action is converted from Kepler’s form

  • [Kepler] W(r,r‘;E)=√μ a Kc [ξ + ε sin(ξ) – ξ’ – ε sin(ξ’)], with eccentricity ε, eccentric anomaly ξ to
  • [Lambert] W(r,r‘;E)=√μ a Kc[γ + sin(γ) – δ – sin(δ)], with
    sin2(γ/2)=(r+r’+ |r‘-r|)/(4a) and sin2(δ/2)=(r+r’- |r‘-r|)/(4a).

The derivation is also discussed in detail in our paper [Kanellopoulos, Kleber, Kramer: Use of Lambert’s Theorem for the n-Dimensional Coulomb Problem Phys. Rev. A, 80, 012101 (2009), free arxiv version here]. The Coulomb problem of the hydrogen atom is equivalent to the gravitational Kepler problem, since both are subject to a 1/r potential. Some readers might have seen the equation for the action in Gutzwiller’s nice book Chaos in classical and quantum mechanics, eq. (1.14). It is worthwhile to point out that the series solution given by Lambert (and Gutzwiller) for the time of flight can be summed up easily and is denoted today by an inverse sine function (for hyperbolic motion a hyperbolic sine, a function later introduced by Riccati and Lambert). Again, the key-point is the introduction of the linear ficticious ellipse by Lambert which avoids integrating along elliptical arcs.

The surprising conclusion: the nodal pattern of the hydrogen atom can be viewed as resulting from a double-slit interference along two principal ellipses. The interference determines the eigenenergies and the eigenstates. Even the notorious difficult-to-calculate van Vleck-Pauli-Morette (VVPM) determinant can be expressed in short closed form with the help of Lambert’s theorem and our result works even in higher dimensions. The analytic form of the action and the VVPM determinant becomes essential for our continuation of the classical action into the forbidden region, which corresponds to a tunneling process, see the last part of our paper.

Lambert is definitely a very fascinating person. Wouldn’t it be nice to discuss with him about philosophy, life, and science?


2 thoughts on “Trilobites revived: fragile Rydberg molecules, Coulomb Green’s function, Lambert’s theorem”

  1. I’m doubly impressed by the fact that you not only managed to evaluate the Van Vleck determinant analytically, but also that you did it using an idea from the 17th century! I remember feeling very frustrated by the complexity of the Gutzwiller method when we attempted this. I also appreciate you extension of the semiclassical method into the classically forbidden region, which always felt like an ugly aspect of the numerical approach. I’m reassured that the conclusion of the analytical work confirms that the wavefunction can still be understood in terms of a simple two-path interference, something we always felt must reflect some simplicity of the Coulomb potential that we didn’t explore in our paper.

    1. Your work in Phys. Rev. A 64, 042508 (2001) Quantum and semiclassical analysis of long-range Rydberg molecules together with B Granger and C Greene was a great inspiration for us! It is indeed a relief to see that just two paths are enough to generate the complex trilobite pattern of the Coulomb Green’s function at all locations. The Coulomb Green’s function is also relevant for the photoionization of atoms.

      Interestingly the superposition of Coulomb and homogeneous electric field is also integrable and related to another astronautical problem already attacked by Lagrange, see the beautiful book “Essays on the Motion of Celestial Bodies” by Vladimir V. Beletsky, Birkhäuser, Basel, 2001 and the application to the semiclassical theory of photoionization in Phys. Rev. A 68, 012709 (2003) by Ch. Bordas and F. Lépine, C. Nicole, and M. J. J. Vrakking.

Leave a Reply

Fill in your details below or click an icon to log in: Logo

You are commenting using your account. Log Out / Change )

Twitter picture

You are commenting using your Twitter account. Log Out / Change )

Facebook photo

You are commenting using your Facebook account. Log Out / Change )

Google+ photo

You are commenting using your Google+ account. Log Out / Change )

Connecting to %s