Direct sampling of electric-field vacuum fluctuations

See allHide authors and affiliations

Science  23 Oct 2015:
Vol. 350, Issue 6259, pp. 420-423
DOI: 10.1126/science.aac9788

Probing the fluctuating vacuum

According to quantum mechanics, a vacuum is not empty space. A consequence of the uncertainly principle is that particles or energy can come into existence for a fleeting moment. Such vacuum or quantum fluctuations are known to exist, but evidence for them has been indirect. Riek et al. present an ultrafast optical based technique that probes the vacuum fluctuation of electromagnetic radiation directly.

Science, this issue p. 420


The ground state of quantum systems is characterized by zero-point motion. This motion, in the form of vacuum fluctuations, is generally considered to be an elusive phenomenon that manifests itself only indirectly. Here, we report direct detection of the vacuum fluctuations of electromagnetic radiation in free space. The ground-state electric-field variance is inversely proportional to the four-dimensional space-time volume, which we sampled electro-optically with tightly focused laser pulses lasting a few femtoseconds. Subcycle temporal readout and nonlinear coupling far from resonance provide signals from purely virtual photons without amplification. Our findings enable an extreme time-domain approach to quantum physics, with nondestructive access to the quantum state of light. Operating at multiterahertz frequencies, such techniques might also allow time-resolved studies of intrinsic fluctuations of elementary excitations in condensed matter.

Vacuum fluctuations give rise to a variety of phenomena, from spontaneous photon emission (1, 2) and the Lamb shift (3) via the Casimir force (4) to cosmological perturbations (5, 6). Representing the ground state, the quantum vacuum does not possess intensity. However, finite noise amplitudes of electric and magnetic fields should exist because of Heisenberg’s uncertainty principle. These fluctuations are best explained by analogy with a harmonic oscillator of mass m, resonance angular frequency Ω, and total energy

Embedded Image (1)

Quantization results in noncommuting operators for momentum p and displacement x. The Gaussian wave function of the ground state exhibits a root-mean-square (RMS) standard deviation of Δx = (ħ/2Ωm)1/2 (7, 8), where ħ is the reduced Planck constant. The total energy of a radiation field of wavevector k in free space, with electric and magnetic amplitudes E and B (respectively), and vector potential A in the Coulomb gauge is (9)

Embedded Image(2)

Considering one polarization direction and the transverse character of electromagnetic waves, Eq. 1 maps onto Eq. 2 by replacing x with A (amplitude of vector potential A), m with ε0V 0, vacuum permittivity; V, spatial volume), and Ω with ck ≡ Ω (c, speed of light; k= |kk|). Instead of x and p, an uncertainty product now links E and B or the amplitudes and phases of E, B, or A. An RMS amplitude of vacuum fluctuations ΔA = (ħ/2Ωε0V)1/2 results. In contrast to the mechanical case where Δx is known, understanding ΔA is less straightforward: Outside any cavities, there are no obvious boundaries that define a normalization volume V. This situation raises the question of whether direct measurement of the vacuum field amplitude in free space is physically meaningful and feasible.

The quantum properties of light (10) are typically analyzed either by photon correlation (1114), homodyning (1518), or hybrid measurements (19). In those approaches, information is averaged over multiple cycles, and accessing the vacuum state requires amplification. Femtosecond studies still rely on pulse envelopes that vary slowly relative to the carrier frequency (2023). In our work, we directly probed the vacuum noise of the electric field on a subcycle time scale using laser pulses lasting a few femtoseconds. In ultrabroadband electro-optic sampling (2427), a horizontally polarized electric-field waveform (red in Fig. 1A) propagates through an electro-optic crystal (EOX), inducing a change Δn of the linear refractive index n0 that is proportional to its local amplitude ETHz (Fig. 1A and fig. S1). The geometry is adjusted so that a new index ellipsoid emerges under 45° to the polarization of ETHz, with ny' and nx' = n0 ± Δn. An ultrashort optical probe pulse at a much higher carrier frequency νp (green in Fig. 1A; intensity, Ip; electric field, Ep) copropagates with ETHz at a variable delay time td. The envelope of Ip has to be on the order of half a cycle of light at the highest frequencies Ω/2π of ETHz that are detected. We used probe pulses as short as tp = 5.8 fs, corresponding to less than 1.5 optical cycles at νp = 255 THz (fig. S2). Upon passage through the EOX, the x' and y' components of Ep acquire a relative phase delay proportional to Δn and ETHz(td). The final polarization state of the probe is analyzed with ellipsometry. The differential photocurrent Δ I/I is proportional to the electric field ETHz(td). We used a radio-frequency lock-in amplifier (RFLA) for readout.

Fig. 1 Our experimental principle, with a demonstration of ultrahigh bandwidth and sensitivity.

(A) Scheme of electro-optic sampling of an electric-field waveform (red) by an ultrafast probe pulse (green), consisting of an EOX, a quarter-wave plate (λ/4), a Wollaston polarizer (WP), and a differential photocurrent detector (DD). (B) Classical electro-optic signal ΔI/I and corresponding electric-field amplitude versus delay time td (red line). The intensity envelope of the 5.8-fs probe pulse is shown in arbitrary units for comparison (green line). (C) Spectral multiterahertz amplitude (red) and phase (blue) obtained by Fourier transform.

We adjusted for optimal conditions to measure the vacuum signal by studying classical multiterahertz transients, which were synchronized to the probe (8). In Fig. 1B, Δ I/I is plotted in red against delay time td. Figure 1C shows the amplitude spectrum (red) and phase deviations (blue) within ±π, corroborating calculations (8) of an effective sampling bandwidth of Δν = ΔΩ/2π = 66 THz (figs. S3 and S4) around a center frequency of νc = Ωc/2π = 67.5 THz (free-space wavelength λc = 4.4 μm). The electric-field amplitude ĒTHz(td) is calibrated using (2830)

Embedded Image (3)

r41 denotes the electro-optic coefficient, and l is the thickness of the EOX. The amplitude response |Rc)| includes the pulse duration of the probe and velocity matching to the multiterahertz phase (8). The classical field transient in Fig. 1B was sampled with a signal-to-noise ratio better than 103 at a RFLA detection bandwidth set to 94 Hz. From the confocal amplitude trace and cross section, we estimated a mean photon number below 900 per pulse. This result proves the capability of our approach to characterize ultrabroadband coherent wave packets containing less than 10−3 photons, on average, within 1 s.

But can we directly access the ground state Φ0 of the radiation field? With the pump branch switched off, electro-optic phase shifts might still be caused by vacuum fluctuations copropagating with the probe. This effect should lead to a statistical distribution of the signal around the average of vac> = 0. The ground-state expectation value of the squared operator for the electric field in free space (31) yields the RMS amplitude

Embedded Image (4)

âΩ and âΩ are the operators for annihilation and creation of a photon with angular frequency Ω, respectively (ν, frequency; h, Planck constant). Because of the commutation relation [âΩ,âΩ] = 1, only âΩâΩ provides a nonvanishing contribution. Summing frequencies over our finite sensitivity interval ensures convergence of Eq. 4. The lateral extension of the volume V is now identified with the effective cross section Aeff, defined by the Gaussian intensity profile of the near-infrared probe beam inside the EOX. Theoretical modeling based on Laguerre-Gaussian modes (30) yields Aeff = w02π, where w0 is the probe spot radius (8). Because V = Aeff L, only the length L remains to be determined. Periodic boundary conditions are applicable when the EOX is short relative to the Rayleigh range of the multiterahertz transverse mode, resulting in a density of free-space modes L/c. Summing over all longitudinal modes within a bandwidth of Δν eliminates L, and we obtain

Embedded Image (5)

A factor of n0–1/2 accounts for dielectric screening inside the EOX (8). Thus, the vacuum amplitude is maximized when averaging over a minimal space-time volume, determined in the transverse directions by w0 = 4.25 μm (fig. S5). The longitudinal cross section cn0/(νcΔν) is defined by the Fourier transform of R(Ω), containing the intensity envelope of the 5.8-fs probe pulse and phase-matching conditions within the EOX (8).

Are such fluctuations discernible on top of the shot noise due to the Poissonian photon statistics of the coherent probe? An average number of Np = 5 × 108 photons detected per pulse causes a relative RMS shot-noise current of ΔISN/I = Np–1/2. With Eq. 3, we obtain the noise-equivalent field

Embedded Image (6)

Because the shot noise of the near-infrared probe, which is centered around νp, and the vacuum fluctuations at multiterahertz frequencies Ω are uncorrelated with each other and lack spectral overlap, the two contributions add up in quadrature. Therefore, the RMS width of the total detected noise distribution is expected to rise by a factor ofEmbedded Image (7)corresponding to a 4.7% increase, due to the multiterahertz vacuum noise.

To experimentally access the statistics of the quantum vacuum, we extended the RFLA bandwidth to 1.6 MHz and sampled the probability distribution of the electric field P(Etotal) every 5 μs. The contribution of the multiterahertz vacuum can be modified to discriminate against the shot-noise baseline by longitudinal or transverse expansion of the probed space-time volume (Eq. 5). In the first approach, we decreased νc and Δν by chirping the probe pulse to 100 fs (fig. S3), via translation of an SF10 prism in the compressor stage (Fig. 2A). A distinct reduction in peak counts around P(Etotal = 0) is observed when comparing the probability distribution obtained with the 5.8-fs probe (green in Fig. 2B) to the measurement with a stretched pulse (black). Also, the probabilities in the wings of the distribution including the multiterahertz vacuum (5.8-fs probe) are consistently higher than the corresponding values in the stretched-pulse distribution. The total change of the normalized noise amplitude amounts to 4%, in good agreement with the theoretical considerations underlying Eqs. 5 to 7. The red histogram in Fig. 2B emerges from a deconvolution algorithm that searches for the best link between distributions of P(Etotal) obtained with and without vacuum noise. This result directly mirrors the ground-state wave function |Ψ0(E)|2 of the electromagnetic field in the polarization plane and space-time volume that we probed. From |Ψ0(E)|2, a RMS standard deviation of ΔĒvac = 18 V/cm is obtained, in good agreement with the theoretical prediction of 20.2 V/cm in Eq. 5.

Fig. 2 Studying vacuum fluctuations via statistic readout and longitudinal modification of the probed space-time volume.

(A) Diagram showing longitudinal expansion of the probe volume: Stretching the sampling pulse from 5.8 fs (green) to 100 fs (black) causes temporal averaging over the vacuum field (red), leading to a reduction of the detected noise amplitude (tp, pulse duration; vgt, group velocity). (B) Normalized counting probability as a function of electro-optic readout by the short pulse (green) and long pulse (stretched to 100 fs; black). The deconvolved wave function |Ψ0|2 of the electric-field ground state is shown in red.

In the transverse option, we kept the short pulse duration and expanded the probe radius w0 by translating the EOX out of the confocal plane (Fig. 3A). Averaging over a larger cross section causes a decrease in fluctuation amplitude, which is projected onto the transverse mode of the gate. The effect of progressive narrowing is emphasized with differential probabilities obtained by subtracting a distribution at w0 = 4.25 μm from P(Etotal) sampled at increasing spot radii (Fig. 3B). When all original histograms are normalized, the maximum change in probability ΔP(Etotal = 0) of 0.04 ≡ 4% directly corresponds to the difference between the relative noise amplitudes measured with and without multiterahertz vacuum fluctuations, in quantitative agreement with Eq. 7. The dependence of the vacuum RMS amplitude on the transverse extension of the probed space-time volume is shown in Fig. 4. The normalized increase of total noise, measured with respect to bare shot noise (right vertical axis), is plotted against the probe spot radius w0 (blue squares). Conversion to the vacuum electric amplitude ΔEvac (left vertical axis) has been carried out analogously to |Ψ0(E)|2 in Fig. 2B. The functional dependence expected from Eq. 5 is shown as a red line. The inset in Fig. 4 illustrates the data recorded at low beam cross sections on a linear scale to highlight the hyperbolic increase of vacuum fluctuations for the smaller space-time volumes that we probed.

Fig. 3 Detection by transverse expansion of the space-time segment.

(A) Sketch showing the lateral increase of the sampling cross section, which leads to averaging over noise patterns within the circled areas. (B) Differential histograms obtained by subtracting the result for the confocal detector position with w0 = 4.25 μm from the results for positions with the beam diameter at 4.25 μm (black), 17 μm (purple), 25 μm (cyan), 50 μm (orange), and 85 μm (red).

Fig. 4 Dependence of the vacuum amplitude on transverse extension of the probed space-time volume.

Relative excess noise of electro-optic signal ΔI/I (right vertical axis) and RMS vacuum amplitude ΔEvac (left vertical axis) versus probe radius w0 (blue squares). Red lines represent a theoretical assessment based on Eq. 5.

In our study, we directly monitored vacuum fluctuations without amplifying them. The only effective part ∑Ω>0 âΩâΩ of the operator that extracts the variance of the field in Eq. 4 indicates that vacuum fluctuations correspond to photons, which spontaneously arise and vanish in the ground state Φ0. Time-energy uncertainty demands that virtual excitations have a limited life-time on the order of their oscillation cycle (32). The subcycle temporal resolution provided by the ultrashort probe ensures that we can directly detect effects originating from purely virtual photons. Phase-matched copropagation of the vacuum field and probe inside the EOX maximizes those signals. But does this measurement influence the quantum vacuum at all? Based on the electro-optic change of the refractive index Δnp~r41ETHz, the local multiterahertz field imprints a phase shift onto the ultrashort probe, which we detected. Because sum- and difference-frequency mixing occur simultaneously in this process (29), it requires no net transfer of energy, momentum, or angular momentum, and it even avoids modulation of the refractive index at frequencies Ω/2π << νp. Our second-order nonlinear element operates far from resonance. Virtual driving of the transitions avoids problems with decoherence, distinguishing our experiment from detection approaches in quantum optics or circuit quantum electrodynamics in which resonant two-level systems are involved (33). In consequence, our approach may be used to study the multiterahertz ground state while imposing negligible influence on it. Back-action might arise only in third order: The nonlinear refractive index n2 generates a local anomaly of phase velocity copropagating with the intensity envelope of the probe, because ΔnΩ ~ n2|Ep|2. When Np/w02tp suffices to induce phase shifts of the multiterahertz field during passage through the EOX, depletion of the vacuum amplitude in the sampled space-time volume and enhanced fluctuations in an adjacent interval are expected.

Supplementary Materials

Materials and Methods

Figs. S1 to S5

References (3440)


  1. Materials and methods are available as supplementary materials on Science Online.
  2. Acknowledgments: Support by the European Research Council (Advanced Grant 290876 “UltraPhase”), by Deutsche Forschungsgemeinschaft (SFB767), and by NSF via a postdoctoral fellowship for D.V.S. (award no. 1160764) is gratefully acknowledged.
View Abstract

Navigate This Article