Glauber-Sudarshan P-function

The quantum world of light is notoriously abstract, defying the intuitive frameworks we use to understand our everyday classical reality. How can we create a picture or a map of a quantum state, like the light from a laser or even a single photon, using concepts we are familiar with? This fundamental challenge in quantum optics—bridging the gap between the probabilistic quantum realm and deterministic classical descriptions—led to the development of one of the field's most insightful tools: the Glauber-Sudarshan P-function. This article serves as a guide to this powerful phase-space representation.

Our journey begins in the first chapter, ​​Principles and Mechanisms​​, where we will explore the foundational idea of representing any quantum state as a mixture of "classical-like" coherent states. We will see how this elegant formalism simplifies calculations but also reveals the profound weirdness of the quantum world when the P-function behaves in ways impossible for a classical probability distribution. In the second chapter, ​​Applications and Interdisciplinary Connections​​, we will move from theory to practice, discovering how this mathematical tool is essential for modeling lasers, understanding decoherence, engineering novel quantum states, and even describing the physics of accelerating observers in empty space. By the end, the P-function will be revealed not just as a calculational method, but as a deep conceptual lens into the nature of light itself.

Imagine you are a cartographer tasked with mapping a new, strange land. This land is the quantum world of light. Our classical maps, with their familiar coordinates of position and momentum, don't quite work here. The landscape is fuzzy, uncertain, and ruled by the bizarre laws of quantum mechanics. Yet, we yearn for a picture, a map that can guide our intuition. The Glauber-Sudarshan P-function is one of the most remarkable attempts to create such a map. It's a bold effort to describe the quantum state of light using the language of classical waves, and in its successes and, more importantly, its failures, it reveals the very soul of what makes the quantum world different.

The fundamental idea, conceived independently by Roy Glauber and George Sudarshan, is as elegant as it is powerful. They asked: can we describe any quantum state of light, represented by its density operator $\hat{\rho}$, as a kind of statistical mixture of the most "classical-like" states we know? These classical-like states are the ​​coherent states​​, denoted by $|\alpha\rangle$. A coherent state $|\alpha\rangle$ is the quantum description of an ideal laser beam—a perfect, single-frequency wave with a definite amplitude $|\alpha|$ and phase. The complex number $\alpha$ is essentially the classical field amplitude; its real and imaginary parts encode the wave's quadratures, akin to position and momentum.

The P-representation proposes that any state $\hat{\rho}$ can be written as a weighted sum over all possible coherent states:

Here, the integral is over the entire complex plane, our "phase space" map. The function $P(\alpha)$ is the weighting factor. It tells us "how much" of each coherent state $|\alpha\rangle$ is present in our mixture. If $P(\alpha)$ were a simple, positive probability distribution, our job would be easy. The quantum state would just be a classical statistical ensemble of ideal laser beams.

The true beauty of this construction is how it simplifies calculations. To find the average value of some observable quantity, you often need to wrestle with non-commuting quantum operators. But with the P-function, the process magically transforms. For any ​​normally ordered​​ product of creation operators ($\hat{a}^\dagger$) and annihilation operators ($\hat{a}$), the quantum expectation value becomes a classical-looking average over phase space:

On the right side, the quantum operators $\hat{a}^\dagger$ and $\hat{a}$ have been replaced by simple complex numbers, $\alpha^*$ and $\alpha$. We have built a bridge from the quirky quantum world to the familiar territory of classical probability and integration. The question is, how sturdy is this bridge?

To get a feel for our new map, let's sketch a few portraits.

What is the P-function for a single, perfect coherent state $|\alpha_0\rangle$? Since the state is this one coherent state and nothing else, its P-function must be infinitely peaked at $\alpha_0$ and zero everywhere else. This is precisely the definition of a two-dimensional ​​Dirac delta function​​: $P(\alpha) = \delta^{(2)}(\alpha - \alpha_0)$. Our map shows a single, sharp pinprick at the location $\alpha_0$.

Now, let's consider a more "classical" kind of uncertainty. Imagine a laser field with a very stable amplitude, $\alpha_0$, but a completely random phase. This is a ​​phase-randomized coherent state​​. In our phase space map, this corresponds to a uniform distribution over a circle of radius $\alpha_0$. The P-function captures this intuition perfectly, becoming a ring-shaped delta function: $P(\alpha) = \frac{1}{2\pi\alpha_0} \delta(|\alpha|-\alpha_0)$. So far, so good. The map is behaving just as we'd expect.

The dynamics also make intuitive sense on this map. If you take a state described by $P(\alpha)$ and apply a ​​displacement operator​​ $\hat{D}(\beta)$—the quantum equivalent of giving the field a "kick" and shifting its amplitude by $\beta$—the new P-function is simply the old one shifted in phase space: $P'(\alpha) = P(\alpha - \beta)$. The entire distribution just slides over without changing its shape. It's as simple as moving a drawing on a piece of paper.

For a while, it seems our classical analogy holds perfectly. The P-function appears to be a genuine probability distribution. But this comfortable illusion shatters the moment we try to map a truly quantum state, one with no classical counterpart.

Consider the simplest case: a single-photon state, or ​​Fock state​​, $|1\rangle$. What is its P-function? A single photon is not a tiny wave; it's a discrete quantum of energy. It cannot be described as a classical mixture of coherent states. When we force it into the P-representation, the mathematics shrieks in protest. The resulting $P(\alpha)$ is not a nice, positive function. For a state that is a mixture of the vacuum $|0\rangle$ and the one-photon state $|1\rangle$, $\hat{\rho} = (1-p) |0\rangle\langle0| + p |1\rangle\langle1|$, the P-function turns out to be:

Look at that second term! It involves derivatives of a delta function. This is not a function in the ordinary sense; it is a highly singular "generalized function" or distribution. It can take on values that are effectively more negative and more sharply peaked than a simple delta function. For a pure Fock state $|n\rangle$, the situation gets even more extreme, involving $2n$-th order derivatives of the delta function.

This is the crucial lesson: ​​The P-function is a quasi-probability distribution​​. The prefix "quasi" is doing a lot of work. It means that $P(\alpha)$ is not bound by the rules of classical probability. It can be negative, or it can be more singular than a Dirac delta function. These wild behaviors are not a flaw in the formalism; they are its most important feature. They are the unambiguous signature of ​​non-classicality​​. A state is non-classical if and only if its P-function cannot be interpreted as a true probability distribution. Our map, by becoming strange and contorted, is telling us that we have left the classical world behind.

These highly singular P-functions, while mathematically precise, are difficult to visualize. It's like having a map with infinitely sharp peaks and valleys. Can we "tame" this wild landscape? Yes, by looking at it through a pair of blurry glasses. This process, known as smoothing or convolution, has both a mathematical and a deep physical meaning.

Mathematically, we can convolve the P-function with a smooth kernel, like a Gaussian, to wash out the sharp features. The ​​Wigner function​​, $W(\alpha)$, another famous phase-space map, is precisely such a smoothed version of the P-function. The two are related by a convolution with a specific Gaussian kernel. The Wigner function can still be negative (a sign of non-classicality), but it is always a well-behaved, regular function.

If we blur the P-function even more, we arrive at the ​​Husimi Q-function​​, $Q(\alpha)$. The Q-function is obtained by convolving the P-function with an even wider Gaussian. The magic of the Q-function is that it is always non-negative. It's a true probability distribution. Let's see what this does to the monstrous P-function of a Fock state $|n\rangle$. After smoothing it to get the Q-function, the infinitely singular object transforms into a beautiful, smooth ring in phase space:

For $n=1$, this is a single ring. For larger $n$, the ring expands. We've tamed the beast and created an intuitive (though blurry) picture of a state with $n$ photons.

This smoothing isn't just a mathematical trick. It corresponds to a real physical process: losing information or adding noise. Imagine a quantum system, like a harmonic oscillator, initially in a pure coherent state (a delta-function P-function). If this system interacts with a warm environment, it starts to decohere. Its P-function, which began as an infinitely sharp spike, will spread out over time, evolving into a broad Gaussian whose width depends on the temperature of the environment. The quantum purity is lost, and the state becomes more classical-like.

In fact, one can ask: how much noise do we need to add to any quantum state to guarantee that its P-function becomes a well-behaved, non-negative distribution? The answer is remarkably simple and profound. If you mix any state, no matter how exotic and non-classical, with a thermal field having an average of just one photon ($\bar{n}_{th} = 1$), the resulting P-function will always be non-negative. This one-photon threshold is a fundamental measure of the noise required to erase all signatures of non-classicality from the P-representation. It's the amount of blurring required to turn our strange quantum map into a familiar classical one, effectively turning the P-function into the always-positive Q-function.

The Glauber-Sudarshan P-function, therefore, is more than a computational tool. It's a deep probe into the quantum nature of light. Where it behaves classically, it offers powerful intuition. Where it misbehaves, it provides a precise and unambiguous flag for the quantum phenomena that defy classical description, revealing the beautiful and strange boundary between the two worlds.

Now that we have wrestled with the peculiar mathematics of the Glauber-Sudarshan P-function in the previous chapter, you might be wondering, quite reasonably: what is it good for? Is it merely a theorist's mathematical plaything, a clever but ultimately esoteric way of rewriting what we already know? The answer, you will be happy to hear, is a resounding no. The P-function is not just a description; it is a powerful and practical tool, a kind of "phase-space lens" through which we can understand, predict, and even engineer the behavior of light and matter in the quantum world. Its utility stretches from the design of everyday technology to the exploration of the most profound questions about the nature of reality.

Perhaps the most immediate and tangible applications of the P-function are found in the field it helped to create: quantum optics. Many of the devices that define our modern technological world, from lasers to optical amplifiers, operate on principles that can only be fully understood through a quantum lens. The P-function provides the perfect language for this.

Consider the laser, the quintessential quantum optical device. How does a jumble of excited atoms in a cavity conspire to produce a beam of perfectly orderly, coherent light? The P-function was the key that unlocked this mystery. In the celebrated Scully-Lamb theory of the laser, the evolution of the light field is described by a Fokker-Planck equation for its P-function. This equation elegantly captures the competition between two opposing forces: the gain from the energized atoms, which tries to build up the field, and the inevitable leakage of photons from the cavity mirrors, which tries to dissipate it. The equation also includes a diffusion term, representing the random "kicks" given to the field by spontaneous emission. The steady-state solution of this equation tells a remarkable story. Below a certain pumping threshold, the P-function is peaked at the origin—there is no light. But above the threshold, the distribution dramatically transforms into a narrow ring in phase space, the signature of a stable, coherent laser beam with a well-defined amplitude but a random phase. This framework is so powerful that it can be extended to include more complex, real-world effects, such as the saturation of the gain medium at high intensities.

The story is just as compelling for optical amplifiers, which are the backbone of our global fiber-optic communications network. If you have a faint optical signal, you want to boost it. But quantum mechanics exacts a price: any process that amplifies a signal must, by necessity, add noise. The P-function provides a beautifully intuitive picture of this fundamental trade-off. The transformation of a state by a quantum-limited amplifier is described as a convolution in phase space. The P-function of the input state is effectively "smeared out" by a Gaussian function representing the injection of vacuum fluctuations, or quantum noise. If you feed a highly non-classical state, like a single photon, into an amplifier, its singular P-function gets smoothed into a broader, more regular distribution by this added noise. While the output may still exhibit non-classical features (like negative values), the amplifier's noisy nature has irrevocably blurred some of its quantum character.

The universe is a noisy place. Quantum states, with their delicate superpositions and correlations, are notoriously fragile. Their interaction with the surrounding environment—a process called decoherence—tends to wash away their "quantumness," making them behave more like classical objects. The P-function formalism provides an exceptionally clear moving picture of this process.

Let's take the simplest model of an open quantum system: a single harmonic oscillator (like a mode of light in a cavity) coupled to a thermal reservoir. The P-function's evolution is governed by a Fokker-Planck equation that contains two essential parts. The first is a "drift" term, which acts like a drag force, pulling the average amplitude of the state towards the origin (the vacuum). This represents energy dissipation. The second is a "diffusion" term, which causes the P-function to spread out, representing the random thermal fluctuations imparted by the reservoir. The strength of this diffusion, $D_P = \gamma N$, is directly proportional to both the damping rate $\gamma$ and the average number of thermal photons $N$ in the environment, elegantly connecting the quantum dynamics to the temperature of the outside world.

To make this more vivid, imagine watching a single photon—a quintessential non-classical state—as it decays. At time zero, its P-function is not a simple function at all, but a highly singular object involving derivatives of a Dirac delta function. This bizarre mathematical form is the P-function's way of shouting "This is not a classical state!" Now, we let it evolve under amplitude damping (leaking out of a cavity at zero temperature). The evolution equation shows how the singular, non-classical part of the P-function shrinks over time, its coefficient decaying exponentially as $e^{-\gamma t}$. Meanwhile, a classical part, a delta function at the origin, grows to take its place. We are literally watching the quantum state decohere, its non-classical essence draining away until all that is left is the vacuum.

As we've seen, the P-function's character is a direct indicator of the quantum nature of a state. If $P(\alpha)$ is positive and well-behaved everywhere, it could, in principle, describe a classical statistical mixture of light waves. But if it ever dips into negative values, or if it becomes more singular than a Dirac delta function, there is no classical explanation. The state is certifiably quantum. This makes the P-function a prime tool for identifying and characterizing non-classicality.

This property is not just for diagnostics; it's a guide for creation. Quantum engineers use this formalism as a kind of recipe book. The rules that translate operator actions (like adding or subtracting a photon) into differential operations on the P-function allow one to predict the non-classical character of engineered states. For example, by taking a simple thermal state (with a Gaussian P-function) and applying the "recipe" for adding a single photon, one can derive the P-function for the resulting state—a single-photon-added thermal state (SPATS). The resulting P-function is a "Mexican hat" shape, dipping into negative values around the origin, providing a clear signature of its engineered non-classicality.

The power of the P-function extends beyond single modes of light into the strange and wonderful realm of entanglement. For a system with multiple modes, the P-function becomes a distribution over a higher-dimensional phase space, $P(\alpha_1, \alpha_2, \dots)$. Its structure now encodes not only the properties of the individual modes but, crucially, the quantum correlations between them. Consider the tripartite W-state, a fundamental entangled state of three particles (or modes) in quantum information science. Its P-function is a fascinating object. It contains terms that are products of derivatives with respect to the different modes' phase-space coordinates, all acting on a delta function locking the fields to the vacuum. These cross-derivative terms are the mathematical embodiment of entanglement; they tell us that what happens in the phase space of one mode is inextricably linked to what happens in the others.

At this point, you might be feeling that this is all well and good for the theorists, but how does an experimentalist ever get their hands on this P-function? We can't measure it directly. The beautiful answer lies in a technique called quantum tomography. It is conceptually similar to a medical CAT scan. In a CAT scan, a series of 2D X-ray images (projections) are taken from different angles, and a computer algorithm reconstructs a full 3D image of the inside of the body. In quantum tomography, the experimentalist measures the probability distributions of the field's quadratures at all different phases $\theta$. These are the "projections." A mathematical procedure known as a filtered back-projection, a cousin of the Radon transform, can then be used to reconstruct a phase-space function from this complete set of measured data. Because the raw P-function can be highly singular and ill-behaved, this reconstruction often yields a "regularized" or smoothed version, but it nonetheless provides a direct experimental window into the phase-space landscape of the quantum state.

Finally, we arrive at the most astonishing and profound application of the P-function—one that bridges quantum optics with general relativity and quantum field theory. It is a cornerstone of modern physics that the vacuum is not truly empty but is seething with quantum fluctuations. The Unruh effect is the prediction that an observer undergoing constant acceleration through this vacuum will not perceive it as empty space. Instead, they will find themselves immersed in a bath of thermal radiation, as if the vacuum itself were glowing with a specific temperature proportional to their acceleration.

What state does the accelerating observer see? When we trace out the parts of the universe they cannot see, the state they are left with for a given mode of the field is precisely a thermal state. And we know exactly what the P-function for a thermal state is: a simple, elegant Gaussian distribution centered at the origin. By calculating the properties of the Unruh effect, one can derive the mean photon number of this thermal bath and, from that, write down its P-function explicitly in terms of the observer's acceleration $a$ and the mode's frequency $\omega$. Think about that for a moment. The P-function, a tool we developed to understand the light from a laser, turns out to also describe a fundamental property of spacetime as perceived by an accelerating observer. It reveals a deep, hidden unity in the fabric of physics, connecting the laboratory bench to the cosmos itself. This, in the end, is the true beauty of a powerful physical idea: its ability not only to solve the problems it was designed for but also to illuminate unexpected corners of the universe and reveal the interconnectedness of it all.