ASTRODREWS✦

ASTRODREWS / RESEARCH · SUMMER 2025

Simulating quantum noise
at the University of Tokyo.

Exploring Kerr squeezing in optical fibers with stochastic simulations.

During summer 2025, I worked with Zicong Xu and Professor Takuro Ideguchi at the University of Tokyo on a quantum-optics project: simulating how the quantum noise of ultrashort light pulses evolves as they travel through nonlinear optical fibers. Our goal was to develop a computational tool that could help identify promising conditions for Kerr-squeezing experiments and future quantum spectroscopy techniques.

The project became my poster, Simulating Quantum Noise and Anti-Squeezing in Kerr Media Using Stochastic Methods, presented at the 2025 MPS Research Fair in Berkeley.

Why squeeze light?

Quantum noise limits how precisely we can measure light. A squeezed state redistributes this uncertainty: fluctuations are reduced in one component of the optical field, called a quadrature, while they increase in another. That increase is known as anti-squeezing.

A conceptual comparison of a circular quantum-noise distribution and a squeezed distribution with reduced uncertainty in one quadrature.
Figure 1 · Conceptual illustration of quantum squeezing. From poster Figure 1, credited to Wang et al. (2024), reference 1 in the full poster. Select the figure to enlarge.
Photocurrent noise over time, showing a narrower distribution in the squeezed region.
Figure 2 · Illustration of a narrower photocurrent-noise distribution. From poster Figure 2; see reference 2 in the full poster. Select the figure to enlarge.

In an optical fiber, the Kerr effect makes the refractive index depend on the light’s intensity. This nonlinear interaction can change the shape and orientation of the field’s noise distribution. Understanding that evolution is important for designing experiments that aim to reduce noise in the quantity being measured.

A noise ellipse rotating and stretching in phase space during propagation through optical fiber.
Figure 3 · Kerr evolution changes the shape and orientation of the noise distribution. From poster Figure 3; see reference 3 in the full poster. Select the figure to enlarge.

Building the simulation

We started with PyNLO, a Python package for modeling nonlinear light–matter interactions. Its classical pulse-propagation solver uses the split-step Fourier method to solve the nonlinear Schrödinger equation. We extended that framework to include stochastic linear and nonlinear noise terms so we could investigate quantum-noise evolution.

Classical nonlinear Schrödinger equation and its stochastic extension, with linear and nonlinear noise terms highlighted.
Figure 4 · Extending the classical nonlinear Schrödinger equation with stochastic noise terms. From poster Figure 4; see reference 6 in the full poster. Select the figure to enlarge.

We also replaced the simplified Raman-response model with a 10-Lorentzian approximation to better represent Raman effects in the quantum-noise calculation. The simulations included Raman scattering and self-steepening.

One simulation configuration followed 200-femtosecond pulses at a wavelength of 1550 nanometers through 50 meters of fiber, averaging over 500 noise realizations. By comparing input and output distributions in phase space, we could examine how propagation changed the noise ellipse and its squeezing and anti-squeezing directions.

Centered input and output phase-space samples: the output distribution is elongated, with arrows marking the analyzed quadrature directions.
Figure 5 · Simulated input and output noise distributions after 50 meters of fiber. Both distributions are centered at the origin; the output shows strong anti-squeezing along its elongated direction. Select the figure to enlarge.

Exploring experimental conditions

We then varied pulse duration, pulse energy, and fiber length to map the simulated noise behavior. For the parameter maps shown in the poster, each 16 × 16 grid used 50 noise realizations per point; interpolation made the trends easier to visualize. These maps helped us explore how the regions of stronger anti-squeezing shifted as the fiber length changed.

Six anti-squeezing maps versus pulse energy and duration for fiber lengths from 5 to 50 meters.
Figure 6 · Simulated anti-squeezing versus pulse energy and duration for fiber lengths of 5, 10, 20, 30, 40, and 50 meters. The 16 × 16 simulation grids were interpolated for display. Select the figure to enlarge.

Comparing the model with measurements

For a separate experimental comparison, we examined anti-squeezing as a function of fiber length at four pulse energies. Both the simulations and the measurements showed the same broad trend: higher pulse energies and longer fibers generally produced greater anti-squeezing. The experimental curves also followed the simulated ordering by pulse energy.

Simulated and experimental anti-squeezing versus fiber length, with four curves representing different pulse energies.
Figure 7 · Simulation (top) and experiment (bottom) at four pulse energies. The ordering and broad length dependence agree qualitatively, while their numerical values differ. Select the figure to enlarge.

This was qualitative agreement, rather than an exact numerical match. Propagation loss and phase-noise effects that were not fully captured in the model could contribute to the remaining offsets. Further work was needed to predict squeezing levels accurately and refine the conditions for maximizing noise reduction.

What I took from the project

This work connected computational modeling with experimental quantum optics. It gave me experience extending a scientific simulation package, analyzing ensembles of noisy trajectories, and comparing model predictions with laboratory measurements. The result was a tool for exploring experimental conditions, alongside a clearer picture of the model’s strengths and limitations.

I’m grateful to Zicong Xu for his continued support and quantum-optics insights, and to Professor Takuro Ideguchi for supporting the work through his laboratory and experimental apparatus.

RESEARCH POSTER · AUGUST 2025

Explore the full project.

Open full-size poster PDF

Andrew McHaty · Zicong Xu · Takuro Ideguchi
University of California, Berkeley · University of Tokyo

Research poster presenting the stochastic simulation method, quantum-noise distributions, parameter maps, and comparison of simulated and experimental anti-squeezing. The article above summarizes the findings.

Select the poster to zoom into its figures and read the references in the original PDF.

My mentor Zicong Xu (left) and me, Andrew McHaty (right), at the University of Tokyo
With my mentor, Zicong Xu (left), at the University of Tokyo during summer 2025.