Parabolic Coordinates

In the world of mathematics and physics, the Cartesian coordinate system offers a simple and powerful grid of straight lines to describe space. However, many physical phenomena, from the trajectory of a projectile to the fields around a charged object, do not conform to this rigid rectilinear structure. Forcing problems with inherent curves into a square framework can lead to unnecessary complexity, obscuring the underlying elegance of the physics. This raises a fundamental question: how can we choose a descriptive language that matches the natural geometry of a problem?

This article explores a powerful alternative: parabolic coordinates. By redrawing our map of space using intersecting parabolas, we can unlock simpler solutions to otherwise challenging problems. We will begin by dissecting the core principles and mathematical machinery of this elegant system, understanding how to measure and describe motion in a world without straight grid lines. Following this, we will journey through its diverse applications, revealing how parabolic coordinates provide crucial insights in fields ranging from classical mechanics and electromagnetism to the quantum mechanical analysis of the Stark effect.

Imagine you have a perfectly flat sheet of paper. The most straightforward way to tell someone where a point is on that sheet is to use a simple rectangular grid, the familiar Cartesian coordinates $(x, y)$. You go $x$ units across and $y$ units up. The grid lines are straight, parallel, and evenly spaced. It’s wonderfully simple. But what if the problem you're trying to solve—say, the path of a comet or the shape of an electric field—doesn't respect this rigid grid? What if the problem has a natural "parabolic" character? Then, forcing it into a square box is like trying to fit a round peg in a square hole. It’s possible, but it’s awkward.

This is where the power of different coordinate systems comes in. We can redraw our map of the plane using new lines that are better suited to the problem at hand. Parabolic coordinates are one such redrawing.

Let's throw away the rigid grid of $x$ and $y$ and define a new set of coordinates, which we'll call $\sigma$ and $\tau$. We relate them to our old Cartesian friends with a pair of transformation equations:

What does this grid look like? If we hold $\sigma$ constant (say, $\sigma=1$) and let $\tau$ vary, the equations trace out a parabola that opens upwards. If we pick a different constant $\sigma$ (say, $\sigma=2$), we get another, wider parabola. On the other hand, if we hold $\tau$ constant and let $\sigma$ vary, we trace out parabolas that open downwards.

The result is a beautiful and elegant grid composed entirely of two families of intersecting parabolas. Our flat plane is now tessellated not by squares, but by a graceful web of curves. The key question is: how do we do physics in this new world? How do we measure distances, describe motion, and calculate forces?

In the Cartesian world, the "direction" of $x$ and the "direction" of $y$ are the same everywhere. Not so in our parabolic grid. The "direction of $\sigma$" is the direction you'd move if you walked along a constant-$\tau$ parabola. This direction changes at every single point!

To navigate, we need to define local basis vectors. These are vectors, $\vec{e}_\sigma$ and $\vec{e}_\tau$, that are tangent to our coordinate curves at any given point. We can find them by seeing how our position vector, $\vec{r} = x\hat{i} + y\hat{j}$, changes as we tweak $\sigma$ and $\tau$ slightly. A bit of calculus gives us:

These vectors are our "local rulers." They point along the grid lines at every point, telling us "this way for $\sigma$" and "that way for $\tau$." Now, let's ask a crucial question: are these two directions perpendicular? We can check this with the dot product:

They are perfectly orthogonal, everywhere! This is a wonderfully convenient property, making calculations much simpler. But we must appreciate that this is a special feature of this particular transformation. If we had defined a slightly "skewed" system, for instance by adding a small term like $y = \frac{1}{2}(\tau^2 - \sigma^2) + \alpha \sigma \tau$, this orthogonality would be lost, and the dot product would no longer be zero. The fact that our standard parabolic coordinates are ​​orthogonal​​ is a gift.

What about the length of our local rulers?

This is fascinating! Not only are the rulers' lengths not equal to one, but they change depending on where we are on the grid. As we move away from the origin (where $\sigma=0$ and $\tau=0$), our grid lines spread out, and the local rulers get longer. These lengths, $h_\sigma$ and $h_\tau$, are known as ​​scale factors​​.

All of this geometric information—the lengths of the basis vectors and the angles between them—is neatly packaged into a single object called the ​​metric tensor​​, $g_{ij}$. Its components are simply all the possible dot products of the basis vectors: $g_{ij} = \vec{e}_i \cdot \vec{e}_j$. For our system, this gives a beautifully simple matrix:

The zeros on the off-diagonal confirm the orthogonality we discovered. The diagonal entries are the squares of our scale factors, telling us how the grid stretches from point to point. This metric tensor is the fundamental tool for doing geometry in any coordinate system.

There's another way to think about this stretching. Imagine a tiny rectangle in an abstract $(\sigma, \tau)$ plane with sides $d\sigma$ and $d\tau$. When we map this rectangle onto our real $(x, y)$ plane using our transformation rules, it becomes a small, curved quadrilateral. How much has its area changed?

This question is answered by the ​​Jacobian matrix​​, which describes how the output coordinates $(x, y)$ change in response to infinitesimal changes in the input coordinates $(\sigma, \tau)$. It's a matrix of all the partial derivatives:

The determinant of this matrix, $\det(J) = \tau^2 - (\sigma)(-\sigma) = \sigma^2 + \tau^2$, gives the local area-stretching factor. A patch of area $d\sigma d\tau$ in the "parameter space" becomes a patch of physical area $dA = (\sigma^2 + \tau^2) d\sigma d\tau$ in our plane.

Notice something wonderful? The determinant of the metric tensor, $g$, is $g = \det(g_{ij}) = (\sigma^2 + \tau^2)^2$. The area element in general relativity and differential geometry is given by $dA = \sqrt{g} \, d\sigma d\tau$. In our case, this is $\sqrt{(\sigma^2+\tau^2)^2} \, d\sigma d\tau = (\sigma^2+\tau^2) d\sigma d\tau$. It's the exact same result! The Jacobian from multivariable calculus and the metric tensor from geometry are two different languages telling the same beautiful story about the fabric of our coordinate system.

Now for a puzzle that reveals a deep truth about coordinates. Imagine a steady, uniform river flowing purely to the right. In Cartesian coordinates, this is described by a simple vector field, $\vec{V} = V_0 \hat{e}_x$. Its components are $(V_x, V_y) = (V_0, 0)$. They are constant, unchanging, everywhere.

What would an observer using our parabolic grid measure for the components of this same vector field? We might intuitively think they should still be constant. But they are not. After applying the correct transformation rules, the components in the parabolic basis, $(V^\sigma, V^\tau)$, turn out to be:

The components are no longer constant! They depend explicitly on the position $(\sigma, \tau)$. This is a profound lesson. The physical reality—the river's flow—is unchanged. But its description depends entirely on the language we use to measure it. The "constancy" of the Cartesian components was an artifact of choosing a grid that perfectly aligned with the flow. When we look at the same flow through the "lens" of our curved parabolic grid, we see its components changing from point to point because our very basis vectors (our local rulers) are changing direction and length everywhere. A vector is a geometric object, but its components are just shadows it casts onto a chosen set of axes.

Let's take this one step further. If you were an ant, condemned to walk only along the curved grid lines, how would you perceive "straight-line motion"? An object moving in a true straight line in the flat plane (like a coasting spaceship) would appear to you, the grid-bound observer, to be constantly turning. You would have to invoke some kind of "force" to explain why its path is bending relative to your grid lines.

This apparent force, which arises purely from the curvilinearity of the coordinates, is mathematically captured by objects called ​​Christoffel symbols​​. In our flat Cartesian grid, all Christoffel symbols are zero, which is the mathematical way of saying that straight lines look straight.

But what if we calculate them for our parabolic system? Even though the underlying plane is perfectly flat, the Christoffel symbols are ​​not​​ zero. For example, one of them is:

This non-zero value is the "ghost of curvature." It doesn't mean the space itself is curved. It means our coordinate system is curved. The Christoffel symbol is precisely the correction term needed to account for the fact that our basis vectors are changing from point to point. It's the mathematical equivalent of the fictitious Coriolis and centrifugal forces you experience on a spinning merry-go-round. The ground beneath the merry-go-round is flat, but your rotating frame of reference introduces apparent forces that deflect moving objects. In the same way, the Christoffel symbols tell us how to correctly write down the laws of motion, like Newton's second law, in a world viewed through the elegant but ever-changing lens of parabolic coordinates.

Now that we have acquainted ourselves with the machinery of parabolic coordinates, we might be tempted to ask, "What is all this for?" It is a fair question. Why should we bother with these graceful curves when we have the sturdy, reliable grid of Cartesian coordinates, which has served us so well? The answer, as is so often the case in physics, lies in the profound connection between symmetry and simplicity. The world is not always made of straight lines and right angles. By choosing a language that matches the geometry of a problem, we can often transform a forbiddingly complex calculation into something manageable, and in doing so, reveal the inherent beauty of the physical laws at play.

Let's start with the most fundamental aspect of physics: motion. In our familiar Cartesian world, the kinetic energy of a particle is a simple sum of squares, $T = \frac{1}{2}m(\dot{x}^2 + \dot{y}^2)$. What happens when we describe this same motion using parabolic coordinates? The transformation is a straightforward, if slightly tedious, application of the chain rule. Yet, the result is wonderfully illuminating. The kinetic energy becomes $T = \frac{1}{2}m(\sigma^2 + \tau^2)(\dot{\sigma}^2 + \dot{\tau}^2)$. Notice this! The expression is no longer a simple sum. The particle's speed depends not only on how fast the coordinates $(\sigma, \tau)$ are changing, but also on where the particle is located, through the factor $(\sigma^2 + \tau^2)$. This is the essence of curvilinear coordinates: the very fabric of our coordinate grid stretches and warps from point to point.

This complexity becomes even more apparent when we consider acceleration. A particle moving with constant $\dot{\sigma}$ and $\dot{\tau}$ is not moving in a straight line with constant velocity. The expression for acceleration in parabolic coordinates is filled with terms that look like corrections, accounting for the curving and stretching of the coordinate lines. These are not "fictitious" forces in the colloquial sense; they are real consequences of describing motion from a non-rectilinear viewpoint.

This leads us to a deeper, more geometric question. What does a "straight line"—the simplest of all paths—look like in this new system? In the language of geometry, a straight line is a geodesic, the shortest path between two points. If we ask what equations govern a geodesic path on a flat plane, but expressed in parabolic coordinates, we are led through the powerful formalism of Lagrangian mechanics to a set of equations that describe the "acceleration" needed to stay on this straight path. The resulting expressions look complicated, yet they describe the simplest possible motion. It is a beautiful paradox: the complexity is not in the physics, but in our choice of description. This is a central lesson of differential geometry and, in a grander sense, Einstein's theory of general relativity, where gravity itself is reinterpreted as the geometry of spacetime.

The utility of parabolic coordinates extends far beyond the motion of single particles. Consider the physics of continuous materials. When an engineer analyzes the stress within a steel plate, they use a mathematical object called the strain tensor to describe how the material is being stretched and sheared at every point. If this plate has a crack or a notch with a parabolic tip, describing the stress field in Cartesian coordinates becomes a nightmare. However, by transforming the strain tensor into parabolic coordinates, the boundary of the crack aligns perfectly with a coordinate line, dramatically simplifying the problem.

A more advanced tool in this domain is the Airy stress function, a scalar potential $\Phi$ whose derivatives give the stress components. For a body in equilibrium, this function must satisfy a challenging fourth-order differential equation known as the biharmonic equation, $\Delta^2 \Phi = 0$. Solving this equation is key to understanding how stress concentrates, for instance, at the tip of a crack—a critical issue in preventing material failure. For a parabolic crack, the natural choice is, of course, parabolic coordinates. The Laplacian operator $\Delta$, and by extension the biharmonic operator $\Delta^2$, takes on a specific form in these coordinates. While the conformality of the coordinates (the fact that the scale factors are equal) simplifies the Laplacian, the biharmonic operator remains complex, highlighting the subtle interplay between coordinate choice and the mathematical structure of our physical laws.

These ideas of transforming fields are not limited to mechanics. In electromagnetism and fluid dynamics, we often work with vector fields and their corresponding potentials. A covector field (or 1-form), such as $\alpha = y \, dx - x \, dy$, which is intimately related to angular momentum, can be transformed into the parabolic system. The result reveals how the components of the field are perceived from the new geometric perspective. This principle of general covariance—that the physics must be independent of the coordinates we use to describe it—is a powerful guiding light, forcing us to express our laws in the language of tensors and differential forms that behave predictably under such transformations.

Perhaps the most celebrated application of parabolic coordinates is in quantum mechanics, in solving the Schrödinger equation for a hydrogen atom placed in a uniform external electric field. This phenomenon is known as the Stark effect. The external field introduces a potential term of the form $-eEz$, which breaks the perfect spherical symmetry of the isolated atom. The problem no longer separates cleanly in spherical coordinates.

However, a miracle occurs. The problem possesses a hidden parabolic symmetry. By transforming the Schrödinger equation into three-dimensional parabolic coordinates, the equation once again becomes separable. This allows for an exact solution, something that would be nearly impossible otherwise. In this quantum world, the operators for physical observables also take on new forms. For example, one might wonder about the relationship between a component of angular momentum, say $L_x$, and one of the parabolic coordinates, $\xi = r+z$. Calculating the commutator $[L_x, \xi]$ reveals a fundamental incompatibility: one cannot simultaneously know with perfect precision both the $x$-component of angular momentum and the particle's position on a given paraboloid of revolution. This is not just a mathematical curiosity; it is a direct statement about the quantum nature of reality, revealed through the lens of parabolic coordinates.

Finally, let us not forget the workhorse of theoretical physics: integration. Whenever we need to calculate a total quantity—be it the mass of a nebula, the charge on a conductor, or the probability of finding a particle in a certain region—we must perform an integral. To do this in a new coordinate system, we must know how the element of area or volume transforms. This transformation factor is the famous Jacobian determinant. For 2D parabolic coordinates, the area element is $dA = (\sigma^2 + \tau^2) \, d\sigma \, d\tau$, and for their 3D rotational counterparts, the volume element becomes $dV = \frac{1}{4}(\xi+\eta)\,d\xi\,d\eta\,d\phi$. The Jacobian tells us precisely how much our coordinate grid is stretching or shrinking at each point.

This tool finds surprising applications in fields as seemingly distant as statistics. Imagine trying to calculate the probability that a random variable, drawn from a standard bell curve (a Gaussian distribution), falls within a region bounded by parabolas. In Cartesian coordinates, this integral is dreadful. But by a clever change to parabolic coordinates, the integrand and the domain of integration can simplify dramatically, turning an intractable problem into a solvable one.

From the trajectory of a planet to the electron cloud of an atom, from the stresses in a bridge to the mathematics of probability, parabolic coordinates demonstrate a beautiful principle: nature rewards us for speaking its language. The physicist's art is not just in discovering the laws, but in finding the right perspective from which those laws appear in their simplest and most elegant form.