Efficient Description of Parametric Amplification of Quantum Pulses

This paper presents an efficient, analytical method for determining the quantum states of output modes in parametric amplification by applying the amplification to the vacuum and then transforming the input creation operator, a technique demonstrated on various input states including coherent, Schrödinger cat, and single-photon pulses.

The Gist

Imagine you have a magical machine called a Parametric Amplifier. Think of this machine like a very sophisticated, high-speed photocopier for light pulses, but instead of just copying a picture, it stretches, squeezes, and multiplies the "quantum information" inside the light.

The problem scientists usually face with this machine is that it's incredibly messy. If you send one specific pulse of light in, the math says the output should be a tangled mess spread across infinite possible frequencies and time slots. Calculating what the final light looks like is usually like trying to solve a puzzle with billions of pieces—it's computationally impossible for complex inputs.

The Big Discovery
The authors of this paper found a shortcut. They proved that no matter what kind of "quantum pulse" you send in (even a single photon or a complex "Schrödinger cat" state), the machine only ever sends the information out through at most two specific channels (or "modes"). All the other infinite channels just carry empty, amplified background noise.

This is like realizing that even though a radio station broadcasts on a million frequencies, your specific song only ever plays on two specific stations. The rest of the dial is just static.

The New "Recipe" (The Method)
The paper introduces a new, efficient way to calculate exactly what comes out of those two channels. Instead of trying to track the complicated input pulse through the machine step-by-step, they flip the script:

1. Start with Nothing: First, they calculate what happens when you send nothing (a vacuum) into the machine. This gives them a "squeezed vacuum" state—a specific, predictable pattern of background noise.
2. Transform the Recipe: Next, they take the "recipe" used to create the original input pulse (the mathematical instructions that turned nothing into a photon or a cat state) and translate it into the language of the machine's output.
3. Mix and Match: Finally, they apply this translated recipe to the "squeezed vacuum" they calculated in step one.

The Analogy:
Imagine you want to know what a cake looks like after it's been baked in a weird, shape-shifting oven.

• The Old Way: You try to track every single grain of flour and sugar as they twist and turn through the oven's chaotic heat. It's a nightmare.
• The New Way: You first bake an empty cake pan (the vacuum) to see how the oven warps the pan itself. Then, you take the instructions for your specific cake batter (the input state), rewrite them to fit the warped pan, and apply those instructions to the warped pan. You get the final cake instantly, without tracking the individual grains.

What They Tested
To prove this works, they ran three specific examples through their new method:

• A Coherent State: Like a standard, steady laser beam.
• A Schrödinger Cat State: A weird quantum state that is like a cat being both alive and dead at the same time (a superposition).
• A Single Photon: Just one single particle of light.

They showed that for the laser and the "cat," the output stays in just one channel. But for the single photon, the information splits into two entangled channels.

The "Heralding" Trick
The paper also describes a cool trick called "heralding." Imagine you have two output channels, but one is much emptier than the other. If you put a detector on the empty channel and it says, "Hey, I detected absolutely nothing (vacuum)," you can be sure that the other channel now contains a much cleaner, higher-quality version of your quantum state.

It's like having two buckets of water. If you check the smaller bucket and find it perfectly dry, you know the water in the big bucket is now pure and undiluted. This process "purifies" the quantum state, making it more useful for future tasks.

Why It Matters
This method is fast and analytical. It doesn't require supercomputers to solve complex equations for every new input. It allows scientists to quickly predict how quantum information will behave when amplified, which is crucial for building quantum networks, quantum computers, and ultra-sensitive sensors. The authors also noted this method could be applied to other systems like optical parametric amplifiers and even different types of light waves (like those carrying orbital angular momentum).

Technical Summary

Technical Summary: Efficient Description of Parametric Amplification of Quantum Pulses

Problem Statement
Describing the evolution of bosonic quantum states traveling through open quantum networks is challenging, particularly when the radiation explores a continuum of frequency modes requiring a multi-mode description. While single-mode operations are tractable, the output of a parametric amplifier driven by a single input pulse is generally distributed in an entangled manner over modes experiencing varying degrees of squeezing. Existing methods for calculating these output states, such as Bogoliubov eigenmode decomposition, can suffer from computational complexity comparable to the boson sampling problem when the input state does not occupy only a few eigenmodes. Furthermore, traditional approaches often rely on ancillary modes to calculate the output, which increases the dimensionality of the Hilbert space unnecessarily.

Methodology
The authors present an efficient, analytical method to determine the quantum state of output modes resulting from a single quantum pulse undergoing parametric amplification. The core of the method relies on a specific decomposition of the input-output relation derived from the work of Christiansen et al. [30], which establishes that a single input pulse $u(\omega)$ contributes to at most two output modes, $v_1(\omega)$ and $v_2(\omega)$, while all other output modes contain only vacuum-derived (squeezed) contributions.

The proposed approach diverges from standard techniques by decoupling the calculation of the vacuum transformation from the input state transformation. The method proceeds in three steps:

1. Vacuum Transformation: The parametric amplifier's effect is first calculated on the vacuum state. This yields a squeezed multi-mode vacuum state $|\xi\rangle = \hat{U}|0\rangle$, characterized by a covariance matrix and displacement vector (which is zero for vacuum). This state is described using multidimensional Hermite polynomials in the Fock basis.
2. Operator Decomposition: The input annihilation operator $\hat{a}_u$ is expanded in terms of the relevant output operators $\hat{b}_{v_1}, \hat{b}_{v_1}^\dagger, \hat{b}_{v_2}, \hat{b}_{v_2}^\dagger$. The coefficients of this expansion depend on the input mode function $u(\omega)$ and the system's transfer kernels $F$ and $G$.
3. State Reconstruction: For an arbitrary input state defined as $|\psi_u\rangle = f(\hat{a}_u, \hat{a}_u^\dagger)|0\rangle$, the method applies the unitary transformation $\hat{U}$ by transforming the creation/annihilation operators within the function $f$. Specifically, $\hat{U}f(\hat{a}_u, \hat{a}_u^\dagger)\hat{U}^\dagger$ becomes $\tilde{f}(\hat{b}_{v_1}, \hat{b}_{v_1}^\dagger, \hat{b}_{v_2}, \hat{b}_{v_2}^\dagger)$. The final output state is obtained by applying this transformed operator $\tilde{f}$ to the pre-calculated squeezed vacuum state $|\xi\rangle$.

This formalism avoids the need for ancillary modes, restricting the calculation to the physical output modes occupied by the input signal, thereby ensuring computational efficiency regardless of the input state's complexity (e.g., Fock state superpositions or mixed states).

Key Contributions and Results
The paper provides a general formalism applicable to any quadratic Hamiltonian, demonstrated through an optical parametric oscillator (OPO) model. The authors validate the method by computing output states for three distinct input scenarios:

• Coherent States: The method confirms that a coherent input state occupies a single output mode, resulting in a squeezed coherent state. The displacement vector is non-zero, while the covariance matrix matches that of the transformed vacuum.
• Schrödinger Cat States: For cat states (e.g., $|\alpha\rangle \pm i|-\alpha\rangle$), the method successfully calculates the output density matrix in the single occupied mode, preserving the non-Gaussian features of the input.
• Fock States: For a single-photon input (or general Fock states), the method correctly identifies that the output populates two entangled modes ($v_1$ and $v_2$). The authors numerically compute the Wigner functions for these joint modes, showing the distribution of the non-Gaussian character across the two modes.

A significant result presented is the heralded purification of the output state. By detecting vacuum in the less populated output mode ($v_2$), the authors demonstrate that the quantum state in the primary mode ($v_1$) becomes significantly purer.

• For a single-photon input, the purity of the output in $v_1$ increases from $\approx 0.75$ (without heralding) to $\approx 0.97$ upon detecting vacuum in $v_2$ (with a success probability of $\approx 52\%$).
• For an odd Schrödinger cat state input, the purity increases from $\approx 0.75$ to $\approx 0.93$ (with a success probability of $\approx 43\%$).
  In both cases, the Wigner negativity is also enhanced, indicating a more useful non-classical resource.

Significance and Scope
The paper claims that this method offers a computationally efficient alternative to existing techniques by leveraging the specific mode decomposition of the input-output relation. By focusing only on the physically relevant output modes and separating the vacuum squeezing calculation from the input state manipulation, the method prevents the computational complexity issues associated with mixing many modes.

The authors emphasize that while the current implementation focuses on temporal/spectral modes in an OPO, the formalism is general. It can be applied to other quadratic systems such as optical parametric amplifiers and traveling wave parametric amplifiers. Furthermore, the framework is extensible to other mode bases, including transverse spatial modes in 2D/3D systems, Laguerre-Gaussian modes carrying orbital angular momentum, and relativistic spatio-temporal modes. The work also suggests potential future extensions to include non-Gaussian operations and higher-order nonlinearities via cumulant or mean-field approximations, though these are not implemented in the current study.