A canonical approach to quantum fluctuations

This paper introduces a canonical formalism for analytically computing quantum fluctuations in the positions, velocities, norms, and phases of soliton and breather solutions within integrable partial differential equations, specifically applying it to nonlinear Schrödinger equation models under coupling-constant quenches with both white-noise and correlated-noise vacuum states.

The Gist

The Big Picture: The "Perfect" Wave vs. The "Fuzzy" Reality

Imagine you are watching a perfect, solitary wave travel across a calm ocean. In physics, this is called a soliton. It's a special kind of wave that holds its shape together as it moves, unlike normal waves that spread out and fade away.

In the world of classical physics (the physics of big, visible things), we can describe this wave with a perfect, smooth equation. If you know the wave's speed, position, and size right now, you can predict exactly where it will be tomorrow. It's like a perfectly choreographed dance.

But, we live in a quantum world (the world of atoms and tiny particles). In this world, nothing is ever perfectly smooth or predictable. Everything has a tiny bit of "fuzziness" or "jitter" called quantum fluctuations. Even if the wave looks perfect to the naked eye, at the atomic level, it's actually shaking and vibrating slightly.

The Problem:
Scientists wanted to know: If we create a complex wave made of two or three of these solitons stuck together (called a "breather"), how much does this quantum jitter mess up their positions, speeds, and sizes?

Previously, calculating this was like trying to solve a massive, 10,000-piece puzzle by hand. It took supercomputers days to crunch the numbers, and sometimes the results were just guesses based on messy calculations.

The Solution:
The authors of this paper invented a new "shortcut" or a new set of tools. They found a way to translate the messy, fuzzy quantum world into a clean, organized mathematical language (called canonical formalism). This allowed them to solve the puzzle with a simple pen and paper (or a laptop) in a few hours, getting exact answers instead of messy approximations.

---

The Analogy: The "Magic Map"

To understand how they did it, imagine you have a Magic Map of a city.

1. The Old Way (The Hard Way):
   Imagine you want to know how a tiny earthquake (quantum fluctuation) affects the traffic in a specific neighborhood. The old method was to look at every single car, every single driver, and every single road condition in the whole city. You'd have to calculate how a bump in the road in the north affects a driver in the south. It's overwhelming, slow, and prone to errors.

2. The New Way (The Canonical Approach):
   The authors realized that the city (the soliton) has a few "control knobs" that define its entire state:

   • Where is it? (Position)
   • How fast is it going? (Velocity)
   • How big is it? (Size/Number of particles)
   • What is its "phase" (like the timing of a heartbeat)?

   Instead of tracking every single atom (every car), they realized they could just track how the control knobs move when the earthquake hits. They built a "Magic Map" (a mathematical transformation) that connects the messy, fuzzy atomic world directly to these few control knobs.

   Because this map is built on the rules of symmetry (mathematical balance), they could flip the map around. Instead of struggling to figure out how the atoms move the knobs, they could easily figure out how the knobs move the atoms. This turned a super-hard math problem into a simple one.

---

The Experiment: The "Snap" of the Coupling Constant

The paper focuses on a specific experiment involving Bose-Einstein Condensates (BECs). Think of a BEC as a super-cold cloud of atoms that all act like a single giant atom.

1. The Setup: Scientists create a single, perfect soliton (a "Mother Soliton") in this cloud.
2. The "Quench" (The Snap): Suddenly, they change a dial on the machine (the "coupling constant"). Imagine you have a rubber band holding the atoms together. Suddenly, you snap the tension of the rubber band to be four times stronger (or nine times stronger).
3. The Result: Because of this sudden change, the single "Mother Soliton" splits apart into two or three "Daughter Solitons" that are stuck together in a breathing motion (expanding and contracting).

The Question:
According to the "perfect" classical theory, these daughter solitons should be born perfectly still relative to each other. They should have zero speed and zero distance between them.

The Reality:
Because of quantum fluctuations, they aren't perfectly still. They have a tiny, random jitter. They might drift apart slightly or move at a tiny speed. The paper calculates exactly how much they jitter.

---

The Two Types of "Noise"

The authors tested their new method using two different ideas about what the "fuzziness" looks like:

1. White Noise (The Static): Imagine the static on an old TV. It's random, and every point is independent of its neighbor. This is a simple, idealized model.
2. Correlated (Colored) Noise (The Ripple): Imagine a ripple in a pond. If you push one spot, the water nearby moves too. The "fuzziness" here is connected; it has a pattern. This is a more realistic model for how atoms actually behave.

The Surprise:
The authors found that for most of their calculations, it didn't matter much which model they used. The "simple" white noise gave almost the same answer as the "complex" correlated noise. This is a huge relief because it means scientists can use the simpler model for future experiments without losing accuracy.

---

Why This Matters

1. Speed and Simplicity: Before this paper, calculating these quantum effects for a 3-soliton system would have taken a supercomputer days (or might have been impossible). Now, a researcher can do it on a laptop in a few hours.
2. New Discoveries: Because the math is now so easy, they could solve a problem for a 3-soliton system that was previously too hard to tackle. They found the exact "jitter" values for these complex systems.
3. Experimental Verification: This gives experimentalists a precise target. If they build a machine to split solitons, they can now say, "We expect the solitons to drift apart by exactly this much due to quantum mechanics." If they see that amount, they have proven that quantum mechanics is at work even in these large, macroscopic objects.

Summary in a Sentence

The authors built a new mathematical "translator" that turns the impossible-to-calculate jitter of trillions of atoms into a simple calculation of a few control knobs, allowing them to predict exactly how quantum fuzziness makes giant atomic waves wiggle and drift.

Technical Summary

1. Problem Statement

The paper addresses the computation of quantum fluctuations in the macroscopic parameters (position, velocity, norm, and phase) of soliton breathers within integrable classical field theories, specifically the Nonlinear Schrödinger Equation (NLSE).

• Context: The NLSE serves as a mean-field (classical) approximation for underlying many-body quantum systems, such as Bose-Einstein Condensates (BECs) in highly anisotropic traps. While mean-field theory predicts that soliton breathers (nonlinear superpositions of solitons) created via a sudden quench of the coupling constant should have zero relative velocity and displacement between constituent solitons, the underlying quantum many-body theory predicts non-zero variances (fluctuations) in these parameters.
• The Challenge: Previous methods (e.g., in Phys. Rev. Lett. 125, 050405 (2020)) relied on "quasi-inner products" and intermediate numerical calculations. These approaches were computationally expensive, often requiring numerical integration for days even for 2-soliton cases, and became intractable for 3-soliton systems. Furthermore, they lacked a rigorous canonical framework, relying on unproven assumptions regarding the orthogonality of continuum fluctuations.
• Goal: To develop a streamlined, analytic method to compute these fluctuations for $n$-soliton breathers (specifically 2- and 3-soliton cases) under different noise models (white and correlated noise).

2. Methodology: Canonical Formalism

The authors propose a canonical formalism that exploits the integrability of the underlying classical field theory to transform the problem from a field-theoretic calculation into a finite-dimensional canonical transformation problem.

Core Principles

1. Canonical Transformation: The method relies on the existence of a canonical transformation between:
   • "Old" variables: The continuous field degrees of freedom (real and imaginary parts of the field $\Psi(x,t)$).
   • "New" variables: A finite set of discrete canonical coordinates ($Q_i, P_i$) characterizing the solitons (e.g., position, momentum, norm, phase).
2. Inverse Derivative Strategy:
   • Directly computing derivatives of soliton parameters with respect to the field ($\partial Q / \partial \Psi$) is difficult because it requires solving the inverse scattering problem.
   • However, expressing the field $\Psi$ in terms of soliton parameters ($\Psi(Q, P)$) is analytically tractable for integrable systems (via the inverse scattering method).
   • The authors utilize Lagrange brackets and the direct conditions of canonical transformations to invert the derivatives. Specifically, they use the relation:
     $$\frac{\partial Q_i}{\partial q_j} = \frac{\partial p_j}{\partial P_i}$$
     This allows them to compute the required derivatives of the "new" variables with respect to the "old" field variables by differentiating the known field solution with respect to the soliton parameters.
3. Fluctuation Calculation:
   • The fluctuations in the new variables ($\delta Q, \delta P$) are expressed as linear functionals of the field fluctuations ($\delta \Psi$).
   • The variance is computed by integrating the product of these derivatives against the noise correlators (vacuum expectation values):
     $$\langle \delta Q_i^2 \rangle = \int \int \frac{\partial \Psi}{\partial Q_i} \frac{\partial \Psi^*}{\partial Q_i} \langle \delta \Psi(x) \delta \Psi^*(y) \rangle dx dy$$
   • Due to the specific structure of the NLSE solutions (exponential/algebraic combinations), these integrals can be evaluated analytically.

Noise Models

The paper evaluates fluctuations under two distinct vacuum states:

1. White-Noise Vacuum: Assumes uncorrelated fluctuations ($\langle \delta \Psi(x) \delta \Psi^\dagger(y) \rangle = \delta(x-y)$). This is a standard approximation in quantum optics.
2. Correlated-Noise (Colored) Vacuum: Uses particle-number-conserving (U(1)-symmetry-conserving) Bogoliubov modes. This model corrects the standard Bogoliubov approach by projecting out the zero-mode component to preserve the total particle number, providing a more physically accurate description for BECs.

3. Key Contributions

• Analytic Solutions: The primary contribution is the derivation of closed-form analytic solutions for quantum fluctuations in 2- and 3-soliton breathers. This replaces the previous reliance on heavy numerical integration.
• Computational Efficiency: The method reduces the computational complexity significantly. Calculations that previously took days (or were impossible for 3-soliton cases) are now solvable on a laptop in hours.
• Conceptual Rigor: The formalism provides a rigorous canonical structure, eliminating the unproven assumptions about orthogonality found in previous works. All degrees of freedom (including continuum radiation) are treated within a unified canonical framework.
• Validation of Corrections: The authors investigate the impact of particle-number-conserving corrections (correlated noise). Surprisingly, they find that for most macroscopic parameters (relative position, velocity, norm), these corrections do not alter the final results compared to the simpler white-noise model, though they do affect the phase fluctuations.

4. Results

The authors applied their formalism to Odd-Norm-Ratio (ONR) breathers created by quenching the coupling constant of an attractive 1D BEC.

• 2-Soliton Breather:

  • Created by a 4-fold quench of the coupling constant ($g \to g/4$), resulting in two daughter solitons with norm ratios 1:3.
  • The authors computed the initial variances and covariances for relative position ($b$), velocity ($v$), norm difference ($n$), and phase difference ($\theta$).
  • White Noise: Results matched previous numerical findings in [38] exactly.
  • Correlated Noise: Results matched [38] to the reported significant figures.
  • Scaling: Fluctuations involving canonical coordinates (position/phase) scale as $1/N$, while those involving momenta (velocity/norm) scale as $N$.

• 3-Soliton Breather:

  • Created by a 9-fold quench ($g \to g/9$), resulting in three solitons with norm ratios 1:3:5.
  • This case was previously intractable analytically. The authors successfully derived analytic expressions for all variances and covariances.
  • The results confirm that the canonical formalism scales effectively to higher $n$.

• Physical Insight:

  • The non-zero variances in relative velocity and displacement imply that even if a breather is created with perfect symmetry in the mean-field limit, quantum fluctuations will cause the constituent solitons to separate immediately after the quench.
  • The particle-number-conserving corrections (correlated noise) were found to have a negligible effect on the relative position and velocity variances, validating the use of simpler models for these specific observables in many experimental contexts.

5. Significance

• Experimental Relevance: The results provide precise theoretical predictions for experiments with ultracold atomic gases (e.g., Lithium BECs) where soliton breathers are created via Feshbach resonance quenches. The predicted fluctuations are within the detection limits of modern experiments, offering a way to observe quantum many-body effects in macroscopic solitons.
• Methodological Advancement: The paper establishes a general "canonical approach" applicable to any integrable PDE with a known Hamiltonian structure and soliton solutions. This opens the door to studying quantum fluctuations in other integrable systems (e.g., KdV, Sine-Gordon) without resorting to brute-force numerical methods.
• Bridging Scales: It successfully bridges the gap between the microscopic quantum many-body description (Lieb-Liniger model) and the macroscopic classical field description (NLSE), quantifying the "beyond-mean-field" corrections in a rigorous, analytic manner.