Fieldoscopy at the quantum limit

TL;DR

This paper introduces a new method to measure electric fields of single photons with extreme sensitivity and precision, enabling new quantum technologies.

Contribution

A novel technique for measuring petahertz-scale electric field transients at the single-photon level with yoctojoule sensitivity.

Findings

Observed a classical regime breakdown consistent with Monte Carlo simulations.

Achieved yoctojoule-level sensitivity and over 90 dB dynamic range.

Measured intrapulse light coherence inaccessible to traditional spectroscopy.

Abstract

We demonstrate a novel concept for measuring time-varying electric field transients of petahertz-scale photons down to a single-photon regime. We observe a clear breakdown of the classical regime consistent with our Monte Carlo model. We reach unprecedented yoctojoule-level (10⁻²⁴ J) sensitivity and a dynamic range exceeding 90 decibels. We utilize this capability to measure intrapulse light coherence - a regime inaccessible to conventional, time-averaged spectroscopy. This opens new avenues for quantum information, cryptography, and quantum light-matter interactions on sub-cycle time scales with attosecond precision.

Figures

Experimental concept.a Schematic of the optical setup. $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\rm{E}}}_{{\rm{Sampilng}}}$$\end{document}$ and $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\rm{E}}}_{{\rm{Test}}}$$\end{do

Measurements of strong coherent states.a Measured electric fields of coherent states consisting of a large number of photons. b Central part of the pulse ( $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left\langle n\right\rangle$$\end{document}$ = 6.42 x $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \

Measurements of weak coherent states.a Measured electric fields of several weak coherent states consisting of a small number of photons. b Measured and simulated breakdown of the mean signal ( $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\bar{{\rm{S}}}$$\end{document}$ ) from the linear scaling. c Measured and simulated breakdown of the linear scaling of the measured standard deviation (σ). Poisson (d) and Bose

Measurements of the intrapulse coherence.a Measured electric field of a weak coherent state ( $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left\langle n\right\rangle$$\end{document}$ = 4.38) with corresponding Gaussian windows used for the evaluation of intrapulse photon statistics in (b) Measured breakdown of the linear scaling of the mean signal ( $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepacka

Full text

Introduction

Light coherence describes the relationship between the fluctuations of photons in a light source and their degree of distribution in time and space^1^. This property is fundamental for distinguishing laser light^2^ from incoherent sources and has enabled widespread application across disciplines, from fundamental physics^3,4^, precision measurements^5,6^, industrial manufacturing^7^, data storage^8^ or medical diagnostics^9–11^. Coherence enables probing of ultrafast processes in matter with femtosecond lasers, offering insights into chemical reactions^12,13^ and material dynamics^14^. By measuring the electric field of light, coherent processes can be studied on sub-cycle time scales with attosecond precision and unprecedented sensitivity^11,15^. This concept, known as fieldoscopy^16^, has recently been utilized to measure the formation of the optical response^15^, ultrafast magnetism^17^, and the light-matter energy transfer^18^ in photoexcited solids.

Coherence is directly linked to the quantum properties of light, i.e. photon statistics, making it a crucial quantity in quantum optics and quantum information science. For example, weak coherent states of a laser with a low mean photon number are widely utilized in quantum cryptography to create single-photon pulses. High-frequency photons are particularly advantageous for generating single-photon pulses, as the substantial amount of energy can be concentrated in a single quantum. Light coherence and photon statistics however can be affected by nonlinear effects^19^. In a laser, self-phase modulation, gain saturation and recovery^20^, amplified spontaneous emission^21^, intracavity soliton dynamics^22^ or mode-lock instabilities^23^ are some of the examples. In intense pulsed lasers (e.g. mode-locked lasers) the electric field strength varies across the laser pulse which inevitably leads to the temporal dependence of coherence, i.e. intrapulse coherence. While pulse-averaged coherence can be measured^24^, the coherence variation within a pulse (intrapulse coherence) requires a direct access to the electric field in the regime dominated by the quantum nature of light.

Fieldoscopy enables direct access to the electric field of light and has seen rapid developments in recent years^25,26^. Measuring electric field of petahertz-scale photons however typically requires intense, few-cycle carrier-envelope phase (CEP) stable pulses, complex experimental infrastructure, and highly non-linear light matter interaction^27–29^. This complexity limits the sensitivity required to access the single-photon regime, where the quantum nature of light is dominant.

Here we perform the first electric field-resolved measurements in the single-photon regime. We measure 1.2 eV (193 zJ, 1030 nm) photons with the average pulse energy down to 876 yoctojoules – a regime where 99.55% of pulses are in a vacuum state and out of the remaining 0.45% pulses, 99.78% are in a single photon state. Our approach does not rely on complex laser infrastructure, highly nonlinear light–matter interactions, or carrier-envelope phase (CEP)-stable pulses. We observe a distinct breakdown of the classical regime, which is well captured by our Monte Carlo model (more in Methods). Furthermore, we employ this concept to probe the time-dependent intrapulse coherence of a mode-locked laser, revealing how coherence varies across the pulse. We find that the more intense regions of the pulse exhibit reduced coherence compared to the weaker ones. We attribute this to the presence of an intensity-dependent decoherence process which admixes a thermal component to the ideal coherent state.

Our work extends petahertz fieldoscopy into the quantum regime opening new possibilities for sub-cycle quantum information, cryptography, and precision metrology. The ability to directly access the electric field of weak states of light represents a paradigm shift in quantum optics, providing a new tool for investigating coherence dynamics at sub-cycle timescales and paving the way towards quantum lightwave electronics.

Results

The concept, schematically shown in Fig. 1a, is based on the recently demonstrated generalized heterodyne optical sampling technique (GHOST)^30^. This technique is similar to the well-established electro-optic sampling^31,32^ which is extensively utilized^33,34^ to sample the electric field of light in the terahertz domain. In these schemes, two pulses, the sampling pulse and the test pulse (whose electric field is of interest) with a temporal delay τ,

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${E}_{S}(t)={\bar{E}}_{S}{e}^{-i({\omega }_{S}t+{\varphi }_{S})}$$\end{document} \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${E}_{T}(t,\tau )={\bar{E}}_{T}{e}^{-i({\omega }_{T}t+{\varphi }_{T}+{\omega }_{T}\tau )}$$\end{document}

interact in a nonlinear crystal to produce a nonlinearly generated pulse (NP) that encodes the test-field’s instantaneous value:

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${E}_{{NP}}\left(t,\tau \right)\approx {\bar{E}}_{S}^{n-1}{\bar{E}}_{T}{e}^{-i((n{{{-}}}1)(\omega_{S}t+{\varphi }_{S})\pm {\omega }_{T}t\pm \varphi _{T}+{\omega }_{T}\tau )}$$\end{document}

Fig. 1. Experimental concept.a Schematic of the optical setup. \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\rm{E}}}_{{\rm{Sampilng}}}$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\rm{E}}}_{{\rm{Test}}}$$\end{document} denote the time-varying electric fields of the sampling and test pulses, respectively. \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\rm{E}}}_{{\rm{LO}}}$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\rm{E}}}_{{\rm{NP}}}$$\end{document} represent the electric fields of the nonlinearly generated local oscillator and the nonlinearly mixed pulse. NC, WG, BP and τ refer to the nonlinear crystal, wire-grid polarizer, bandpass filter and the time delay between sampling and test pulses respectively. b, Benchmarking of the measured GHOST signals with CEP-stabilized and CEP-unstabilized pulses. c Fourier-transform spectra of the signals in (b) illustrating the frequency content with and without CEP stabilization. d Measured GHOST traces of the test pulse transmitted through glass with three slightly different thicknesses. The measurement highlights the signal’s sensitivity to the carrier-envelope phase (CEP) of the test pulse. The medium’s thickness was adjusted using a pair of fused silica wedges. The shaded areas in (b−d) correspond to one standard deviation of the measurement

Here, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\bar{E}}_{S}$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\bar{E}}_{T}$$\end{document} are field strengths, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\omega }_{S}$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\omega }_{T}$$\end{document} are frequencies, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varphi }_{S}$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\varphi }_{T}$$\end{document} are global phases of the sampling and test fields. The NP can be any nonlinear mixing signal^30^, such as sum-frequency (SFG), difference-frequency (DFG) generation, four-wave mixing (FWM), etc. By interfering the nonlinearly generated field with a “stationary” local oscillator (LO):

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${E}_{{LO}}\left(t\right)\approx {\bar{E}}_{S}^{m}{e}^{-{im}\left({\omega }_{S}t+{\varphi }_{S}\right)}$$\end{document}

one can retrieve the electric field of the test pulse at each delay. The \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m$$\end{document} (∈ 1, 2, ...) and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n$$\end{document} (∈ 2, 3, ...) factors correspond to the light-matter interaction orders used for LO and NP, respectively. The LO can be simply the sampling field itself (m = 1) or a nonlinearly generated light, e.g. second (m = 2), third (m = 3), or higher harmonics^30^. The intensity of the interference I(τ, ω_D_) at the detection frequency ω_D_ is measured with a photodiode and is proportional to the electric field of the test pulse at time τ. Experimentally, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{\rm{\omega }}}_{{\rm{D}}}={{\rm{m}}{\rm{\omega }}}_{{\rm{S}}}=\left(n-1\right){{\rm{\omega }}}_{{\rm{S}}}\pm {{\rm{\omega }}}_{{\rm{T}}}$$\end{document} , and is set by the bandpass filter (BP),

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${I}_{{LO}+{NP}}\left(\tau ,{\omega }_{D}\right)\approx {\bar{E}}_{S}^{m+n-1}{\bar{E}}_{T}\cos \left({\omega }_{T}\tau -(m-{n+1)\varphi }_{S}\pm {\varphi }_{T}\right)$$\end{document}

As we show further, unlike electro-optic sampling where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m=1$$\end{document} , \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n=2$$\end{document} , and previous implementations of GHOST (m = 3, n = 2^30^, or m = 2, n = 3^35^) balanced \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m=n$$\end{document} combination provides a previously unexplored opportunity. In the unbalanced sum-frequency-like case, at each time delay \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\tau$$\end{document} , \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${I}_{{LO}+{NP}}\approx {-\left(m-n+1\right)}({\varphi }_{S}+\Delta {\varphi }_{S})+{\varphi }_{T}+\Delta {\varphi }_{T}$$\end{document} (Eq. 5). If the detection is not performed in a single-shot, unevenly amplified global phase instabilities \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\triangle {\varphi }_{S}$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\triangle {\varphi }_{T}$$\end{document} , i.e. CEP instability, would average \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${I}_{{LO}+{NP}}$$\end{document} to zero. On the other hand, in the balanced case, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${I}_{{LO}+{NP}}\approx {\varphi }_{T}-{\varphi }_{S}+\triangle {\varphi }_{T}-\triangle {\varphi }_{S}$$\end{document} . When CEPs of the sampling and test pulses are locked, i.e. when pulses originated from the same laser, the CEP fluctuations are locked ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\triangle {\varphi }_{S}=\triangle {\varphi }_{T}$$\end{document} ) and cancel each other out enabling extremely sensitive detection without the need for CEP-stable pulses.

In this scheme, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varphi_{S}=0$$\end{document} , fieldoscopic measurement with CEP-stabilized and unstabilized pulses is expected to yield identical results. To confirm this, we set up a GHOST detection based on the balanced SHG/SFG channel ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m=n=2$$\end{document} ) utilizing a state-of-the-art octave-spanning (450–1000 nm), CEP-stable, sub-3 fs pulse^36^, which is split into two orthogonally polarized replicas—sampling and test—with a controlled time delay τ between them. After setting the field strength and the CEP of each pulse, the pulses are recombined and collinearly focused on the nonlinear medium. The test pulse was kept weak while the field strength of the sampling pulse was ~1 V/Å. We use a ∼ 12 μm thick z-cut α-quartz crystal as the nonlinear medium. The spectral response of the SHG/SFG channel covers the entire spectrum of the test pulse^30^, while the CEP stability of the pulses can be switched off on demand. We set the CEP of the sampling pulse to zero and perform sampling of the test electric field with CEP-stabilized and CEP-unstabilized pulses. Figure 1b, c confirms that the balanced ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m=n$$\end{document} ) detection yields identical field traces whether or not the pulses are CEP-stabilized, as predicted.

To confirm that the measurement remains sensitive to the CEP of the test pulse, we vary its CEP by inserting a small amount of material using fused silica wedges (more in the Supplementary Information) and repeat the measurement. As can be seen in Fig. 1d, the change in the CEP of the test pulse is clearly captured by the measurement. We also note the enormous sensitivity of the measurement, which is typical for the heterodyne detection.

Fieldoscopy of strong (classical) coherent states

After confirming that the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m=n$$\end{document} channel provides a unique opportunity to access the electric field of light with CEP-unstabilized pulses, we test the concept directly with a laser oscillator. In this experiment we utilize commercial tabletop mode-locked Yb-based laser (Pharos, Light Conversion). In contrast to the previous measurements, which require sophisticated infrastructure for the generation of short (sub-3 fs) and intense (~0.4 mJ) few-cycle CEP-stable laser pulses at 1 kHz repetition rate, we employ a simple mode-locked laser delivering ~150 fs, CEP-unstable, ~8 nJ pulses centered at 1030 nm wavelength and with a 75 MHz repetition rate (more in the Supplementary Information). Figure 2 shows the results of the experiment demonstrating that the fieldoscopic measurement previously feasible only with sophisticated infrastructure is now possible with simple table-top lasers.Fig. 2. Measurements of strong coherent states.a Measured electric fields of coherent states consisting of a large number of photons. b Central part of the pulse ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left\langle n\right\rangle$$\end{document}  = 6.42 x \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${10}^{6}$$\end{document} ) from the panel (a) showing clear electric field oscillations with high fidelity. c Benchmarking of the measured (blue dots) spectrum obtained by a Fourier transformation against a grating spectrometer (red curve). d Dependence of the measured signal ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\bar{{\rm{S}}}$$\end{document} ) on the electric field strength of the test pulse. e Dependence of the measured standard deviation (σ) on the electric field strength of the test pulse. f Dependence of the measured signal ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\bar{{\rm{S}}}$$\end{document} ) on the electric field strength of the sampling pulse

We study our concept in detail by measuring the electric field of the laser oscillator at various test and sampling pulse energies (Fig. 2a). For each measurement, we record a mean value of the signal and its standard deviation with a lock-in amplifier. Since our signal is expected to be proportional to the electric field, we expect a linear scaling of the mean signal with the field strength of the test pulse – that is confirmed by the experiment (Fig. 2d). On the other hand, three photons of the sampling pulse contribute to the heterodyne detection in the SHG/SFG ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m=n=2$$\end{document} ) channel (Eq. 5). We conduct a similar power scaling experiment by varying the field strength of the sampling pulse while keeping the test field unchanged (Fig. 2f).

In ideal environment with only quantum noise (due to the probabilistic nature of photons) and in the absence of classical noise (eg. jitter, electrical or environmental instabilities) one expects the standard deviation of the electric field amplitude to remain constant, regardless of the number of photons in the pulse^37^. In our experiment however, we observe that the standard deviation of the measurement scales linearly with the test field strength (Fig. 2e). This is expected since any practical system is a subject of the classical noise. We attribute this scaling to the dominance of the classical noise over quantum in the regime of strong (classical) coherent states consisting of many photons. This suggests that in this regime, the standard deviation is primarily limited by experimental stability, rather than the probabilistic nature of the photons. However, the contribution of the classical noise is expected to decrease as the test field (i.e., the number of photons) is decreased. Conversely, quantum noise—which is inherently associated with photon statistics – becomes more pronounced as the number of photons in the pulse is decreased. Therefore, we expect the quantum noise to be dominant in the regime of weak coherent states.

Fieldoscopy of weak (quantum) coherent states

Our approach not only allows us to utilize stable table-top lasers at high 75 MHz repetition rate, but it also removes the noise associated with the CEP fluctuations and highly non-linear light-matter interaction. Additionally, in contrast to broadband few cycle pulses (Fig. 1), narrow bandwidth laser oscillator pulses allow us to utilize thick nonlinear crystal (for higher signal strength) since the phase matching restriction is drastically eased. Having these advantages, we perform electric field measurement of very weak laser pulses to enter the regime of weak coherent states where quantum properties dominate over classical. Since our laser photon energy corresponds to about 1.2 eV (1030 nm wavelength) one enters this regime at pulse energies on the order of 193 zepto joules. We attenuate the test pulse down to the zepto joule level and perform the measurement. We measure the test field scaling as in Fig. 2d but in the zepto joule regime. Surprisingly, despite the fact that the pulse energy is of zepto joules only, the sensitivity of the detection remains similar to the strong (classical) pJ-scale regime. As discussed earlier, we attribute this to the absence of the CEP noise and low-order non-linear light-matter interaction. However, stability of the table-top laser, MHz repetition rate, thick nonlinear medium and eased phase matching condition are beneficial as well.

We note that in our classical regime, the measured signal and a corresponding standard deviation scale linearly with the test field strength (Fig. 2d, e). Since we are primarily interested in the deviation from the classical regime, in Fig. 3b, c we divide the mean and the standard deviation by the test field strength (~ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{\left\langle n\right\rangle }$$\end{document} ), so that in the classical regime these ratios are constant (which we further normalize to 1). Interestingly, when the mean photon number of the test pulse approaches 1 (~193 zJ), we observe a clear breakdown of the classical regime, both in the measured signal and the standard deviation. In particular, in the regime of weak coherent states, the scaling of the signal with the electric field becomes nonlinear. We note that in the classical description of electromagnetic waves, the scaling remains strictly linear regardless of the light field strength. The deviation from the linear scaling is a clear signature of the quantum nature of light, indicating that quantum properties (associated with coherence and photon statistics) become dominant and captured by our experiment.Fig. 3. Measurements of weak coherent states.a Measured electric fields of several weak coherent states consisting of a small number of photons. b Measured and simulated breakdown of the mean signal ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\bar{{\rm{S}}}$$\end{document} ) from the linear scaling. c Measured and simulated breakdown of the linear scaling of the measured standard deviation (σ). Poisson (d) and Bose-Einstein (e) probability distributions for mean photon numbers in (a)

Photon statistics analysis

Different quantum states of light are associated with unique photon statistics. In particular, coherent states follow the Poisson probability distribution, while non-coherent (thermal) light can be described by the Bose-Einstein probability distribution. We evaluate these distributions for the experimental mean photon numbers (Fig. 3d, e), with further details provided in Methods. As can be clearly seen, when the mean photon number \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left\langle n\right\rangle$$\end{document} approaches 1, the probability of the test pulse containing zero photons exponentially increases in both, Poisson, and Bose-Einstein cases.

In contrast to intensity-based measurements, direct access to the electric field serves as an observable even in cases where the test pulse contains no photons (i.e., the vacuum state). This can be intuitively understood as follows. In intensity-based measurements, reducing the photon number by half (or blocking every other pulse) would halve the detected signal. When measuring the field however, reducing photon number by half reduces the field amplitude by only ~ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{2}$$\end{document} (since the signal ∝ electric field, which scales with the square root of intensity), which is not equivalent to blocking every other pulse (more in the Supplementary Information).

Since our measurements are conducted with a laser, we expect the experiment to mainly follow the photon statistics of coherent light described by the Poisson probability distribution. We first model our experiment (more in Methods) based on the Poisson probability distribution. Indeed, our model—where the electric field serves as an observable—predicts a deviation from the linear classical regime when the mean photon number approaches 1 and below (Fig. 3b, c). For the scaling of the mean value of the signal (Fig. 3b), our model agrees with the experiment. Surprisingly, however, the scaling of the experimental standard deviation significantly disagrees with the model. In particular, as the mean photon number approaches 1, the standard deviation is expected to increase (Fig. 3c, red curve) reaching its maximum when the mean photon number is 1. This is because for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left\langle n\right\rangle$$\end{document}  = 1 the probability of a pulses in the vacuum state is close to 50% (Fig. 3d). Overall, the power scaling of the standard deviation based on Poisson statistics exhibits a distinct peak at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left\langle n\right\rangle$$\end{document}  = 1. In our experiment however, this feature is much less pronounced, which indicates a deviation from a perfectly coherent state.

To verify that the experimental observation is indeed a signature of the deviation from coherent properties, we re-model our experiment using the Bose-Einstein probability distribution (Fig. 3b, c, yellow curve). Interestingly, we find that the scaling of the mean value of the signal follows a similar trend to that predicted by Poisson statistics. This occurs because, although the Poisson and Bose-Einstein probability distributions describe different states of light, they appear similar in the regime of weak coherent states with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left\langle n\right\rangle$$\end{document} \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\approx$$\end{document} 1 (Fig. 3b, c, green and yellow bars). However, the scaling of the standard deviation according to the Bose-Einstein probability distribution differs significantly from that of Poisson statistics (Fig. 3c, yellow curve). In particular, the standard deviation monotonically decreases as the mean photon number decreases, which is similar to our experimental observation. Overall, our experiment appears to lie between of these two extremes—perfectly coherent and perfectly thermal states.

Intrapulse coherence

The observed deviation from perfectly coherent light could potentially be attributed to the variations in the coherence across the pulse. The direct access to the electric field in our experiment allows us to perform the same analysis for different regions of the pulse—the front, center, and tail. We utilize a Gabor method where a portion of the measured mean values and standard deviations of the pulse is multiplied by a Gaussian window as shown in Fig. 4a. The results of this analysis are presented in Fig. 4b, c. Regarding the mean value of the signal, we observe a trend similar to Fig. 3b, with no significant difference in scaling between different pulse regions. However, the standard deviation analysis reveals that different parts of the pulse exhibit different photon statistics. Notably, a characteristic feature of coherent Poisson statistics is clearly observed at the regions of the pulse with low intensity. As predicted by our model (Fig. 3c), the standard deviation peaks around \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left\langle n\right\rangle$$\end{document}  = 1, a hallmark signature of a coherent state described by Poisson statistics. However, this feature is much less pronounced for intense parts of the pulse. Since most of the pulse energy is concentrated in this region, the overall pulse coherence is dominated by this region. This explains why the overall pulse coherence (Fig. 3b, c) looks similar to the coherence of the most intense part of the pulse. As previously discussed, variations in intrapulse coherence are generally expected in intense mode-locked lasers.Fig. 4. Measurements of the intrapulse coherence.a Measured electric field of a weak coherent state ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left\langle n\right\rangle$$\end{document}  = 4.38) with corresponding Gaussian windows used for the evaluation of intrapulse photon statistics in (b) Measured breakdown of the linear scaling of the mean signal ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\bar{{\rm{S}}}$$\end{document} ). c Measured breakdown of the linear scaling of the standard deviation (σ). d Simulated breakdown of the linear scaling of the mean signal ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\bar{{\rm{S}}}$$\end{document} ). e Simulated breakdown of the linear scaling of the standard deviation (σ). The colors of the curves in (b, c) correspond to the Gaussian windows in (a)

Finally, we model intrapulse coherence as a time-varying combination of coherent and thermal states of light, described by the Poisson and Bose-Einstein probability distributions (Fig. 4d, e). We note a deviation from purely Poissonian statistics towards a mixed contribution of Bose–Einstein distribution, which increases with the number of photons in the pulse.

Discussion

The temporal dependence of the intrapulse coherence might have multiple origins. Although identification of the actual mechanisms requires further investigation, we attribute this deviation to an intensity-dependent decoherence mechanism (e.g., an elastic dephasing in the laser medium)—for example, due to fluctuations in the lasing transition or nonlinear coupling that effectively admixture a thermal component into the light field. In other words, the more intense parts of the pulse suffer slight decoherence. Such an effect could emerge from an elastic dephasing process \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\propto {\hat{a}}^{\dagger }\hat{a}$$\end{document} in the medium.

In summary, we presented a novel approach for sampling of ultrafast light fields without the need for CEP-stabilized pulses. We benchmarked our method against a state-of-the-art CEP-stable detection, finding identical results, thereby validating its efficacy. The technique enables direct measurement of the light electric field with unprecedented sensitivity, reaching about 880 yoctojoules of average pulse energy. For our ~150 fs pulses, this corresponds to about 20 yoctojoules of average energy per optical cycle. We experimentally demonstrate the concept for pulse energies up to 1.24 pJ, achieving a dynamic range of >90 decibels. The high sensitivity allows measuring a time varying electric field of both, strong (classical) and weak (quantum) states of light. We have shown that the direct access to the electric field in the weak coherent state regime enables the study of intrapulse coherence and time-dependent photon statistics with sub-cycle temporal resolution.

By extending petahertz fieldoscopy into the quantum regime, our study provides a powerful new tool for investigating coherence at the most fundamental level. Freed from CEP-stabilization, our approach can be implemented with simple tabletop lasers, greatly broadening access to field-sensitive measurements. By bridging ultrafast science with quantum optics, our work provides a new platform for exploring quantum coherence and light-matter interactions at sub-cycle timescales. Future experiments could enable the direct measurement of nonclassical states of light, such as entangled photons and squeezed vacuum in a single photon regime, to gain deeper insight into quantum decoherence mechanisms. This could advance quantum-enhanced metrology, sensing, and optical quantum information processing, where control over coherence and photon statistics is essential. Beyond practical applications, our concept represents a new paradigm for exploring light–matter interactions at the quantum level, enabling tests of fundamental quantum optics with unprecedented sub-cycle precision.

Materials and methods

Theoretical model

A coherent state of light is a quantum state of the electromagnetic field that can be described as the closest quantum analog to a classical light wave. It can be expressed in terms of Fock states (also known as photon number states), which are the eigenstates of the photon number operator:

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hat{n}={\hat{a}}^{\dagger }\hat{a}$$\end{document}

A coherent state \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$|{\alpha }$$\end{document} is defined as an eigenstate of the annihilation operator \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hat{a}$$\end{document} :

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hat{a}|{\alpha }{{\rangle }}=\alpha |{\alpha }{{\rangle }}$$\end{document}

where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha$$\end{document} is a complex number representing the amplitude and phase of the state. The magnitude \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${|\alpha |}^{2}$$\end{document} is the mean photon number of the coherent state. The coherent state \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$|{\alpha }{{\rangle }}$$\end{document} can be expanded as a superposition of Fock states \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${|n}{{\rangle }}$$\end{document} where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${n}$$\end{document} is the number of photons in the state:

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$|{\alpha }{{\rangle }}={e}^{-{\left|{\alpha }\right|}^{2}/2}\mathop{\sum }\limits_{n=0}^{\infty }\frac{{\alpha }^{n}}{\sqrt{n!}}{|n}{{\rangle }}$$\end{document}

The probability of finding \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n$$\end{document} photons in the coherent state is governed by the Poisson distribution

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${P}_{{Poisson}}\left(n\right)={\left\langle {n|}\alpha \right\rangle }^{2}=\frac{{|\alpha |}^{2n}}{n!}{e}^{-{|\alpha |}^{2}}$$\end{document}

with mean photon number \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left\langle n\right\rangle ={|\alpha |}^{2}$$\end{document} . The complex amplitude \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha =|\alpha |{e}^{i\varphi }$$\end{document} encodes both the intensity \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$({|\alpha |}^{2})$$\end{document} and phase ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varphi$$\end{document} ) of the light. Similarly, the photon statistics of the non-coherent thermal light obeys Bose-Einstein probability distribution:

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${P}_{{Bose}-{Einstein}}\left(n\right)=\frac{1}{(\left\langle n\right\rangle +1)}{\left(\frac{\left\langle n\right\rangle }{\left\langle n\right\rangle +1}\right)}^{n}$$\end{document}

In the Poisson distribution, the variance in the photon number equals the mean \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta {n}^{2}=\left\langle n\right\rangle ={|\alpha |}^{2}$$\end{document} , while the standard deviation \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma =\sqrt{n}$$\end{document} . In contrast, a thermal (Bose-Einstein) distribution has a variance larger than its mean (i.e., > \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left\langle n\right\rangle$$\end{document} ), indicating larger number fluctuations than a coherent state.

In our experiment, the strong sampling pulse of 63 pJ pulse energy and with a photon energy of 1.2 eV (1.9^-19^ Joules, 1030 nm wavelength) corresponds to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left\langle n\right\rangle =$$\end{document} 3.27∗10^9^. The test pulse energies and corresponding mean photon numbers \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left\langle n\right\rangle$$\end{document} utilized in the experiments on weak coherent states are summarized in the table below.Energy3.196.3713.5426.2853.35105.91212.61422.83845.6517043376.2zJ \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left\langle n\right\rangle$$\end{document} 0.01650.0330.07020.13630.27660.54911.10242.19244.38488.835717.5062

We statistically model the sampling pulse based on Eq. 4 with a mean photon number \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left\langle n\right\rangle =1000000$$\end{document} , while the test pulse is modeled with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left\langle n\right\rangle$$\end{document} from the table above.

We assume that the photon statistics of the sampling pulse follows the Poisson distribution \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{P}_{{sampling}}\left(n\right)=P}_{{Poisson}}\left(n\right)$$\end{document} , while the test pulse is a linear combination of the Poisson and Bose-Einstein distributions: \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${{P}_{{test}}\left(n\right)=A* P}_{{Poisson}}\left(n\right)+\left(1-A\right)* {P}_{{Bose}-{Einstein}}\left(n\right)$$\end{document} .

The experiment involves three stages: (1) Second harmonic generation from the sampling pulse; (2) Sum-frequency generation between the sampling and test pulses; (3) Interference of the second harmonic and the sum-frequency fields.

For modeling of the second harmonic, two numbers of photons ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${n}_{1}$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${n}_{2}$$\end{document} ) are randomly chosen from \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${P}_{{sampling}}\left(n\right)$$\end{document} . Since in the second harmonic generation two fundamental photons combine to create a photon with higher energy, the number of second harmonic photons is limited either by \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${n}_{1}$$\end{document} or \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${n}_{2}$$\end{document} . For instance, when \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${n}_{1}=10$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${n}_{2}=20$$\end{document} , \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${n}_{{SHG}}$$\end{document}  = 10.

In the sum-frequency generation two numbers of photons ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${n}_{1}$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${n}_{2}$$\end{document} ) are randomly chosen from \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${P}_{{sampling}}\left(n\right)$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${P}_{{test}}\left(n\right)$$\end{document} . Similarly, to the second harmonic generation, in the sum-frequency generation, a photon from the sampling pulse and a photon from the test pulse combine to create a sum-frequency photon. The number of second harmonic photons is limited either by \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${n}_{1}$$\end{document} or \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${n}_{2}$$\end{document} . For instance, when \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${n}_{1}=10$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${n}_{2}=20$$\end{document} , \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${n}_{{SFG}}$$\end{document} is 10.

Finally, the measured signal due to the SHG/SFG interference is modeled as

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$S\left({n}_{{SHG}},\,{n}_{{SFG}}\right)\approx 2* \sqrt{{n}_{{SHG}}}* \,\sqrt{{n}_{{SFG}}}$$\end{document}

where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$S$$\end{document} is a measured signal. With this we conduct a Monte Carlo simulation where we repeat this procedure 10000 times to obtain a statistical distributing of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$S$$\end{document} , from which we evaluate the mean ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\bar{{\rm{S}}}$$\end{document} ) and the standard deviation ( \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma$$\end{document} ).

Benchmarking measurements with few-cycle CEP-stabilized pulses

The laser beamline (more in the Supplementary Information) comprises a Ti:Sa oscillator (Rainbow 2, Spectra Physics), followed by chirped pulse amplification to 1 mJ pulse energy at 3 kHz repetition rate, further spectral broadening in a hollow-core fiber and a chirped mirror compressor. The experimental data acquisition was performed with a dual-phase lock-in amplifier (SR-830, Stanford Research Systems).

Measurements with CEP-unstabilized oscillator pulses

The laser beamline was a commercial tabletop mode-locked Yb-based laser (Pharos, Light Conversion).

Supplementary information

Supplementary materials for fieldoscopy at the quantum limit